結果
問題 |
No.1173 Endangered Species
|
ユーザー |
![]() |
提出日時 | 2020-08-14 22:31:18 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 105 ms / 2,000 ms |
コード長 | 1,152 bytes |
コンパイル時間 | 2,388 ms |
コンパイル使用メモリ | 195,736 KB |
最終ジャッジ日時 | 2025-01-13 00:14:09 |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 |
ソースコード
#include <bits/stdc++.h> using namespace std; using Int = long long; const char newl = '\n'; template<typename T1,typename T2> inline void chmin(T1 &a,T2 b){if(a>b) a=b;} template<typename T1,typename T2> inline void chmax(T1 &a,T2 b){if(a<b) a=b;} template<typename T> void drop(const T &x){cout<<x<<endl;exit(0);} template<typename T=int> vector<T> read(size_t n){ vector<T> ts(n); for(size_t i=0;i<n;i++) cin>>ts[i]; return ts; } struct Precision{ Precision(){ cout<<fixed<<setprecision(12); } }precision_beet; //INSERT ABOVE HERE signed main(){ cin.tie(0); ios::sync_with_stdio(0); int n; cin>>n; auto ps=read<double>(n); auto qs=read<double>(n); auto as=read(n); using D = double; auto calc=[&](D w)->D{ D res=0; for(int i=0;i<n;i++){ D tmp=(1.0-qs[i])/(1.0-qs[i]*w); res+=ps[i]*tmp; } // cout<<w<<':'<<res<<endl; return res; }; D L=0,R=1; // calc(L) > L // calc(R) <= R for(int k=0;k<100;k++){ D M=(L+R)/2; if(calc(M)<M) R=M; else L=M; } D ans=0; for(int i=0;i<n;i++) ans+=as[i]*log((1.0-qs[i])/(1.0-qs[i]*L)); cout<<ans<<newl; return 0; }