栈 / 卡特兰数：火车进出栈问题

传送门

题意：一列火车n节车厢，依次编号为1,2,3,…,n。
每节车厢有两种运动方式，进栈与出栈，问n节车厢出栈的可能排列方式有多少种。
数据范围
1≤n≤60000
输入样例：
3
输出样例：
5
思路：进出栈问题–>序列问题–>路径问题。
这显然是上一道火车进栈的升级版，数据太大不能再进行递归处理了。问题转换：寻找满足条件的 ‘+’,’-‘序列，即把进栈用’+‘标记，出栈用’-‘标记；满足前缀和中’-‘的数量小于等于’+‘的数量就是合法的。例如：231转换成序列便是+ + - + - -。
就相当于在2n个位置找n个位置放+的方案数减去不合法的方案数。

而这个不合法数怎么求呢？这就很神奇了，我们可以把它转化成路径问题：在一个坐标走上，从(0,0)走到(n,n)的位置，往右表示’+’，往上表示’-’；那么所有在y=x以下的路径都是合法路径(满足x>=y，即前缀和的’+’ > ‘-’)；而所有与y=x+1有交点的路径都是不合法的路径，那么我们找到每条不合法路径与y=x+1的第一个交点，将交点以后的部分关于y=x+1对称后终点都会是（n-1,n+1），如此每一天不合法路径都和其改变后的路径相对应，由此所有从（0，0）走到（n-1,n+1）的路径就是不合法的路径。

显然的用高精度来解决这个问题，但是如果我们先算乘法再来做除法的话就会超时，这样我们可以想办法只用乘法：找出分子分母共有的质因子数，将它们相互抵消（如分子含义2 ^ x,分母含义2 ^ y，那么我们就可以相互抵消掉后分母乘以2^(y-x)）；如此找到所有的质因子数并进行抵消最后得到答案。
代码实现：

#include <iostream>
#include <vector>
#include<stdio.h>
using namespace std;
typedef long long ll;
const int N = 6000010, M = 120010;


ll res[N], tt;
int q[M];
bool s[M];
//高精度乘法
void multi(int b)
{
    ll t = 0;
    for (int i = 0; i <= tt; i ++ ){
        res[i] = res[i] * b + t;
        t = res[i] / 1000000000;
        res[i] %= 1000000000;
    }
    while (t){
        res[++tt] = t % 1000000000;
        t /= 1000000000;
    }
}
//压位输出，可以提高速度
void out()
{
    printf("%lld", res[tt]);
    for (int i = tt - 1; i >= 0; i -- ) 
         printf("%09lld", res[i]);
    cout << endl;
}
//返回n!里面有多少个因子p，即n!里面p的次数是多少
int get(int n, int p)
{
    int s = 0;
    while (n){
    s += n / p;
     n /= p;
    }
    return s;
}

int main()
{
    int n;
    cin >> n;
    //筛素数
    for (int i = 2; i <= 2 * n; i ++ )
        for (int j = 2 * i; j <= 2 * n; j += i)
            s[j] = 1;

    for (int i = 2; i <= n * 2; i ++ )
        if (!s[i]){
            q[i] = get(n * 2, i) - get(n * 2 - n, i) * 2;
        }

    int k = n + 1;
    for (int i = 2; i <= k; i ++ )
        while (k % i == 0){
            k /= i;
            q[i] -- ;
        }

    res[0] = 1;
    for (int i = 2; i <= n * 2; i ++ )
        while (q[i] -- )
            multi(i);

    out();
    return 0;
}