結果
問題 | No.456 Millions of Submits! |
ユーザー | 🐬hec |
提出日時 | 2016-12-08 11:37:42 |
言語 | C++14 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 4,385 ms / 4,500 ms |
コード長 | 1,605 bytes |
コンパイル時間 | 2,007 ms |
コンパイル使用メモリ | 163,608 KB |
実行使用メモリ | 19,824 KB |
最終ジャッジ日時 | 2023-08-19 15:09:39 |
合計ジャッジ時間 | 11,419 ms |
ジャッジサーバーID (参考情報) |
judge10 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 3 ms
8,084 KB |
testcase_01 | AC | 2 ms
8,164 KB |
testcase_02 | AC | 3 ms
8,068 KB |
testcase_03 | AC | 3 ms
8,160 KB |
testcase_04 | AC | 2 ms
8,156 KB |
testcase_05 | AC | 3 ms
8,072 KB |
testcase_06 | AC | 3 ms
8,164 KB |
testcase_07 | AC | 7 ms
8,044 KB |
testcase_08 | AC | 7 ms
8,040 KB |
testcase_09 | AC | 47 ms
8,088 KB |
testcase_10 | AC | 47 ms
7,960 KB |
testcase_11 | AC | 439 ms
10,888 KB |
testcase_12 | AC | 4,385 ms
19,824 KB |
ソースコード
#include <bits/stdc++.h> #define _overload(_1,_2,_3,name,...) name #define _rep(i,n) _range(i,0,n) #define _range(i,a,b) for(int i=int(a);i<int(b);++i) #define rep(...) _overload(__VA_ARGS__,_range,_rep,)(__VA_ARGS__) #define _rrep(i,n) _rrange(i,n,0) #define _rrange(i,a,b) for(int i=int(a)-1;i>=int(b);--i) #define rrep(...) _overload(__VA_ARGS__,_rrange,_rrep,)(__VA_ARGS__) #define _all(arg) begin(arg),end(arg) #define uniq(arg) sort(_all(arg)),(arg).erase(unique(_all(arg)),end(arg)) #define getidx(ary,key) lower_bound(_all(ary),key)-begin(ary) #define clr(a,b) memset((a),(b),sizeof(a)) #define bit(n) (1LL<<(n)) #define popcount(n) (__builtin_popcountll(n)) using namespace std; template<class T>bool chmax(T &a, const T &b) { return (a<b)?(a=b,1):0;} template<class T>bool chmin(T &a, const T &b) { return (b<a)?(a=b,1):0;} using R=double; const R eps=1e-9; // Problem Specific Parameter: int a[1000010]; int b[1000010]; R t[1000010]; inline R f(R x,int i){return pow(x,a[i])*pow(log(x),b[i])-t[i];} inline R g(R x,int i){return pow(x,a[i]-1.0)*pow(log(x),b[i]-1.0)*(a[i]*log(x)+b[i]);} int main(void){ int m; scanf("%d",&m); rep(i,m) scanf("%d %d %lf",a+i,b+i,t+i); rep(i,m){ R ans=10.0; if(a[i]==0){ ans=exp(pow(t[i],1.0/b[i])); }else if(b[i]==0){ ans=pow(t[i],1.0/a[i]); }else{ R l=1.0,r=7.0; rep(hoge,10){ R m=(l+r)/2.0; if(f(m,i)<0.0) l=m; else r=m; } ans=(l+r)/2.0; //printf("%.12f\n",ans); rep(hoge,10){ R nxt=ans-f(ans,i)/g(ans,i); ans=nxt; } } printf("%.12f\n",ans); } return 0; } // n=e(t/n^a)^1/b