QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#749005#9576. Ordainer of Inexorable JudgmentPlentyOfPenalty#WA 0ms11784kbC++202.7kb2024-11-14 22:20:062024-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:10]
  • 评测
  • 测评结果:WA
  • 用时:0ms
  • 内存:11784kb
  • [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'