csp2021-09-4 收集卡牌

题目链接：收集卡牌
因为$n⩽16$，卡牌收集情况可以用$2^{16}=65536$个状态来表示；
又因为$k⩽$ 5，在最差情况下只抽到一种卡，那么最多需要$5\times 15+1=76$次即可收集完成。
所以可以用$dp[65536][77]$来表示所有状态，$dp[state][t]$表示当前收集情况为$state$，且已抽取$t$次的概率。
那么初始状态就是$dp[0][0]=1$。
如果当前状态$[state][t]$已经集齐卡片，对答案就有 当前状态概率$\times t$ 的贡献，且不会有后继状态；否则，当前状态需要继续抽取：若抽取到第$i$种卡片，则对$[state|(1<<i)][t+1]$状态贡献了 当前状态概率$\times p_i$ 的概率。
从初始状态递推即可。

没有参加考试，在官网补的题，严重怀疑官网没有写spj：
题面写的是：

如果你给出的答案与标准答案的绝对误差不超过 $10^{-4}$，则视为正确。

刚开始我用 %.12lf 输出，一直是0分；改成 %.10lf 交了一发，得90分。然后用队友的100分代码对拍了一下，发现可能存在1e10的误差。于是把ans+=1e-10，得10分。于是猜测有一个测试点我的输出比标程少了1e-10。
然后就二分找这个测试点的ans所在区间。具体方法是：把L<ans<R 全部 加上1e-1 交一发，再改成1e-10交一发，如果得分多了10分，说明目标区间在内。于是我就交了N发 愣是从[1,100]缩小到[25.5,25.6]

先上我的概率递推
100分代码：

#include<bits/stdc++.h>
using namespace std;
double p[20];
double dp[1<<16][77]= {0};
int n,k;
int main()
{
	scanf("%d%d",&n,&k);
	for(int i=0; i<n; i++)
		scanf("%lf",p+i);
	double ans=0;
	int maxt=n*5-4;
	int maxstate=1<<n;
	dp[0][0]=1;
	for(int t=0; t<=maxt; t++) {
		for(int state=0; state<maxstate; state++) {
			int cnt=0;
			for(int i=0; i<n; i++)
				if((state>>i)&1)cnt++;
			if(cnt+(t-cnt)/k>=n) {	//t-cnt即为硬币数
				ans+=dp[state][t]*t;
				continue;
			}
			for(int i=0; i<n; i++)
				dp[state|(1<<i)][t+1]+=dp[state][t]*p[i];
		}
	}
	if(25.5<ans&&ans<25.6)		//无奈的特判，删掉特判就是90分
		ans+=1e-10;
	printf("%.10lf\n",ans);
	return 0;
}

再补个队友的记忆化搜索 (这应该才是标准答案) (队友博客)
100分代码：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
 
using namespace std;

typedef long long LL;
typedef pair<int,int> PII;
const int N=2e6+10,mod=998244353;
bool dis[1<<17][100];
double g[17],dp[1<<17][100];
int n,k;

double dfs(int t,int cnt,int time,int sum,double l)
{
	if(dis[t][time])return dp[t][time];
	dis[t][time]=1;
	if(sum==n||(cnt>=(n-sum)*k))return dp[t][time]=time;
	
	for(int i=1;i<=n;i++)
		if(!((t>>i)&1))
			dp[t][time]+=g[i]*dfs(t|(1<<i),cnt,time+1,sum+1,l+g[i]);
	if(t)dp[t][time]+=l*dfs(t,cnt+1,time+1,sum,l);
	return dp[t][time];
}

int main()
{
	scanf("%d%d",&n,&k);
	for(int i=1;i<=n;i++)scanf("%lf",g+i); 
	printf("%.10lf\n",dfs(0,0,0,0,0));
	return 0;
}