hdu 4373 数论

额，首先要推个公式。。就是一个T1循环后面跟着m个TYPE2循环的复杂度（显示语句次数）

这公式刚开始也不会推了，后来随便整，发现了个规律

ANS = C(m,0)+C(m,1)+...+C(m+n-1,n-1)  

组合数学当中有个公式，C(M,0)+C(M,1)+....+C(m+n-1,n-1) = C(m+n-1,n) 

接下来就是求C(M+N-1,N)%MOD

出题人太坑爹了，模数竟然能够分解成97和3761599。。。。实在缺德

这样的话就是求组合数模两个素数的乘积，我们令这个结果为x，组合数为C

设 C% p1 =m1 ,C%p2 =m2

则 x = m1 (mod p1),x= m2(mod p2)

用中国剩余定理求x即可，求模97的余数用LUCAS定理，求大素数的余数直接搞逆元

LUCAS定理一直搞错了，是这样 L(N,M) = C(N%P,M%P) * L(N/P,M/P) ......附代码：

#include <iostream>
#include <cstdio>

#define ll long long
#define ss(a) scanf("%d",&a)
#define N 1100005
#define p1 97
#define p2 3761599
#define mod (p1*p2)
 
using namespace std;

int a[20],c[100][100];
ll y1,y2,mul[N];

void extend(ll a0,ll b0,ll &x,ll &y)  
{  
    if (!b0)  
    {  
        x=1;y=0;  
    }  
    else  
    {  
        extend(b0,a0%b0,x,y);  
        ll temp=x;x=y;  
        y=temp-a0/b0*y;  
    }  
}  

ll niyuan(ll t,ll p)
{
             ll x,y;
             extend(t,p,x,y);
             return (x+p)%p;
}

ll lucas(ll x,ll y)  
{  
    if(y>x) return 0;  
    ll re=1,dx,dy;  
    dx = x%p1;  
    dy = y%p1;  
    re*=c[dx][dy];  
    dx = x/p1;  
    dy = y/p1;  
    if(dx || dy)  
    {  
        re*=lucas(dx,dy);  
    }  
    re%=p1;  
    return re;  
}  

ll C(int n,int m)
{
         ll r1,r2,re;
         r1=lucas(n,m);
         r2=mul[n]*niyuan(mul[m]*mul[n-m]%p2,p2)%p2; 
         re=(r1*y1%mod+r2*y2%mod)%mod;
         return re;
}

void init()
{
     y1=p2*niyuan(p2,p1)%mod;
     y2=p1*niyuan(p1,p2)%mod;
     int i,j;
     for (i=1;i<100;i++)
         for (j=0;j<=i;j++)
         {
             if (j==0) c[i][j]=1;
             else if (j>i-j) c[i][j]=c[i][i-j];
             else c[i][j]=(c[i-1][j-1]+c[i-1][j])%p1;
         }
     mul[0]=1;
     for (i=1;i<N;i++) mul[i]=(mul[i-1]*i)%p2;
}

int main()
{
    int i,j,T,cas,n,k,m;
    ss(T);
    init(); 
    for (cas=1;cas<=T;cas++)
    {
        ss(n);ss(m);ss(k); 
        for (i=1;i<=k;i++) ss(a[i]);
        a[k+1]=m;
        ll res=1;
        for (i=1;i<=k;i++) 
        {
            int t=a[i+1]-a[i];
            if (t==1) res=(res*n)%mod;
            else res=res*C(n+t-1,t)%mod;
        }
        printf("Case #%d: %d\n",cas,res);
    }   
    return 0; 
}