【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=0,sign=1; char ch=getchar();
    while(ch<'0'||ch>'9') {if (ch=='-') sign=-1; ch=getchar();}
    while(ch>='0'&&ch<='9') {v=v*10+ch-'0'; ch=getchar();}
    return v*=sign;
}
/*******************tamplate********************/
const int N=100086;
typedef long long LL;
LL n,m,P,fac[N],inv[N],ans;
LL pow(LL a,LL b){
    LL r=1,base=a;
    for(;b;b>>=1,base=base*base%P)
        if (b&1) r=r*base%P;
    return r;
}
inline LL c(LL n,LL m){
    if (n<m) return 0;
    if (n<P && m<P){
        LL a=1,b=1;
        for(;m;m--,n--){
            a=a*n%P;
            b=b*m%P;
        }
        return a*pow(b,P-2)%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
",c(n+m,m));
    }
    return 0;
}