EpicSevenAlgorithm/dp.cpp at main · zhaoyifan0528/EpicSevenAlgorithm · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
#include<iostream>
using namespace std;

long double f[40], g[40], fp[40]; // Ored 1pur
long double speedR[4] = {0.33223, 0.33223, 0.33223, 0.00332};
int hasV[5] = {2, 6, 10, 14, 14};
long double speedRp[4] = {0.03833, 0.34843, 0.34843, 0.26481};

void generatePur() {
    double long vp = (0.3/0.8) * 0.333333333 + (3.0/10.0) * 0.5;
    for (int i = 0; i < 40; ++i) fp[i] = 0;
    for (int i = 3; i <= 4; ++i) fp[i] = 1.0 * vp * speedRp[i - 1];

    for (int k = 1; k <= 4; ++k) {
        for (int i = 30; i >= 2; --i) {
            fp[i] *= (k == 4 ? 0.75 : 0.666666666666666667);
            for (int j = 5; j >= 2; --j) {
                if (i - j < 0) continue;
                fp[i] += (fp[i - j] * speedRp[j - 2]) *
                         (k == 4 ? 0.25 : 0.3333333333333);
            }
            if (i < hasV[k]) fp[i] = 0;
        }
    }

    for (int i = 34; i >= 17; --i) {
        fp[i] += fp[i + 1];
    }
    for (int i = 17; i <= 35; ++i) {
        cout << i << " " << 1.0 / fp[i] << endl;
    }
}

void generateRed() {
    double long vp = (0.4/0.8) * 0.3333333333 + (4.0/10.0) * 0.5;
    for (int i = 0; i < 40; ++i) f[i] = 0;
    for (int i = 2; i <= 5; ++i) f[i] = 1.0 * vp * speedR[i - 2];

    for (int k = 1; k <= 5; ++k) {
        for (int i = 35; i >= 2; --i) {
            f[i] = f[i] * 0.75;
            for (int j = 6; j >= 3; --j) {
                if (i - j < 0) continue;
                f[i] += (f[i - j] * speedR[j - 3]) * 0.25;
            }
            if (i < hasV[k - 1]) f[i] = 0;
        }
    }

    for (int i = 34; i >= 17; --i) {
        f[i] += f[i + 1];
    }
    for (int i = 17; i <= 35; ++i) {
        cout << i << " " << 1.0 / f[i] << endl;
    }
}

int main() {
    generateRed();
    generatePur();
    return 0;
}