QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5255|回复: 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 编辑 3 r( `: E' {' `0 M7 q# q" x9 G7 v
    + F+ z" y+ r5 {; [, H6 ]
    n:=12;n;
    0 L( v/ h0 S2 r0 a0 ?: a5 u7 _0 |Factorial(n);求阶乘
    0 d: [$ P' S) CFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));+ J4 p. N( f7 c0 N# g1 g. O# B0 {
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    2 Q: e) W) D* A( f  K' BNumberOfPermutations(n, 4);
    + e, u& v. u" n6 h! MNumberOfPermutations(n, 11);$ }. a. b$ @: F8 z
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;. r% s) y  ]& j2 `9 F0 n" T
    Binomial(n, 3) ;* G5 E5 @& _& v! I7 N9 h; Q/ t
    Binomial(n, 9) ;
    2 D5 J- Y3 U& l8 }% M- q. e! bBinomial(n, 10) ;
    ' u0 _# l& a. E- k2 a$ |4 GBinomial(n, 11) ;
    ) O- x4 G4 H/ G, d! `# [Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=8316004 G9 k* J! s% n' I, Y

    ! g- j# L  p0 zFibonacci(n);斐波数Fibonacci(n-1);
    5 c# d5 @% U" L9 Y. SFibonacci(n+1);
    $ }% \: S, E5 w4 c# C% j) p( g8 uGeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;8 Y! J4 Q: v2 w& F
    GeneralizedFibonacciNumber(0, 1, n) ;) `, p; W' ~& |9 D, a
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))
      j( ?+ M- }4 ~& G% q) z9 _k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    ( O  \, `! T! y$ t& T! ]! \Catalan(1);
    : d/ F2 E: H# @3 V, q# TCatalan(2);( b. Y6 a3 s# s% F. Q: e
    Catalan(3);Catalan(4);- r0 Q! l8 Q) }- E
    Catalan(12);$ x$ x- _1 c* W4 l5 @/ v' \- ~! E7 J
    7 \& R7 ~! @% e4 A5 {9 D* d8 C
    Lucas(n);卢卡斯数
    ; U7 T3 {: U2 b12
    7 q' O2 a% B5 N$ |" o479001600
      m+ E$ g# w, C- G. p, _831600& t9 R8 c# G# C8 t/ G
    128 S; j3 Q: u2 m
    1323 K/ }$ A! H" \* t  W. S
    11880
    5 |0 W' m; c* X479001600. s5 X6 U3 [+ `0 u$ o: ^' }
    12/ ~+ Y- F( C! J. G' U6 `0 `
    66$ e! b- e) O: P( m6 V
    220
    ) t+ m6 i2 V( ?2201 ]. V1 H# i! B9 Y: U
    664 {4 t1 A2 a* q- W: s, u
    12
      @* w# d1 H4 V! x* S831600
    7 J" r& x! [% H' V6 M1449 `- D% q1 i! N" |# Q
    89
    * D+ W+ K& [8 T' {233
      ]2 a7 v9 `- Y! z3 n: ^233
    1 s' y% P/ A9 q9 e/ T4 v( d) H610+ R5 R2 i' e' `% V' k: {
    1449 }( c. u; s' C8 O. Z8 \
    208012
    ( F( n6 f2 Q. y% Y; P$ P3 L208012! A5 }. f6 d8 K/ [
    1
    ) N& v2 a, D% \5 E2
    : f' c- {: [0 z5
    3 l, {$ N" g5 r7 b0 d14
    & V' e# B  R' Y5 m& i3 B- Y208012
    ( T* B9 ]+ l9 O322
    . u  H" B9 p, z  b& `, {0 y0 b# Z8 w
    0 L2 h- u# Z+ u: `3 s& E/ k: w7 [
    卡特兰数=(2n)!/((n)!*(n+1)!))( f1 n$ w1 E( s2 x" w- a
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: ' b  D# z" }+ g( G& z1 U; g2 K
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY1 J$ _! M, M* j5 u$ x8 w
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数:
    # J: C+ q9 ^6 N- h  X% r+ Q((())) ()(()) ()()() (())() (()())' |* d9 }) g8 R( P1 G4 \& K
    Cn表示有n+1个叶子的二叉树的个数。
    ! u$ f" W0 V. P8 k
    % S. x6 ]# l1 X/ {( P( O+ fCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) 1 i4 n7 j0 e- ?8 r5 \/ k
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况: ; Q) \6 S1 ~5 k7 J5 N

    : F4 r; ^3 _9 b4 [, @0 f' ?Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    4 ~5 H$ O  f9 K/ yCn表示集合{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. ; k  X. B+ o' m' l1 A
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    ( O7 v0 {- j+ _: b* Q8 f) R  @/ ?# g  N

    : s" ]8 \. L  _$ n
    ' j+ N- Y; c  v卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。
    + u/ Q/ Y  d( c/ ]  G% O* n7 M2 s8 c' w, T3 X" Y7 X
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    # ]& x: g+ z; f6 l+ z7 K/ d6 W! j9 g: \- ~. _/ x
    / N3 @7 O- ?) X6 T
    & C6 l; B% }% i. F) A: l" b
    n:=100;n;# i7 |4 t# ]! k" D0 t
    a:=Lucas(n);a;
    9 m) P& L1 d2 p+ A. \b:=Fibonacci(n+1)+Fibonacci(n-1);b;
    * W( I$ w" G# u3 P) ELucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    % n% p7 h/ O8 s! o% W% c
    + [" A( X& @& J100
    & q" Q; X3 P; |# L: o792070839848372253127% p) k; v* P0 C3 V2 I- |6 P6 x- P5 J
    792070839848372253127
    , B/ [; V$ f; j0 v1771124240896309575375  `" O: V: U$ \" w0 J
    1771124240896309575375

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

    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 编辑
    + N/ F3 T: }! J5 f1 g, A
    - i  K6 j/ j" t/ `  L5 Q6 R反费波那西数列反费波那西数列的递归公式如下:
    7 l2 @* m  C  I' t1 F8 E. M6 B
      s$ O" f# i3 k) u5 XGn + 2 = Gn − Gn + 1
    4 g$ q2 l8 \. y$ D* T' \如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...
    3 H+ }/ C$ L1 V+ X$ W4 P: N/ R9 K$ d/ u: y' \
    即是F2n + 1 = G2n + 1,F2n = − G2n。
    + k  N$ h( U" |
    0 \4 f' ]) B1 C6 X% ABell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);% g- G  x0 @+ m+ t, O, S
    StirlingFirst(4,3);
    - ^6 L2 P' D! X" i) ]) O! JStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);: c) S. t+ _. i  B
    StirlingSecond(4, 3);" ~3 ~$ u1 @6 \0 m2 a
    2; t" ]5 n4 A2 ]; m* F: \
    52# V6 Z# N4 U7 J0 V8 x
    58 |/ {4 S* V* P1 Y/ T
    15+ _: g8 _1 E; Q
    -6
    3 j$ ^3 k3 r" p$ T. m/ x11
    * p7 @8 G* C& E; C4 v-6: L( X4 ~% p( E% |) \) j
    1
    4 F1 a4 ~8 j7 A2 h1 Y4 e5 h6 H7
    ( j4 s8 p. m, h0 B6
    ; I' h7 T2 W) x" y5 l, `# {8 s  ]/ {
    - _: H0 s  @0 Q: O* F1 cBn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    ! W4 M1 }, [' S! V- b" V% F& ?  z. A' }
    {{a}, {b}, {c}} ' A2 ~# x1 s/ ?' P% U9 S: U' l
    {{a}, {b, c}} & p8 Y/ t$ i7 b, k
    {{b}, {a, c}}
    ( K- W7 @) ^  J! D* t0 n: U) b{{c}, {a, b}} $ O( E2 s% F  C$ O
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
    . n4 D9 w7 Y" b6 R) K4 ms(n,k)是递降阶乘多项式的系数
    0 V$ |3 x( ^) {3 D2 V/ d有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)
    / l: u+ D+ k1 _+ W8 d
    + k6 p! a0 v) n6 k* }5 r8 r* N换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):
    3 t9 E8 g) i- m( D6 T- N5 J* Z+ A
    {A,B},{C,D} ' x1 ^8 _+ A0 b6 O) V5 ^$ ^' b- V
    {A,C},{B,D}
    - D! W  H* M7 {7 Q2 F% ?' y0 T{A,D},{B,C} - A' }' c( j; J/ v
    {A},{B,C,D} ) l% A5 M4 ^' _8 a( g( V
    {A},{B,D,C}
    0 h# E! _2 n5 D/ J{B},{A,C,D}
      d) C0 l" O# N% F0 L{B},{A,D,C}
    : e9 n6 p/ n+ P% V{C},{A,B,D}
    ; O" E$ H/ a' n2 W7 Z{C},{A,D,B}
    / h! ?, U) b* C: J{D},{A,B,C}
    * S+ J  g* A" E1 ?" c8 @' Q' [0 E% X{D},{A,C,B}
    3 q$ D6 s( }( i5 K' h- {3 W3 e; T
    第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S6 G& g: H, H& ^2 }7 z$ n5 K) v
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) : u0 z. J, s) |" _
    S(n,n − 1) = C(n,2) = n(n − 1) / 2
    * H$ b7 e/ D( ~S(n,2) = 2n − 1 − 1 0 ~+ `& \! h) P

    ' A( y& Y$ k* L$ c换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
    $ C; r- j- r9 Z* Q+ x7 Z. N/ E5 o' M. z! _
    {A,B},{C,D}
    ' K- D( E& y& R{A,C},{B,D}
    : O9 G! _' _2 `& c7 S{A,D},{B,C}
    3 F" {' t! o! b4 o+ M* `1 T{A},{B,C,D}   z+ W/ A% w; \! U( G$ D6 b
    {B},{A,C,D}
    9 k8 k& Z! w- C: @* w{C},{A,B,D} ; F5 d  O- x3 q% V4 P" h/ Q: I
    {D},{A,B,C} ' b7 ~1 t9 s6 c; n" V' b+ U; V
    因此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 编辑 ' `! n6 F& m6 b& y) K$ E. u
    9 K! q" k, p9 D
    n:=5;r:=3;5 [+ J2 L: f7 I- \& ~
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    7 `# }# l* `( B0 ]; F. LBernoulliPolynomial(n) ;伯努利多项式/ }5 R3 [: O: L7 `
    # F2 _6 y' z0 a5 U* P+ o
    266 S) ]+ o. p$ j, h
    137/60
      e$ u2 h/ G* S3 \8 S0' y# r6 a' r4 X' b) D7 s6 k
    0.0000000000000000000000000000000 v2 l, x! L6 m4 H
    $.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑 1 r( I, o. g% a& N, @2 \

    7 o2 O. I4 |/ u3 z' g$ U% ~, D/ H( F( b' F/ ^' B/ Q

    0 A. V2 M6 k  B* v$ O: I伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    + c- ]8 Z5 e6 D3 n2 p2 k
    9 \* c0 z6 ^8 n* E: H. j2 D- k伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关; |# u7 M' V+ _9 |! @4 Z9 ^
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑
    $ ~! @% e! i. Q3 Y! B. E4 b% l3 g0 [$ Z+ Z9 S' L
    拆分 。。。。强!2 D8 |2 g; Q+ }, i
    ) V, `$ W  D. d3 D
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;  `  |2 B- a. g' U3 H

    8 A  q: T+ U& h: \$ \70 p9 M0 M4 O; e  n+ l2 }
    190569292$ o% @$ L% D7 }$ L; l8 L; g
    [
    ! T) j! t" j# m4 D/ N    [ 10 ],- j2 N% B7 J" s5 j
        [ 9, 1 ],0 u. _- i: k5 ^) W6 Z* ?
        [ 8, 2 ],! q1 J5 Z: \8 i* X9 N- G
        [ 8, 1, 1 ],  s* `8 z7 t9 A( N3 _5 @7 j* L
        [ 7, 3 ],# a; z" F/ T/ O6 m: Z# J7 s* t. `! [
        [ 7, 2, 1 ],7 l  B- E' \1 r# C9 ?: `2 l( q
        [ 7, 1, 1, 1 ],
    . T* r/ i# a7 {  [    [ 6, 4 ],
    9 m7 t% e% |5 V! w    [ 6, 3, 1 ],  L7 \6 X8 H' D1 O5 ^2 ]- n
        [ 6, 2, 2 ],
    . J8 F6 t) n. E    [ 6, 2, 1, 1 ],% i+ K: f- N, R+ y, U% x
        [ 6, 1, 1, 1, 1 ],0 G8 L( J, E: O8 S5 I. H6 S+ y
        [ 5, 5 ],
    - L! v+ c- L$ W( N    [ 5, 4, 1 ],# s' [7 m1 U1 j  U) G. f; f
        [ 5, 3, 2 ],
    ! M( S7 `; p- y4 E; w4 t    [ 5, 3, 1, 1 ],
    % n+ Y8 }, z" [0 E& Q% [* N    [ 5, 2, 2, 1 ],
    9 }1 F3 x' A8 ]% C1 m% K    [ 5, 2, 1, 1, 1 ],. ?! L2 O" r: ], k
        [ 5, 1, 1, 1, 1, 1 ],  e5 ~$ b' y0 m6 u0 ?
        [ 4, 4, 2 ],! `; c7 b" K$ g& m2 u3 u' n
        [ 4, 4, 1, 1 ],9 u, T2 t' e# ^) E6 `
        [ 4, 3, 3 ],4 E' K$ e# K0 K
        [ 4, 3, 2, 1 ],6 W% s: f+ i) I. g
        [ 4, 3, 1, 1, 1 ]," Y, v: b" h: t( x0 b2 F, K$ R/ }
        [ 4, 2, 2, 2 ],
    $ ?. ^9 d2 h5 @9 D% w0 [  S2 d8 x    [ 4, 2, 2, 1, 1 ],
    " O* n. l  a/ ?/ V3 c    [ 4, 2, 1, 1, 1, 1 ],
    $ B: v2 i% ?$ f- p  w, T% x    [ 4, 1, 1, 1, 1, 1, 1 ],- o; m' c: K, J2 R0 [  i8 ~  t6 b
        [ 3, 3, 3, 1 ],
    0 A- h4 p% F; ]  Q9 G    [ 3, 3, 2, 2 ],. i1 {1 F4 ]* |
        [ 3, 3, 2, 1, 1 ],
    * m( W2 m" y: `8 f    [ 3, 3, 1, 1, 1, 1 ],2 R8 T! W3 F/ g2 Z
        [ 3, 2, 2, 2, 1 ],% I% X$ ]9 U4 |- U  w
        [ 3, 2, 2, 1, 1, 1 ],) Q; Q  t) d' S3 g3 ?) Q3 w) V
        [ 3, 2, 1, 1, 1, 1, 1 ],9 M/ ~: x0 S. W; m6 @- c; k8 Q
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    , r3 T: s" K' X0 j    [ 2, 2, 2, 2, 2 ],
    ( P2 \0 D" t# _) y    [ 2, 2, 2, 2, 1, 1 ],
    , E/ B/ W4 z% F  @' z( O6 w/ e    [ 2, 2, 2, 1, 1, 1, 1 ],
    ; ^+ i3 h& [+ j( R    [ 2, 2, 1, 1, 1, 1, 1, 1 ],, ^; D6 M4 b" }
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],+ k2 D8 P" |( Y9 }
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]/ w& A  ?5 o2 m
    ]
    回复

    使用道具 举报

    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-10 01:20 , Processed in 0.513152 second(s), 90 queries .

    回顶部