bzoj3029 守望者的挑战

传送门

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](1p[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);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值