【HDOJ】【3037】Saving Beans

排列组合

　　啊……这题是要求c(n-1,0)+c(n,1)+c(n+1,2)+......+c(n+m-1,m)

　　这个玩意……其实就等于c(n+m,m)

　　好吧然后就是模P……Lucas大法好= =

　　我SB地去预处理<P的所有fac和inv了……果断TLE

　　事实上Lucas时对于<P的部分直接暴力算就好了

 //HDOJ 3037
 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 using namespace std;

 int getint(){
     int v=,sign=; char ch=getchar();
     while(ch<''||ch>'') {if (ch=='-') sign=-; ch=getchar();}
     while(ch>=''&&ch<='') {v=v*+ch-''; ch=getchar();}
     return v*=sign;
 }
 /*******************tamplate********************/
 const int N=;
 typedef long long LL;
 LL n,m,P,fac[N],inv[N],ans;
 LL pow(LL a,LL b){
     LL r=,base=a;
     for(;b;b>>=,base=base*base%P)
         if (b&) r=r*base%P;
     return r;
 }
 inline LL c(LL n,LL m){
     if (n<m) return ;
     if (n<P && m<P){
         LL a=,b=;
         for(;m;m--,n--){
             a=a*n%P;
             b=b*m%P;
         }
         return a*pow(b,P-)%P;
     }
     return c(n%P,m%P)*c(n/P,m/P)%P;
 }
 int main(){
 //    freopen("input.txt","r",stdin);
     int T=getint();
     while(T--){
         n=getint(); m=getint(); P=getint();
         printf("%lld\n",c(n+m,m));
     }
     return ;
 }