HDU 3398 String

题目大意：一个长为n的01字符串，使前缀任意0的数量不大于1的数量，求方案数……

题解：高一模拟赛时做过，是卡特兰数的几何意义，将字符串变为矩阵寻路，不可越过对角线，那么就是卡特兰数了，C(n+m, n)-C(n+m,n+1)=(n+1-m)(n+m)!/(n+1)!m!。需要注意的是取模的问题，如果用高精度最后取模会太慢了，会超时，所以直接用power定理分解素数，对于每个素数分别算幂，取模相乘即可。

#include <cstdio>  
#include <cstring>   
const int N=1000001;  
using namespace std;  
typedef long long  LL;   
int pri[N*2+1],p[N*2+1],tot;
int cal(int pr,int n){int rs=0;while(n)n/=pr,rs+=n;return rs;}  
void initp(){  
    memset(pri,0,sizeof pri); tot=0;  
    for(int i=2;i<=N*2;i++){  
        if(pri[i])continue; p[tot++]=i; 
        for(int j=i*2;j<=N*2;j+=i)pri[j]=1;  
    }  
}   
int main(){  
    int cas,n,m; initp();  
    scanf("%d",&cas);  
    while(cas--){  
        scanf("%d%d",&n,&m);         
        LL rs=1; int nm=n-m+1;  
        for(int i=0;i<tot&&p[i]<=n+m;i++){  
            int cnt=0;  
            while(nm%p[i]==0)nm/=p[i],cnt++;  
            int ipow=cnt+cal(p[i],n+m)-cal(p[i],n+1)-cal(p[i],m);  
            for(int j=1;j<=ipow;j++){  
                rs=(rs*p[i])%20100501;  
            }  
        }  
        printf("%lld
",rs);  
    }  
    return 0;  
}