[Sdoi 2016] 排列计数

题意描述：

求有多少种长度为 n 的序列 A，满足以下条件：

1 ~ n 这 n 个数在序列中各出现了一次

若第 i 个数 A[i] 的值为 i，则称 i 是稳定的。序列恰好有 m 个数是稳定的

满足条件的序列可能很多，序列数对 10^9+7 取模。

题目分析：

错排公式+组合数
答案为 $mar[n-m]*C(n,m)$

题目链接：

Luogu 4071
BZOJ 4517
COGS 2224

Ac code :

#include <cstdio>
#include <iostream>
#define ll long long
#define il inline
const int mod=1e9+7;
const int maxm=1000010;
ll mar[maxm]={1,0,1};
ll mul[maxm];
il void make()
{
    for(int i=3;i<=maxm;i++)
     mar[i]=(1ll*(i-1)*(mar[i-1]+mar[i-2])%mod)%mod;
    mul[0]=1;
    for(int i=1;i<=maxm;i++)
     mul[i]=(mul[i-1]*1ll*i)%mod;
}
il ll fast_pow(ll x,int y)
{
    ll ans=1;
    while(y)
    {
        if(y%2) ans=(ans*x)%mod;
        y/=2;
        x=(x*x)%mod;
    }
    return ans%mod;
}
ll C(int n,int m)
{
    ll x=mul[n]%mod;
    ll y=mul[m]*mul[n-m]%mod;
    ll ret=(x*fast_pow(y,mod-2))%mod;
    return ret;
}
int main()
{
    freopen("menci_permutation.in","r",stdin);
    freopen("menci_permutation.out","w",stdout);
    make();
    int q;
    scanf("%d",&q);
    while(q--)
    {
        int n,m;
        scanf("%d%d",&n,&m);
        printf("%lld\n",(mar[n-m]*C(n,m))%mod);
    }
    return 0;
}