Codeforces 623D [Amazing概率题]
Posted RDCのACM奇幻之旅
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Codeforces 623D [Amazing概率题]相关的知识,希望对你有一定的参考价值。
很有趣的一道题吖!
做法:贪心+迭代
- Sigma(i*(pr[i]-pr[i-1])))=n-sigma(pr[i]), 所以我们贪心地是pr[i]尽可能大。
- 也就是让pr[i]/pr[i-1]尽可能大。
这种类型的贪心,和17EC-Final的那个最小化方差的题挺相似的。
#include <iostream> #include <algorithm> #include <cstdio> #include <vector> #include <cmath> #include <cstring> #include <map> #include <set> using namespace std; typedef long long LL; #define rd(x) scanf("%d",&x) #define prt(x) printf("%d\n", x); #define prtvec(v) for(int i=0;i<v.size();i++) printf("%d%c", v[i], i==(v.size()-1)?‘\n‘:‘ ‘); #define sz(x) (int)x.size() #define pb(x) push_back(x) #define rep(i,x,y) for(int i=x;i<=y;i++) #define per(i,y,x) for(int i=y;i>=x;i--) const int N=300000+10; const double EPS = 1e-8; int n; double f[N][102],g[N],p[N]; int main(){ rd(n); rep(i,1,n) scanf("%lf",&p[i]), p[i]/=100.0; rep(i,1,n) f[n][i] = p[i]; rep(i,n+1,N-1) { double mx = 0; int bst = -1; rep(j,1,n) { double tmp = (1-f[i-1][j])*p[j]/f[i-1][j]; if (tmp > mx) { mx=tmp, bst = j; } } rep(j,1,n) f[i][j]=f[i-1][j]; f[i][bst]=(1-f[i-1][bst])*p[bst] + f[i-1][bst]; } double res = 0; rep(i,n,N-1) { g[i]=1; rep(j,1,n) g[i]*=f[i][j]; res=res+i*(g[i]-g[i-1]); } printf("%.5f\n", res); }
以上是关于Codeforces 623D [Amazing概率题]的主要内容,如果未能解决你的问题,请参考以下文章
Python伊始——From Boring to Amazing
Why The Yeezy Boost is amazing and popular ?
EOJ Monthly 2021.9 Sponsored by TuSimple——A.Amazing.Discovery(分治or二次剩余)