catalan---卡特兰数(小结)

 

(关于卡特兰数的详细介绍)http://baike.baidu.com/view/2499752.htm 

下面有练习的题目:

经过测试,_int64/long long 最大只能表示到33位,超过这个范围就要用大数来表示。。。

有几个重要的公式是要记得的

F(n)=f(n-1)*(4*n-2)/(n+1);

还有 f(n)=C(2n,n)/(n+1);

或者 f(n)=c(2n,n)-c(2n,n+1);

http://acm.hdu.edu.cn/showproblem.php?pid=1023

代码:

 1 #include<stdio.h>
 2 #define maxn 100
 3 int arr[101][maxn+1]={{1},{1},{2}};
 4 //利用f(n)=f(n-1)*(4n-2)/(n+1);
 5 int main()
 6 {     
 7       int s, c=0,i,j,k;
 8     for(i=3;i<=100;i++)
 9     {
10         for(j=0;j<=maxn;j++)   //先计算乘法部分
11         {
12            arr[i][j]=arr[i-1][j]*(4*i-2)+c;
13            c=arr[i][j]/10;
14            arr[i][j]%=10;
15         }
16         s=c=0;
17         for(k=maxn;arr[i][k]==0;k--);      //计算除法部分
18         for(j=k;j>=0;j--)
19         {
20             s=arr[i][j]+10*c;
21             c=s%(i+1);
22             arr[i][j]=(s-c)/(i+1);
23         }
24     }
25 
26         int n;
27     while(scanf("%d",&n)!=EOF)
28     {
29         for(i=maxn;arr[n][i]==0;i--);
30         for(j=i;j>=0;j--)
31             printf("%d",arr[n][j]);
32            puts("");
33     }
34     return 0;
35 }
36 //后记,由于测试的时候,发现_int64/long long,都wa了,所以果断的用大数来
37 //做了...
38 //f(n)=f(n-1)*(4*n-2)/(n+1)......f(n)=c(2n,n)-c(2n,n+1)...
39 //f(n)=c(2n,n)/(n+1)...卡特兰数..很重要的...再次说明..希望不要嫌啰嗦、
40  
View Code

http://acm.hdu.edu.cn/showproblem.php?pid=4165

代码

 1 #include<iostream>
 2 using namespace std;
 3 int main()
 4 {
 5     long long  arr[31]={1,1,2,5};
 6     for(int i=4;i<=30;i++)
 7         arr[i]=arr[i-1]*(4*i-2)/(i+1);
 8     int n;
 9     while(cin>>n,n)
10         cout<<arr[n]<<endl;
11 return 0;
12 }
View Code

 

http://acm.hdu.edu.cn/showproblem.php?pid=1130

代码

 1 #include<iostream>
 2 #define maxn 100
 3 using namespace std;
 4 int arr[maxn+1][maxn+1]={{1},{1},{2}};
 5 int main()
 6 {
 7     int i,j,k,s,c,n;
 8     s=c=0;
 9     for(i=3;i<=maxn;i++)
10     {
11         for(j=0;j<=maxn;j++)
12         {
13             s=arr[i-1][j]*(4*i-2)+c;
14             c=s/10;
15             arr[i][j]=s%10;
16         }
17         s=c=0;
18         for(k=maxn;arr[i][k]==0;k--);
19         for(j=k;j>=0;j--)
20         {
21           s=arr[i][j]+10*c;
22           c=s%(i+1);
23           arr[i][j]=s/(i+1);
24         }
25 
26     }
27     while(cin>>n)
28     {
29         for(i=maxn;arr[n][i]==0;i--);
30         for(j=i;j>=0;j--)
31             cout<<arr[n][j];
32         cout<<endl;
33     }
34     return 0;
35 }
View Code

http://poj.org/problem?id=2084

代码:

 1  1 #include<iostream>
 2  2 #define maxn 80
 3  3 using namespace std;
 4  4 int arr[101][maxn+1]={{1},{1},{2}};
 5  5 int main()
 6  6 {
 7  7     int i,j,k,s,c,n;
 8  8     for(i=3;i<=100;i++)
 9  9     {
10 10       for(j=0,s=c=0;j<=maxn;j++)
11 11       {
12 12           s=(arr[i-1][j]*(4*i-2)+c);
13 13           arr[i][j]=s%10;
14 14           c=(s-arr[i][j])/10;
15 15       }
16 16       for(k=maxn;arr[i][k]==0;k--);
17 17 
18 18       for(j=k,s=c=0;j>=0;j--)
19 19       {
20 20         s=(arr[i][j]+10*c);
21 21         c=s%(i+1);
22 22         arr[i][j]=(s-c)/(i+1);    
23 23       }
24 24     }
25 25     while(cin>>n,n>=0)
26 26     {
27 27         for(i=maxn;arr[n][i]==0;i--);
28 28         for(j=i;j>=0;j--)
29 29             cout<<arr[n][j];
30 30         cout<<endl;
31 31     }
32 32     return 0;
33 33 }View Code 
View Code

 

