1. 程式人生 > >[AHOI2012]樹屋階梯,洛谷P2532,Catalan+高精度

[AHOI2012]樹屋階梯,洛谷P2532,Catalan+高精度

正題

      一看到輸入3,輸出5,我們就可以想到Catalan數,但是要詳細說明一下。

      因為每次可以從n個角中選取一個角,將他填充掉,然後把左右兩邊分成兩部分,所以很明顯:

      f(n)=\sum_{i=0}^{n-1}f(i)*f(n-i-1)=Catalan(n)

      又因為它的n是500,所以高精度即可。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;

int n;
struct node{
	int l,d[10010];
	node operator*(int a)const{
		node q;
		for(int i=1;i<=l;i++)
			q.d[i]=d[i]*a;
		q.l=l;
		for(int i=1;i<=q.l;i++)
			if(q.d[i]>9){
				if(i==q.l) q.l++,q.d[i+1]=0;
				q.d[i+1]+=q.d[i]/10;
				q.d[i]%=10;
			}
		return q;
	}1
	node operator/(int a)const{
		node q;
		int tot=0;
		for(int i=l;i>=1;i--){
			tot=tot*10+d[i];
			q.d[i]=tot/a;
			tot%=a;
		}
		q.l=l;
		while(q.d[q.l]==0 && q.l>1) q.l--;
		return q;
	}
	
}h[510];

int main(){
	scanf("%d",&n);
	h[0].l=h[1].l=1;
	h[0].d[1]=h[1].d[1]=1;
	for(int i=2;i<=n;i++) h[i]=h[i-1]*(4*i-2)/(i+1);
	for(int i=h[n].l;i>=1;i--) printf("%d",h[n].d[i]);
}