ZOJ 3557 How Many Sets II (组合计数，Lucas定理，费马小定理)

题意： 从1~n里选出m个数，要求这m个数不相邻，问有几种求法，对p取模。n<=1e9,m<=1e4,p<=1e9。

解法： 对于一个数论盲来说，首先想到的是设状态f(n,m)为n个数选m个不相邻的方法数，然后通过讨论第一个数的选和不选，可以容易得到方程：

f(n,m) = f(n-1,m) + f(n-2,m-1)，这个方程对于巨大的n无力，于是变成这样：

f(n,m) = sigma(f(i,m-1)),1<=i<=n-2，然后这个方程和上面那个其实具有相同的性质，即：对解题没有帮助- -。

然后换个思路，n个数，选出来m个，等价于不选n-m个，然后由于m个数不能相邻，所以，m个数只能在n-m+1个空隙里选，到这里了，高中数学：ans = C(n-m+1,m)，轻巧的组合数，比递推式上档次n+1倍。。

因为本题中n比m大5个数量级，因此尝试化简约掉n!：

C(n,m) = n! / ( m! * (n-m)! )

             = n * (n-1) * (n-2) * ... * (n-m+1) / m!

m数量级1e4，可递推解决。

这样没完，因为模数p不是一个确定的巨大素数(like 1e9+7)，而可能只是一个小素数，那么在求a对模p的逆的过程中可能发生(a,p) = p （即 a = k*p）的情况，使得a^-1不能正确求出，这里引入Lucas定理：

C(n,m) % p = Lucas(n,m,p) = C(n%p,m%p,p) * Lucas(n/p,m/p,p)

保证了m与p的互素，逆元可以正确求出。

btw，继昨天exgcd之后，又学了一种求乘法逆元的新姿势：

a^-1 = a^(p-2) % p

简单推一下：

前提：                a 与p互素

费马小定理：     a^(p-1) = 1 % p

-->                        a * a^(p-2) = 1 % p

而乘法逆元：     a * a^-1 = 1 % p

可知：                a^-1 = a^(p-2) % p

/* Created Time: Wednesday, November 06, 2013 PM10:51:00 CST */
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
typedef long long lld;
lld Pow(lld x,lld k,lld p) {
    lld ret = 1;
    while (k) {
        if(k&1) ret = ret*x%p;
        x = x*x%p;
        k >>= 1;
    }
    return ret;
}
lld exgcd(lld a,lld b,lld &x,lld &y) {
    if (!b) {
        x = 1,y = 0;
        return a;
    } else {
        lld c = exgcd(b,a%b,y,x);
        y -= a/b*x;
        return c;
    }
}
lld inv(lld m,lld p) {
    lld x,y;
    exgcd(m,p,x,y);
    x = (x%p+p)%p;
    return x;
}
lld C(lld n,lld m,lld p) {
    lld ret = 1;
    for (lld i = 1; i <= m; i ++) {
        ret = ret*(n-i+1)%p*inv(i,p)%p;
        // inv(i,p) == Pow(i,p-2,p) % p
    }
    return ret;
}
lld Lucas(lld n,lld m,lld p) {
    if (n<m) return 0;
    if (m==0) return 1;
    // 用Lucas定理来保证a与p互质关系
    return C(n%p,m%p,p)*Lucas(n/p,m/p,p)%p;
}
int main() {
    lld n,m,p;
    while (~scanf("%lld%lld%lld",&n,&m,&p)) {
        printf("%lld\n",Lucas(n-m+1,m,p));
    }
    return 0;
}