QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5174|回复: 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 编辑 8 c1 \: a) F5 b
    7 D; g) Y6 s" k4 o2 j5 r
    n:=12;n;5 {3 W9 j0 K+ ^6 Z4 ]- f
    Factorial(n);求阶乘, n% P: a" P- w( U7 T, A  M
    Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));/ Y' A$ C. e- Q( w; u# t+ ~" ^
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);% K3 J. [7 D6 |2 }  F% y7 E  z
    NumberOfPermutations(n, 4);
    * X* A1 a7 R( L" p& t3 `1 lNumberOfPermutations(n, 11);
    " |1 d* J& n/ l; d& r$ m# F0 |Binomial(n, 1) ;二项式系数Binomial(n, 2) ;& J- J, n& m: ]$ E9 H
    Binomial(n, 3) ;
    - \! {  |/ m& z( h" R3 s$ K1 `* s' k4 PBinomial(n, 9) ;( Q1 }" R$ L# z; q& a2 d
    Binomial(n, 10) ;
    * c$ T7 L1 W  c2 M( o( P& z  BBinomial(n, 11) ;
    $ i5 A& \& V5 V- y1 M: v! V8 X, `Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=8316005 k  j& R9 E& N0 J- \* H

    $ G4 D4 T& c5 |$ ]# SFibonacci(n);斐波数Fibonacci(n-1);( R7 E2 h; g3 Y2 ^4 {
    Fibonacci(n+1);/ e2 e2 Q2 g9 a
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;
    ! F2 K/ S/ H  y4 X- j, E) P% ~GeneralizedFibonacciNumber(0, 1, n) ;
    1 P% N/ T0 ~: [- Q5 k( [Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))# L- W1 f& y5 ^9 X
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    9 q2 U1 D: Q: H' t) \8 Z( M+ h- yCatalan(1);, M/ G! n7 Z  @! ~
    Catalan(2);
    ) m9 I! e* T# tCatalan(3);Catalan(4);
    ! A8 ?" R* C" x: }& tCatalan(12);0 [) f$ J8 P8 C( }
    , G6 G$ p% T( _; \% r" ]5 Y
    Lucas(n);卢卡斯数0 i/ ]; o7 z' ]4 f1 \3 ]: o4 |, H
    12
    # ^+ s7 e, v- i479001600- Y% Q& x5 a0 Z2 [- t6 `
    831600- S$ W3 ]6 z3 c- `! d
    12# v! O. G# C5 A) |2 P
    132) d( H: H) z) [: Z- O: T: ^3 K
    118802 Z& o8 S' D! y
    479001600
    5 \' e) C/ o. R6 v12
    " p; y: t0 E+ ]) x) m# M66
      d% s. f- b, w3 h5 X' e220% a% @4 J: X1 u
    220+ s2 p) X+ Z) j  P. l
    66
    4 M! o0 s# v! t9 E6 \12
    , P# r+ z* o  y8 n1 m" t+ |) g831600
    ' `) ]/ u( D$ z144+ k6 O0 B1 N9 U- \; M9 v$ r* ~6 W, x
    892 I* V  j6 j% M6 q5 V
    233
      E0 k: X; K" D* V5 R% C9 Z: |233  i  F8 \4 {1 i; F
    610" Y( A7 }) B' h# s  b# O
    144  {/ C; t; y/ L% A) e; F
    208012
    ( a# }9 |1 |* e& ~! q2080125 ]$ @) z- g. Q+ k# ~. t2 T0 f. t8 B
    1
    : Q" ?5 S- q7 c- x' Q20 ?2 {7 g- c" I* M2 z- L
    5+ M, C0 u, \& I( M
    14
      J" f# j! R: b0 m5 y  ]4 w8 s* r208012
    : J$ {' z. |! d322
    % V% q7 F& X& [# ~; G0 p) }5 q  b' ^* ~

    6 t" W1 u8 u/ z! d. b1 E( ?卡特兰数=(2n)!/((n)!*(n+1)!))" y, `1 E. b# W; J  B: Y
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    , \, J3 Z- P+ Z' L' a' k9 T" o**YYY XYXXYY XYXYXY XXYYXY XXYXYY
    1 N$ ~+ g: P: e/ ?将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: ; ]  _4 G, }) r# U4 S7 `! ^
    ((())) ()(()) ()()() (())() (()())& h( x5 m0 u' U8 K. s
    Cn表示有n+1个叶子的二叉树的个数。 # b. q9 Y' l5 J4 s# ]( L

    : S- f5 Z- R7 g: Y  o; gCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) - p, D( j0 P3 N# w) V$ Z
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:   {+ X& j  j9 D& F8 P

    ' `! _5 }' J% Z$ {4 OCn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    & v5 U: w0 P  ]. ^1 W: ~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.
    # z" r( z/ r' ?8 n5 i/ V/ ~; `/ c2 fCn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:
    ! x! ?$ z; R, R( @# p. }

    & ^5 g, A& V3 G3 p: F  d6 a& J1 ~) t8 P+ m

    0 O3 Z/ Q4 S, K" i卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。* T& L  G) \% q

      z7 d6 Z! c7 o  F1 H- l( V' Q  c但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

      |0 x) C3 _/ g% `# h) b
    - |, d+ x9 |2 ^0 M, [, J1 [# ^6 \' p* ~' V4 h, N3 C
    ) x. ?* z0 {7 L, o' ^; f) b$ T9 \% |
    n:=100;n;5 o& y, @& S+ y0 K9 y' {  b- i: |
    a:=Lucas(n);a;
    ! ]. }% Q+ S7 r& ]  b! pb:=Fibonacci(n+1)+Fibonacci(n-1);b;' s. \! T7 Y. i2 c- s
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    4 K+ |4 P+ n) }( K
    & o6 F$ T6 D: r. e100
      y  b1 L* v' A: d7 E  [' h+ y792070839848372253127
    2 l7 P. n) ]* {4 v4 L5 d5 o7920708398483722531278 V0 w! O/ u  Z/ w' h  W1 V
    1771124240896309575375- e7 \+ V% U3 [' ]
    1771124240896309575375

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

    12.JPG

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

    74

    主题

    6

    听众

    3302

    积分

    升级  43.4%

  • 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 编辑 ) A6 u* U, u) U) z. ~' K( h7 k
    $ S' O+ c3 W3 P
    反费波那西数列反费波那西数列的递归公式如下:
    8 W9 h/ [1 l( _! O, \5 Q4 v% ?
    3 x' B1 n$ Z% ^Gn + 2 = Gn − Gn + 1
    3 V0 v6 z6 }: K! \$ _: m! ^- I; L如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...
    7 n) U2 j# Z! }, v1 p
    3 i  m9 n' S- y7 i6 a. o即是F2n + 1 = G2n + 1,F2n = − G2n。  p7 n8 ?& a; b  T$ f' |

    - }* [: l' n. a: F! t. [4 ]Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);' {  t- b! ^% T" k) t0 K9 j
    StirlingFirst(4,3);! r! H3 V! Y  N( B) ^
    StirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);  F% p" E8 d) {) M8 m6 U
    StirlingSecond(4, 3);
    " t- f/ K6 o& ^7 ~  `* T! R" W1 T$ N2
    ! G1 d% ^  ]$ \. [2 `& u2 H2 g521 S- Y- q! _- q3 m0 x5 ]
    5
    ( l4 s. d& R! ~7 H8 z6 j9 H; j15
    + l+ r0 _9 M" g-6+ T( O! ]! o. y3 c9 n
    117 R0 o6 l+ D0 S( ~! o, r4 U* n' V
    -6% M6 s5 G: A/ l5 X* c; K
    15 h; d. B  [0 K- [. `2 d
    7
    3 @0 Z4 L# s7 P/ x6
    ' a* W$ C& N( |* f
    " E$ l8 K2 w) e( U$ P, YBn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:3 Z5 L# K& K0 O1 \* L
    8 Z5 e: v6 _5 n" Z0 K
    {{a}, {b}, {c}} ; H5 c% |# ~" Y2 e6 G
    {{a}, {b, c}} ! ?: s. v, B- G7 B
    {{b}, {a, c}} 5 b) x% e5 F+ R9 K+ p# [# O
    {{c}, {a, b}}
    ' z/ I$ U" C. v5 [* q  C% C{{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s# x! j% E# i5 l
    s(n,k)是递降阶乘多项式的系数
    4 [2 M# |( N, g有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)
    2 v- B# D# F+ K) i/ m, F: @- E( Q7 G# w% Y3 ~( d
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):9 t' [' l) ^6 N
    6 ~5 W9 r+ A( \5 A$ ], r3 Y" o
    {A,B},{C,D}
    * W; k1 s% F8 Q$ z4 B. [) N$ c{A,C},{B,D} 1 }: C/ X6 c/ l) [; Z3 c
    {A,D},{B,C}   f/ z0 {1 D' k* N
    {A},{B,C,D}
    # q& n0 q3 b! x/ n0 n. o{A},{B,D,C}
    + D5 c8 x3 R6 P/ k7 c2 N* F{B},{A,C,D} 4 v  \4 p) P& w, g1 J7 }
    {B},{A,D,C} + }( |  h2 X7 A5 q, T
    {C},{A,B,D}
    4 u, A8 C5 a+ t  i{C},{A,D,B} ; {3 n" i1 c+ m" N% s! v
    {D},{A,B,C} ) D8 C' J2 e, Z" a0 U8 L1 C: p
    {D},{A,C,B}

    2 e9 h" v. l8 O$ f# C0 Z' I/ l/ W第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    4 n; |* E2 ~1 s1 \给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k)
    & ]; K: N5 k3 S% M% p( w- KS(n,n − 1) = C(n,2) = n(n − 1) / 2 1 N/ w2 T5 `8 U5 s5 I  E% z
    S(n,2) = 2n − 1 − 1
    8 K" I9 e+ L& f' k3 O8 e
    / y; ]5 O5 d5 {" m# m* H/ O! @$ y换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:( ]4 J7 p* ?# A( Z( a, f
    : o! o1 ~6 [, a5 F+ t
    {A,B},{C,D} + h9 p0 t1 }$ d
    {A,C},{B,D}
    $ B1 t" t7 J+ L5 r7 P{A,D},{B,C} $ O0 Y& X5 M3 n* G6 \
    {A},{B,C,D}
    * C% W+ M" A0 o4 f0 r- J: b# C{B},{A,C,D} 2 ]$ `% ?1 m* p4 q: v7 J
    {C},{A,B,D} / |8 b1 r$ o3 \2 D$ R' D
    {D},{A,B,C}
    % I& M6 Y  _7 ^因此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 编辑 4 U9 K+ ^5 }  W* z$ i

    ( v0 n8 H9 y: E: I8 ^9 |$ |5 Yn:=5;r:=3;& _# ~# P2 u4 [4 P& }& V$ K
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    ( x7 j1 k' V. U# _9 X% eBernoulliPolynomial(n) ;伯努利多项式
    * H8 S& S/ O0 x0 `6 A( y: ?, G0 Y9 e" Z, j! G3 y
    26. Q; j/ N& p0 s7 y
    137/60
    + L8 J  I0 H# ]' w) u* H0
    * H  j  Y- I' w0.000000000000000000000000000000
      g- G) ^" C4 ~" M, X$.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑 ! e! p8 k* H) u8 M/ y/ H
      z9 h, n. W+ o4 i, w& I; h

    % n& l- y  a' I' n6 p% ~7 C8 g
    , F; C& `: E/ `5 r' s5 [* O- p9 u伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。& {& a8 w. S* j& b7 }
    5 ]0 W  |/ {6 r  R: a
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关: e% r; A. `& ^/ o
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑
    4 E; o' Z( j, x
    4 b2 y" D7 p8 E. }$ o, f拆分 。。。。强!
    , y$ I- Q8 k7 H& o) b4 A9 N5 o* y' O) {' x, y7 u
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;
    ! i% F/ j9 U; \8 e$ g' r+ F) c9 W! W, k) N! w6 l
    7  L5 ]: o5 A/ }3 P* R
    1905692926 d) E' Z. x1 L# J: v# \
    [
    1 u7 p& k1 A8 A/ C1 w    [ 10 ],
    / r$ o1 Z- {8 s: ?; u8 q5 N8 ?3 G    [ 9, 1 ],8 q2 B$ I0 s* J$ t/ I! U
        [ 8, 2 ],
    & H/ @- I$ |4 @$ B' G  C8 a* t    [ 8, 1, 1 ],7 P" d2 @& j/ P3 C: `+ _- g/ J
        [ 7, 3 ],/ ]) f# x) W* D: p- N8 E
        [ 7, 2, 1 ],3 N+ W/ o' b$ U* h
        [ 7, 1, 1, 1 ],' q8 V# h& Y) X1 ^9 C! F1 j  B5 F
        [ 6, 4 ],
    % Y, I! U; k! G% V" `    [ 6, 3, 1 ],5 q5 F! B8 K2 y& T- }% ^
        [ 6, 2, 2 ],0 T  W) ]! R& p, G5 i; L; K% F
        [ 6, 2, 1, 1 ],
    / d2 z* ^3 |: s( k7 S. Z9 p' P1 k    [ 6, 1, 1, 1, 1 ],$ @  p$ q/ N3 M& K/ v
        [ 5, 5 ],( o. I, p0 P+ K8 H4 Z6 u* p6 U
        [ 5, 4, 1 ],+ L2 f9 b/ v+ x% p0 C  d+ i; t, M0 `
        [ 5, 3, 2 ],
    5 Y, m/ e$ {0 ^' T    [ 5, 3, 1, 1 ],% V- U3 H& y% U' ?2 F: C
        [ 5, 2, 2, 1 ],0 I/ ]' v' z1 m  |3 A/ F. ?% z
        [ 5, 2, 1, 1, 1 ],7 _( [  |, [- p" D- X8 d
        [ 5, 1, 1, 1, 1, 1 ],
      n3 F9 @# S  m9 P& p9 ~+ m3 X9 Z    [ 4, 4, 2 ],. X3 C; r3 j$ L; z8 n3 o, {$ V
        [ 4, 4, 1, 1 ],* C1 s# d# e  y% T' w' {" Q
        [ 4, 3, 3 ],
    7 r- @- L' M# {/ H  S& H: b. {0 }    [ 4, 3, 2, 1 ],6 M) C: J2 C) z& ]% A
        [ 4, 3, 1, 1, 1 ],% y+ R8 A4 L) q, p7 x- b
        [ 4, 2, 2, 2 ],; ~/ D* D3 w0 v) w
        [ 4, 2, 2, 1, 1 ],5 r, r/ {# R+ P/ F  e
        [ 4, 2, 1, 1, 1, 1 ],- I$ O  j. Q) {' O. E9 {
        [ 4, 1, 1, 1, 1, 1, 1 ],) D2 Y" _3 r) o. [( y
        [ 3, 3, 3, 1 ],  v/ [$ E3 k7 h& |
        [ 3, 3, 2, 2 ],
    7 V# S, a" [0 M! J2 h1 K4 |2 i& S; c    [ 3, 3, 2, 1, 1 ],
    8 l6 c; o0 I' N9 T0 K# W: w, Q    [ 3, 3, 1, 1, 1, 1 ],
    ; \7 K4 m' J: y8 ]+ W: K* W# h6 @    [ 3, 2, 2, 2, 1 ],2 V, I. G8 p- P' l
        [ 3, 2, 2, 1, 1, 1 ],6 I# `' n3 y1 v& z% H6 ]
        [ 3, 2, 1, 1, 1, 1, 1 ],9 q& A3 f% l9 z6 w$ R: v- A
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    ( M1 E2 ]. C1 v( U& C% L  H    [ 2, 2, 2, 2, 2 ],
    & G( c6 `* j$ t. a    [ 2, 2, 2, 2, 1, 1 ],# @1 r0 b/ i% q& l
        [ 2, 2, 2, 1, 1, 1, 1 ],6 ?) L, ?, ?8 b* _
        [ 2, 2, 1, 1, 1, 1, 1, 1 ],0 c, g7 V7 u* C% S
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],0 }* p# c/ |) F8 [9 a; a
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    0 b1 {) S3 C" r]
    回复

    使用道具 举报

    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-5-25 13:11 , Processed in 0.504640 second(s), 92 queries .

    回顶部