結果
問題 | No.456 Millions of Submits! |
ユーザー | koba-e964 |
提出日時 | 2016-12-08 01:06:55 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 3,815 ms / 4,500 ms |
コード長 | 1,120 bytes |
コンパイル時間 | 600 ms |
コンパイル使用メモリ | 42,912 KB |
実行使用メモリ | 16,264 KB |
最終ジャッジ日時 | 2023-08-27 04:35:08 |
合計ジャッジ時間 | 44,438 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge11 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2,348 ms
16,160 KB |
testcase_01 | AC | 2,339 ms
16,192 KB |
testcase_02 | AC | 2,350 ms
16,192 KB |
testcase_03 | AC | 2,352 ms
16,264 KB |
testcase_04 | AC | 2,351 ms
16,224 KB |
testcase_05 | AC | 2,342 ms
16,132 KB |
testcase_06 | AC | 2,343 ms
16,264 KB |
testcase_07 | AC | 2,350 ms
16,140 KB |
testcase_08 | AC | 2,341 ms
16,136 KB |
testcase_09 | AC | 2,361 ms
16,132 KB |
testcase_10 | AC | 2,368 ms
16,200 KB |
testcase_11 | AC | 2,497 ms
16,196 KB |
testcase_12 | AC | 3,815 ms
16,192 KB |
ソースコード
#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 T = 100 * 16384 + 1; double memo_l[T]; /* 0.01 <= z <= 100 */ double lambert_W(double z, bool advised = true) { double lo = 0.0, hi = 100; if (advised) { int idx = z * 16384; lo = memo_l[idx], hi = memo_l[idx + 1]; } REP(i, 0, 45) { if (hi - lo <= 5e-12) { break; } double mid = (lo + hi) / 2; if (mid * exp(mid) <= z) { lo = mid; } else { hi = mid; } } if (not advised) { int idx = z * 16384; memo_l[idx] = lo; } return lo; } double solve(int a, int b, double t) { int idx = floor(t * 10000 + 0.5); 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 / 16384.0, false); } REP(i, 0, m) { int a, b; double t; scanf("%d%d%lf", &a, &b, &t); printf("%.15f\n", solve(a, b, t)); } }