AtCoder Regular Contest 096 E

题意：

有n个数，$2^n$个集合，从中选择一些集合，问有多少种方案使每个数至少出现两次。

题解：

直接计数不好计，考虑容斥，枚举不合法。
$$ans=\sum _{i=0}^n{2}^{2^{n-i}}{C}_n^i\sum _j^i{C}_i^j\sum _k^{j}S(j,k)(2^{n-i}{)}^k$$
$S(n,m)$是第二类斯特林数，$j$为枚举有多少个数只出现一次，再将他们分到若干集合中。
考虑优化：
$$ans=\sum _{i=0}^n{2}^{2^{n-i}}{C}_n^i\sum _k^i(2^{n-i}{)}^k\sum _{j=k}^i{C}_i^{j}S(j,k)$$

我们有：
$$\sum _{j=k}^i{C}_i^{j}S(j,k)=S(i+1,k+1)$$

这个很好理解，相当于枚举有多少不和$i+1$一个集合。
$$ans=\sum _{i=0}^n{2}^{2^{n-i}}{C}_n^i\sum _k^i(2^{n-i}{)}^{k}S(i+1k+1)$$

$O(n^{2}logm)$

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
LL C[3010][3010],S[3010][3010];
LL n,p;
void pre()
{
	for(LL i=0;i<=3001;i++)
	{
		C[i][0]=1;
		for(LL j=1;j<=i;j++) C[i][j]=(C[i-1][j]+C[i-1][j-1])%p;
	}
	S[0][0]=1;
	for(LL i=1;i<=3001;i++)
	{
		S[i][1]=1;
		for(LL j=2;j<=i;j++) S[i][j]=(S[i-1][j]*j%p+S[i-1][j-1])%p;
	}
}
LL pow(LL a,LL b,LL mod)
{
	LL ans=1;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		a=a*a%mod;b>>=1;
	}
	return ans;
}
int main()
{
	LL ans=0;scanf("%lld %lld",&n,&p);
	pre();
	for(LL i=0;i<=n;i++)
	{
		LL s=C[n][i]*pow(2,pow(2,n-i,p-1),p)%p;
		LL tot=0;
		for(LL j=0;j<=i;j++) tot=(tot+pow(2,(n-i)*j%(p-1),p)*S[i+1][j+1]%p)%p;
		ans=(ans+(i&1?-1:1)*s%p*tot%p)%p;
	}
	printf("%lld",(ans+p)%p);
}