bzoj2982: combination 组合数 卢卡斯定理

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;
}