bzoj2822 [AHOI2012]树屋阶梯（高精度+Catalan）

Description

暑假期间，小龙报名了一个模拟野外生存作战训练班来锻炼体魄，训练的第一个晚上，教官就给他们出了个难题。由于地上露营湿气重，必须选择在高处的树屋露营。小龙分配的树屋建立在一颗高度为N+1尺（N为正整数）的大树上，正当他发愁怎么爬上去的时候，发现旁边堆满了一些空心四方钢材（如图1.1），经过观察和测量，这些钢材截面的宽和高大小不一，但都是1尺的整数倍，教官命令队员们每人选取N个空心钢材来搭建一个总高度为N尺的阶梯来进入树屋，该阶梯每一步台阶的高度为1尺，宽度也为1尺。如果这些钢材有各种尺寸，且每种尺寸数量充足，那么小龙可以有多少种搭建方法？（注：为了避免夜里踏空，钢材空心的一面绝对不可以向上）

以树屋高度为4尺、阶梯高度N=3尺为例，小龙一共有如图1.2所示的5种
搭建方法：

Input

一个正整数 N(1≤N≤500)，表示阶梯的高度

Output

一个正整数，表示搭建方法的个数。（注：搭建方法个数可能很大。）

Sample Input
3

Sample Output
5

HINT
1≤N≤500

分析：
这个问题就是我们在学习Catalan的时候总结的典型问题
那写一个Catalan的板子不就好了吗

然而当我在写KSM的时候
发现没有模数了，很慌啊
我们只能选择高精度了（或者来个NTT？）

tip

注意初始化
s和len都要清零
板子一定要打对
高精度一定要写对！！！

//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int N=2003;
int n;
int sshu[N],tot,num[N],mp[N];
bool no[N];
struct node{
    int len;
    int s[2002];
    void clear()
    {
        len=0;
        memset(s,0,sizeof(s));
    }
};
node ans;

node operator *(const node &a,const node &b)
{
    node an;
    an.clear();
    int n=a.len+b.len;
    int len=max(a.len,b.len);
    for (int i=1;i<=len;i++)
    {
        int x=0;
        for (int j=1;j<=len;j++)
        {
            an.s[i+j-1]+=a.s[i]*b.s[j]+x;
            x=an.s[i+j-1]/10;
            an.s[i+j-1]%=10;
        }
        an.s[len+i]=x;
    }
    while (an.s[n]==0) n--;
    an.len=n;
    return an;   
}

void cl()
{
    memset(no,0,sizeof(no));
    for (int i=2;i<=N;i++)
    {
        if (!no[i]) sshu[++tot]=i,mp[i]=tot;
        for (int j=1;j<=tot&&sshu[j]*i<=N;j++)
        {
            no[sshu[j]*i]=1;
            if (i%sshu[j]==0) break;
        }
    }
}

void calc(int x,int v)
{
    for (int i=1;i<=tot&&sshu[i]*sshu[i]<=x;i++)
        if (x%sshu[i]==0)
            while (x%sshu[i]==0)
            {
                num[i]+=v;
                x/=sshu[i];
            }
    if (x>1) num[mp[x]]+=v;  //
}

node KSM(int a,int b)
{
    node aa;
    node t;
    t.clear();  aa.clear();
    t.len=1;t.s[1]=1;
    while (a)
    {
        aa.len++;
        aa.s[aa.len]=a%10;
        a/=10;
    }
    while (b)
    {
        if (b&1)
           t=t*aa;
        b>>=1;
        aa=aa*aa;
    }
    return t;
}

int main()
{
    scanf("%d",&n);
    cl();
    for (int i=n+1;i<=n*2;i++) calc(i,1);
    for (int i=1;i<=n;i++) calc(i,-1);
    calc(n+1,-1);
    ans.clear();
    ans.len=1; ans.s[1]=1;
    for (int i=1;i<=tot;i++)
        ans=ans*KSM(sshu[i],num[i]);
    for (int i=ans.len;i>=1;i--)
        printf("%d",ans.s[i]);
    return 0;
}