hdu4373 Mysterious For（剩余定理+Lucus定理）

题意：给你m重循环每一层循环从上一层或者0开始，到n-1结束，求循环的总次数。
解法：$a1<=a2$<=>$a1<a2+1$
因此由插空法可得计算公式为$C(n+d_{k+1}-d_k-1,d_{k+1}-d_k)$的连乘形式
但是问题的关键在于结果要取模，且取模对象不是素数，因此需要对mod进行分解，分别Lucus定理求解后再用中国剩余定理合并求解。
用Lucus定理的原因是C(n,m)中的n和m大于p，直接求逆元为0，不能使用，因此需要减小n和m的取值范围。
用中国剩余定理是合并模的方法
注意：Lucus定理求解时，需要判断a和b的大小关系！

#include<iostream>
#include<cstdio>
#include<string.h>
#include<queue>
#include<cmath>
#include<vector>
#include<algorithm>
#include<map>
#include<stdlib.h>
#include<set>
#define ll long long
using namespace std;

const ll maxn = 3761599+10;
const ll mod = 364875103;
ll f[2][maxn],ans;
ll T,n,k,m,d[30];
ll mo[2];

ll quick(ll a,ll b,ll m){
    ll ans=1;
    while(b){
        if(b&1) ans=(ans*a)%m;
        a=(a*a)%m; b>>=1;
    }
    return ans;
}

void init(){
    f[0][0]=f[0][1]=f[1][0]=f[1][1]=1;
    mo[0]=97,mo[1]=3761599;
    for(ll i=2;i<mo[0];i++){
        f[0][i]=(f[0][i-1]*i)%mo[0];
    }
    for(ll i=2;i<mo[1];i++){
        f[1][i]=(f[1][i-1]*i)%mo[1];
    }
}

ll C2(ll id,ll a,ll b,ll m){
    if(a<b) return 0;//!!!
    return f[id][a]*quick(f[id][b]*f[id][a-b]%m,m-2,m)%m;
}

ll Lucus(ll a,ll b,ll m,ll id){
    ll res=1;
    while(a&&b&&res){
        res=(res*C2(id,a%m,b%m,m))%m;
        a/=m,b/=m;
    }
    return res;
}

ll china(ll n,ll a[],ll m[]){
    ll M=1,ans=0;
    for(ll i=0;i<n;i++) M*=m[i];
    for(int i=0;i<n;i++){
        ll p=quick(M/m[i],m[i]-2,m[i]);
        ans=(ans+p*M/m[i]%M*a[i]%M)%M;
    }
    return ans;
}

ll C(ll a,ll b){
    ll ans[2];
    ans[0]=Lucus(a,b,mo[0],0);
    ans[1]=Lucus(a,b,mo[1],1);
    return china(2,ans,mo);
}

int main(){
    init();
    //freopen("a.txt","r",stdin);
    scanf("%I64d",&T);
    for(ll cas=1;cas<=T;cas++){
        scanf("%I64d%I64d%I64d",&n,&m,&k);
        for(ll i=0;i<k;i++) scanf("%I64d",&d[i]);
        d[k]=m,ans=1;
        for(ll i=0;i<k;i++) ans=ans*C(n+d[i+1]-d[i]-1,d[i+1]-d[i])%mod;
        printf("Case #%I64d: %I64d\n",cas,ans);
    }
    return 0;
}