1. 程式人生 > >卡特蘭數 性質及應用

卡特蘭數 性質及應用

卡塔蘭數組合數學中一個常在各種計數問題中出現的數列。以比利時的數學家歐仁·查理·卡塔蘭18141894)命名。歷史上,清代數學家明安圖(1692年-1763年)在其《割圜密率捷法》最早用到“卡塔蘭數”,遠遠早於卡塔蘭[1][2][3]。有中國學者建議將此數命名為“明安圖數”或“明安圖-卡塔蘭數”[4]

卡塔蘭數的一般項公式為C_n = \frac{1}{n+1}{2n \choose n} = \frac{(2n)!}{(n+1)!n!}

前20項為(OEIS中的數列A000108):1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900, 2674440, 9694845, 35357670, 129644790, 477638700, 1767263190

性質[編輯]

Cn的另一個表達形式為C_n = {2n\choose n} - {2n\choose n+1} \quad\mbox{ for }n\ge 1所以,Cn是一個自然數;這一點在先前的通項公式中並不顯而易見。這個表達形式也是André對前一公式證明的基礎。(見下文的第二個證明。)

遞推關係

C_0 = 1 \quad \mbox{and} \quad C_{n+1}=\sum_{i=0}^{n}C_i\,C_{n-i}\quad\mbox{for }n\ge 0.

它也滿足

C_0 = 1 \quad \mbox{and} \quad C_{n+1}=\frac{2(2n+1)}{n+2}C_n,

這提供了一個更快速的方法來計算卡塔蘭數。

卡塔蘭數的漸近增長為

C_n \sim \frac{4^n}{n^{3/2}\sqrt{\pi}}

它的含義是當n → ∞時,左式除以右式的商趨向於1。(這可以用n!的斯特靈公式來證明。)

所有的奇卡塔蘭數Cn都滿足n=2^k-1。所有其他的卡塔蘭數都是偶數。

應用[編輯]

組合數學中有非常多的組合結構可以用卡塔蘭數來計數。在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:
XXXYYY XYXXYY XYXYXY XXYYXY XXYXYY
  • 將上例的X換成左括號,Y換成右括號,Cn表示所有包含n組括號的合法運算式的個數:
((())) ()(()) ()()() (())() (()())
  • Cn表示有n個節點組成不同構二叉樹的方案數。下圖中,n等於3,圓形表示節點,月牙形表示什麼都沒有。
  • Cn表示有2n+1個節點組成不同構滿二叉樹(full binary tree)的方案數。下圖中,n
    等於3,圓形表示內部節點,月牙形表示外部節點。本質同上。
Catalan number binary tree example.png

證明:

令1表示進棧,0表示出棧,則可轉化為求一個2n位、含n個1、n個0的二進位制數,滿足從左往右掃描到任意一位時,經過的0數不多於1數。顯然含n個1、n個0的2n位二進位制數共有{2n \choose n}個,下面考慮不滿足要求的數目。

考慮一個含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的二進位制數。反之亦然(相似的思路證明兩者一一對應)。

從而C_n = {2n \choose n} - {2n \choose n + 1} = \frac{1}{n+1}{2n \choose n}。證畢。

  • Cn表示所有在n×n格點中不越過對角線的單調路徑的個數。一個單調路徑從格點左下角出發,在格點右上角結束,每一步均為向上或向右。計算這種路徑的個數等價於計算Dyck word的個數:X代表“向右”,Y代表“向上”。下圖為n= 4的情況:
Catalan number 4x4 grid example.svg
  • Cn表示通過連結頂點而將n + 2邊的凸多邊形分成三角形的方法個數。下圖中為n= 4的情況:
Catalan-Hexagons-example.svg
  • Cn表示對{1, ...,n}依序進出置換個數。一個置換w是依序進出棧的當S(w) = (1, ...,n),其中Sw)遞迴定義如下:令w=unv,其中nw的最大元素,uv為更短的數列;再令S(w) =S(u)S(v)n,其中S為所有含一個元素的數列的單位元。
  • Cn表示用n個長方形填充一個高度為n的階梯狀圖形的方法個數。下圖為n = 4的情況:
Catalan stairsteps 4.svg
  • 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 .... 這道題隨便弄幾個數就發現是卡特蘭數然而為什麽是呢? 我們發現我們在增加一列時,如果這一個東西(那一列)他就一格,那麽就是上一次的方案數,並沒有任何改變,他占滿了也是,然後他