QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5194|回复: 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 编辑
    ! n4 v4 `$ g& G" k: K. O0 E7 ?3 H  @/ \; r; L9 n
    n:=12;n;3 X- Q% x! J8 a  H7 s
    Factorial(n);求阶乘
    # w# d7 a' `% ~6 \9 B0 {, aFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
      e- v( h" N- w2 u; `2 VNumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);6 z* K& F0 V  \% b
    NumberOfPermutations(n, 4);4 Q- }- p. j5 }
    NumberOfPermutations(n, 11);) {, V; W% |! T" \) M' U7 R
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;" J& E! k1 h, ~7 ~2 r1 }# |, G
    Binomial(n, 3) ;
    . `. D; Z" Y& p, W% F- Z0 Z) _* FBinomial(n, 9) ;
    1 O6 r/ m6 i$ K/ _1 [4 IBinomial(n, 10) ;6 W2 W; A# ~- e7 R# J
    Binomial(n, 11) ;
    / ]* J3 D9 A1 i% U0 f! m) d; QMultinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    9 V( N3 U: G" R5 s7 n, Q( u7 W) X+ j. w9 B
    Fibonacci(n);斐波数Fibonacci(n-1);
    * I# |! J0 v8 \8 U1 kFibonacci(n+1);
    3 W- P2 b  \4 j* l) TGeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;5 f7 Z7 `7 K2 M$ y  K7 X
    GeneralizedFibonacciNumber(0, 1, n) ;, P& w6 E! P3 b; I1 w
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))9 o. V( e4 P! N# r. ]
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    * C8 p. k  n# x5 R+ C" f- `# WCatalan(1);7 d, D% X, F/ D, W8 X; p
    Catalan(2);! x0 m+ Y# D7 X* Z
    Catalan(3);Catalan(4);: Z6 \, R7 o0 j" u% P  v
    Catalan(12);
    % H  Z: z* y% h
    * i1 [# I9 l" }! x' @! dLucas(n);卢卡斯数
    % e8 \; R! N2 k4 a. Z5 C4 J+ I12
    $ n7 l+ r3 P. D: ]479001600- z; Q3 }' f( Q* r% [! x/ I
    831600# p  |# R1 `3 q% R
    12) h% [. {" E) K* _8 t- {7 }% u
    132
    7 t% M; L7 \4 e  y11880, ~. a1 i# J3 K2 ^' O
    479001600
    $ z, d# N: \5 G1 U, t12# y5 z+ ?1 n% U/ f7 c/ t6 w
    66! s" w8 c/ p' B, R/ w7 R
    2209 D' e: u, ^( n$ h( ^: }" h- w
    220* O, _3 C. D4 }5 r  @# f" D
    66$ x% h5 ^3 A7 b8 \! p, ^
    12
    9 e, L3 o' k. C831600
    ' U- t9 q6 g$ |8 w144- v' f' [  ^2 |: X: w* d# J
    89
    & e9 d- J9 |; v: s+ X233
    $ y& q/ J5 u) r- k) u8 m" l233/ A# s9 s# C3 A
    6105 }, a5 U- Y6 w% C8 t, C$ K
    1449 C0 ?- E6 W' x$ t0 w3 v5 o8 ~
    208012  N5 _2 L3 @7 F9 E  ~, L
    208012
    5 Z3 d# B! F5 L% v1
    3 k% D5 `4 a( B- N2
    + \6 ]- x. c- g: n& l  i59 E( s- G2 f3 O0 P
    14
    # M* n, Q  f- ^2080120 x/ H$ R$ h2 v/ T5 }( c
    322& _% M3 ~  \) {" Z
    / c* N- O7 n  S. ]+ r9 X

      r) s- M6 M1 M. n0 x: T卡特兰数=(2n)!/((n)!*(n+1)!))
    / z, j* F8 m9 V1 h& nCn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    $ G/ p0 b1 h& y$ d/ S) U# R  }**YYY XYXXYY XYXYXY XXYYXY XXYXYY, z* O* ?( S, x
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: . F) R5 b( I  i' l
    ((())) ()(()) ()()() (())() (()())
    & g8 J" P7 L: f* E" WCn表示有n+1个叶子的二叉树的个数。
    & W+ \" |& i) z# b; {" ^& Q- P/ t# H& U
    Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    ) q5 ^9 P+ G9 c8 i: T) l$ p" BCn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况: ) m( A4 f* B" |
    : `: d2 P$ S6 z2 ~; \
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    . ^% I$ e) Z& \1 L% A) G% z9 K. ACn表示集合{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.
    + Y, K, t* l3 E/ }# ^5 b* NCn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    7 L: {* o- a$ o  Y% U# L. K* H7 y) F" ^

    # [5 V$ o; H5 X- }/ z- v
    1 @! y' d+ c8 M0 n0 f6 S卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。
    ( _3 b, p/ d8 c( s; O% Z- E$ u" g, G, c7 |( \! I
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同
    : Z5 M$ a9 c# }

    ! _" B' Q5 p; i& U' v4 p+ t* _' w: J' e2 r

    2 w( I; Z  t% Y: V. z, ^n:=100;n;
    " L$ L$ Q! s/ `a:=Lucas(n);a;' Q6 s/ q2 |+ d1 q' N8 w3 u3 p
    b:=Fibonacci(n+1)+Fibonacci(n-1);b;2 h1 U2 [( ?' I0 S) G
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    8 c) z1 h* \8 G  A4 b& y. ]! L. B7 {
    100
    : D/ I8 S8 U. }8 }, h792070839848372253127; r0 l* O! w- H8 J4 k" j, @
    7920708398483722531277 R2 e, ?* {  c, s6 d: v, _8 B5 @
    1771124240896309575375
    $ B! j2 A+ k" ~3 G% s+ r+ B1771124240896309575375

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

    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 编辑
    : j+ l( _- l1 X1 M) f+ _/ j" |5 }) {2 S5 M3 w
    反费波那西数列反费波那西数列的递归公式如下:/ x6 ]* \" c" Y7 G; s/ r

    8 A( x& G$ p. n. f$ t1 p  TGn + 2 = Gn − Gn + 1 ! u. q. w$ d) r  l* {* I+ q! g
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, .../ q: }7 P/ y' O. R( U+ @! L( \1 [

    1 k2 e$ G, o$ s: m2 w: {即是F2n + 1 = G2n + 1,F2n = − G2n。
    ' \: f# R" B8 i2 n: J+ Z' d2 B) c* H! |# ]5 L) @
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);. @/ e. w$ V% y
    StirlingFirst(4,3);' J' n1 T$ G- t4 T4 J
    StirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);" G2 `$ z- y6 D8 a2 s3 _0 N0 f
    StirlingSecond(4, 3);- ]( H7 f$ g4 A# v
    28 f% ~2 @( y9 l( n  w- B- i8 r
    52, _. v5 U8 w2 M0 q% V: z
    52 I5 n9 @6 Q4 B' G. v
    15" ~9 }' |; o6 L" e& q( l; v
    -62 m$ x" }( H# \7 O0 z1 Q) ]7 G3 s
    11
    $ b1 ], ^# a- m$ j$ D-66 p" v' `/ }" R* ]6 H5 _* X
    1
    ; b+ D" j4 d# w3 g1 U6 w) A% r7
    . E# O2 a/ k5 G5 v& \( r6' M, T" n! z" k" ~) J" G
    3 L, H( {& Z8 C9 v  _
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:2 a0 w5 b% t3 _2 z
    * t# d. G* `; u! t6 A
    {{a}, {b}, {c}} ; m' F- {! [4 [% @* j
    {{a}, {b, c}} * E0 M. l, h5 C0 H) V6 g
    {{b}, {a, c}}
    , ~! c! |9 a, j8 A' n# C" W1 `{{c}, {a, b}}
    . n6 d5 y' O/ e{{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s2 M; l# ]4 H' e6 P1 s8 e; g
    s(n,k)是递降阶乘多项式的系数
      C2 z6 v3 |5 ^有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)7 E" G+ \- y/ d8 _3 {

    5 P! X" @: F. c) U" R8 T换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):
    2 U7 J+ G( J: M1 t8 q) C' t5 e/ ~! u5 i' m0 g  P% l
    {A,B},{C,D}
    & B) Z4 {0 ~% U' J5 h  v1 B{A,C},{B,D} 5 L; d6 {7 N. N5 \4 o2 t9 ~
    {A,D},{B,C}
    0 [5 E" B. f' U3 A. I9 v{A},{B,C,D}
    8 ?3 ?7 \* [. R5 z) f5 J& s/ G7 U{A},{B,D,C}   z; |2 M$ Q3 O- Y
    {B},{A,C,D} 0 }- E* v: X* _. ~8 K8 O7 ^
    {B},{A,D,C} ! N$ l1 L. w- z6 C: f
    {C},{A,B,D}
    & |& c1 T$ @: m+ O$ w) V{C},{A,D,B} 0 U7 o4 L& {1 [
    {D},{A,B,C}
    5 s& M$ P: v9 S) f8 |# v+ k{D},{A,C,B}

    6 X6 E) _' q' l5 S5 y5 E: ]第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S7 W: U6 d6 J7 V& W0 m2 z* f
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k)
    : @+ Y) S+ d7 U9 a6 x9 C, FS(n,n − 1) = C(n,2) = n(n − 1) / 2
    : m& I3 T& j$ O3 h8 p3 B2 P3 ?( VS(n,2) = 2n − 1 − 1 8 C1 M, ]+ ?9 _

    " I" E% C5 V5 W* y4 g+ k" D换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
    ' x1 o4 N3 G0 q! S3 C, J
    ; a2 J0 c: Y: l4 V) i0 F. g{A,B},{C,D} ! T: Y! g. _8 g6 S! ~
    {A,C},{B,D}
      C) \- w* I, S8 w& d{A,D},{B,C}
    8 S' W; ^, x0 p9 F$ {* n& R; J{A},{B,C,D} 8 ^& X' u+ |' w6 j+ W9 \
    {B},{A,C,D} ) p# Z/ L% S! k! v
    {C},{A,B,D} , S9 v* h8 x& `6 u
    {D},{A,B,C} , l9 K8 e$ Z& H- J6 R; 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 编辑 % ~' k+ F. Q  P! a& s1 r2 n

    5 N& G) U7 H8 ]* T" Zn:=5;r:=3;
    1 ^% e$ h7 m3 }+ H! `/ DEulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;% L; b$ D) L# u. k" h! U
    BernoulliPolynomial(n) ;伯努利多项式5 W3 b) m7 X' g% T8 |: l

    9 D. o6 R6 ]  J7 g  x1 v; q26
    ! |0 u: g+ W6 B& q! K, P137/60
    9 g1 w: v! A. i3 W0) q& y$ F% `6 a7 q7 D3 g
    0.000000000000000000000000000000
    ( c" W/ m( t7 h1 k: R1 u' ]  y$.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 19:56 编辑
    ! T4 R, M  P# Q2 ?! B. H6 J7 s' u; [. Q6 Y7 i: g1 x
    " [+ F4 k$ o! v/ N. V/ i4 L

    & ?9 j$ ^- `, ?! a伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    " l$ e! \, [  J  A5 @3 S1 W$ j
    9 J4 t% S8 {1 S. b% ?4 c伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关2 M# v0 u0 ^/ M0 n/ D% W7 u
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑
    ) D$ I# F; K. c) `9 t
    4 [. q) ~' f( E) A6 e: z% F拆分 。。。。强!! K8 i' X0 \4 K8 W7 ~2 P$ `, n
    . f5 J8 L3 S# _  m! k
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;: j& ?6 b( r8 ?# \# Z
    " T) d+ M8 l2 Q7 z0 c
    7# `' v" ~( |  U& U8 _! A
    1905692923 c- T8 L+ ^1 S1 k5 M" r+ h( `
    [
    8 Y! Q4 `: |  H0 Z- M    [ 10 ],
    4 {7 g5 N0 ~% X% ?% ?! F, G    [ 9, 1 ],
    7 t- I& y. S5 F3 k: r4 \    [ 8, 2 ],
    * R  a6 V4 P) \* l% X) @' y- [    [ 8, 1, 1 ],
    ; b) h1 t1 Q+ g& _    [ 7, 3 ],- p4 p+ B$ s* C. X
        [ 7, 2, 1 ],6 C1 p! R2 y# w" I2 T# F- m
        [ 7, 1, 1, 1 ],
    / i+ r, B* k6 `& Q    [ 6, 4 ],
    9 n4 Z8 b6 ?4 w. \7 ]8 z3 W% s    [ 6, 3, 1 ],6 Z6 v( G9 y/ L, F! T
        [ 6, 2, 2 ],; o' v! x! j9 y) t+ W
        [ 6, 2, 1, 1 ],4 b; B/ m9 w; t/ h% f% o( ?
        [ 6, 1, 1, 1, 1 ],5 I; Y5 p( B+ `+ H  D& C
        [ 5, 5 ],
    5 H" L9 e; J  `3 K- F1 Y: O, j    [ 5, 4, 1 ],
    " `4 Y! h$ b9 N9 p% `    [ 5, 3, 2 ],, S7 d, [* D$ J
        [ 5, 3, 1, 1 ],/ A1 k1 S6 F7 e
        [ 5, 2, 2, 1 ],
    0 N! i4 b* v8 ^' C0 G1 A5 G; X    [ 5, 2, 1, 1, 1 ],2 X( d6 f7 |6 l+ T4 s9 C  T
        [ 5, 1, 1, 1, 1, 1 ],. L; R4 y, x3 ^; R( ?5 E
        [ 4, 4, 2 ],: ]7 h: u, ~8 U/ g6 w- F, e+ ]6 k
        [ 4, 4, 1, 1 ],$ _+ y! J7 Y; g
        [ 4, 3, 3 ],5 y  X4 |, e; u( H/ Q$ I% I5 Q, j1 P
        [ 4, 3, 2, 1 ],0 e8 l. q0 F& ?9 s8 u7 n
        [ 4, 3, 1, 1, 1 ],
    , T. p! A2 o4 H# e+ K5 o, ^% x    [ 4, 2, 2, 2 ],! C! e% ?" \9 P: X( k+ k
        [ 4, 2, 2, 1, 1 ],
    & v3 b3 o. z0 q4 `: B" y    [ 4, 2, 1, 1, 1, 1 ],) e- d2 Q$ C& f0 U
        [ 4, 1, 1, 1, 1, 1, 1 ],
    7 ?% h) v: P( H    [ 3, 3, 3, 1 ],# J# I8 k! ]# B
        [ 3, 3, 2, 2 ],& u' o9 l" X% o5 ^2 t' d7 |* w
        [ 3, 3, 2, 1, 1 ],( b, I" j5 |9 F9 N3 B
        [ 3, 3, 1, 1, 1, 1 ],2 p. q8 ^6 K& B# c$ k* P
        [ 3, 2, 2, 2, 1 ],
    ! |6 X7 K6 O. I& c    [ 3, 2, 2, 1, 1, 1 ],
    8 q" _$ K1 p) G, l    [ 3, 2, 1, 1, 1, 1, 1 ],( Z0 |7 q. ~: \. W' |
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    1 g/ X% C, G# K5 ]5 o7 Y    [ 2, 2, 2, 2, 2 ],) {" F+ g, S4 x3 w: P3 B
        [ 2, 2, 2, 2, 1, 1 ],
    $ J5 F! c% Q% v9 l0 w) {* G    [ 2, 2, 2, 1, 1, 1, 1 ],- q, L  ?- W- b: h# h$ h
        [ 2, 2, 1, 1, 1, 1, 1, 1 ],% S  U( Q. d' ^, H7 [7 X+ Z
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],
    " {4 W8 E0 [2 `% F! W8 B$ \    [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]1 }  K! r9 G/ m- Z9 q0 y
    ]
    回复

    使用道具 举报

    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-6-14 05:26 , Processed in 0.548334 second(s), 92 queries .

    回顶部