Train Problem II
不得不说,这是一道特别好的题目!我从一开始的找不到思绪,再到了解卡特兰数,了解到大数的乘除法,每一步都收获良多!
先理清一下学到的知识点吧!
1.卡特兰数!又称
卡塔兰数
卡塔兰数是组合数学中一个常在各种计数问题中出现的数列。以比利时的数学家欧仁·查理·卡塔兰
卡塔兰数的一般项公式为
前几项为 (OEIS中的数列A000108): 1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796,
58786, 208012, 742900,2674440, 9694845, 35357670, 129644790, 477638700, 1767263190,
6564120420, 24466267020, 91482563640,343059613650, 1289904147324, 4861946401452, ...
性质
Cn的另一个表达形式为所以,Cn是一个自然数;这一点
在先前的通项公式中并不显而易见。这个表达形式也是André对前一公式证明的基础。
它也满足
-
//这个在本题中用来推导卡特兰数!
- 很重要的一个推导式!!
这提供了一个更快速的方法来计算卡塔兰数。
卡塔兰数的渐近增长为
它的含义是左式除以右式的商趋向于1当n→∞。(这可以用n!的斯特灵公式来证明。)
所有的奇卡塔兰数Cn都满足。所有其他的卡塔兰数都是偶数。
卡特兰数的应用:
出栈次序
一个 栈 (无穷大)的 进栈 序列为1,2,3,…,n,有多少个不同的 出栈 序列?
2.在c以及c++中,EOF被定义成-1!
3.如何处理大数的乘除法!很妙的方法!参考代码!
代码如下:
//输出卡特兰数
//首先需要肯定,程序是正确的
//这算是大数乘除法!记住他们是如何处理的!由于数据很大,用基本数据类型根本无法满足要求,只能用数组来表示!
#include <iostream>
#include<cstdio>
#include<memory.h>
using namespace std;
#define MAX 101
#define BASE 10000//base只是一个基度,对最终取值并没有影响,相反,base取值愈大,计算量愈小
//base发生改变的时候,下面的输出也要相应地做出调整,否则也会输出错误答案!除非当base取10!
void multiply(int a[],int len,int b)//乘法
{
for(int i=len-1,carry=0;i>=0;--i)//从最后一位开始相乘,依次向前与每一位相乘
{//问题在于,为什么BASE=10000?
carry+=b*a[i];
a[i]=carry%BASE;
carry/=BASE;
//cout<<"carry="<<carry<<" "<<"a["<<i<<"]="<<a[i]<<endl;//以4个0为一组
}
}
void divide(int a[],int len,int b)//除法,很妙的!这种除法可能想不到,仔细体会!
{//应当如何除呢?
for(int i=0,div=0;i<len;++i)//从高位除起
{
div=div*BASE+a[i];
a[i]=div/b;//b为除数
div%=b;
}
}
int main()
{
int i,j,h[101][MAX];
memset(h[1],0,MAX*sizeof(int));//赋值,每一个都置为0
for(i=2,h[1][MAX-1]=1;i<=100;++i)//运用递归,并且h[1]=1;
{
memcpy(h[i],h[i-1],MAX*sizeof(int));//h[i]=h[i-1];按字节拷贝,保证了h[i]和h[i-1]指向数组的一致性
multiply(h[i],MAX,4*i-2);//h[i]*=(4*i-2);
divide(h[i],MAX,i+1);//h[i]/=(i+1);
}//递归得到前100项的卡特兰数!
while(cin>>i && i>=1 && i<=100)//输入i的值
{
// for(i=1;i<=100;i++)
// {
for(j=0;j<MAX && h[i][j]==0;++j);//从0位开始搜索,找到不为0的第一个数
//printf("%d\n",EOF);在c语言中,EOF=-1;
printf("%d",h[i][j++]);//像是这个输出,就很妙了,第一位可能不足四位,就地输出!
for(;j<MAX;++j)
{
// if(h[i][j]==0)
printf("%04d",h[i][j]);//处在中间的值也可能没有四位,这时候要注意了,往左边加0,凑足4位,不然答案会出错!
// else
// printf("%d",h[i][j]);//不断输出值
}
printf("\n");
}
system("pause");
return 0;
}