记录编号 376763 评测结果 AAAAAAAAAAAAAAAAAAAA
题目名称 [SDOI 2016] 排列计数 最终得分 100
用户昵称 Gravatarsxysxy 是否通过 通过
代码语言 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;
}