BZOJ 4402: Claris的剑|组合数学

由于存在集合等价的情况，我们就只考虑字典学最小的方案
假设最大值是$max$,那么构造出的序列有两种情况
$1,2,1,2,1,2....3,2,3,2,3....max-2,max-1,max-2,max-1....max$
$1,2,1,2,1,2....3,2,3,2,3....max-2,max-1,max-2,max-1....max,max-1$
$1..max$的所有数都至少出现一次然后剩余的空间必须成对的插入.
那么我们就先构造出了$1,2,3..max$的数列然后再插入其他的数
比如我们可以在$1$和$2$之间插入任意个数对$(2,1)$
更一般的来说可以在$i$和$i+1$之间插入任意个数对$(i+1,i)$
假设剩余的空间有$n$个，可以选择的数对有$m$种，那么方案数就是

$$C(\lfloor\frac{n}{2}\rfloor+m,m)$$
第一次做组合数学的题……被惨虐QaQ

#include<bits/stdc++.h>
#define N 2000200
#define R 1000000007
using namespace std;
long long fac[N],inv[N],ans;
int n,m;
long long cal(int n,int m)
{
    if(n<0)return 0;
    n>>=1;if(n==0)return 1;
    return fac[n+m]*inv[n]%R*inv[m]%R;
}
int main()
{
    cin>>n>>m;
    fac[0]=fac[1]=inv[0]=inv[1]=1;
    for(int i=2;i<=n;i++)fac[i]=fac[i-1]*i%R;
    for(int i=2;i<=n;i++)inv[i]=inv[R%i]*(R-R/i)%R;
    for(int i=2;i<=n;i++)inv[i]=inv[i]*inv[i-1]%R;
    if(n&&m)ans=1;
    for(int i=2;i<=m;i++)
        (ans+=cal(n-i,i-1))%=R,
        (ans+=cal(n-i-1,i-1))%=R;
    cout<<ans;
    return 0;
}