题目描述
传送门
题意:有一个大小为k的缓存区,每次从n种物品中按照一定的概率选取一种物品尝试放进去.同一个物品每一次选取的概率都是相同的.如果这种物品已经放进去过就不再放进去.如果缓存区满了就把放进去的时间离现在最远的物品拿出来.问10^100次后每个物品在缓冲区中的概率.
题解
10100
趋近于
+∞
,也就是说,做很多很多次之后,缓存区一定是满的,缓存区不满的概率趋近于0
因为操作的次数太多,实际上刚开始的操作都是没有用的,我们可以倒着考虑,问题就变成了每一次选一种加入,加满就停止的概率
有了这个转化,就可以dp了
f(i)表示缓存区的状态为i的概率,那么 f(i|2j−1)=∑j∉if(i)∗p(j)+f(i|2j−1)∗∑j∈i|2j−1p(j) ,移项整理得转移方程。
代码
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 25
double eps=1e-12;
int n,m,k;
double p[N],ans[N],cnt[1<<21],f[1<<21];
int main()
{
scanf("%d%d",&n,&k);m=n;
for (int i=0;i<n;++i)
{
scanf("%lf",&p[i]);
if (p[i]<eps) --m;
}k=min(k,m);
f[0]=1;
for (int i=0;i<1<<n;++i)
{
int cnt=0;double now=0;
for (int j=0;j<n;++j)
if ((i>>j)&1) ++cnt,now+=p[j];
if (cnt==k)
{
for (int j=0;j<n;++j)
if ((i>>j)&1) ans[j]+=f[i];
continue;
}
if (cnt>k) continue;
f[i]/=1-now;
for (int j=0;j<n;++j)
if (!((i>>j)&1)) f[i|(1<<j)]+=f[i]*p[j];
}
for (int i=0;i<n;++i) printf("%.16lf%c",ans[i]," \n"[i==n-1]);
}