守卫者的挑战(据说在bzoj有但我没找到)

芒果君:一看就是概率dp(可是我不会啊,就算再裸也不会啊)。然后先从最后想,能够满足题意的状态是 挑战次数>=L,获得价值>=0,那一定有f[总挑战数i][挑战成功数j][价值k]。转移很好写,无非两种情况,成功和不成功。不成功,i+1,其他不变,值*(1-p[i]);成功,i+1,j+1,k+w,值是 ∑各种情况*p[i]。题目中有个关键信息是n<=200,也就是说极限情况是200* ±1,我们只需要开200*2的数组,将负数移个位。空间炸了怎么办,滚一滚啊。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 using namespace std;
 5 double ans,p[210],f[2][210][410];
 6 int n,t,m,w;
 7 int main()
 8 {
 9     scanf("%d%d%d",&n,&t,&m);
10     m=min(m,n);
11     for(int i=1;i<=n;++i) scanf("%lf",&p[i]),p[i]/=100.0;
12     f[0][0][n+m]=1;
13     for(int i=1;i<=n;++i){
14         scanf("%d",&w);
15         for(int j=0;j<=n;++j) for(int k=0;k<=(n<<1);++k) f[i&1][j][k]=f[i&1^1][j][k]*(1-p[i]);
16         for(int j=0;j<n;++j) for(int k=1;k<=(n<<1);++k) f[i&1][j+1][min(k+w,n<<1)]+=f[i&1^1][j][k]*p[i]; 
17     }
18     for(int i=t;i<=n;++i) for(int j=n;j<=(n<<1);++j) ans+=f[n&1][i][j];
19     printf("%.6lf
",ans);
20     return 0;
21 }
原文地址:https://www.cnblogs.com/12mango/p/7711524.html