守衛者的挑戰

_君の名は發表於2024-05-27

狀態轉移

我們假設 \(f{_i}{_,j}{_,k}\),表示前 \(i\) 場,贏了 \(j\)場,
目前揹包容量為 \(k\) 的機率,每一項挑戰有兩種狀態,勝或失敗,兩種情況答案不同,所以要分開計算,

失敗狀態:

\[f_{i,j,k}+=f_{i-1,j,k}*(1-p[i]) \]

成功狀態:

\[f_{i,j+1,k+a[i]}+=f_{i-1,j,k}*p[i],(k+a[i]>=0) \]

實現注意

在實現時我們需要注意某些地方,

  • 1:實現有正有負,要判非負再轉移
  • 2:揹包容量最多為 \(n\),(因為最多 \(n\) 個碎片)
  • 3:要按揹包容量排序,先拿揹包,再拿碎片,這樣一定能保證狀態轉移全,但不排序則會導致揹包容量不夠的情況(還有一些是把 \(n\) 加200來保證空間,兩種都可以)

注意3證明

實際上可以看一下題目,它只是要求 \(n\) 場完了之後容量夠即可,所以實際上挑戰中的順序是不重要的,排序保證了容量是夠的,如果揹包拿完了還裝不下碎片那就是不能了

#include<bits/stdc++.h>
using namespace std;
const int maxn=3e2+10; 
double f[210][210][210];
int n,l,k;
struct node
{
	int a;
	double p; 
}m[maxn];
bool cmp(node a,node b)
{
	return a.a>b.a;
}

int main()
{
	scanf("%d%d%d",&n,&l,&k);
	k=min(k,n);
	for(int i=1;i<=n;i++)scanf("%lf",&m[i].p),m[i].p/=(double)100;
	for(int i=1;i<=n;i++)scanf("%d",&m[i].a);
	sort(m+1,m+n+1,cmp);
	f[0][0][k]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<=i;j++)
		{
			for(int k=0;k<=n;k++)
			{
				f[i][j][k]+=f[i-1][j][k]*(1-m[i].p);
				int x=k+m[i].a;
				if(x<0)continue;
				x=min(x,n);
				f[i][j+1][x]+=f[i-1][j][k]*m[i].p;
			}
		}
	}
	double ans=0;
	for(int i=l;i<=n;i++)
		for(int j=0;j<=n;j++)ans+=f[n][i][j];
	printf("%.6lf\n",ans);
	
	return 0;
}

相關文章