守卫者的挑战(据说在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\n",ans);
20     return 0;
21 }

 

转载于:https://www.cnblogs.com/12mango/p/7711524.html