QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#162218#7110. Kuririn MIRACLEucup-team1198#WA 555ms238248kbC++201.9kb2023-09-03 06:46:452023-09-03 06:46:46

Judging History

你现在查看的是最新测评结果

  • [2023-09-03 06:46:46]
  • 评测
  • 测评结果:WA
  • 用时:555ms
  • 内存:238248kb
  • [2023-09-03 06:46:45]
  • 提交

answer

#include <bits/stdc++.h>

#pragma GCC optimize("fast-math")

using namespace std;
#define ll long long
#define pii pair<int, int>
#define ld long double
#define all(a) (a).begin(), (a).end()

const ld pi = atan2l(1, 0) * 2;
const int T = 15000000;


ld get_dphi(ld phi, ld u, ld r) {
    ld sinp = sinl(phi);
    return (- u * sinp + sqrtl(u * u * (sinp * sinp + 3))) / (2 * r);
}

long double tphi[T + 2];

void precalc() {
    ld cur_t = 0;
    for (int i = 0; i < T; ++i) {
        ld phi = pi * i / T;
    
        ld dphi = get_dphi(phi, 1, 1);  // at further point better
        cur_t += pi / T / dphi;
        tphi[i + 1] = cur_t;
    }
}

ld get(ld phi, ld r) {
    ld prop = phi / pi * T;
    int i = floor(prop);

    return r * (tphi[i] * (1 - (prop - i)) + tphi[i + 1] * (prop - i));
}

long double rush(ld phi, ld r, ld d, ld t) {
    d -= 2 * r + t;

    ld x = - 2 * r * cosl(phi);
    ld y = 2 * r * sinl(phi);

    ld x1 = d - x;
    ld y1 = y;

    ld l = sqrtl(x1 * x1 + y1 * y1);
    ld fac = 2 / l;

    x1 *= fac;
    y1 *= fac;

    x1 -= 1;

    if (x1 * x > y * y1 - 1e-7) {
        // cerr << l << '\n';
        return l / 2;
    } else {
        return -1;
    }
}

void solve() {
    ld u, r, d;

    cin >> u >> r >> d;

    if (tphi[T] * r <= d - 4 * r) {
        ld lp = pi / 2;
        ld rp = pi;

        for (int i = 0; i < 300; ++i) {
            ld mp = (lp + rp) / 2;

            ld t = get(mp, r);
            if (rush(mp, r, d, t) >= 0) {
                rp = mp;
            } else {
                lp = mp;
            }
        }

        ld t = get(rp, r);
        cout << (t + max((ld)0, rush(rp, r, d, t))) / u << '\n';

    } else {
        cout << d / u << '\n';        
    }
}

signed main() {
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);

    precalc();

    cout << fixed << setprecision(15);

    // cout << get(pi, 1, 1) << '\n';

    int t = 1000;
    cin >> t;
    while (t--) solve();


    return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Test #1:

score: 100
Accepted
time: 554ms
memory: 238116kb

input:

1
2.00 3 30.0

output:

8.310579893565113

result:

ok found '8.3105799', expected '8.3105799', error '0.0000000'

Test #2:

score: -100
Wrong Answer
time: 555ms
memory: 238248kb

input:

100
8.66 6.05 71.59
8.44 8.67 82.55
4.23 9.68 24.74
6.63 9.95 7.45
1.2 7.27 97.55
2.08 2.94 31.61
8.48 5.56 20.16
7.35 5.55 52.48
1.35 8.15 51.23
3.88 9.02 72.37
7.03 1.97 57.03
1.78 4.44 93.44
5.5 6.11 52.94
1.09 9.64 69.01
1.01 3.54 16.78
1.63 4.55 44.06
5.82 5.51 5.92
5.62 2.93 76.01
4.75 4.43 73...

output:

4.457837097151715
5.488360789590511
5.848699763593381
1.123680241327300
43.266479814427429
5.780351146592027
2.377358490566038
3.611937569814390
37.948148148148148
18.652061855670103
4.155123631331962
27.178962899829640
9.625454545454545
63.311926605504587
16.613861386138614
12.794292498623896
1.017...

result:

wrong answer 6th numbers differ - expected: '8.3054901', found: '5.7803511', error = '0.3040325'