题意: 从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;
}