QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#162218 | #7110. Kuririn MIRACLE | ucup-team1198# | WA | 555ms | 238248kb | C++20 | 1.9kb | 2023-09-03 06:46:45 | 2023-09-03 06:46:46 |
Judging History
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'