传送门:bzoj2339
题解
神仙 D P DP DP。
本质不同有点麻烦,考虑枚举方案之后除
m
!
m!
m!即可(???这一步就没想到)
设 f [ i ] f[i] f[i]表示长度为 i i i的满足条件的音乐个数(排列方式不同算不同)。考虑用容斥的方法计数。
假设已知前 i − 1 i-1 i−1个片段,第 i i i个片段就是唯一确定的,共 A 2 n − 1 i − 1 A_{2^n-1}^{i-1} A2n−1i−1种方案。
第 i i i个片段为空集的方案数为 f [ i − 1 ] f[i-1] f[i−1],与前 i − 1 i-1 i−1个片段中某个重复的概率为 f [ i − 2 ] ∗ ( 2 n − 1 − ( i − 2 ) ) ∗ ( i − 1 ) f[i-2]*(2^n-1-(i-2))*(i-1) f[i−2]∗(2n−1−(i−2))∗(i−1)(删去与它相同的后方案数为 f [ i − 2 ] f[i-2] f[i−2], i i i的方案数为 2 n − 1 − ( i − 2 ) 2^n-1-(i-2) 2n−1−(i−2),任选前面一个片段与它重复( × ( i − 1 ) \times (i-1) ×(i−1))。
所以转移方程为:
f [ i ] = A 2 n − 1 i − 1 − f [ i − 1 ] − f [ i − 2 ] ∗ ( 2 n − 1 − ( i − 2 ) ) ∗ ( i − 1 ) f[i]=A_{2^n-1}^{i-1}-f[i-1]-f[i-2]*(2^n-1-(i-2))*(i-1) f[i]=A2n−1i−1−f[i−1]−f[i−2]∗(2n−1−(i−2))∗(i−1)
初始化 f [ 0 ] = 1 , f [ 1 ] = 0 f[0]=1,f[1]=0 f[0]=1,f[1]=0。
代码
#include<bits/stdc++.h>
using namespace std;
const int N=1e6+100,mod=1e8+7;
typedef long long ll;
int n,m,lim,f[N],a[N],frc;
inline int dc(int x,int y){x-=y;return x<0?x+mod:x;}
inline int ad(int x,int y){x+=y;return x>=mod?x-mod:x;}
inline int fp(int x,int y)
{
int re=1;
for(;y;y>>=1,x=(ll)x*x%mod)
if(y&1) re=(ll)re*x%mod;
return re;
}
int main(){
int i,j;
scanf("%d%d",&n,&m);
a[1]=lim=dc(fp(2,n),1);frc=1;
for(i=1;i<m;++i) a[i+1]=(ll)a[i]*dc(lim,i)%mod;
for(i=2;i<=m;++i) frc=(ll)frc*i%mod;
f[1]=0;f[0]=1;
for(i=2;i<=m;++i)
f[i]=dc(a[i-1],ad(f[i-1],(ll)f[i-2]*(i-1)%mod*(ll)dc(lim,i-2)%mod));
printf("%d",(ll)f[m]*fp(frc,mod-2)%mod);
return 0;
}