O(n) 곱셈 역원

공식은 다음과 같습니다.
inv[i]=(MOD-MOD/i)*inv[MOD%i]%MOD;
전제는 모수이고 MOD는 질수여야 돼요.
코드로 작성된 것은 다음과 같습니다.

inv[1]=1;
for(int i=2;i<MOD;i++)
    inv[i]=(MOD-MOD/i)*inv[MOD%i]%MOD;

시간의 복잡도는 O(n)이다. 즉, 한 번 훑어보면 [1,n] 범위 내의 모든 수의 역원을 구할 수 있다.(n

#include
#include
#define LL long long
#define MOD 
using namespace std;
LL fac[200000],inv[200000],p;
void op(){
    fac[0]=1,inv[1]=1,inv[0]=1;
    for(int i=1;i<=p;i++)fac[i]=(fac[i-1]*i)%p;
    for(int i=2;i
//      i!   
    //           ，          i!   
}
LL C(LL n,LL m){  
    if(m>n)return 0;  
    if(m
return fac[n]*inv[n-m]%p*inv[m]%p;
    return C(n%p,m%p)*C(n/p,m/p)%p;  
}
int main(){
    LL n,m;
    scanf("%lld%lld%lld",&m,&n,&p);
    op();//  
    printf("%lld",(C(m+n,m)+p)%p);
    return 0;
}