結果
問題 | No.456 Millions of Submits! |
ユーザー | koba-e964 |
提出日時 | 2016-12-08 01:29:26 |
言語 | C++11 (gcc 11.4.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,217 bytes |
コンパイル時間 | 2,641 ms |
コンパイル使用メモリ | 42,868 KB |
実行使用メモリ | 5,164 KB |
最終ジャッジ日時 | 2023-09-05 10:19:17 |
合計ジャッジ時間 | 10,183 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 231 ms
5,108 KB |
testcase_01 | AC | 231 ms
5,120 KB |
testcase_02 | AC | 231 ms
4,908 KB |
testcase_03 | AC | 230 ms
4,908 KB |
testcase_04 | WA | - |
testcase_05 | WA | - |
testcase_06 | WA | - |
testcase_07 | WA | - |
testcase_08 | WA | - |
testcase_09 | WA | - |
testcase_10 | WA | - |
testcase_11 | WA | - |
testcase_12 | WA | - |
ソースコード
#include <algorithm> #include <cmath> #include <cstdio> #define REP(i,s,n) for(int i=(int)(s);i<(int)(n);i++) using namespace std; const int RES = 2048; const int T = 100 * RES + 1; double memo_l[T]; /* 0.01 <= z <= 100 */ double lambert_W(double z, bool advised = true) { double lo = 0.0, hi = 100; int idx = z * RES; if (not advised) { if (log(z) >= 1) { lo = memo_l[idx - 1], hi = log(z); } else { lo = z / 2.719, hi = 1; } } if (advised) { lo = memo_l[idx], hi = memo_l[idx + 1]; } REP(i, 0, 45) { if (hi - lo <= 2e-10) { break; } double mid = (lo + hi) / 2; if (mid * exp(mid) <= z) { lo = mid; } else { hi = mid; } } if (not advised) { memo_l[idx] = lo; } return lo; } double solve(int a, int b, double t) { if (b == 0) { return pow(t, 1.0 / a); } if (a == 0) { return exp(pow(t, 1.0 / b)); } double res = lambert_W(pow(t, 1.0 / b) * a / b); return exp(res * b / a); } int main(void){ int m; scanf("%d", &m); REP(i, 0, T) { lambert_W(i / (double)RES, false); } REP(i, 0, m) { int a, b; double t; scanf("%d%d%lf", &a, &b, &t); printf("%.15f\n", solve(a, b, t)); } }