QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5277|回复: 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 编辑
    : C! B0 N* Q) a4 S" j( U6 m, f" A/ b, n. v6 j/ k. c
    n:=12;n;
    7 g3 `$ _, ~7 J4 ^8 l7 Y' SFactorial(n);求阶乘1 G1 J- J5 C9 l% Q
    Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
    9 p. y+ e5 C- U& ?1 J, p% PNumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    # N& p$ ^+ F7 aNumberOfPermutations(n, 4);
    / ?, q9 y( l  x! ?3 I- ~NumberOfPermutations(n, 11);: \8 z) H; F) ?+ `) C4 i: ?5 c
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;; z& u( Y4 N# h/ L' D) ?0 K) P. _
    Binomial(n, 3) ;/ }$ {4 K" z2 K! m2 X3 o
    Binomial(n, 9) ;
    , _4 a# A+ v+ V  D* f0 ?Binomial(n, 10) ;' W! \2 M; `) U1 J2 S
    Binomial(n, 11) ;  z! b) p' Q' K7 o# B" h
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    . X5 L9 ^9 e0 O0 U) j2 {- z" I- ]1 |$ _% t. M( V
    Fibonacci(n);斐波数Fibonacci(n-1);* {, a* ^9 Q; |3 N
    Fibonacci(n+1);" R/ O3 M4 Z) o- J4 H  W
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;6 W* S4 n/ o: }5 z  B$ h
    GeneralizedFibonacciNumber(0, 1, n) ;
    ; A6 s6 v# N. J$ a1 c) g$ iCatalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))
    # L3 ~: O# l) O; c5 Fk:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    , e( ^$ ?% o' L, u: k8 o8 KCatalan(1);
    ) {3 \& b" Z1 T+ a: wCatalan(2);
    / t! e7 h& n# A3 r% CCatalan(3);Catalan(4);
    + ^7 Q3 x' I# O6 D2 HCatalan(12);
    2 j) M3 G8 V6 B
      C! y1 I% }( _" Y# hLucas(n);卢卡斯数
    * b7 `% t7 C" V! g4 M12  e1 z1 G2 G% _, [/ a! X
    479001600) K) j4 n% u! f
    831600
    ; S7 P( O/ t4 y6 }8 N124 Y9 D# t) b1 Y
    132
    / [% E7 m; N# |  z11880
    - |' a/ E7 ^* Y479001600
      A% I5 Y' L& {6 x( ]# }1 u3 E$ Q0 P12
    ( B: E, X8 E/ r+ E0 g66" ~- i* }3 c2 q  U. ^3 t- w  e
    220
    0 z- G* C" G! [220; G- [2 Z( K+ u; r0 E
    66" i) t  x6 |* s5 w9 o' z# K+ X( Z/ l
    12
    ( i/ v! Y, Y. i9 p831600  c$ l& E$ K+ F
    144
    / e$ w6 z& m9 y89
    " k- D6 b* E7 z; {233: l+ \% C0 X: g, b* w4 a* \
    2338 J6 p0 R, |' o( _  K. f1 F
    610
    / R$ q% f( A) @' |5 Y8 G144
    + x. c0 x9 f5 z* t8 R208012
    1 m$ W6 s# ?" @) |: L8 g  d2080128 i7 V5 e9 w) V( n+ e& A
    1
    + V5 S: c+ R' c  F( |9 ^2
    $ X9 G  Z. z- Y. L6 c' @+ a4 z5- b$ }  [! k* G0 @% F
    14& t6 C& W7 j, C
    2080128 Y  h- l4 U* t; Q% @6 \
    3229 i) v% A3 B! a8 r" E
    # k8 S2 E% R7 E0 l& |6 o; z
    * Q4 h: I2 }( q
    卡特兰数=(2n)!/((n)!*(n+1)!))/ r$ i/ m8 y2 @2 w  L1 b
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    + P/ s, Q" o: H; `' _+ z, d! T**YYY XYXXYY XYXYXY XXYYXY XXYXYY
    ; E. W! |0 [+ s# W2 h5 r; s1 A+ v将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: , o! M' K+ w" B
    ((())) ()(()) ()()() (())() (()())# {" g, C( G5 v4 \
    Cn表示有n+1个叶子的二叉树的个数。 6 v( V6 X- t9 S: H

    + I. p6 Z5 M" q- a1 w$ wCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) 6 S# |8 x& k& A( e, {5 Z4 f
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    3 Q, h- r1 j4 ^/ v* b' ]. N3 l$ N$ v
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    ) o- K& D) W# g6 h* wCn表示集合{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. ) n: `" k# K8 a9 I6 y& k
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:
    4 ~- A, q: T8 h( h* n0 X7 l! _
    3 c* `4 b5 g" B: y$ W5 [
    3 ~$ p, k% P2 ^4 G: d" {! Z

      u9 U9 b, ^4 j3 [; E3 r卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。
    ) y3 v- e9 a( S0 z
    - Q# S# F) h9 a) ]! h1 F2 ~5 T3 `但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    # i2 F/ m* T3 R7 ~3 O3 o
    1 B. w1 d4 d# g% P) z5 ]" w" A
    % y0 M5 t3 u$ M/ e
    n:=100;n;: c/ r9 R, J1 H0 _( y8 |0 Y+ I: N
    a:=Lucas(n);a;- J/ b# ]7 y) p1 ~2 X
    b:=Fibonacci(n+1)+Fibonacci(n-1);b;/ ]9 w8 E8 D- u. j9 r
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    : g' E- f2 ^. m' I: Y, {
    2 Y- F/ L- B$ Q5 W: v; N; g& ~; D8 F100
    - G; \; T5 q2 ~$ r# A' i! G6 m792070839848372253127
    5 D2 i3 c- l$ F4 K6 x" [792070839848372253127
    / g. l) A) {# i2 j8 V+ n; C1771124240896309575375% V: g. \, }; [$ \  \) j* G; v) h2 X) N
    1771124240896309575375

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

    12.JPG

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

    74

    主题

    6

    听众

    3303

    积分

    升级  43.43%

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

    [LV.9]以坛为家II

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

    群组数学建摸协会

    群组Matlab讨论组

    群组小草的客厅

    群组数学建模

    群组LINGO

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 18:44 编辑 5 @' q$ @- ]+ k4 l) \/ s

    4 l; Q5 |7 {! O3 P8 q) s# u* J& p反费波那西数列反费波那西数列的递归公式如下:  n  B( B6 p) X1 S# O9 p0 s9 [

    / J; e- j" v: W7 r% j5 N& e1 bGn + 2 = Gn − Gn + 1
    ! @1 F  u, {7 x; u$ X如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...6 o6 v) p. D# E. w  D6 t

    ! G3 o6 M6 d3 _5 d1 ], a% Y即是F2n + 1 = G2n + 1,F2n = − G2n。) y) c0 v9 I, M: U- V+ r
    * Q0 j( w7 x* H. V! M, Z
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    3 ?3 h8 L& B1 ~! A+ s- d+ Z7 J& }StirlingFirst(4,3);
    " q! N8 u0 _7 T. D" YStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
    2 Z# W1 Y% n! S. s4 BStirlingSecond(4, 3);
    8 w, h" `6 u" m* b2
    ! ^4 S! b: Q1 _1 w' q52
    ( D0 V3 Y, f: d7 A; p5 }5) C% \) e0 ]3 }( J# s! R
    150 m/ K9 A2 ~: Z( `. d
    -6
    $ w& u$ c9 m* L& d8 L) P112 y0 e2 C( f5 U3 r" L2 m
    -62 o7 _3 W5 a/ ]1 w, i: x9 q, w
    1/ t5 E* a# Y4 ]. R2 Y
    7" e, L9 [# {3 f7 c3 B% }. k
    67 X, q$ v4 M; b5 n7 S' w% h8 ~
    4 f0 `4 e. ~7 L' F7 ?7 @+ x
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    ) @0 t0 k- D3 D  i+ K7 h0 O. m: l( f% I( ~/ t
    {{a}, {b}, {c}} 4 D7 C; A6 P: d  _4 J9 O# F9 F) T
    {{a}, {b, c}} 4 ^. I# H2 E5 Y+ J
    {{b}, {a, c}} ' E1 q  N( d  @2 Y
    {{c}, {a, b}}
    - \* _. G: R/ }{{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s' ]7 I6 W" Z; ~% L" [2 \
    s(n,k)是递降阶乘多项式的系数
    5 n1 F5 ?, n/ u5 k# \5 O4 m有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)
    ' B6 c: j* i8 z. S0 A8 ~' z7 j8 ~8 q: `
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):, w; G: \" u% ~5 D( ^" o

    ) x" a! m! Q+ F6 `" c{A,B},{C,D}
    ( b, [5 a$ }# I* f{A,C},{B,D}
    ( K2 v6 p* K8 i+ J/ B{A,D},{B,C}
    * t5 @" G1 Y/ a: r8 w/ l% V{A},{B,C,D} 2 b* Z4 X* S( W: R5 F
    {A},{B,D,C}
    / t) f; U  Z/ V/ I! K  Y{B},{A,C,D} 7 n6 D' f: D- x& D: t" P0 `& U
    {B},{A,D,C}
    . J& h, j, c$ p{C},{A,B,D}
    . I1 X4 B: U- V! z{C},{A,D,B}
    . [2 o: R1 O6 S9 \" \; }{D},{A,B,C} % Y& i# m# K! ~4 @
    {D},{A,C,B}
    6 B) w/ a. j0 j- r6 F4 J
    第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S+ @+ K0 Z7 v* U6 @. O5 p7 i
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) ; w* h+ Z7 Z- C6 R6 p. }! Y! @4 [
    S(n,n − 1) = C(n,2) = n(n − 1) / 2 + z' Y' {8 [  t) C7 j+ k9 U
    S(n,2) = 2n − 1 − 1
    / ]( ]5 m  {' a8 K6 y, w* s) u2 m5 x2 H$ a5 J
    换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:2 X. K7 c3 d- s# I  [4 u3 f

    . P& k) `* X5 O$ @{A,B},{C,D}
    ) |( j& X' R+ r5 |{A,C},{B,D} / r1 j" T0 U3 C' P: S+ S4 o9 Y
    {A,D},{B,C} ! L( J3 e$ T$ X0 a- b9 d$ @) m  J
    {A},{B,C,D} $ S# _3 u9 C: |0 f; u7 g
    {B},{A,C,D} 4 y3 S, X- \. N7 o5 r7 h
    {C},{A,B,D}
    0 Y& D7 M" ?  |0 J{D},{A,B,C}
    & \( @% x' t, E) E因此S(4,2) = 7。
    回复

    使用道具 举报

    15

    主题

    4

    听众

    113

    积分

    升级  6.5%

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

    [LV.3]偶尔看看II

    本帖最后由 lilianjie1 于 2012-1-12 20:06 编辑
    $ s. ~( s& W% k( D0 ~' B/ K2 i0 J* G5 N* w# Q6 [% D, |) {
    n:=5;r:=3;* m6 q: I3 m, ^5 Q/ `" u: d
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
      U4 H) E0 |  ~8 b% ^: [: V# z1 wBernoulliPolynomial(n) ;伯努利多项式
    2 R; J* f/ R+ d& t8 g. P0 Y
    # r; v- ~) p: ]1 ~26/ D' [8 }" [- c: N. x
    137/60# E5 c7 E' i+ Z
    0/ C. g5 Q) ^( l; y# U& C# C
    0.000000000000000000000000000000
    " |1 U/ Q, [6 l$.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑 + ~6 Z* {7 D: E/ g' K: o( S" D6 t

    / E, a5 u% O- V8 ]( p' G  \" }+ G; q% V6 e! T/ P

    % F& |! J3 P: n* M; ]6 Y伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。+ x) J5 X; G& U1 I% c0 N; O9 S
    ; ^0 {/ ], z. O* R% |3 J3 v8 ~4 g
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
      I. G. I/ X  f' z+ C  s1 P
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 . ]! G1 [, ]7 U) i. f7 m. p* m
    $ u" X. G# d6 w* s6 H7 o3 D, C
    拆分 。。。。强!8 i0 q$ k5 T( p4 Y

    6 k8 v- Z7 S& B: {- JNumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;( {; O6 Z) Y3 W1 i3 h, q
    ) I+ [+ K- Z9 @" C
    7
    9 K; u" k& ], _3 v# K190569292
    0 Q/ o8 J3 v+ C! `0 s[
    ! P1 `6 Y; d% E; ?' _% _" O    [ 10 ],7 \* m  t4 @2 m- g
        [ 9, 1 ],  u; S5 n8 }4 h# I2 d* Y
        [ 8, 2 ],$ y& M& a& y" x2 h
        [ 8, 1, 1 ],
    4 U3 V$ N0 G9 v/ }% @    [ 7, 3 ],* f( }4 C# y- F5 Q# p4 g
        [ 7, 2, 1 ],. h# M% ]+ U% k1 M& ]2 P
        [ 7, 1, 1, 1 ],$ S# }9 v5 {5 x5 w7 N
        [ 6, 4 ],
    $ |  e; y* X* _8 ^2 b    [ 6, 3, 1 ],/ p8 u3 u5 r1 v4 b
        [ 6, 2, 2 ],
    0 t* D8 L' {8 \0 q1 ~    [ 6, 2, 1, 1 ],1 |* K) U: A+ K* H
        [ 6, 1, 1, 1, 1 ],7 y  ]; \$ y% P: |+ O; |
        [ 5, 5 ],
    7 a1 l8 n  f/ Y$ K% z    [ 5, 4, 1 ],, m, t" Q8 V! C( l3 l" s
        [ 5, 3, 2 ],
    # {, _% t, M8 k" y7 M# s    [ 5, 3, 1, 1 ],: J8 Y# \. Z- {+ @
        [ 5, 2, 2, 1 ],2 u% o& t9 l: _8 L
        [ 5, 2, 1, 1, 1 ],7 \" S, w) n. u" F( r' o* D
        [ 5, 1, 1, 1, 1, 1 ],8 s1 G% ?3 c  _, ]
        [ 4, 4, 2 ],) v% F/ h# _& t5 T
        [ 4, 4, 1, 1 ],: ?$ m3 D5 ^" j. n
        [ 4, 3, 3 ],
    2 u+ I9 S3 y4 u! w% U/ t  |+ E/ y    [ 4, 3, 2, 1 ],
    3 a; L7 j6 ^2 P9 `+ ?    [ 4, 3, 1, 1, 1 ],
    ) k( s& ~- X5 e( T" m$ I! p    [ 4, 2, 2, 2 ]," v4 p+ R3 o6 V' f+ d: Z- d, D
        [ 4, 2, 2, 1, 1 ],# ]3 B- n0 p7 j" y8 _4 N% n
        [ 4, 2, 1, 1, 1, 1 ],1 m# O  [# n0 p- @+ m+ \. W
        [ 4, 1, 1, 1, 1, 1, 1 ],& k/ g% x3 \4 m9 v+ p7 O3 g
        [ 3, 3, 3, 1 ],7 p6 u. s' b: n+ ]4 O6 a
        [ 3, 3, 2, 2 ],3 l1 f- G+ ^, o$ l9 Z8 o
        [ 3, 3, 2, 1, 1 ],
    % s$ c3 f& f) t% l/ |# m$ ~    [ 3, 3, 1, 1, 1, 1 ],, u  g7 [; H# h8 L+ I; n! L
        [ 3, 2, 2, 2, 1 ],
    " _8 Y0 d" J) [    [ 3, 2, 2, 1, 1, 1 ],
    2 W5 c( R4 Q% d- o, Q7 M    [ 3, 2, 1, 1, 1, 1, 1 ],& {) {+ O$ y& A1 _" ~, _: L# X
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],7 a2 z4 x/ U3 u* J$ ]
        [ 2, 2, 2, 2, 2 ],- z2 d/ f" C( `* ]; o7 m  `
        [ 2, 2, 2, 2, 1, 1 ],1 H8 G# v0 d3 ~  o# L
        [ 2, 2, 2, 1, 1, 1, 1 ],+ _+ ^1 ~& f* E
        [ 2, 2, 1, 1, 1, 1, 1, 1 ],6 X( k3 x$ \, y0 ^
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],, d- v! {2 _4 x6 T; U6 H
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]' B  _8 F' p4 c; B: ]) h0 `7 c
    ]
    回复

    使用道具 举报

    xxgzftj        

    4

    主题

    4

    听众

    278

    积分

    升级  89%

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

    [LV.5]常住居民I

    回复

    使用道具 举报

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

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2026-8-23 20:31 , Processed in 0.634212 second(s), 92 queries .

    回顶部