记录编号 |
376763 |
评测结果 |
AAAAAAAAAAAAAAAAAAAA |
题目名称 |
[SDOI 2016] 排列计数 |
最终得分 |
100 |
用户昵称 |
sxysxy |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
10.508 s |
提交时间 |
2017-02-28 08:02:56 |
内存使用 |
12.44 MiB |
显示代码纯文本
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cstdarg>
#include <algorithm>
#include <string>
#define MAXN 1000001
typedef long long LL;
const LL MOD = 1e9+7;
using namespace std;
LL pow_mod(LL a, LL x){
LL ans = 1;
for(a %= MOD; x; x>>=1, a=a*a%MOD)
if(x&1)ans = ans*a%MOD;
return ans;
}
inline LL inverse(LL x){
return pow_mod(x, MOD-2);
}
LL fact[MAXN];
LL D[MAXN];
int main(){
freopen("menci_permutation.in", "r", stdin);
freopen("menci_permutation.out", "w", stdout);
fact[0] = 1;
for(int i = 1; i <= MAXN; i++)
fact[i] = fact[i-1]*i%MOD;
D[1] = 0;
D[2] = 1;
for(int i = 3; i <= MAXN; i++)
D[i] = (i-1)*(D[i-1]+D[i-2])%MOD;
int T; scanf("%d", &T);
while(T--){
int n, m; scanf("%d %d", &n, &m);
if(n <= m){
if(n < m)puts("0");
else puts("1");
continue;
}
LL ans = fact[n];
LL tmp = fact[m]*fact[n-m];
ans = ans*inverse(tmp)%MOD;
ans = ans*D[n-m]%MOD;
printf("%lld\n", ans);
}
return 0;
}