結果

問題 No.456 Millions of Submits!
ユーザー しらっ亭しらっ亭
提出日時 2016-12-07 23:45:25
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 3,306 bytes
コンパイル時間 2,356 ms
コンパイル使用メモリ 183,336 KB
実行使用メモリ 36,192 KB
最終ジャッジ日時 2023-09-05 08:10:08
合計ジャッジ時間 11,294 ms
ジャッジサーバーID
(参考情報)
judge12 / judge11
このコードへのチャレンジ
(要ログイン)

テストケース

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

ソースコード

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];
const int DEPTH = 100;

double calc(double x, int a, int b) {
  return pow(x, a) * pow(log(x), b);
}

void dfs(int li, int ri, double mi, double ma, int a, int b, const vector<double> &T, const vector<int> &I, int depth) {
  //_p("dfs(%d,%d,%f,%f)\n",li,ri,mi,ma);
  if (li >= ri) return;
  if (li + 1 == ri) {
    rep(i, depth, DEPTH) {
      double mid = (mi + ma) / 2;
      bool c = calc(mid, a, b) <= T[li];
      if (c) mi = mid;
      else ma = mid;
    }
    ans[I[li]] = mi;
    return;
  }

  double mid = (mi + ma) / 2;
  double val = calc(mid, a, b);

  if (depth == DEPTH) {
    rep(i, li, ri) {
      ans[I[i]] = val;
    }
    return;
  }

  int th = lower_bound(all(T), val) - T.begin();
  //_p("li: %d, ri: %d, th: %d\n", li, ri, th);

  if (th >= ri) {
    dfs(li, ri, mi, mid, a, b, T, I, depth + 1);
  }
  else if (th <= li) {
    dfs(li, ri, mid, ma , a, b, T, I, depth + 1);
  }
  else {
    dfs(li, th, mi , mid, a, b, T, I, depth + 1);
    dfs(th, ri, mid, ma , a, b, T, I, depth + 1);
  }
}

void Main() {
  int T;
  scanf("%d\n", &T);
  map<pii, vector<pair<double, int>>> Q;
  rep(i, T) {
    int a, b;
    double t;
    scanf("%d%d%lf\n", &a, &b, &t);
    Q[make_pair(a, b)].emplace_back(t, i);
  }

  forr(query, Q) {
    int a = query.fst.fst;
    int b = query.fst.snd;

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

    int n = sz(query.snd);

    sort(all(query.snd));
    vd T(n);
    vi I(n);
    rep(i, n) {
      T[i] = query.snd[i].fst;
      I[i] = query.snd[i].snd;
    }
    //cout << "T:"; rep(ii,sz(T)) cout << ' ' << T[ii]; cout << endl;
    //cout << "I:"; rep(ii,sz(I)) cout << ' ' << I[ii]; cout << endl;

    double mi = 0;
    double ma = 10;

    //[li, ri)
    int li = 0;
    int ri = n;

    dfs(li, ri, mi, ma, a, b, T, I, 0);
  }

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

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