QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#209976 | #7110. Kuririn MIRACLE | PetroTarnavskyi | WA | 418ms | 4208kb | C++17 | 1.5kb | 2023-10-10 20:38:52 | 2023-10-10 20:38:53 |
Judging History
answer
#include <bits/stdc++.h>
using namespace std;
#define FOR(i, a, b) for(int i = (a); i < (b); i++)
#define RFOR(i, a, b) for(int i = (a) - 1; i >= (b); i--)
#define SZ(a) int(a.size())
#define ALL(a) a.begin(), a.end()
#define PB push_back
#define MP make_pair
#define F first
#define S second
typedef long long LL;
typedef vector<int> VI;
typedef pair<int, int> PII;
typedef double db;
const int M = 200000;
const db PI = acos(-1.0);
const db EPS = 1e-7;
db v, r, d;
db f(db phi)
{
db si = sin(phi);
return -v / (2 * r) * (sqrt(3 + si * si) - si);
}
void solve()
{
cin >> v >> r >> d;
db ans = d / v;
db phi = PI;
const db dt = ans / M;
FOR(i, 1, M + 1)
{
db t = i * dt;
phi += dt * f(phi);
db si = sin(phi), co = cos(phi);
db a = d - 2 * r - v * t, b = a - 2 * r * co;
if (abs(b) < EPS && phi < EPS)
{
ans = t;
break;
}
if (b < -EPS)
break;
if (phi < PI / 2)
{
db l = sqrt(4 * r * r + a * a - 4 * r * a * cos(phi));
db siBeta = 2 * r * si / l;
assert(0 <= siBeta && siBeta < 1 + EPS);
db beta = asin(min(1.0, siBeta));
if (phi + beta <= PI / 2 && 2 * cos(phi + beta) >= cos(phi))
{
ans = min(ans, t + l / (2 * v));
break;
}
}
else if (phi < EPS)
{
ans = t + a / (2 * v);
break;
}
}
cout << ans << "\n";
//cerr << ans - 8.310579933902352 << "\n";
}
int main()
{
ios::sync_with_stdio(0);
cin.tie(0);
cout << fixed << setprecision(15);
int t;
cin >> t;
while (t--)
solve();
return 0;
}
详细
Test #1:
score: 100
Accepted
time: 5ms
memory: 4208kb
input:
1 2.00 3 30.0
output:
8.310572894044750
result:
ok found '8.3105729', expected '8.3105799', error '0.0000008'
Test #2:
score: -100
Wrong Answer
time: 418ms
memory: 4024kb
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.457832178995005 5.488357405882677 5.848699763593380 1.123680241327300 43.266428662029696 8.305481819823621 2.377358490566038 4.015483707940933 37.948148148148142 18.652061855670105 4.155118123835539 27.178927696135439 9.625454545454545 63.311926605504588 16.613861386138616 15.093983558138465 1.017...
result:
wrong answer 1st numbers differ - expected: '4.4578371', found: '4.4578322', error = '0.0000011'