用
f
[
i
]
[
j
]
[
k
]
f[i][j][k]
f[i][j][k]表示:进行了
i
i
i轮战斗,赢了
j
j
j轮,当前背包容量为
k
k
k的概率。
当背包容量大于
N
N
N时,注意到一定能装下所有的地图。
所以容量一维最高只用
N
m
a
x
=
200
N_{max}=200
Nmax=200
转移很简单——
f [ i + 1 ] [ j + 1 ] [ k + v a l [ i + 1 ] ] + = f [ i ] [ j ] [ k ] ∗ p [ i ] f[i+1][j+1][k+val[i+1]]+=f[i][j][k]*p[i] f[i+1][j+1][k+val[i+1]]+=f[i][j][k]∗p[i]
f [ i + 1 ] [ j ] [ k ] + = f [ i ] [ j ] [ k ] ∗ ( 1 − p [ i ] ) f[i+1][j][k]+=f[i][j][k]*(1-p[i]) f[i+1][j][k]+=f[i][j][k]∗(1−p[i])
v a l [ i ] val[i] val[i]表示第 i i i场战斗胜利的贡献, p [ i ] p[i] p[i]表示第 i i i场胜利的概率。
最后要求的就是
∑
j
∈
[
L
,
N
]
,
k
∈
[
0
,
N
]
f
[
N
]
[
j
]
[
k
]
\sum_{j∈[L,N],k∈[0,N]}f[N][j][k]
j∈[L,N],k∈[0,N]∑f[N][j][k]
这里的递推可以用滚动数组优化。对于负数下标,平移一下就行了。就是给所有下标加上一个值让他们都变成正的。
#include<bits/stdc++.h>
using namespace std;
const int maxn=205;
const int mx=201;
double f[2][maxn][maxn+mx],p[maxn],ans=0;
int N,L,K,val[maxn],nxt=1,now=0;
int main(){
scanf("%d%d%d",&N,&L,&K);
for(int i=1;i<=N;++i) scanf("%lf",&p[i]),p[i]*=0.01;
for(int i=1;i<=N;++i) scanf("%d",&val[i]);
K=min(K,N),f[0][0][K+mx]=1;
for(int i=0;i<N;++i){
memset(f[nxt],0,sizeof(f[nxt]));
for(int j=0;j<=i;++j)
for(int k=0;k<=N+mx;++k){
if(k+val[i+1]>=0)
f[nxt][j+1][min(N+mx,k+val[i+1])]+=f[now][j][k]*p[i+1];
f[nxt][j][k]+=f[now][j][k]*(1.00-p[i+1]);
}
nxt^=1,now^=1;
}
for(int i=L;i<=N;++i)
for(int j=0;j<=N;++j)
ans+=f[now][i][j+mx];
printf("%.6lf",ans);
}