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\n",rs);

    }

    return 0;

}

巴特西

HDU 3398 String

最新文章

热门文章