lucas定理学习小记

lucas是数论定理，运用到OI上时，就会发现非常便利，它可以在极快的时间内算出C（n,m）mod p，p为质数。时间复杂度为O（log（以p为底）(n)*p)）。

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

我们令n=sp+q , m=tp+r .（q ，r ≤p）
那么：

（在编程时你只要继续对 调用Lucas定理即可。
代码可以递归的去完成这个过程，其中递归终点为t = 0 ；
时间O（log（以p为底）(n)*p)）
模板：（SDOI2010古代猪文部分代码）

function lucas(a,b,x:int64):int64;
begin
if b=0 then exit(1);
exit( c(a mod t[x],b mod t[x],x)*lucas(a div t[x],b div t[x],x)mod t[x] );
end;

证明自行百度，，，