2018.09.27 bzoj3029: 守卫者的挑战(概率dp)

183 篇文章 0 订阅
11 篇文章 0 订阅

传送门
概率dp经典题目。
直接 f [ i ] [ j ] [ k ] f[i][j][k] f[i][j][k]表示当前是第i次挑战,已经胜利了j次,目前的背包剩余空间是k。
然后用前面的转移后面的就行了。
注意第三维可能是负数,需要用一些技巧转化一下(比如把整个数组的下标向右平移)
代码:

#include<bits/stdc++.h>
#define N 205
using namespace std;
int n,l,K,a[N],tmp;
double p[N],f[2][205][605],ans;
inline int calc(int x){
	if(x<0)return x=max(-n,x),N+x;
	return x=min(n,x),N+x;
}
inline int read(){
	int ans=0,w=1;
	char ch=getchar();
	while(!isdigit(ch)){if(ch=='-')w=-1;ch=getchar();}
	while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
	return ans*w;
}
int main(){
	n=read(),l=read(),K=read();
	for(int i=1;i<=n;++i)p[i]=read()/100.0;
	for(int i=1;i<=n;++i)a[i]=read();
	f[0][0][calc(K)]=1;
	for(int i=0;i<n;++i){
		memset(f[tmp^1],0,sizeof(f[tmp^1]));
		for(int j=0;j<=n;++j)
			for(int k=-n;k<=n;++k){
				f[tmp^1][j+1][calc(k+a[i+1])]+=f[tmp][j][calc(k)]*p[i+1];
				f[tmp^1][j][calc(k)]+=f[tmp][j][calc(k)]*(1-p[i+1]);
			}
		tmp^=1;
	}
	for(int i=l;i<=n;++i)for(int j=0;j<=n;++j)ans+=f[tmp][i][calc(j)];
	printf("%.6lf",ans);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值