BZOJ2822: [AHOI2012]树屋阶梯

Catalan数，质因数分解后高精乘单精

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 1e5+10;

int s[maxn];
int p[maxn],pri,minp[maxn];
bool v[maxn];

void pre()
{
    for(int i=2;i<maxn;i++)
    {
        if(!v[i]) minp[i]=i,p[++pri]=i;
        for(int j=1;j<=pri;j++)
        {
            int k=i*p[j];
            if(k>maxn) break;
            v[k]=true; minp[k]=p[j];
            if(i%p[j]==0) break;
        }
    }
}

int f[maxn],fn;

int main()
{
    pre();
    int n; scanf("%d",&n);
    for(int i=2*n;i>n;i--)
        for(int k=i;k!=1;k/=minp[k]) s[minp[k]]++;
    for(int i=2;i<=n;i++)
        for(int k=i;k!=1;k/=minp[k]) s[minp[k]]--;
    for(int k=n+1;k!=1;k/=minp[k]) s[minp[k]]--;

    fn=0;
    for(int i=1;i<=2*n;i++)
    {
        while(s[i]--)
        {
            if(!fn) 
            {
                for(int k=i;k;k/=10) f[++fn]=k%10;
            }
            else
            {
                for(int j=1;j<=fn;j++) f[j]*=i;
                for(int j=1;j<=fn;j++) if(f[j]>9)
                {
                    f[j+1]+=f[j]/10;
                    f[j]%=10;
                }
                while(f[fn+1])
                {
                    fn++;
                    if(f[fn]>9) f[fn+1]+=f[fn]/10;
                    f[fn]%=10;
                }
            }
        }
    }
    while(fn) printf("%d",f[fn--]);
    printf("\n");

    return 0;
}