题目
思路
错排:
给定n元素集合X,它的每一个元素都有一个特定的位置,而现在要求求出集合X的排列中没有一个元素在它指定位置上的排列的数目.
递推公式:
D
1
=
0
D
2
=
1
D_1=0 \ D_2=1
D1=0 D2=1
D
i
=
(
i
−
1
)
∗
(
D
i
−
1
−
D
i
−
2
)
D_i=(i-1)*(D_{i-1}-D_{i-2})
Di=(i−1)∗(Di−1−Di−2)
这道题要保证 m m m个数等于它的下标,即 n − m n-m n−m个数不等于下标 答案为 C m n ∗ D ( n − m ) C_m^n*D(n-m) Cmn∗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;
}