关于卡特兰公式的扩展问题....

关于扩展的卡特兰数:
1.(n-m+1)/(n+1)*c(n+m,n)
2.c[n+m][n]-c[n+m][m-1]
Catalan,Eugene,Charles,卡特兰(1814~1894)比利时数学家,生于布鲁日(Brugge),早年在巴黎综合工科学校就读。1856年任列日(Liege)大学数学教授,并被选为比利时布鲁塞尔科学院院士。

卡特兰一生共发表200多种数学各领域的论著。在微分几何中,他证明了下述所谓的卡特兰定理:当一个直纹曲线是平面和一般的螺旋面时,他只能是实的极小曲面。他还和雅可比(Jacobi,C·G·J)同时解决了多重积分的变量替换问题,建立了有关的公式。

1842年,他提出了一种猜想:方程xz-yt=1没有大于1的正整数解,除非平凡情形32-23=1。这一问题至今尚未解决。

(mathoe注:即除了8、9这两个连续正整数都是正整数的方幂外,没有其他。1962年我国数学家柯召以极其精湛的方法证明了不存在三个连续正整数,它们都是正整数的方幂,以及方程x2-yn=1,n>1,xy≠0无正整数解。并且还证明了如果卡特兰猜想不成立,其最小的反例也得大于1016。)

此外,卡特兰还在函数论、伯努利数和其他领域也做出了一定的贡献。

卡特兰通过解决凸n边形的剖分得到了数列Cn

凸n+2边形用其n-1条对角线把此凸n+2边形分割为互不重叠的三角形,这种分法的总数为Cn

为纪念卡特兰,人们使用“卡特兰数”来命名这一数列。

据说有几十种看上去毫不相干的组合计数问题的最终表达式都是卡特兰数的形式。

卡特兰数在数学竞赛、信息学竞赛、组合数学、计算机编程等都会有其不同侧面的介绍。

前几个卡特兰数:规定C0=1,而

C1=1,C2=2,C3=5,C4=14,C5=42,

C6=132,C7=429,C8=1430,C9=4862,C10=16796,

C11=58786,C12=208012,C13=742900,C14=2674440,C15=9694845。

递推公式

圆周上有标号为1,2,3,4,……,2n的共计2n个点,这2n个点配对可连成n条弦,且这些弦两两不相交的方式数为卡特兰数Cn

2003年浙江省小学数学夏令营竞赛考了这个题:圆周上10个点可以连成既不相交,也没有公共端点的5条线段,不同的连法共有_____种。

答:方法的种数是卡特兰数C5=42,此题被收录进单墫主编的知识出版社出版的《华数奥赛强化训练》小学六年级册的“计数问题”专题。


共六种类型,第1类有5种连法,第2类有2种连法,第3类有10种连法,第4类有10种连法,第5类有10种连法,第6类有5种连法。共有42种连法。

1994年《小学数学》有奖征答竞赛:游乐园门票1元一张,每人限购一张。现在有10个小朋友排队购票,其中5个小朋友每人只有1元的钞票一张,另5个小朋友每人只有2元的钞票一张,售票员没有准备零钱。问:有多少种排队方法,使售票员总能找的开零钱?

(此题也被许多奥数资料收录为例题或习题,《华罗庚学校数学课本》小学六年级册的思维训练也收有此题)

答:现把拿1元的5个小朋友看成是相同的,把拿2元的5个小朋友也看成是相同的,使用我们常用的“逐点累加法”:

图中每条小横段表示拿1元的小朋友,每条小竖段表示拿2元的小朋友,要求从A走到B的过程中网格中任何点均有横段数不小于竖段数:拿1元的要先,且人数不能少于拿2元的,即不能越过对角线AB:每个点所标的数即为从A走到此点的方法数。求从A到B的走法的方法数。逐点累加可求出为42,即卡特兰数C5=42。



