2018.09.28【JoyOI1864】守卫者的挑战（期望DP）（背包问题）

传送门

---

解析：

本来这种题应该很快A掉，结果我读优打错调了半个多小时。。。

思路：

显然这是一个贪心的策略，那么我们就只考虑达到每个状态的概率。
记为$f_{i,j,k}$。表示当前已经进行了$i$次挑战，赢了$j$次，背包剩余容量为$k$的概率。然而这道题题目说的最不清楚的一点：拿到的碎片不一定要直接放入背包，也就是说，背包剩余容量可能为负。

当然你可以排序处理一下，我选择下标平移。

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

inline
int getint(){
	re int num;
	re char c;
	re bool f=0;
	while(!isdigit(c=gc()))f^=c=='-';num=c^48;
	while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
	return f?-num:num;
}

cs int N=205;
int n,L,K;
double f[2][205][610];
short p[205];
int a[205];

inline
int sign(int x){
	if(x<0)return max(-n,x)+N;
	return min(n,x)+N;
}

signed main(){
	n=getint(),L=getint(),K=getint();
	f[0][0][sign(K)]=1.0;
	for(int re i=1;i<=n;++i)p[i]=getint();
	for(int re i=1;i<=n;++i)a[i]=getint();
	
	for(int re i=1;i<=n;++i){
		int now=i&1,pre=now^1;
		memset(f[now],0,sizeof f[now]);
		for(int re j=0;j<=n;++j)
		for(int re k=-n;k<=n;++k){
			f[now][j][sign(k)]+=(f[pre][j][sign(k)])*(1.0-p[i]/100.0);
			f[now][j+1][sign(k+a[i])]+=f[pre][j][sign(k)]*p[i]/100.0;
		}
	}
	double ans;
	for(int re i=L;i<=n;++i)
	for(int re j=0;j<=n;++j)
	ans+=f[n&1][i][sign(j)];
	printf("%.6lf",ans);
	return 0;
}