QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#162199#7110. Kuririn MIRACLEucup-team1198#WA 741ms159904kbC++202.0kb2023-09-03 06:30:532023-09-03 06:30:55

Judging History

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

  • [2023-09-03 06:30:55]
  • 评测
  • 测评结果:WA
  • 用时:741ms
  • 内存:159904kb
  • [2023-09-03 06:30:53]
  • 提交

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 = 20000000;


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

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 v, ld r) {
    ld prop = phi / pi * T;
    int i = floor(prop);

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

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

    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 * v / l;

    x1 *= fac;
    y1 *= fac;

    x1 -= v;

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

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

    cin >> u >> r >> d;

    if (5.246565891571807L * 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, u, r);
            if (rush(mp, r, u, d, t) >= 0) {
                rp = mp;
            } else {
                lp = mp;
            }
        }

        ld t = get(rp, u, r);
        cout << t + max((ld)0, rush(rp, r, u, d, t)) << '\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;
}

详细

Test #1:

score: 100
Accepted
time: 741ms
memory: 159848kb

input:

1
2.00 3 30.0

output:

8.310579903650451

result:

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

Test #2:

score: -100
Wrong Answer
time: 736ms
memory: 159904kb

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.457837103548061
5.488360794442070
5.848699763593381
1.123680241327300
43.266479874204259
8.305490103625902
2.377358490566038
4.015485959893929
37.948148148148148
18.652061855670103
4.155123634393019
27.178962926642404
9.625454545454545
63.311926605504587
16.613861386138614
15.093994414213707
1.017...

result:

wrong answer 39th numbers differ - expected: '5.3370511', found: '3.7511454', error = '0.2971502'