额,同一个代码,贴了三道题 ,都不用改的,卡特兰数的应用,还有hdu1130,hdu1134
出栈次序问题,其实就是卡特兰数的应用
卡特兰公式:
令h(1)=1,h(0)=1,catalan数满足递归式:
h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)h(0) (其中n>=2)
另类递归式:
h(n)=((4*n-2)/(n+1))*h(n-1);
该递推关系的解为:
h(n)=C(2n,n)/(n+1) (n=1,2,3,...)
用递归式来做,接下来就是一个大数问题了,(4*n-2)/(n+1) 不能直接算,未必整除呀ORZ
先算(4*n-2)*h(n-1) ,再除以(n+1)
#include <iostream> using namespace std; #define SIZE 100 int main() { int a[101][SIZE]={0}; //数组用来存放结果 a[1][0]=1; //初始化第一项 int i,j,r=0,temp=0,len=1; //len 表示当前最长的有效位数,初始为1 for(i=2;i<=100;i++) //从第二项到第100项,用公式计算 { for(j=0;j<len;j++) //---------------乘法部分------------------ { a[i][j]=a[i-1][j]*(4*i-2); //乘法从低位到高位 } for(j=0;j<len;j++) //对乘出的结果进行处理,不包括最高位 { temp=a[i][j]+r; a[i][j]=temp%10; r=temp/10; } while(r) //对最高位进位处理 { a[i][len]=r%10; r/=10; len++; } //-----------------除法部分----------------- for(j=len-1,r=0;j>=0;j--) { //除法从高位到低位 temp=r*10+a[i][j]; a[i][j]=temp/(i+1); r=temp%(i+1); } while(!a[i][len-1]) //处理高位的零位 len--; } //------------------------------------------- int n; while(cin>>n&&n!=-1) { for(i=SIZE-1;!a[n][i];i--); for(i;i>=0;i--) cout<<a[n][i]; cout<<endl; } return 0; }