QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5136|回复: 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 编辑 7 y7 [5 {9 o% o2 Y

    ! C/ L( @. F& T' R+ y3 G. \& \n:=12;n;
    % d; w1 A0 G2 A- mFactorial(n);求阶乘
    0 ?3 I1 z2 t% p. B+ u9 F- jFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));, l8 s$ e' L* n+ q
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);2 V4 Z- l& u: t# O6 ]! _- X8 e6 F2 |4 O
    NumberOfPermutations(n, 4);
    ( f; e  I) M. ^- d8 h. TNumberOfPermutations(n, 11);' i: a$ r4 g& o6 i8 s5 D  M
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;
    ; {) f& m$ H3 g/ j% Y4 sBinomial(n, 3) ;
      [3 ^$ t5 }1 {9 I% Q$ _# }Binomial(n, 9) ;0 M; _6 \% r" i7 `+ h  t
    Binomial(n, 10) ;! h+ O2 G9 ?: n* d$ s/ R
    Binomial(n, 11) ;5 |2 G- z7 Y/ m  H1 v: B! D7 u+ h( S
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600# }; C8 d7 Q1 V& ~+ T2 D
    ( E4 y0 j" [, R2 D/ \! v
    Fibonacci(n);斐波数Fibonacci(n-1);
    * {. [2 ^' n5 dFibonacci(n+1);" ^! C) ^: b- G/ q
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;5 y6 z; U1 D4 R3 U
    GeneralizedFibonacciNumber(0, 1, n) ;& {8 m! }, x7 D# `  b6 t3 k* [
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))) _8 ?. d# l$ q9 n
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    2 f: E( e/ @  r' o6 P( @' M) kCatalan(1);
    ; y% b( Q. p9 O& w% qCatalan(2);/ D6 |- |2 Z1 {( B2 v4 Z* p) j
    Catalan(3);Catalan(4);- K; Z7 i" A' |" }
    Catalan(12);/ c  c* G" I! i" _9 `6 \
    1 I$ L3 }  B2 P7 o7 D; y
    Lucas(n);卢卡斯数
    ! v: Q" w" L" I12
    : P/ e1 y$ C5 e; Z# M" @9 A479001600
    + C7 P: W$ Z% z- }) G' _831600* |0 f1 c9 A1 y! x
    12  l4 ~- ?7 {2 q# t" u# s  d2 y
    1329 _0 U8 O' e  m' I
    11880
    3 `) P1 o2 h6 B/ v+ \9 E" q3 A. M8 ]479001600
    ! {, |( I& I2 p& }; V12& W. ]- {' z2 R+ q7 [2 V, e
    663 R1 ^) Z5 Z7 c3 B
    2208 R1 t- {3 o2 L  ]4 `
    220
    & E: z7 B8 M1 @% P, C* p0 a7 O66. k' j* X# H& Z
    12' e5 r( K+ L" P
    831600
    ' L0 h  F& L& [! y7 T) d+ F; r) t' N144
    ' Z! ~" o/ O* B- Z8 x89* c, a7 _, w8 F+ \" K
    233( _; c! n( R$ Q* ^9 n* A
    233
    & [. |- V4 z  n! t1 Q. ?7 ~610
    8 {8 ^5 ^9 n2 f# W5 }% E3 F144
    4 U7 |* t4 g  C& {. K, t, ]9 q208012
    5 |' A( m. [/ h2 R208012. M! a1 C( K. d1 U/ j& C: G
    18 p0 X7 h) q$ h1 Q
    2
    6 h8 S; E' O% I8 l7 `- P5; s, I5 f9 `9 \; L* f! f
    14; X5 x) V" S+ u" }: j
    208012) q/ B# I# e7 N) Y
    322
    - {0 r8 K/ z# i! z+ n+ N0 t2 _" @1 d* N9 z
    5 \* |" Q0 Q# y  s+ {/ V
    卡特兰数=(2n)!/((n)!*(n+1)!))% L! K: L/ ]; L/ y! U% l
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: 7 [4 ?$ q6 ^$ m7 z% W3 B
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY
      E# n" Q. [# G; G( c6 n+ h将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: : N' }/ w6 w$ U
    ((())) ()(()) ()()() (())() (()())
    4 I9 `) r% H; _. \) tCn表示有n+1个叶子的二叉树的个数。 + J9 o: z& \- Y; V) y% g3 H
    6 x  N2 Y+ c! F; n8 Q9 O5 G& V
    Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
      R( [7 `! V1 A" n9 h" GCn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况: 0 i4 W: A# p4 }; _' {

    8 R; q, p4 f/ Z& K2 QCn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    % q& ^( j1 q# o. O. |0 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. $ M7 s5 q* N& }
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    * I4 n$ E6 r- u" T8 L. A) D% e) k6 V
    4 ]2 E& {" K: {0 w" t0 x

    - L) G" @! D; Y& [" K卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。1 ~  _. ~! B2 ~
    $ }- p! M0 A/ r8 c! ~# V8 z# O/ I
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    . B1 t- _9 P% K9 ^6 [6 y3 G( {# a# @9 [6 b

    " C+ B, o3 [9 w9 Y5 W1 E( P# S
    n:=100;n;
    2 }& J. H- Q, u4 e5 wa:=Lucas(n);a;
    , e( [7 F; S7 ^" kb:=Fibonacci(n+1)+Fibonacci(n-1);b;
    9 e3 B& k5 S2 ALucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    . O2 J* e7 T; k+ ^( T% O/ h) k( e( I
    8 S  H6 U- E0 I6 Z100
      ~$ V& U. C% M" e792070839848372253127
    ) t1 |& t& P4 H: i7 C% e5 s792070839848372253127
    , E& A/ m! U, L1771124240896309575375
    8 f3 ?$ R, o" a4 r; j* y1771124240896309575375

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

    12.JPG

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

    74

    主题

    6

    听众

    3299

    积分

    升级  43.3%

  • 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 编辑
    4 h0 C- W' J/ `. m; ?" r% z/ F$ U5 @5 L' f+ {
    反费波那西数列反费波那西数列的递归公式如下:# ]; E& i' R, k# n. _
    6 {( t/ b) S* X# a
    Gn + 2 = Gn − Gn + 1 ; w8 Q+ u& _& V0 D5 T$ V/ L- y
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ..." l: p& b2 s# W. P; i8 {

    8 P( {' b' L' {( B! Y% K6 q9 h! m% H) @& [即是F2n + 1 = G2n + 1,F2n = − G2n。
    ! ^. Z4 t7 Y* `6 k& S( N
    ; h6 [: K3 ~, N4 L8 x! kBell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);2 d0 h* @% o8 _+ O
    StirlingFirst(4,3);
    8 V1 c, d: V4 I) j3 Q$ n" g: oStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
    , M# m& I. ]* N  z! p. ]StirlingSecond(4, 3);( l6 P7 K. Q; }9 P
    2
    * r( Z4 }3 W$ W$ j& D6 Q52' R' [7 n4 G( ?* ~
    5
    , H4 C  [9 l9 J; e1 N1 _( R15, Z, |6 o9 V/ C% G: i! |
    -6; T# e6 x" A5 S' U+ E; x5 s
    11; m! y! P- M" J+ e, H3 R+ S
    -6
    0 S/ b5 u) a$ f% I& @& ]1" \" D$ }9 \* J0 e5 L7 E5 |- I
    7- s! G, u  _3 Y4 Z+ Y
    6) }, c# L  ~) \0 R  T& w7 y

    . Q. I! v, G6 n1 oBn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:* k* o, s! a0 y8 C3 d+ r
    ! m9 m3 c2 |) A) B/ q9 ~
    {{a}, {b}, {c}}
    & v6 Y: _- w3 X, j( B- C{{a}, {b, c}} 5 L0 u- ^# o) N7 d
    {{b}, {a, c}} " a& ?3 a. `# }. Q2 {0 M* [8 Y
    {{c}, {a, b}} 0 a; d8 Z9 }- h, E' u
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
    / I' j3 H9 Y/ _s(n,k)是递降阶乘多项式的系数
    $ \! T7 N+ r/ m9 n& `有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)# x, ]1 r3 _+ x- x

    " ?9 y- {6 t: ?, E4 J换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):# @1 q% y4 U1 J* z) l: t

    6 e9 P& h: R+ F8 o& u, V% V& z{A,B},{C,D} ! V- C4 f0 t! Z
    {A,C},{B,D}
    6 n$ `- R; t1 v7 X% l{A,D},{B,C} # Y' _: n8 S5 ~* p# D
    {A},{B,C,D}
    $ P' ^- [- u2 u' D' |) L% ]{A},{B,D,C} . }2 p% [3 y$ \3 }$ E
    {B},{A,C,D} # q) A# x. f: Z0 ^* m6 w
    {B},{A,D,C}   |* C. ^& R% ^% S' E- H/ X) {
    {C},{A,B,D}
    " c6 A2 M  J; e4 n3 p{C},{A,D,B} & B& ~& j  c- n3 L7 J' Q' b
    {D},{A,B,C} 5 b3 `/ P6 D" I, w/ l$ ~5 n
    {D},{A,C,B}
    4 @& w, `0 V; K: c/ y- n
    第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    ' R# M& H/ ?: N* ?给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) ( C1 Q7 a' ]% n" T' k
    S(n,n − 1) = C(n,2) = n(n − 1) / 2
    1 l2 K5 ]7 y. f: F& wS(n,2) = 2n − 1 − 1 ( F" @# }2 z# e# Z; h% P
    . G, |0 z9 m9 l* s$ G2 `7 Y/ l; f9 o
    换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:. `" S4 r+ a& {% Q' ?; I& Q' B) w

    . Y, r; V2 Q1 R8 y9 X7 U& t{A,B},{C,D} 2 ?2 g, f$ a2 s! d  X
    {A,C},{B,D} 8 D8 H( k9 R3 H' Y9 I, O+ f
    {A,D},{B,C} 0 P* b+ X' J# c; ?& w! x
    {A},{B,C,D} : Q3 K, o7 X1 T& I; F: Z3 h
    {B},{A,C,D}
    " h& l% m, x  }{C},{A,B,D} 4 B( P8 N: o1 l' m6 n1 k
    {D},{A,B,C} : S: p8 Y, N. u
    因此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 编辑 5 q$ n  l* }3 ]8 B: M2 a) c+ o
    $ L. _2 H. o' I7 g( A
    n:=5;r:=3;. w: W/ ~* n& \  I' u; U
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;- A: Y" a6 F$ T! I5 J% U3 i
    BernoulliPolynomial(n) ;伯努利多项式
    5 c# t; P6 R; ]8 C8 O* M: b0 q) K6 i  _) b  s- V5 J" Z
    26
      y# K) t4 \2 i: [0 C. D- \% O137/60% v& v) p  N; I% u5 y6 a& ]
    0) I8 C5 Q, f. @) \! n7 h
    0.000000000000000000000000000000
    & u& {! Z( T+ e: S$.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑 5 O. P8 w3 G6 Z; s

    $ r4 ]" T% t9 r; {0 g% Q: _# @8 O
    7 @. j% U* U) I' h; ?4 w4 u* n$ V* ]$ E$ G9 p( S
    伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。: Q; R" N! u9 W' [  Z; `- N" u5 \6 L
    + _6 O. K& w( p% l  i
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
    0 `' \5 W! ]6 W5 ^4 `) l$ U
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑
    # W% j6 m. K0 L$ ]- B7 I6 n8 a6 P$ P$ B  v0 B2 M/ `8 W- g
    拆分 。。。。强!6 l% R$ u3 P1 x" i9 @0 G) o- W! s

      `$ N. Z; R# p8 }  D2 j: oNumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;7 }$ e1 L/ P# K6 {) s; K

    8 m% v+ N4 b  p# _2 m7+ Z. t$ D! p0 g" X
    1905692928 U# l2 [3 ~' {; R. o  |
    [
    % x+ U2 g$ c7 O1 Z    [ 10 ],
    # X5 O2 M0 _, X" b* ^3 N    [ 9, 1 ],
    3 i  V% a/ I& w; G    [ 8, 2 ],8 O# s% p( U; h, O' z0 [6 L- `
        [ 8, 1, 1 ],
    / ?  a3 R0 g% N" [8 w) ^    [ 7, 3 ],( t7 f0 v  f# l$ n  @  r
        [ 7, 2, 1 ],5 i) j/ }3 w1 N- P/ `: f+ y9 X
        [ 7, 1, 1, 1 ],& e4 \  S- l4 p
        [ 6, 4 ],0 }2 u' \6 r6 H* \
        [ 6, 3, 1 ],, o+ z" U; g5 c
        [ 6, 2, 2 ],  e7 Y- d4 ?$ ^5 U5 y, }
        [ 6, 2, 1, 1 ],- |; A7 k& k, t; B  v
        [ 6, 1, 1, 1, 1 ],
    $ o) k. P# G+ u) z1 H  P    [ 5, 5 ],
    1 p) d* {0 i" [( M% W+ D2 d    [ 5, 4, 1 ],% K9 @4 U1 M& ]1 u$ X, T
        [ 5, 3, 2 ],: W% ?' |+ A$ D4 ?
        [ 5, 3, 1, 1 ],' s) T2 t$ Y% f$ z, E
        [ 5, 2, 2, 1 ],
    2 S/ e- W9 ]+ F$ a) E& S; d# o* e# N2 k    [ 5, 2, 1, 1, 1 ],4 M. d  Y0 ?% C1 v* n
        [ 5, 1, 1, 1, 1, 1 ],
    1 @& H! M0 m$ R" P  I  a1 z8 i    [ 4, 4, 2 ],
    8 p0 u0 ?( p9 s  r2 s    [ 4, 4, 1, 1 ],9 X# J  q- z. u5 U( L2 x% _3 H
        [ 4, 3, 3 ],( g" i9 E+ }9 o9 d: I
        [ 4, 3, 2, 1 ],
    9 u) ?, z. S5 C0 q  Q    [ 4, 3, 1, 1, 1 ],
    & O9 a, ~- M. s    [ 4, 2, 2, 2 ],
    - D& H9 U, S9 d- h    [ 4, 2, 2, 1, 1 ],/ }: W, X* X5 @; J
        [ 4, 2, 1, 1, 1, 1 ],
    3 W: ], m+ q  s- ]" N% U) w3 U    [ 4, 1, 1, 1, 1, 1, 1 ],0 ~, [8 K. `% ]! A, N" [% A) V% @8 L
        [ 3, 3, 3, 1 ],
    7 r  Q* A6 t2 a! }1 o. t! M! D    [ 3, 3, 2, 2 ],& M$ P6 ^* w+ m
        [ 3, 3, 2, 1, 1 ],* I( U( h% ~$ L. j; y
        [ 3, 3, 1, 1, 1, 1 ],% b, [$ y! J! t5 R4 g1 j9 X  ~
        [ 3, 2, 2, 2, 1 ],& i  c" _4 D3 H9 V
        [ 3, 2, 2, 1, 1, 1 ],
    : g. V: o  Q  s1 |( X% ^) C9 q    [ 3, 2, 1, 1, 1, 1, 1 ],
    5 P* V4 n0 E6 f* Y+ @8 |( @; T    [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    0 Y4 h+ }  V! d, _6 b; o    [ 2, 2, 2, 2, 2 ],
    9 y/ r5 b: d  a    [ 2, 2, 2, 2, 1, 1 ],
    / [  m! A5 i/ I    [ 2, 2, 2, 1, 1, 1, 1 ],8 Q) Z/ j) T8 b2 N( K
        [ 2, 2, 1, 1, 1, 1, 1, 1 ],
    7 R* _- ?0 f& q# Q: B5 B+ V    [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],
    / e& X8 o! G% k+ i1 g! W4 m    [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    2 T; [0 ]  Q* X1 L7 J7 Z  k  ^]
    回复

    使用道具 举报

    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-4-14 11:48 , Processed in 0.508226 second(s), 92 queries .

    回顶部