JoyOI火车进出栈问题(水水版+普通版+强化版) 高精度压位+筛素数+统计约数+快速幂=卡特兰数

题目

一列火车n节车厢，依次编号为1,2,3,…,n。每节车厢有两种运动方式，进栈与出栈，问n节车厢出栈的可能排列方式有多少种。
水水版：n<=100
普通版：n<=30000
加强版：n<=50000

题解

$Catalan$数列
给定n个0和n个1，它们按照某种顺序排成长度为2n的序列，满足任意前缀中0的个数都不少于1的个数的序列的数量为：$$Catala{n}_n=C_{2n}^n/(n+1)$$
$Catalan$数还与以下问题有关：
1.n个左括号和n个右括号组成的合法括号序列的数量为$Catala{n}_n$。
2.1,2,…n经过一个栈，形成的合法出栈序列的数量为$Catala{n}_n$。
3.n个节点构成的不同二叉树的数量为$Catala{n}_n$。
4.在平面直角坐标系上，每一步只能向上或向右走，从(0,0)走到(n,m)并且除两个端点外不接触直线y=x的路线数量为$2Catala{n}_{n-1}$。

那么很显然这题就是求$Catala{n}_n$了
然后发现数据好大。。。QwQ。。。emmmm…
优化+优化+优化+优化+…!!!
首先，上一个高精度
然后，压一下高精度的位（4~8位）
接着，筛一下计算质因数所需的质数
下一步，统计每种质因数出现的数量（除法转化为质因数数量的加减，而质因数数量的计算又要用奇奇怪怪的方法：1—n中共有n/i个约数i）
紧接着，上一个快速幂（直接把每个约数的次方计算出来）
最后，输出！！

这个算法效率还是很高的
代码

#include <cstdio>
#include <cmath>

using namespace std;

const int max=5002;
int n,l,ss;
unsigned long long a[max];
int bz[100005],z[100005],b[100005];

unsigned long long ksm(int a,int b){
	int g=b;
	unsigned long long x=1,c=a;
	while (g){
		if (g&1) x=x*c;
		c*=c;
		g>>=1;
	}
	return x;
}

int solve(int v,int l){
	int z=0,cc=l;
	while (l<=v){
		z+=v/l;
		l*=cc;
	}
	return z;
}

int main(){
	scanf("%d",&n);
	if (n==1){printf("1");return 0;} 	
	
	for (int i=2;i<=2*n;i++)
		if (!bz[i]){
			int k=2;
			z[++ss]=i;
			while (k*i<=2*n){
				bz[k*i]=1;
				k++;
			}
		}
	
	for (int i=1;i<=ss;i++)
		b[z[i]]-=solve(n-1,z[i]);
	for (int i=1;i<=ss;i++)
		b[z[i]]-=solve(n+1,z[i]);
	for (int i=1;i<=ss;i++)
		b[z[i]]+=solve(2*n-1,z[i]);
	b[2]++;		
	
	a[max-1]=1;l=max-1;
	for (int i=1;i<=ss;i++)
	if (b[z[i]]){
		unsigned long long k=ksm(z[i],b[z[i]]);
		unsigned long long x=0;
		for (int j=max-1;j>=l;j--){
			a[j]=a[j]*k+x;
			x=a[j]/100000000;
			a[j]%=100000000;
		}
		while (x){
			l--;
			a[l]=x%100000000;
			x/=100000000;
		}
	}
	printf("%d",a[l]);
	for (int i=l+1;i<=max-1;i++){
		if (a[i]<1) {printf("00000000");continue;}
		if (a[i]<10) {printf("0000000%d",a[i]);continue;}
		if (a[i]<100) {printf("000000%d",a[i]);continue;}
		if (a[i]<1000) {printf("00000%d",a[i]);continue;}
		if (a[i]<10000) {printf("0000%d",a[i]);continue;}
		if (a[i]<100000) {printf("000%d",a[i]);continue;}
		if (a[i]<1000000) {printf("00%d",a[i]);continue;}
		if (a[i]<10000000) {printf("0%d",a[i]);continue;}
		printf("%d",a[i]);
	}
}