排列 [HNOI2011]卡农

问题 A: [HNOI2011]卡农
时间限制: 1 Sec 内存限制: 128 MB
题目描述

共有2^n-1种片段组合方式。设数字i出现次数ai,选i个片段是合法总方案数为fi。
这里我们先按不同排列方式为不同情况。
那么对于A（2^n-1,i-1）就是i个位置上有i-1个位置上放上了片段，有两种情况，
1，这i-1个片段已满足要求，不合法。
2，这i-1个片段不满足要求，那么有唯一一个片段使之合法。
但是考虑放i-2个片段时是满足的，加了任意一个其他片段后是i-1个，之后又加了一个相同的片段使之满足，那他就重复了，不合法。所以要把这两种情况去掉。
第一种减去fi-1，第二种减fi-2*(i-1)*(2^n-1-(i-2))。
那我解释一下，对于每种放i-2个时的满足情况，可以放剩余2^n-1-(i-2)个片段中的任一个，并且有i-1个位置可以放。那现在再结合上面的看。
最后除以m!（因为不同排列算一个）

#include<iostream>#include<cstdlib>#include<cstdio>#include<cstring>#include<algorithm>#define ll long long#define mod 100000007#define N 1000005using namespace std;ll n,m,A[N],f[N],xp=1,ans=1;inline ll cheng(ll x,ll n){    ll s=1;x%=mod;    while(n)    {        if(n&1)s=s*x%mod;        x=x*x%mod;        n/=2;    }    return s;}int main(){   // freopen("canon.in","r",stdin);//  freopen("canon.out","w",stdout);    cin>>n>>m;A[0]=1;f[0]=1;    xp=cheng(2ll,n);xp--;if(xp<0)xp+=mod;    for(int i=1;i<=m;i++)A[i]=A[i-1]*(xp-i+1)%mod;    for(int i=3;i<=m;i++)f[i]=(A[i-1]-f[i-1])%mod-(f[i-2]*(xp-i+2)%mod*(i-1))%mod,f[i]%=mod;    for(int i=2;i<=m;i++)ans=(ans*i)%mod;    ans=cheng(ans,mod-2);//求逆元    ans=(ans*f[m]%mod+mod)%mod;    cout<<ans;}