排列计数(数论）

排列计数

Description
求有多少种长度为n的序列A，满足以下条件：
1 n这n个数在序列中各出现了一次
若第i个数A[i]的值为i，则称i是稳定的。序列恰好有m个数是稳定的满足条件的序列可能很多，序列数对109+7取模。

Input
第一行一个数 T，表示有 T 组数据。
接下来 T 行，每行两个整数 n,m。
T=500000，n≤1000000，m≤1000000

output
输出T行，每行一个数，表示求出的序列数

Samples
Input 复制
5
1 0
1 1
5 2
100 50
10000 5000
Output
0
1
20
578028887
60695423

题意：一共n个数，m个数是符合第i个数是在i上的，其他都是乱序**（其他数都不在自己位置上）**

思路：解题分成两步：1.即n个数中选m个数 2.求n-m乱序的可能排序的数量，用f[n-m]表示，那么答案ans=c(n,m)f[n-m];

第一步：n个数中选m个数
c(n,m)=(n!)/(m!*(n-m)!)，在这里n太大了，所以我们用另外一种方式求!!!(俺以前都不知道！！！我太菜了）
c(n,m)=n!（m!的逆元）（（n-m）！的逆元）
求逆元也有一个很牛的式子！！
n的逆元=qpow(n,mod-2) mod必须是素数！！
另外一个式子
n的逆元=（1的逆）（2的逆）……（n的逆）

//这里我们用infact[]表示逆元
//两个式子的写法
//first
void jj(){
   infact[0]=1;
   for(int i=1;i<=n;i++)
   infact[i]=infact[i-1]*qpow(i,mod-2);
}
//second
void jj2(){
 infact[n]=qpow(n,mod-2);
 for(int i=n-1;i>=1;i--)
 infact[i]=infact[i+1]*(i+1);
}

第二步：是一个dp
转移公式：f[i]=(f[i-1]+f[i-2])*(i-1);
解释：已经排完了i-1个数，现在又添加一个新的，他有两种添加方式，一种是和排完序的一个互换其他不变，即f[i-2]，另一种是排完序的一个到他的位置，但添加的数未进队伍，即f[i-1]，一共（i-1）种选择。

下面是完整代码：

int n,m;
ll f[maxn],h[maxn]={1,1},a[maxn];
ll qpow(ll a,ll b){
	ll res=1;
	while(b)
	{
		if(b&1)
		res=res*a%mod;
		b>>=1;
		
		a=(a*a)%mod;
	}
	return res;
}
void intal(){
	int maxx=1000000;
	for(int i=2;i<=maxx;i++) h[i]=h[i-1]*i,h[i]%=mod;
	a[maxx]=qpow(h[maxx],mod-2);
	for(int i=maxx-1;i>=1;i--) a[i]=a[i+1]*(i+1)%mod;
	a[0]=1;f[0]=1,f[1]=0,f[2]=1;
	for(int i=3;i<=maxx;i++) f[i]=(i-1)*(f[i-1]+f[i-2]),f[i]%=mod; 
}
ll c(int x,int y){
	ll ans=h[n]*a[m]%mod*a[n-m]%mod;
	return ans;
}

int main(){
	int T = read();
	intal();
	while(T--){
		n = read(),m = read();
		ll ans=c(n,m)*f[n-m];
		ans%=mod;
		cout<<ans<<"\n";
	}
}