求大组合数mod p，（p不一定为质数）

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define N 2000005
ll p;
ll n,m;
vector<ll>pri;
bool prime[N];

void seive()
{
    for(int i=;i<=N;i++) prime[i]=;
    for(int i=;i<N;i++)
    {
        if(prime[i])
        {
            pri.push_back(i);
            for(int j=i;j<N;j+=i) prime[j]=;
        }
    }
}
ll P(ll a,ll b)
{
    ll res=;
    while(b)
    {
        if(b&) res=res*a%p;
        a=a*a%p;
        b>>=;
    }
    return res;
}
ll work(ll x,ll k)
{
    ll res=;
    while(x)
    {
        res+=x/k;
        x=x/k;
    }
    return res;
}
ll c(ll a,ll b)
{
    ll ans=;
    for(int i=;pri[i]<=a;i++)
    {
        ll aa=work(a,pri[i]);
        ll bb=work(b,pri[i]);
        ll cc=work(a-b,pri[i]);
        ans*=P(pri[i],aa-(bb+cc));
        ans%=p;
    }
    return (ans+p)%p;
}

int main()
{
    seive();
    while(~scanf("%lld%lld%lld",&n,&m,&p))
    {
        ll res1=c(n*,n);
        ll res2=;
        if(m<=n) res2=c(n*,n+m);
        printf("%lld\n",((res1-res2)+p)%p);
    }
    return ;
}

https://ac.nowcoder.com/acm/contest/317/H