結果

問題 No.456 Millions of Submits!
ユーザー しらっ亭しらっ亭
提出日時 2016-12-08 02:47:36
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 934 ms / 4,500 ms
コード長 3,114 bytes
コンパイル時間 1,644 ms
コンパイル使用メモリ 175,156 KB
実行使用メモリ 31,968 KB
最終ジャッジ日時 2023-09-05 11:16:54
合計ジャッジ時間 5,283 ms
ジャッジサーバーID
(参考情報)
judge12 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 1 ms
4,380 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 2 ms
4,384 KB
testcase_07 AC 2 ms
4,376 KB
testcase_08 AC 3 ms
4,380 KB
testcase_09 AC 11 ms
4,376 KB
testcase_10 AC 11 ms
4,376 KB
testcase_11 AC 93 ms
6,428 KB
testcase_12 AC 934 ms
31,968 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;

#define _p(...) (void)printf(__VA_ARGS__)
#define forr(x,arr) for(auto&& x:arr)
#define _overload3(_1,_2,_3,name,...) name
#define _rep2(i,n) _rep3(i,0,n)
#define _rep3(i,a,b) for(int i=int(a);i<int(b);++i)
#define rep(...) _overload3(__VA_ARGS__,_rep3,_rep2,)(__VA_ARGS__)
#define _rrep2(i,n) _rrep3(i,0,n)
#define _rrep3(i,a,b) for(int i=int(b)-1;i>=int(a);i--)
#define rrep(...) _overload3(__VA_ARGS__,_rrep3,_rrep2,)(__VA_ARGS__)
#define all(x) (x).begin(), (x).end()
#define bit(n) (1LL<<(n))
#define sz(x) ((int)(x).size())
#define fst first
#define snd second
using ll=long long;using pii=pair<int,int>;
using vb=vector<bool>;using vs=vector<string>;
using vi=vector<int>;using vvi=vector<vi>;using vvvi=vector<vvi>;
using vl=vector<ll>;using vvl=vector<vl>;using vvvl=vector<vvl>;
using vd=vector<double>;using vvd=vector<vd>;using vvvd=vector<vvd>;
using vpii=vector<pii>;using vvpii=vector<vpii>;using vvvpii=vector<vvpii>;
template<class T,class U>ostream&operator<<(ostream&o,const pair<T,U>&p){o<<'('<<p.fst<<", "<<p.snd<<')';return o;}

double ans[1000001];

inline double POW(double x, int a) {
  if (a == 0) return 1;
  if (a == 1) return x;
  if (a == 2) return x * x;
  if (a == 3) return x * x * x;
  double x2 = x * x;
  if (a == 4) return x2 * x2;
  double x4 = x2 * x2;
  if (a == 5) return x4 * x;
  if (a == 6) return x4 * x2;
  if (a == 7) return x4 * x2 * x;
  if (a == 8) return x4 * x4;
  double x8 = x4 * x4;
  if (a == 9) return x8 * x;
  if (a == 10) return x8 * x2;
  if (a == -1) return 1 / x;
  cout << "a: "<<a << endl;
  assert(0);
}

void Main() {
  int M_;
  scanf("%d", &M_);
  vector<vector<pair<double, int>>> Q(200);
  rep(i, M_) {
    int a_, b_;
    double t;
    scanf("%d%d%lf", &a_, &b_, &t);
    Q[a_ * 16 +  b_].emplace_back(t, i);
  }

  rep(a, 11) rep(b, 11) {
    int k_ = a * 16 + b;
    if (sz(Q[k_]) == 0) continue;
    vector<pair<double, int>> &vq = Q[k_];

    if (a == 0 || b == 0) {
      forr(ti, vq) {
        double t = ti.fst;
        int idx = ti.snd;
        if (a == 0) ans[idx] = pow(M_E, pow(t, 1.0/b));
        else ans[idx] = pow(t, 1.0/a);
      }
      continue;
    }

    int n = sz(vq);

    sort(all(vq), greater<pair<double, int>>());

    //cout << "vq:"; rep(ii,sz(vq)) cout << ' ' << vq[ii]; cout << endl;
    double x = vq[0].fst + 10;

    const double EPS = 1e-10;

    rep(i, n) {
      double t = vq[i].fst;
      int idx = vq[i].snd;
      if (i && t == vq[i-1].fst) {
        ans[idx] = ans[vq[i-1].snd];
        continue;
      }

      rep(it, 500) {
        double lgx = log(x);
        double fx = POW(x, a) * POW(lgx, b) - t;
        double gx = POW(x, a-1) * POW(lgx, b-1) * (a * lgx + b);
        double d = fx / gx;
        x -= d;
        if (abs(d) < EPS) {
          //_p("%d: %d\n", i, it);
          break;
        }
        //if (it == 499) puts("!");
      }
      ans[idx] = x;
    }
  }

  rep(i, M_) {
    _p("%.12f\n", ans[i]);
  }
}

int main() {
  cin.tie(nullptr);
  ios::sync_with_stdio(false);
  cout << fixed << setprecision(10);
  Main();
  return 0;
}
0