BZOJ 4517 组合数+错排
思路:
预处理错排
然后C(n,m)*s[n-m-1]就是答案了
特判n-m-1<0
//By SiriusRen #include <cstdio> using namespace std; #define int long long const int mod=1000000007,N=1000050; int cases,n,m,fac[N],s[N]; int pow(int x,int y){ int res=1; while(y){ if(y&1)res=res*x%mod; x=x*x%mod,y>>=1; }return res; } int C(int x,int y){return fac[x]*pow(fac[x-y],mod-2)%mod*pow(fac[y],mod-2)%mod;} signed main(){ scanf("%lld",&cases);fac[0]=s[1]=1; for(int i=1;i<N;i++)fac[i]=fac[i-1]*i%mod; for(int i=2;i<N;i++)s[i]=(s[i-1]+s[i-2])*i%mod; while(cases--){ scanf("%lld%lld",&n,&m); printf("%lld\n",C(n,m)*s[n-m-1<0?1:n-m-1]%mod); } }