QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5191|回复: 6
打印 上一主题 下一主题

一些组合函数

[复制链接]
字体大小: 正常 放大

15

主题

4

听众

113

积分

升级  6.5%

  • TA的每日心情
    开心
    2012-1-13 11:49
  • 签到天数: 9 天

    [LV.3]偶尔看看II

    跳转到指定楼层
    1#
    发表于 2012-1-12 15:56 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    本帖最后由 lilianjie1 于 2012-1-12 18:09 编辑
    ! M/ s  o6 a6 ^0 D  Y, y
    3 p# u- O  q) K3 L) ^$ Y& Rn:=12;n;, K+ Q6 `4 U% A; l2 M" O% K/ I
    Factorial(n);求阶乘
    3 k2 M# e( Z5 v1 X5 ~% @Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
    2 c4 h, w7 h; c: M- p( \( F  ?NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);5 ]" ?2 [4 y& Y( H* ?8 H3 r
    NumberOfPermutations(n, 4);
    " h: q; A0 Y/ m* N0 x% kNumberOfPermutations(n, 11);- ]; ^3 w6 s. Z  s
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;2 d( q2 [/ s; p
    Binomial(n, 3) ;
    : _1 i! _4 b/ B1 t* \  M& JBinomial(n, 9) ;
    8 {% z8 W: W1 a' L2 Y; |" kBinomial(n, 10) ;+ ^& r9 y2 \# T5 v( `3 q/ [
    Binomial(n, 11) ;" l! a# J; T$ t
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    , O- s6 F3 N4 d
    - ]( }, p" E/ I, V& P( LFibonacci(n);斐波数Fibonacci(n-1);
    9 M* [* X- y7 Q' [' J8 R- MFibonacci(n+1);2 R. P! O- d: q2 S, {% A1 t" i
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;# u  U' u% Y. s- S7 \: V' E3 W
    GeneralizedFibonacciNumber(0, 1, n) ;. q/ ]# B6 f& N( @; B
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))
    ( g& i& B: ~; ]* _k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;6 q* j( y' G: {% T0 _
    Catalan(1);
    3 H  ?6 U6 o# `8 T5 Q9 J& cCatalan(2);
    & f; B! U' g% ?  M" i5 gCatalan(3);Catalan(4);
    ) j0 J/ A! b5 \1 kCatalan(12);
    # J- I# K* u4 U8 R# V
    , h) n! @/ U6 l. |* z& w. ULucas(n);卢卡斯数
    1 B0 m4 V  F; i; w2 w+ [128 ~( v# o1 h' F& y
    4790016008 \- d3 e) ^/ x: T# {  ~7 ]& s+ Z7 o
    831600% Z# o0 {. ?5 d9 p( P
    12
    , \: @9 r3 {/ Z1 M$ a1320 H  v# w: b: B9 M
    11880
    9 P& @% P( W0 l  T479001600/ h/ k5 u( K1 J% h) j
    12, C5 c( _3 y3 r3 i( `# e8 A
    66
    2 z" ^9 I0 p1 M* C$ R% B" D) `2206 F- I3 B: s) N, ~
    2207 Y# C9 z3 E) `0 l! S
    66
    3 n% {. E  y$ H' D7 t) P. S3 d9 @4 z12
    * I5 R$ @& {# y; d* S  e% ?  h7 O831600
    $ ^( r3 I- u( [0 ^' Z144" P' z! |, v# @& M8 H
    89+ K5 i4 G' j) F
    233+ F' R1 Q" Y9 d/ V9 X7 k8 D3 f. ]
    233* `8 h8 _& d: v! C) X
    610, _* {) P8 V+ r; b- V
    144
    7 p# S! \; j+ f208012% \1 b. `1 @7 ?+ i  P6 e
    208012) `0 I' _9 O: w9 K
    1% m8 V) ?+ e: |) v
    2" V9 P2 u! \; J+ ?
    57 ~, @4 C! n; n4 V/ ~' J; y
    14# r4 e  I9 k- {) G
    2080121 v# ~6 E$ o! b/ n2 ~7 y3 S4 m$ A
    322) j" I0 H9 ~: v4 m, R& N

    " h( ]5 N  f$ `- Y' L0 `4 {( N! b+ k5 ^! B
    卡特兰数=(2n)!/((n)!*(n+1)!))0 y$ m/ W6 [+ M% |( E4 C
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: / d3 Q; g! p* J2 N7 n7 P
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY* x# E. E/ e$ M% K
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: & h# {: W% o' @0 ]) S
    ((())) ()(()) ()()() (())() (()())$ {/ ?5 g0 p  _2 a; Y
    Cn表示有n+1个叶子的二叉树的个数。 ; n& R7 i+ i$ o! [7 \8 j" R9 c
    $ n1 k  X. Q5 j; _4 x1 ^
    Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    ! W5 G3 M$ N9 ^/ \* B: }2 ?0 ]Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    6 l& N0 a  z) [) o: s! M) H6 T9 ^( g- h7 ]. a5 [% `
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。 , M( w$ _5 ^9 d& x
    Cn表示集合{1, ..., n}的不交叉划分的个数. 那么, Cn 永远不大于第n项贝尔数. Cn也表示集合{1, ..., 2n}的不交叉划分的个数,其中每个段落的长度为2。综合这两个结论,可以用数学归纳法证明 that all of the free cumulants of degree more than 2 of the Wigner semicircle law are zero. This law is important in free probability theory and the theory of random matrices. . h2 A8 A' M( h6 p1 i
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:
    5 z  \- k, f0 E' h" y  y8 e! B) u

    " C" \+ K9 E; F. a2 }' H! T* |
    4 g3 D0 {/ `9 h, v( C8 `, Z
    7 q$ M' Q. l$ g( ^# b% V卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。' P. w6 P5 ~# O
    7 c! }* r2 i! I5 y$ q/ U) ]4 N2 N
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    8 E1 o, G  o; s) C$ C4 ]
    4 O: {7 Q  }# J6 L
    1 g  p- p8 P" U8 s# h: ~2 s* ~# f
    n:=100;n;* A% A! K2 h) j$ f* P, H
    a:=Lucas(n);a;
    3 f2 v0 g. e$ l: p) d0 ]b:=Fibonacci(n+1)+Fibonacci(n-1);b;. l7 G3 C7 m( w) p
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    ) r% n* W" L  ^' H# ]3 `- F( p7 ]4 P0 Y, ~  x% E3 s
    100
    2 L  R1 `' x  X3 r! I. ?4 u792070839848372253127" g8 Q: D+ e* h6 N7 V/ x
    792070839848372253127. @: ]/ E2 m9 P. m$ a
    1771124240896309575375& s" x8 x1 G7 S+ b. s' l
    1771124240896309575375

    12.JPG (43 KB, 下载次数: 428)

    12.JPG

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    xxgzftj        

    4

    主题

    4

    听众

    278

    积分

    升级  89%

  • TA的每日心情
    奋斗
    2012-11-16 20:15
  • 签到天数: 57 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

  • TA的每日心情
    开心
    2012-1-13 11:05
  • 签到天数: 15 天

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑   O; k2 P. A1 d

    - B1 i$ _% S' b, m拆分 。。。。强!
    1 s7 S  S+ }! I. l( M% Q7 v0 R1 T% M" t3 v' a
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;+ i9 `3 \9 @7 F% Z3 g, L) F

    * R- i6 ^8 g7 ^% p# @# o7- I1 \) W9 H* O
    190569292
    + q3 |2 }+ J0 A. ~! \. x. a/ Q[. g4 G1 L9 f* k9 F: j' d
        [ 10 ],
    ; }7 y! L* G5 B1 M& G    [ 9, 1 ],  y  ^) L$ M2 j; V
        [ 8, 2 ],5 K2 o3 [; N/ I  y
        [ 8, 1, 1 ],
    , L. u$ X1 o; c) {    [ 7, 3 ],1 _: T7 L6 S) e/ j. Z" k
        [ 7, 2, 1 ],( [, u, H; i* j1 a" h0 p
        [ 7, 1, 1, 1 ],
    4 p8 V: W8 B, ~/ `% ?9 ^    [ 6, 4 ],8 p, U" `7 s" E
        [ 6, 3, 1 ],; r0 T" ?# x% s. W8 S3 S* D" D
        [ 6, 2, 2 ],
    ; Y& I* @! W8 {1 L0 K' k6 @7 ~# }% v    [ 6, 2, 1, 1 ],
    - F6 x1 G0 ]2 k, b, T+ m    [ 6, 1, 1, 1, 1 ],: B% Y3 _, j& Q* }9 n' j  q4 T
        [ 5, 5 ],
    * o6 f, y3 v/ L& U1 M    [ 5, 4, 1 ],( y) |/ T, ~: l( _- Q6 z- x
        [ 5, 3, 2 ],
    0 O' l- V/ x& M# j2 }5 ]    [ 5, 3, 1, 1 ],% l0 v+ S# F1 C6 B& V. H  K/ z
        [ 5, 2, 2, 1 ],  T9 J. E/ {6 E# P- P- D3 H
        [ 5, 2, 1, 1, 1 ],
    : y/ V% ?3 P$ n3 Q# i. H# O    [ 5, 1, 1, 1, 1, 1 ],# x# R$ J6 o8 ?; c; [- S+ j
        [ 4, 4, 2 ],
    - T3 K2 w  ^8 Q2 L, }    [ 4, 4, 1, 1 ],
    ( A+ @6 W9 S: {: ?  Z    [ 4, 3, 3 ],
    % p+ ~" V2 A. m" K3 s. N    [ 4, 3, 2, 1 ],
    4 W0 X) G* W: j) p7 [    [ 4, 3, 1, 1, 1 ],: M8 v0 |6 r) b  ~+ c  ]
        [ 4, 2, 2, 2 ],
    ) p: `7 B: t3 s4 m# u/ Y' M    [ 4, 2, 2, 1, 1 ],2 ]+ x% T9 ~+ u+ e, V# e( M' n
        [ 4, 2, 1, 1, 1, 1 ],: E0 L- G/ b! W2 |* J: K+ _  o
        [ 4, 1, 1, 1, 1, 1, 1 ],
    # J  k! F1 N/ h    [ 3, 3, 3, 1 ],
    / d+ k' F: w6 Y: a    [ 3, 3, 2, 2 ],
    $ m/ b1 H# {- t    [ 3, 3, 2, 1, 1 ],
    5 O  n/ y" ?( v9 X5 `7 o3 F" j    [ 3, 3, 1, 1, 1, 1 ],
    8 U; g. l2 g) y8 w! s/ h    [ 3, 2, 2, 2, 1 ],5 Y" k: g* O2 R1 A2 j
        [ 3, 2, 2, 1, 1, 1 ],7 D+ ?' x- Y+ k0 Y8 q& ~6 R- G( M
        [ 3, 2, 1, 1, 1, 1, 1 ],1 [, @! O3 L- x6 ?
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    ! j& F1 q/ n6 D" o7 k- Y; V# {    [ 2, 2, 2, 2, 2 ],! z: M/ S( O( Z' F2 r
        [ 2, 2, 2, 2, 1, 1 ],
    8 ?! m1 e$ v; C! D" O    [ 2, 2, 2, 1, 1, 1, 1 ],
    8 M; n! x: g9 U% z/ w6 Y    [ 2, 2, 1, 1, 1, 1, 1, 1 ],
    4 q; _6 D+ C; ]+ H1 Z+ Q0 d    [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],$ O2 W4 B% y  e/ z. t+ R) e' Z
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]9 c  l+ k6 u8 i4 N
    ]
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

  • TA的每日心情
    开心
    2012-1-13 11:05
  • 签到天数: 15 天

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑
    6 l/ u1 ]( k) @* i7 |
    1 n6 V; Y5 F4 [: B1 c( b9 c' x/ k: F

    0 G( D; U- W# X  a* D/ F4 {# a伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    1 w$ h- d! B. E' h, A$ F: T- M' t$ J* x
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关" K. n- K6 p/ L3 s
    回复

    使用道具 举报

    15

    主题

    4

    听众

    113

    积分

    升级  6.5%

  • TA的每日心情
    开心
    2012-1-13 11:49
  • 签到天数: 9 天

    [LV.3]偶尔看看II

    本帖最后由 lilianjie1 于 2012-1-12 20:06 编辑
    & P5 i- L; {1 H8 N5 X5 _- w# G
    ) N. p9 g! e" u- C6 X: f4 in:=5;r:=3;
    3 F, u' G, r3 w* c+ H% mEulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    $ h: Z! n. U( H. E" _  F8 p. n5 iBernoulliPolynomial(n) ;伯努利多项式
    9 d7 j2 A+ q" O: F2 X- _
    ; C, o4 T* k# x& f26( Y+ g7 Z) w! T) b8 X0 l/ k" c
    137/60
    7 v1 S' \7 K' a: D7 A+ D0  g% U- j0 U  v) U
    0.000000000000000000000000000000
    5 U; V; R3 |4 ~0 y  S0 a$.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

    22.JPG (56.17 KB, 下载次数: 382)

    22.JPG

    33.JPG (50.55 KB, 下载次数: 408)

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

  • TA的每日心情
    开心
    2012-1-13 11:05
  • 签到天数: 15 天

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 18:44 编辑
    $ g6 N. H2 r) G4 L8 K% j% d1 c
    7 h- `6 ]! F* N9 ]8 o反费波那西数列反费波那西数列的递归公式如下:
    7 h6 E4 B8 f: M$ z- O2 f
    + K; F$ Q* Q6 _, _  PGn + 2 = Gn − Gn + 1 3 U4 o$ n( u4 t1 l' D" L
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...
    - B/ [4 w9 C* ~7 {& Z1 e
    & |$ d4 _2 e7 X; n9 |5 X6 L即是F2n + 1 = G2n + 1,F2n = − G2n。
    % n1 c4 ?$ ?; j; `+ F+ e7 R& k& A  r  y
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    # L) M  X  m' Z0 `: xStirlingFirst(4,3);! p5 q- f* l' u* t
    StirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);3 q: t& k1 m/ k
    StirlingSecond(4, 3);
    % e* z, c& N, w$ M% [# D2
      {5 X) `- d' y4 X$ s3 ~0 }- m4 d52
    " j( ~% V; y0 Z- e' Y) O: |' T5
    1 s3 a; B0 q+ q: A15# d3 W! B( s$ T6 V1 ]8 x
    -6
    & l8 u- l) S3 f- }11
    " b  h5 Q5 r) H7 k3 n-66 V% u2 x" s2 [# o# {7 ]
    1
    1 k8 W) e8 y/ \% G9 i. B( `" D7
    + c( R3 h. A; a/ {' t+ s6' [' D9 I6 P# w) n
    9 T* E; L, Z5 [, ~0 ?! X
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    6 `+ l5 D5 {6 e% {  A+ N' P; S( e0 ^4 N+ f
    {{a}, {b}, {c}} ' ]9 A" {: X9 u& w6 C
    {{a}, {b, c}} 8 h2 h, C" f/ u( C2 _' i, I
    {{b}, {a, c}} 3 ?0 }& a* S& H0 A9 M. k
    {{c}, {a, b}} $ k  D% T- g( W+ h
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s8 \' _; G. q! L
    s(n,k)是递降阶乘多项式的系数
    % i. e; G" g. Q4 `6 f9 J有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)  P- C2 A$ M1 \  F: n% W+ ]

    ! ?0 j' }4 X3 L3 O- S. ?换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):
    8 j9 e0 p; v* d9 y) S7 h) t
    / M+ l8 S1 }/ v- A$ K- J{A,B},{C,D}
    1 @; C2 F* B' l8 U5 w* s( C# L/ u{A,C},{B,D} ; ]$ h# W) u. I& e7 f3 r3 V3 t
    {A,D},{B,C}
    % M4 O/ U0 ~$ R! I+ [# U{A},{B,C,D}
    8 L+ u, e3 {$ E: m) ^1 ~{A},{B,D,C} 7 w6 s, c1 J7 L) `0 B
    {B},{A,C,D} % R, Y9 d- S+ ~! X: G3 D
    {B},{A,D,C} / A# R% s& t- c. O: |# O
    {C},{A,B,D}
    ' o* ]) z6 U! [3 E{C},{A,D,B}   z' J1 k9 V! A$ a: j3 @- h6 D" L) M
    {D},{A,B,C}
    % s- A2 S- V' @* o% |{D},{A,C,B}

    0 V! X6 c, `( U0 E0 b: o7 D/ @$ {第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S$ e* u% Y/ M& j' y- W1 n+ O
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) " h4 Z9 T2 O: I% \. b% M/ G# _
    S(n,n − 1) = C(n,2) = n(n − 1) / 2 ; q3 `: k* P$ _
    S(n,2) = 2n − 1 − 1 , c  t- V2 ?* `6 B% ]9 K
    # c/ |* q3 r( o" x1 }, g2 e/ c7 J
    换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
    , O: v7 R2 S+ A  \/ x8 m) K$ g( k) Q
    {A,B},{C,D} ' {! a3 k4 y$ k
    {A,C},{B,D}
    ' r6 Q! ~+ Z6 P7 ]{A,D},{B,C}
    : W6 o' P7 Y" d7 u& r) m+ s{A},{B,C,D} ! j$ M8 |  f" Z' f
    {B},{A,C,D}
    1 R  o$ r: o* `; P7 J{C},{A,B,D} + r+ @) N. m6 F0 j# n* A) i
    {D},{A,B,C} 3 A) ]( c! w1 P  z+ N# i  i4 ?
    因此S(4,2) = 7。
    回复

    使用道具 举报

    74

    主题

    6

    听众

    3302

    积分

    升级  43.4%

  • TA的每日心情
    无聊
    2015-9-4 00:52
  • 签到天数: 374 天

    [LV.9]以坛为家II

    社区QQ达人 邮箱绑定达人 发帖功臣 最具活力勋章

    群组数学建摸协会

    群组Matlab讨论组

    群组小草的客厅

    群组数学建模

    群组LINGO

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-12 21:47 , Processed in 0.645337 second(s), 94 queries .

    回顶部