这题主要就是用错排公式 + 组合数学:
错排公式: D[n] = (n-1) * ( D[n-1] + D[n-2] )
求组合数: 我们注意到n、m都比较大,所以我们在求C(n, m) = n!/((n-m)!*m!)时候,直接求的话会出现精度问题,所以我们可以通过求逆元来避免这个问题,同时我们又注意到MOD=1000000007,这是个素数,所以我们可以根据费马小定理直接算出逆元。详情见代码
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
const LL MOD = 1000000007;
const int maxn = 10005;
const int maxc = 105;
int T, n, k;
LL f[maxn], D[maxn], inv[maxn];
LL pow_w(LL a, LL b, LL c){
LL ans = 1;
a = a%c;
while(b){
if(b&1) ans = ans*a%c;
a = a*a%c;
b >>= 1;
}
return ans;
}
void deal(){
f[0] = f[1] = 1;
for(int i = 2; i <= 10000; ++i)
f[i] = f[i-1]*i%MOD;
for(int i = 0; i <= 10000; ++i)
inv[i] = pow_w(f[i], MOD-2, MOD);
D[0] = 1; D[1] = 0;
for(int i = 2; i <= 10000; ++i)
D[i] = ((i-1)*(D[i-1] + D[i-2]))%MOD;
}
LL c(LL n, LL m){
//利用逆元将:a/b mod c <=转化为=> a*x mod c ; x为b的逆元。
return (f[n]*inv[n-m]%MOD)*inv[m]%MOD;
}
int main(){
deal();
scanf("%d", &T);
while(T--){
scanf("%d%d", &n, &k);
LL ans = 0;
for(k ; k <= n; ++k){
ans = (ans+c(n, k)*D[n-k])%MOD;
}
printf("%lld\n", ans);
}
return 0;
}