[错排] luogu 4071

题目

思路

错排：

给定n元素集合X,它的每一个元素都有一个特定的位置,而现在要求求出集合X的排列中没有一个元素在它指定位置上的排列的数目.

递推公式：
$D_1=0D_2=1$
$D_i=(i-1)\ast (D_{i-1}-D_{i-2})$

这道题要保证$m$个数等于它的下标，即$n-m$个数不等于下标 答案为$C_m^n\ast D(n-m)$

代码

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<ctime>
#include<iostream>
#include<string>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<vector>
#include<iomanip>
#include<list>
#include<bitset>
#include<sstream>
#include<fstream>
#include<complex>
#include<algorithm>
#if __cplusplus >= 201103L
#include <unordered_map>
#include <unordered_set>
#endif
#define int long long
using namespace std;
const int INF = 0x3f3f3f3f;
int inv[1000010],f[1000010],d[1000010];
const int mod = 1e9+7;
int quick_pow(int a,int b){
	int res=1;
	while(b){
		if(b&1){
			res=res*a%mod;
		}
		b>>=1;
		a=a*a%mod;
	}
	return res;
}
void pre(){
	d[0]=1,d[2]=1,d[3]=2;
	for(int i=4;i<=1000010;i++){
		d[i]=(i-1)*(d[i-1]+d[i-2])%mod;
	}
	f[0]=1;
	for(int i=1;i<=1000010;i++){
		f[i]=f[i-1]*i%mod;
	}
	for(int i=1;i<=1000010;i++){
		inv[i]=quick_pow(f[i],mod-2)%mod;
	}
}
signed main(){
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
	int t;
	cin>>t;
	pre();
	while(t--){
		int n,m;
		cin>>n>>m;
		if(n-m==1) cout<<"0"<<endl;
		else if(n==m) cout<<1<<endl;
		else if(m==0) cout<<d[n]<<endl;
		else cout<<f[n]*inv[m]%mod*inv[n-m]%mod*d[n-m]%mod<<endl;
	}
    return 0;
}