BZOJ 2784 时间流逝

BZOJ 2784 时间流逝

古典概率论...

可以发现由于能量圈数量限制，所以所构成的必定为树状结构（即便是转成最小能量圈和能量圈权值和之后存在重复状态，但是每个状态的含义不同，而且不能自身转移自身就没有问题，而且每个重复状态的本质是相同的）。

那么很显然，可以看作树上高斯消元解决。

剩下的就模板了.jpg

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <iostream>
#include <bitset>
using namespace std;
struct Segment
{
	double k,b;
	Segment(){}
	Segment(double x,double y){k=x,b=y;}
	Segment operator + (const Segment &a) const {return Segment(k+a.k,b+a.b);}
	Segment operator * (const double &a) const {return Segment(k*a,b*a);}
}f[55][55];
int n,T,a[55],vis[55][55];double p;
Segment dfs(int s,int x)
{
	if(s>T)return Segment(0,0);if(vis[s][x])return f[s][x];vis[s][x]=1;
	Segment ret=Segment(0,0);
	double np=s?p:0,nq=(1-np)/x;
	for(int i=1;i<=x;i++)ret=ret+dfs(s+a[i],i);
	return f[s][x]=Segment(np/(1-nq*ret.k),(ret.b*nq+1)/(1-nq*ret.k));
}
int main()
{
	while(scanf("%lf%d%d",&p,&T,&n)==3)
	{
		for(int i=1;i<=n;i++)scanf("%d",&a[i]);
		sort(a+1,a+n+1);memset(vis,0,sizeof(vis));
		printf("%.3lf\n",dfs(0,n).b);
	}
}