Time Limit: 1 Sec Memory Limit: 128 MB
Description
LMZ有n个不同的基友,他每天晚上要选m个进行[河蟹],而且要求每天晚上的选择都不一样。那么LMZ能够持续多少个这样的夜晚呢?当然,LMZ的一年有10007天,所以他想知道答案mod 10007的值。(1<=m<=n<=200,000,000)
Input
第一行一个整数t,表示有t组数据。(t<=200)
接下来t行每行两个整数n, m,如题意。
Output
T行,每行一个数,为C(n, m) mod 10007的答案。
Sample Input
4
5 1
5 2
7 3
4 2
Sample Output
5
10
35
6
noip要来了,写点模板题找找自信?
一看题简直裸得不能再裸了,然后直接上lucas。
然而没有预处理阶乘的逆元和阶乘,于是愉快的跑了20MS。
惨啊。
然后研究了一下怎么线性求1~n每个数的逆元,inv[i] ≡ -p/i * inv[p % i] (mod p)
然后再乘一次就得到阶乘的逆了
#include <cstdio>
#include <algorithm>
using namespace std;
const int mod = 10007;
int fac[mod + 3], inv[mod + 3];
void init() {
fac[0] = 1;
for(int i = 1; i < mod; i++) fac[i] = (long long)fac[i - 1] * i % mod;
inv[0] = inv[1] = 1;
for(int i = 2; i < mod; i++) inv[i] = (long long)(mod - mod / i) * inv[mod % i] % mod;
for(int i = 1; i < mod; i++) inv[i] = (long long)inv[i] * inv[i - 1] % mod;
}
int C(int a, int b) {
if(b > a) return 0;
if(a < mod) return (long long)fac[a] * inv[b] % mod * inv[a - b] % mod;
return (long long)C(a / mod, b / mod) * C(a % mod, b % mod) % mod;
}
int main() {
init();
int T, a, b;
scanf("%d", &T);
while(T--) {
scanf("%d%d", &a, &b);
printf("%d\n", C(a, b));
}
return 0;
}