QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#749005 | #9576. Ordainer of Inexorable Judgment | PlentyOfPenalty# | WA | 0ms | 11784kb | C++20 | 2.7kb | 2024-11-14 22:20:06 | 2024-11-14 22:20:10 |
Judging History
你现在查看的是最新测评结果
- [2024-12-23 14:23:26]
- hack成功,自动添加数据
- (/hack/1303)
- [2024-12-06 11:32:56]
- hack成功,自动添加数据
- (/hack/1271)
- [2024-11-14 22:20:06]
- 提交
answer
#include <bits/stdc++.h>
#include <cmath>
#include <iomanip>
#define all(x) begin(x), end(x)
#define sz(x) ((int)(x).size())
#ifdef memset0
#define log(...) fprintf(stderr, __VA_ARGS__)
#else
#define endl '\n'
#define log(...) (void(0))
#endif
using namespace std;
using ll = long long;
using lf = long double;
using ull = unsigned long long;
const int N = 200011;
const lf pi = M_PI;
struct pt {
ll x, y;
pt(ll _x = 0, ll _y = 0) { x = _x, y = _y; }
pt operator-(const pt &you) { return pt(x - you.x, y - you.y); }
pt rotate() { return pt(-y, x); }
ll dot(pt you) { return (ll)x * you.x + (ll)y * you.y; }
ll len2() { return (ll)x * x + (ll)y * y; }
double angle() { return atan2(y, x); }
pt operator-() { return pt(-x, -y); }
} a[N];
lf dist(pt st, pt ed, pt p) {
pt r = p - st, v = ed - st;
return sqrtl(r.len2() - lf(r.dot(v)) * (r.dot(v)) / v.len2());
}
lf fx[N], U[N], L[N], R[N];
lf cover(lf l1, lf r1, lf l2, lf r2) {
if (l1 > l2) {
return cover(l2, r2, l1, r1);
}
return max(lf(0), r1 - l2);
}
int main() {
#ifdef memset0
freopen("M.in", "r", stdin);
#endif
cin.tie(0)->sync_with_stdio(0);
int n, x0, y0, r, _t;
cin >> n >> x0 >> y0 >> r >> _t;
lf t = _t;
pt O(0, 0);
for (int i = 1; i <= n; ++i) {
cin >> a[i].x >> a[i].y;
if (a[i].len2() <= (ll)r * r) {
log("touched!!\n");
printf("%.9lf\n", double(t));
return 0;
}
}
// sort(a + 1, a + n + 1);
for (int i = 1; i <= n; ++i) {
U[i] = atan2l(a[i].y, a[i].x);
if (U[i] < 0) U[i] += 2 * pi;
lf ang = asinl(r / sqrtl(a[i].len2()));
L[i] = U[i] - ang;
R[i] = U[i] + ang;
if (L[i] < 0) {
L[i] += 2 * pi;
R[i] += 2 * pi;
}
}
lf ansl = 1000, ansr = -1000;
for (int i = 1; i <= n; i++) {
int f = i;
int g = i % n + 1;
lf diff = U[g] - U[f];
if (diff < 0) diff += 2 * pi;
lf l, r;
if (diff <= pi) {
l = L[f];
r = R[g];
} else {
l = L[g];
r = R[f];
}
log("f=%d g=%d l=%.5Lf r=%.5Lf\n", f, g, l, r);
if (r < l) {
r += 2 * pi;
}
ansl = min(ansl, l);
ansr = max(ansr, r);
}
log(">> ansl=%.10Lf ansr=%.10Lf\n", ansl, ansr);
if (ansr - ansl >= 2 * pi) {
log("all covered!!\n");
printf("%.9lf\n", double(t));
return 0;
}
lf len = ansr - ansl;
lf cir = floor(t / (2 * pi));
lf ans = len * cir;
lf s = atan2l(y0, x0);
if (s < 0) {
s += 2 * pi;
}
lf e = s + (t - 2 * pi * cir);
// cerr << "cir: " << fixed << setprecision(12) << cir << endl;
ans += cover(ansl, ansr, s, e);
ans += cover(ansl, ansr, s + 2 * pi, e + 2 * pi);
cout << fixed << setprecision(12) << ans << endl;
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 0
Wrong Answer
time: 0ms
memory: 11784kb
input:
3 1 0 1 1 1 2 2 1 2 2
output:
1.570796326795
result:
wrong answer 1st numbers differ - expected: '1.0000000', found: '1.5707963', error = '0.5707963'