线性逆元模板P4071

分析：

线性求逆元：https://blog.csdn.net/qq_34564984/article/details/52292502

code：

#include<cstdio>
using namespace std;
const long long mod=1000000007;
long long ni[1000005],cheng[1000005],dao[1000005],d[1000005];
int main()
{
    d[0]=1;
    d[1]=0;
    d[2]=1;
    for(long long i=3;i<=1000000;i++)
    {
        d[i]=((i-1)*(d[i-1]+d[i-2]))%mod; 
    }//错排递推公式！！！
    ni[1]=1;//1的逆元为1 
    for(long long i=2;i<=1000000;i++)
    {
        ni[i]=(mod-mod/i)*ni[mod%i]%mod;
    }//求出i的逆元(线性求逆元板子
    cheng[0]=1;
    for(long long i=1;i<=1000000;i++)
    {
        cheng[i]=(cheng[i-1]*i)%mod;
    } //求出i的正常阶乘
    dao[0]=1;
    for(long long i=1;i<=1000000;i++)
    {
        dao[i]=(dao[i-1]*ni[i])%mod;
    } 
    long long T;
    scanf("%lld",&T);
    while(T--)
    {
        long long n,m;
        scanf("%lld%lld",&n,&m);
        printf("%lld\n",(cheng[n]%mod*dao[m]%mod*dao[n-m]%mod*d[n-m]%mod)%mod);
    }
    return 0;
}