n<=20种东西,有个大小k<=n的箱子,每次会以固定的概率从所有东西里选一种,若箱子里有空位且这种东西没出现过就丢进去,若箱子满了且这种东西没出现过就把最早访问过的一个丢掉,(只要在每次操作最早的“挑一种物品”选中某个数都算“访问”)问10^100操作后每个数存在箱子里的概率。
由于操作过多,最后箱子可以看成满的。可以发现由于操作跟时间有关,只有最后新进入箱子的k种东西在进入箱子时会对答案产生影响,因此只用看后面几次操作即可。可以把这个过程等同于:从头开始拿,问刚把箱子填满时的概率。
在后面的操作中,虽然只看最后k种东西被访问的过程,但期间也有可能多次访问到同种东西,咋办呢?状压,f(i)表示状态i出现的概率,$f(i)=\sum f(j)*a_{j关于i的补集}+\sum_{k\epsilon i} a_k*f(i)$,错误!最终状态不会再自己转移到自己,所以判一下i是不是终态即有没有k个数,没有才把后面那部分算上去。
还错!数据中有一些是箱子不可能满,所以要把箱子大小和概率>0的种数取个min。
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
1 #include<iostream> 2 #include<cstring> 3 #include<cstdlib> 4 #include<cstdio> 5 //#include<bitset> 6 #include<algorithm> 7 #include<cmath> 8 using namespace std; 9 10 int n,K; 11 #define maxn 1100011 12 double a[33],ans[maxn],eve[33]; 13 int main() 14 { 15 scanf("%d%d",&n,&K); 16 int cnt=n; 17 for (int i=1;i<=n;i++) scanf("%lf",&a[i]),cnt-=a[i]<1e-10?1:0; 18 K=min(cnt,K); 19 ans[0]=1; for (int i=1;i<=n;i++) eve[i]=0; 20 for (int i=1;i<(1<<n);i++) 21 { 22 int one=0; ans[i]=0; double kk=0; 23 for (int j=1;j<=n;j++) if ((i>>(j-1))&1) one++,ans[i]+=ans[i^(1<<(j-1))]*a[j],kk+=a[j]; 24 if (one==K) for (int j=1;j<=n;j++) if ((i>>(j-1))&1) eve[j]+=ans[i]; 25 ans[i]/=(1-kk); 26 } 27 for (int i=1;i<=n;i++) printf("%.8lf ",eve[i]); 28 return 0; 29 }
这里有一种容斥做法。
比如说我最后序列里存在第i种东西,且i是最后拿的,那么之前可以拿很多很多的东西,但只能拿另外的n-1种中的k-1种。比如说某k-1种的概率加起来是x,i的概率是a,那么这k-1种的状态对i的答案的贡献就是:a*P(这k-1种的状态)=a*(1+x+x^2+x^3+……)=a/(1-x),等比数列求和,简单。等会确定是这样?这里面包含了选k-2种,k-3种,……,0种的情况!!!那就挑掉k-2的,然后可能在k-2里面挑多了,那k-3再补回来……这个系数可以用组合数计算:(copy from liurunda)