又由于每个小朋友是不相同的,所以共有42×5!×5!=42×120×120=604800种情况。

若把此题的10个人,拿1元的有5人,拿2元的有5人改为共有2n个人,拿1元的n人,拿2元的n人,则符合要求的排队方法数为:





再一个卡特兰数的例子:

甲乙两人比赛乒乓球,最后结果为20∶20,问比赛过程中甲始终领先乙的计分情形的种数。

即甲在得到1分到19分的过程中始终领先乙,其种数是卡特兰数





再一个卡特兰数的例子

饭后,姐姐洗碗,妹妹把姐姐洗过的碗一个一个放进碗橱摞成一摞。一共有n个不同的碗,洗前也是摞成一摞的,也许因为小妹贪玩而使碗拿进碗橱不及时,姐姐则把洗过的碗摞在旁边,问:小妹摞起的碗有多少种可能的方式?

答:得数是第n个卡特兰数Cn

再一个卡特兰数的例子

一个汽车队在狭窄的路面上行驶,不得超车,但可以进入一个死胡同去加油,然后再插队行驶,共有n辆汽车,问共有多少种不同的方式使得车队开出城去?

答:得数是第n个卡特兰数Cn

卡特兰数

求证:卡特兰数Cn是整数。

证明:

①取整函数不等式:对任意实数x,y有[x+y]≥[x]+[y]。这里[x]表示不大于实数x的最大整数。

解:由定义x≥[x]……(1)

                 y≥[y]……(2)以上两式相加,得:x+y≥[x]+[y],

       把上式再取整,得:[x+y]≥[[x]+[y]]=[x]+[y],即[x+y]≥[x]+[y]。

②1000!的末尾0的个数249个。(现在有的小学奥数书上出现了100!末尾有几个零的题目:24个)

解:1000÷5=200,

        200÷5=40,

        40÷5=8,

        8÷5=1……3

       以上各商相加,即得1000!末尾0的个数=200+40+8+1=249个。

③n!的质因数分解式中质因子p的幂次数:


…………(1)

k!的质因数分解式中质因子p的幂次数

…………(2)

(n-k)!的质因数分解式中质因子p的幂次数


…………(3)

这里写成西格马求和式时使用了无穷的形式,但是从某一确定项之后的每项都是0,为了统一,都写成了“∞”形式。

④组合数是整数

解:

⑤卡特兰数是整数

⑥卡特兰数是整数的另外一个证明

④组合数是整数




⑤卡特兰数是整数



⑥卡特兰数是整数的另一个证明




凸六边形剖分成三角形的14种方法,是卡特兰数C4





从左下角(0,0)走到右上角(4,4),只允许向上、向右走,但不允许穿过对角线的方法数是14种,是卡特兰数C4

1936第40届匈牙利奥林匹克数学竞赛 第1题考了Catalan恒等式的证明。





1979第21届国际数学奥林匹克 第1题考了一个卡特兰恒等式的应用的题目






此题由1989年第1届匈牙利-以色列数学竞赛题改编。

有下面的练习题:

http://acm.hdu.edu.cn/showproblem.php?pid=2067

代码:

 1 #include<iostream>
 2 #define maxn 40
 3 using namespace std;
 4 int arr[36][maxn+1]={{0},{2}};
 5 int main()
 6 {
 7     int i,j,k,c,s;
 8     for(i=2;i<=35;i++)
 9     {
10         for(j=0,c=0,s=0;j<=maxn;j++)
11         {
12           s=(arr[i-1][j]*(4*i-2)+c);
13           arr[i][j]=s%10;
14           c=(s-arr[i][j])/10;
15         }
16         for(j=maxn;arr[i][j]==0;j--);
17         for(k=j,s=c=0;k>=0;k--)
18         { 
19            s=(arr[i][k]+c*10);
20            c=s%(i+1);
21            arr[i][k]=(s-c)/(i+1);
22         }
23     }
24     int n,count=1;
25     while(cin>>n,n>0)
26     {
27         cout<<count++<<" "<<n<<" ";
28         for(i=maxn;arr[n][i]==0;i--);
29         for(j=i;j>=0;j--)
30             cout<<arr[n][j];
31         cout<<endl;
32     }
33     return 0;
34 }
View Code

暂时,就归纳这么多吧!!。。。。

posted @ 2013-08-02 15:25  龚细军  阅读(2398)  评论(0编辑  收藏  举报