結果
問題 | No.1173 Endangered Species |
ユーザー |
![]() |
提出日時 | 2020-04-24 02:32:52 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,001 ms / 2,000 ms |
コード長 | 1,248 bytes |
コンパイル時間 | 2,595 ms |
コンパイル使用メモリ | 196,988 KB |
最終ジャッジ日時 | 2025-01-09 22:56:28 |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define rep2(i, m, n) for(int i=int(m); i<int(n); ++i) #define rep(i, n) rep2(i, 0, n) using ll = long long; using ld = long double; using V = vector<int>; using Vld = vector<ld>; struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_; template<typename T> istream &operator>>(istream &is, vector<T> &v) { for (auto &e : v) is >> e; return is; } template<typename T> ostream &operator<<(ostream &os, const vector<T> &v) { for (auto &e : v) os << e << " "; return os; } const ld EPS = 1e-14; ld log_prob(ld a, ld z) { return log(1-a) - log(1-a*z); } int main() { ll n; cin >> n; Vld p(n), q(n); cin >> p >> q; V x(n); cin >> x; auto fun = [&] (ld z) { ld f = 0; rep(i, n) f += p[i] * exp(log_prob(q[i], z)); // rep(i, n) f += p[i] * (1-q[i]) / (1-q[i]*z); // cout << z << " " << f << endl; return (z < f); }; ld ok = 0; ld ng = 1; while (abs(ok - ng) > EPS) { ld md = (ok + ng) / 2; if (fun(md)) ok = md; else ng = md; // cout << md << endl; } ld ans = 0; rep(i, n) ans += x[i] * log_prob(q[i], ok); cout << ans << endl; // cout << exp(ans) << endl; return 0; }