HZOJ 礼物

其实是比较简单的一道期望状压dp，考试时一直在想数组表示概率，然而最后出的数总是小于一，于是无奈的把第一个点判掉放弃了其他点。

设f[i]为状态为i时到全部买到的期望次数，$f[i]=∑f[j]*p[k]+（1-∑p[k]）+1$,f[（1<<n）-1]=0,倒着推，k为j中的元素，$i|（1<<(k-1)）=j$，然后可以用高斯消元求解（显然会T啊），其实只要移项就可以了，，时间复杂度：O(2^N)。

 

#include<iostream>
#include<cstdio>
#include<bitset>
#define int LL
#define LL long long
using namespace std;
double f[1<<22],sp;
int n,w[21],sum;
double p[21];
signed main()
{
//	freopen("in.txt","r",stdin);

	cin>>n;
	for(int i=1;i<=n;i++)
	{
		cin>>p[i]>>w[i];sp+=p[i];sum+=w[i];
	}		
	f[(1<<n)-1]=0;
	for(int i=(1<<(n))-2;i>=0;i--)
	{
		double sum=0,spp=0;
		for(int k=1;k<=n;k++)
			if(((1<<(k-1))|i)!=i)
				spp+=p[k],sum+=p[k]*f[i|(1<<(k-1))];
		f[i]=(sum+1)/spp;
	}
	printf("%lld\n%0.3lf\n",sum,f[0]);
}

 

 还有一点就是要用longlong,int会wa50，虽然算着只是擦边，但是这种情况还是开longlong比较好。