题意:一列火车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;
}