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 amazing_op.py

Amazing!!CSS 也能实现极光?

Amazing!巧用 CSS 视差实现酷炫交互动效

Python伊始——From Boring to Amazing

Why The Yeezy Boost is amazing and popular ?

EOJ Monthly 2021.9 Sponsored by TuSimple——A.Amazing.Discovery(分治or二次剩余)