結果

問題 No.620 ぐるぐるぐるりん
ユーザー pekempeypekempey
提出日時 2017-12-20 03:03:13
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 2,356 bytes
コンパイル時間 1,147 ms
コンパイル使用メモリ 88,204 KB
実行使用メモリ 10,752 KB
最終ジャッジ日時 2024-05-09 13:31:59
合計ジャッジ時間 6,783 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
9,088 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 WA -
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 WA -
testcase_11 AC 2 ms
5,376 KB
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 WA -
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 AC 2 ms
5,376 KB
testcase_25 AC 2 ms
5,376 KB
testcase_26 WA -
testcase_27 WA -
testcase_28 TLE -
testcase_29 -- -
testcase_30 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
#include <random>

using namespace std;

const double pi = acos(-1);
mt19937 mt(123456);

int t;
double p, w, v, gx, gy;

bool solve2() {
  double gr = hypot(gx, gy);
  double mn = 1e100;

  double x = v+1;
  double y = w;

  double det = (1+v)*(1+v) + w*w;

  double gx = ::gx;
  double gy = ::gy;

  for (int i = 0; i < t - 1; i++) {
    double gx0 = gx;
    double gy0 = gy;

    gx = gx0*(v+1) + gy0*w;
    gy = gy0*(v+1) - gx0*w;
    gx /= det;
    gy /= det;

    double dist = hypot(gx, gy);

    double score = dist*dist + x*x + y*y;
    if (score < p - 0.00001) {
      cerr << "ok" << endl;
      printf("%.16f %.16f\n", -x, -y);
      for (int j = 0; j < t - 3 - i; j++) {
        printf("%.16f %.16f\n", 0.0, 0.0);
      }
      printf("%.16f %.16f\n", gx, gy);
      printf("%.16f %.16f\n", 0.0, 0.0);
      return true;
    }

  }
  cerr << "ng" << endl;
  return false;
}

bool solve() {
  double gr = hypot(gx, gy);
  double mn = 1e100;

  for (int ii = 0; ii < 5000; ii++) {
    double x = 1;
    double y = 0;

    vector<double> vx(t);
    vector<double> vy(t);

    double det = hypot(1 + v, w);

    double theta = uniform_real_distribution<double>(0, 2*pi)(mt);
    double c = uniform_real_distribution<double>(1e-6, 10)(mt);
    vx[0] = cos(theta) - 1;
    vy[0] = sin(theta);

    for (int i = 0; i < t; i++) {
      double x0 = x;
      double y0 = y;

      x = x0 - y0*w + x0*v + vx[i];
      y = y0 + x0*w + y0*v + vy[i];

      double r = hypot(x, y);
      double tr;

      tr = gr*pow((double)i / t, c);
      double d = tr - r;

      vx[i] += d*x / r;
      vy[i] += d*y / r;

      x += d*x / r;
      y += d*y / r;
    }

    vx[t - 1] += gx - x;
    vy[t - 1] += gy - y;

    double s = 0;
    for (int i = 0; i < t; i++) {
      s += vx[i]*vx[i] + vy[i]*vy[i];
    }
    mn = min(mn, s);

    if (s < p - 0.001) {
      for (int i = 0; i < t; i++) {
        printf("%.16f %.16f\n", vx[i], vy[i]);
      }

      fprintf(stderr, "Sum: %.16f\n", s);
      fprintf(stderr, "p:   %.16f\n", p);

      return true;
    }
  }
  cout << "FAIL" << endl;
  cout << endl;
  cout << mn << endl;
  return false;
}

int main() {
  int T;
  cin >> T;
  while (T--) {
    cin >> t >> p >> w >> v >> gx >> gy;
    if (!solve2()) {
      solve();
    }
  }
}
0