火车进出栈问题

一列火车 n 节车厢，依次编号为 1,2,3,…,n。

每节车厢有两种运动方式，进栈与出栈，问 nn 节车厢出栈的可能排列方式有多少种。

输入格式

输入一个整数 n，代表火车的车厢数。

输出格式

输出一个整数 s 表示 n 节车厢出栈的可能排列方式数量。

数据范围

1≤n≤60000

输入样例：

3

输出样例：

5

分析： 这道题n的数值太大，用dfs是过不了的，为此我就想着直接算出有多少种情况，假如规定 1 表示栈中进入一节车厢，0 表示栈中弹出一节车厢，那么每种火车进出栈方案都能用一串 01 序列来表示。
由于每节车厢都要进栈一次出栈一次，所以每种方案所对应的 01 序列的长度都是 2n。
那么火车进出栈方案的总数量就可以看成是所有合法的 01 串的个数。

现在我们来想一下什么样的01串是合格的，我们要知道的是栈中有车厢才能弹出，所以在任意位置，其前面的 1 的个数大于等于其前面 0 的个数。

这里就转化成了求卡特兰数问题（这个大家有兴趣的可以看一下如何推导），我们可以知道答案为

你以为这道题就这样结束了，大错特错！！！如是进入了更恶心的环节，n的阶乘和（2*n)的阶乘都太大了！！！因此我们只能用一种可以在数据不超范围的情况下求阶乘的方法，这里我们是吧（2*n)的阶乘和n的阶乘全部质因式分解，相同的因式全部消掉，这里用到的方法参考我的这道题阶乘分解。

然而！！！这样还不行哭了哭了，因为最后的结果还是很大不能让剩下的数直接相乘，这里我用到了vector来存储每一位数字，每次相乘的时候拿出一位相乘，这是一种高精度相乘的方式，具体可以看代码中的mul（）函数，但尽管我这样写了又报了超时的错！！！！

姐妹们兄弟们，如果你们能看到这里，我敬你们是条汉子，又报错了我该怎么办呢，我原本是用vector数组每个值都只存了一位数字，于是我让它们都存了八位，你以为这样过了？

nonono还有最后一丢丢丢错误，哎！！字字全是血泪，在输出的时候除了数组的最后一位，其他记得不够八位的补0，这个大家知道为什么把，最后一位中的八位数是整个数的最前面的数，前缀的0可以不要，但其他的八位数都是中间的，所以0不能省略！！！

代码来了来了：

#include <bits/stdc++.h>

using namespace std;
const int N=120010,mod=1e8;
int st[N],prim[N],nu[N];
int cnt;
int n;
vector<long long> ans;
void func(int x)
{
    for(int i=2;i<x;i++)
    {
        if(st[i]) continue;
        for(int j=2;j*i<=x;j++) st[j*i]=1;
    }
    for(int i=2;i<=x;i++)
    {
        if(!st[i]) prim[++cnt]=i;
    }
}
vector<long long> mul(vector<long long> an,int x)
{
    vector<long long>te;
    long long t=0;
    for(int i=0;i<an.size();i++)
    {
        t+=an[i]*x;
        te.push_back(t%mod);
        t/=mod;
    }
    while(t)
    {
        te.push_back(t%mod);
        t/=mod;
    }
    return te;
}
int main()
{
    cin>>n;

    func(2*n);
    for(int i=1;i<=cnt;i++)
    {
        int k=2*n;
        int num=0;
        while(k)
        {
            k/=prim[i];
            num+=k;
        }
        nu[prim[i]]=num;
    }

    for(int i=1;i<=cnt&&prim[i]<=n;i++)
    {
        int k=n;
        int num=0;
        while(k)
        {
            k/=prim[i];
            num+=k;
        }
        nu[prim[i]]-=2*num;
    }
    int k=n+1;
    for(int i=1;i<=cnt&&prim[i]<=k;i++)
    {
        int num=0;
        while(k%prim[i]==0)
        {
                num++;
                k/=prim[i];
        }
        nu[prim[i]]-=num;
    }
    ans.push_back(1);
    for(int i=1;i<=2*n;i++)
    {
        while(nu[i]--) ans=mul(ans,i);
    }
    
    printf("%lld",ans.back());
    for(int i=ans.size()-2;i>=0;i--) printf("%08lld",ans[i]);
    return 0;
}