组合数取模模板

组合数取模模板（Pascal（到时候改C++））

function fact(x,y:longint):note;
var
    k:note;
begin
    if x<a[y] then
    begin
        fact.e:=sum[x,y];fact.f:=0;exit;
    end;
    fact.f:=x div a[y];k:=fact(x div a[y],y);
    fact.f:=fact.f+k.f;
    fact.e:=k.e*sum[x mod g[y],y] mod g[y]*mi(sum[g[y]-1,y],x div g[y],g[y]) mod g[y];
end;
function c(x,y:longint):longint;
var
    i,j:longint;
    k1,k2,k3:note;
begin
    if x>y then exit(0);
    for i:=1 to a[0] do
    begin
        k1:=fact(y,i);k2:=fact(x,i);
        k3:=fact(y-x,i);
        an[i]:=mi(a[i],k1.f-k2.f-k3.f,g[i])*k1.e 
        mod g[i]*ni(k2.e,g[i]) mod g[i]*ni(k3.e,g[i]) mod g[i];
    end;c:=0;
    for i:=1 to a[0] do
    begin
        j:=p div g[i];
        c:=(c+j*ni(j,g[i]) mod p*an[i] mod p) mod p;
    end;
end;