卡特蘭數 性質及應用
卡塔蘭數是組合數學中一個常在各種計數問題中出現的數列。以比利時的數學家歐仁·查理·卡塔蘭(1814–1894)命名。歷史上,清代數學家明安圖(1692年-1763年)在其《割圜密率捷法》最早用到“卡塔蘭數”,遠遠早於卡塔蘭[1][2][3]。有中國學者建議將此數命名為“明安圖數”或“明安圖-卡塔蘭數”[4]。
卡塔蘭數的一般項公式為
前20項為(OEIS中的數列A000108):1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900, 2674440, 9694845, 35357670, 129644790, 477638700, 1767263190
性質[編輯]
Cn的另一個表達形式為所以,Cn是一個自然數;這一點在先前的通項公式中並不顯而易見。這個表達形式也是André對前一公式證明的基礎。(見下文的第二個證明。)
它也滿足
這提供了一個更快速的方法來計算卡塔蘭數。
卡塔蘭數的漸近增長為
它的含義是當n → ∞時,左式除以右式的商趨向於1。(這可以用n!的斯特靈公式來證明。)
所有的奇卡塔蘭數Cn都滿足。所有其他的卡塔蘭數都是偶數。
應用[編輯]
組合數學中有非常多的組合結構可以用卡塔蘭數來計數。在Richard P. Stanley的Enumerative Combinatorics: Volume 2一書的習題中包括了66個相異的可由卡塔蘭數表達的組合結構。以下用n=3和n=4舉若干例:
- Cn表示長度2n的dyck word的個數。Dyck word是一個有n個X和n個Y組成的字串,且所有的字首字串皆滿足X的個數大於等於Y的個數。以下為長度為6的dyck words:
- 將上例的X換成左括號,Y換成右括號,Cn表示所有包含n組括號的合法運算式的個數:
- Cn表示有n個節點組成不同構二叉樹的方案數。下圖中,n等於3,圓形表示節點,月牙形表示什麼都沒有。
- Cn表示有2n+1個節點組成不同構滿二叉樹(full
binary tree)的方案數。下圖中,n
證明:
令1表示進棧,0表示出棧,則可轉化為求一個2n位、含n個1、n個0的二進位制數,滿足從左往右掃描到任意一位時,經過的0數不多於1數。顯然含n個1、n個0的2n位二進位制數共有個,下面考慮不滿足要求的數目。
考慮一個含n個1、n個0的2n位二進位制數,掃描到第2m+1位上時有m+1個0和m個1(容易證明一定存在這樣的情況),則後面的0-1排列中必有n-m個1和n-m-1個0。將2m+2及其以後的部分0變成1、1變成0,則對應一個n+1個0和n-1個1的二進位制數。反之亦然(相似的思路證明兩者一一對應)。
從而。證畢。
- Cn表示所有在n×n格點中不越過對角線的單調路徑的個數。一個單調路徑從格點左下角出發,在格點右上角結束,每一步均為向上或向右。計算這種路徑的個數等價於計算Dyck word的個數:X代表“向右”,Y代表“向上”。下圖為n= 4的情況:
- Cn表示對{1, ...,n}依序進出棧的置換個數。一個置換w是依序進出棧的當S(w) = (1, ...,n),其中S(w)遞迴定義如下:令w=unv,其中n為w的最大元素,u和v為更短的數列;再令S(w) =S(u)S(v)n,其中S為所有含一個元素的數列的單位元。
- Cn表示集合{1, ...,n}的不交叉劃分的個數.那麼,Cn永遠不大於第n項貝爾數.Cn也表示集合{1, ..., 2n}的不交叉劃分的個數,其中每個段落的長度為2。綜合這兩個結論,可以用數學歸納法證明:在魏格納半圓分佈定律中度數大於2的情形下,所有自由的累積量s 為零。 該定律在自由概率論和隨機矩陣理論中非常重要。
- Cn表示用n個長方形填充一個高度為n的階梯狀圖形的方法個數。下圖為n = 4的情況:
- Cn表示表為2×n的矩陣的標準楊氏矩陣的數量。 也就是說,它是數字 1, 2, ..., 2n被放置在一個2×n的矩形中並保證每行每列的數字升序排列的方案數。同樣的,該式可由勾長公式的一個特殊情形推導得出。
- Cn表示n個無標號物品的半序的個數。
附上bin神的模板:
#include <iostream>
#include <stdio.h>
#include <algorithm>
#include <string.h>
using namespace std;
/*
* 完全大數模板
* 輸入cin>>a
* 輸出a.print();
* 注意這個輸入不能自動去掉前導0的,可以先讀入到char陣列,去掉前導0,再用建構函式。
*/
#define MAXN 9999
#define MAXSIZE 1010
#define DLEN 4
class BigNum
{
private:
int a[500]; //可以控制大數的位數
int len;
public:
BigNum(){len=1;memset(a,0,sizeof(a));} //建構函式
BigNum(const int); //將一個int型別的變數轉化成大數
BigNum(const char*); //將一個字串型別的變數轉化為大數
BigNum(const BigNum &); //拷貝建構函式
BigNum &operator=(const BigNum &); //過載賦值運算子,大數之間進行賦值運算
friend istream& operator>>(istream&,BigNum&); //過載輸入運算子
friend ostream& operator<<(ostream&,BigNum&); //過載輸出運算子
BigNum operator+(const BigNum &)const; //過載加法運算子,兩個大數之間的相加運算
BigNum operator-(const BigNum &)const; //過載減法運算子,兩個大數之間的相減運算
BigNum operator*(const BigNum &)const; //過載乘法運算子,兩個大數之間的相乘運算
BigNum operator/(const int &)const; //過載除法運算子,大數對一個整數進行相除
BigNum operator^(const int &)const; //大數的n次方運算
int operator%(const int &)const; //大數對一個int型別的變數進行取模運算
bool operator>(const BigNum &T)const; //大數和另一個大數的大小比較
bool operator>(const int &t)const; //大數和一個int型別的變數的大小比較
void print(); //輸出大數
};
BigNum::BigNum(const int b) //將一個int型別的變數轉化為大數
{
int c,d=b;
len=0;
memset(a,0,sizeof(a));
while(d>MAXN)
{
c=d-(d/(MAXN+1))*(MAXN+1);
d=d/(MAXN+1);
a[len++]=c;
}
a[len++]=d;
}
BigNum::BigNum(const char *s) //將一個字串型別的變數轉化為大數
{
int t,k,index,L,i;
memset(a,0,sizeof(a));
L=strlen(s);
len=L/DLEN;
if(L%DLEN)len++;
index=0;
for(i=L-1;i>=0;i-=DLEN)
{
t=0;
k=i-DLEN+1;
if(k<0)k=0;
for(int j=k;j<=i;j++)
t=t*10+s[j]-'0';
a[index++]=t;
}
}
BigNum::BigNum(const BigNum &T):len(T.len) //拷貝建構函式
{
int i;
memset(a,0,sizeof(a));
for(i=0;i<len;i++)
a[i]=T.a[i];
}
BigNum & BigNum::operator=(const BigNum &n) //過載賦值運算子,大數之間賦值運算
{
int i;
len=n.len;
memset(a,0,sizeof(a));
for(i=0;i<len;i++)
a[i]=n.a[i];
return *this;
}
istream& operator>>(istream &in,BigNum &b)
{
char ch[MAXSIZE*4];
int i=-1;
in>>ch;
int L=strlen(ch);
int count=0,sum=0;
for(i=L-1;i>=0;)
{
sum=0;
int t=1;
for(int j=0;j<4&&i>=0;j++,i--,t*=10)
{
sum+=(ch[i]-'0')*t;
}
b.a[count]=sum;
count++;
}
b.len=count++;
return in;
}
ostream& operator<<(ostream& out,BigNum& b) //過載輸出運算子
{
int i;
cout<<b.a[b.len-1];
for(i=b.len-2;i>=0;i--)
{
printf("%04d",b.a[i]);
}
return out;
}
BigNum BigNum::operator+(const BigNum &T)const //兩個大數之間的相加運算
{
BigNum t(*this);
int i,big;
big=T.len>len?T.len:len;
for(i=0;i<big;i++)
{
t.a[i]+=T.a[i];
if(t.a[i]>MAXN)
{
t.a[i+1]++;
t.a[i]-=MAXN+1;
}
}
if(t.a[big]!=0)
t.len=big+1;
else t.len=big;
return t;
}
BigNum BigNum::operator-(const BigNum &T)const //兩個大數之間的相減運算
{
int i,j,big;
bool flag;
BigNum t1,t2;
if(*this>T)
{
t1=*this;
t2=T;
flag=0;
}
else
{
t1=T;
t2=*this;
flag=1;
}
big=t1. len;
for(i=0;i<big;i++)
{
if(t1.a[i]<t2.a[i])
{
j=i+1;
while(t1.a[j]==0)
j++;
t1.a[j--]--;
while(j>i)
t1.a[j--]+=MAXN;
t1.a[i]+=MAXN+1-t2.a[i];
}
else t1.a[i]-=t2.a[i];
}
t1.len=big;
while(t1.a[len-1]==0 && t1.len>1)
{
t1.len--;
big--;
}
if(flag)
t1.a[big-1]=0-t1.a[big-1];
return t1;
}
BigNum BigNum::operator*(const BigNum &T)const //兩個大數之間的相乘
{
BigNum ret;
int i,j,up;
int temp,temp1;
for(i=0;i<len;i++)
{
up=0;
for(j=0;j<T. len;j++)
{
temp=a[i]*T.a[j]+ret. a[i+j]+up;
if(temp>MAXN)
{
temp1=temp-temp/(MAXN+1)*(MAXN+1);
up=temp/(MAXN+1);
ret.a[i+j]=temp1;
}
else
{
up=0;
ret.a[i+j]=temp;
}
}
if(up!=0)
ret.a[i+j]=up;
}
ret.len=i+j;
while(ret.a[ret.len-1]==0 && ret.len>1)ret.len--;
return ret;
}
BigNum BigNum::operator/(const int &b)const //大數對一個整數進行相除運算
{
BigNum ret;
int i,down=0;
for(i=len-1;i>=0;i--)
{
ret.a[i]=(a[i]+down*(MAXN+1))/b;
down=a[i]+down*(MAXN+1)-ret.a[i]*b;
}
ret.len=len;
while(ret.a[ret.len-1]==0 && ret.len>1)
ret.len--;
return ret;
}
int BigNum::operator%(const int &b)const //大數對一個 int型別的變數進行取模
{
int i,d=0;
for(i=len-1;i>=0;i--)
d=((d*(MAXN+1))%b+a[i])%b;
return d;
}
BigNum BigNum::operator^(const int &n)const //大數的n次方運算
{
BigNum t,ret(1);
int i;
if(n<0)exit(-1);
if(n==0)return 1;
if(n==1)return *this;
int m=n;
while(m>1)
{
t=*this;
for(i=1;(i<<1)<=m;i<<=1)
t=t*t;
m-=i;
ret=ret*t;
if(m==1)ret=ret*(*this);
}
return ret;
}
bool BigNum::operator>(const BigNum &T)const //大數和另一個大數的大小比較
{
int ln;
if(len>T.len)return true;
else if(len==T.len)
{
ln=len-1;
while(a[ln]==T.a[ln]&&ln>=0)
ln--;
if(ln>=0 && a[ln]>T.a[ln])
return true;
else
return false;
}
else
return false;
}
bool BigNum::operator>(const int &t)const //大數和一個int型別的變數的大小比較
{
BigNum b(t);
return *this>b;
}
void BigNum::print() //輸出大數
{
int i;
printf("%d",a[len-1]);
for(i=len-2;i>=0;i--)
printf("%04d",a[i]);
printf("\n");
}
BigNum f[110];//卡特蘭數
int main()
{
f[0]=1;
for(int i=1;i<=100;i++)
f[i]=f[i-1]*(4*i-2)/(i+1); //卡特蘭數遞推式
int n;
while(scanf("%d",&n)==1)
{
if(n==-1)break;
f[n].print();
}
return 0;
}
相關推薦
卡特蘭數 性質及應用
卡塔蘭數是組合數學中一個常在各種計數問題中出現的數列。以比利時的數學家歐仁·查理·卡塔蘭(1814–1894)命名。歷史上,清代數學家明安圖(1692年-1763年)在其《割圜密率捷法》最早用到“卡塔蘭數”,遠遠早於卡塔蘭[1][2][3]。有中國學者建議將此數命名為“明安
卡特蘭數(Catalan)及其應用
入棧 一個 其中 無法 數列 選擇 每天 編號 匹配 卡特蘭數 大佬博客https://blog.csdn.net/doc_sgl/article/details/8880468 卡特蘭數是組合數學中一個常出現在各種計數問題中出現的數列。 卡特蘭數前幾項為 : C0=
卡特蘭數catalan證明及應用舉例
卡塔蘭數是組合數學中一個常在各種計數問題中出現的數列。其計算公式是 Cn=Cn2nn+1=(2n)!(n+1)!n!,n為自然數 C n
卡特蘭數的性質及其應用擴充套件
問題描述:卡塔蘭數,是組合數學中一個常出現在各種計數問題中出現的數列。輸入一個整數n,計算h(n)。其遞迴式如下:h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(
Catalan(卡特蘭)數及定理的簡要證明------附上簡要程式碼
Catalan數很重要, 學計算機的, 沒有不知道這個的, 我這個非計算機專業的學生, 也來湊湊熱鬧: catalan數和上述定理的應用非常普遍, 也是很多IT公司筆試面試的常考點之一, 其變換方式層出不窮, 有興趣的朋友可以百
卡特蘭數及括號正確匹配個數問題解釋
(7)圓桌握手問題: 圓桌周圍有 2n個人,他們兩兩握手,但沒有交叉的方案數。 (8)本文的“n對括號正確匹配組成的字串數”問題。 雖然同屬於卡特蘭數問題,但是用同樣的方法解釋似乎行不通(看百度百科上的解釋:關於出棧的問題。還是比較容易理解的,但是我卻不能通過它理解本題) 最後發現是《程式設計
卡特蘭數在多種問題下的應用 組合數學-Catalan數
4)排隊順序問題1】長度為2n的Dyck words有多少種?(Dyck words是由n個X和n個Y組成的字串,其有一個特點:從左往右,對X和Y分別計數,Y的個數始終不大於X的個數。)轉化思路:可以把X看作入棧,Y看作出棧,Y的個數始終不大於X的個數這一性質正好和空棧無法再出棧相一致,所以Dyck word
一道面試題到卡特蘭數及其應用
繪畫展覽門票每張5元,如果有2n個人排隊購票,每人一張,並且其中一半人恰有5元錢,另一半人恰有10元錢,而票房無零錢可找,那麼如何將這2n個人排成一列,順次購票,使得不至於因票房無零錢可找而耽誤時間,應該採用什麼演算法解決呢?(分支限界法) 這個是卡特
卡特蘭數及其應用
一道棧的題目 最近做一道有關棧的題目,題目是這樣的: 若一序列進棧順序為e1,e2,e3,e4,e5,問存在多少種可能的出棧序列? 這道題用遞推可以算出來,假設第k個數進棧,設f(k)表示k個數的總出棧序列數,則它前面的k-1個數有f(k-1)種出棧
關於卡特蘭數及典型例題
關於卡特蘭數:f[0] = 1, f[1] = 1; for(int i = 2; i <= n; i++) for(int j = 0; j < i; j++) f[i] += f[j] * f[i-j-1]; 關於乘法原理和加法原理:1
卡特蘭數,程式實現,遞迴,迴圈,BST和出入棧順序的應用
卡特蘭數是組合數學中的一種數列,它的來歷和重要性可以自行百度,我主要說它的特徵和程式設計實現。 前幾項是1, 1, 2, 5, 14, 42, 132……, 如果令h(0)=h(1)=1,那麼h(n)=h(0)*h(n-1)+h(1)*h(n-2) + .
卡特蘭數應用--n個元素的出棧順序與從(0,0)到(n,n)不穿過對角線的方法數
1.出棧順序方法數: hdoj1023 求出棧序列,比如1,2,3,出棧序列為3 2 1,1 2 3,1 3 2,2 1 3,2 3 1,一共5種 第一種思路: 我們把入棧看做1,出棧看做0,那麼入棧
卡特蘭數(catalan數)總結 (卡特蘭大數、卡特蘭大數取模、卡特蘭數應用)
本文講解卡特蘭數的各種遞推公式,以及卡特蘭數、卡特蘭大數、卡特蘭大數取模的程式碼實現,最後再順帶提一下卡特蘭數的幾個應用。 什麼是卡特蘭數呢?卡特蘭數無非是一組有著某種規律的序列。重要的是它的應用。
卡特蘭數的應用(易懂版)
不符合要求的數的特徵是由左而右掃描時,必然在某一奇數位2m+1位上首先出現m+1個0的累計數和m個1的累計數,此後的2(n-m)-1位上有n-m個 1和n-m-1個0。如若把後面這2(n-m)-1位上的0和1互換,使之成為n-m個0和n-m-1個1,結果得1個由n+1個0和n-1個1組成的2n位數,即一個不
(轉載)Catalan數——卡特蘭數
出現 註意 城市 ads 大於 編號 只有一個 導致 一個點 Catalan數——卡特蘭數 今天阿裏淘寶筆試中碰到兩道組合數學題,感覺非常親切,但是筆試中失蹤推導不出來後來查了下,原來是Catalan數。悲劇啊,現在整理一下一、Catalan數的定義令h(1)=1,Cata
卡特蘭數
ini bits clas cnblogs 操作 div esp class 序列 卡特蘭數是組合數學 常見的數列 主要有4中形式: 1: h(n)= C 2n n /(n+1) 2: h(n)= C 2n n - C 2n n-1 3: h(n)= h(n
HDU 1133 Buy the Ticket 卡特蘭數
i++ ava () pos str mat bre == ann 設50元的人為+1 100元的人為-1 滿足前隨意k個人的和大於等於0 卡特蘭數 C(n+m, m)-C(n+m, m+1)*n!*m! import java.math.*; import java
HDU 1134 Game of Connections(卡特蘭數)
cut res ras sam eof cpp ont des tel 題目代號:HDU 1134 題目鏈接:http://acm.hdu.edu.cn/showproblem.php?pid=1134 Game of Connections Time Limit: 200
bzoj2822[AHOI2012]樹屋階梯(卡特蘭數)
n+1 amp nbsp put mat pan cat limit 一個 2822: [AHOI2012]樹屋階梯 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 879 Solved: 513[Submit][Sta
【BZOJ 2822】[AHOI2012]樹屋階梯 卡特蘭數+高精
div cnblogs operator line code clu while pan .... 這道題隨便弄幾個數就發現是卡特蘭數然而為什麽是呢? 我們發現我們在增加一列時,如果這一個東西(那一列)他就一格,那麽就是上一次的方案數,並沒有任何改變,他占滿了也是,然後他