QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5190|回复: 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 编辑
    ! a4 s) p, }, P
    2 |2 R# k" @3 c2 P6 r  yn:=12;n;
    2 ^. D# O: z& r- vFactorial(n);求阶乘4 c) W: n# K0 U1 H3 x. @
    Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
    ( h4 l9 p% B! f3 gNumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    ! D( U& s" |7 I7 N& ENumberOfPermutations(n, 4);
    , J! B: N  l' l2 jNumberOfPermutations(n, 11);+ V# Q" H( r) ?9 Z( Y( n) |( Y
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;
    ( r" p5 Z1 e: x5 lBinomial(n, 3) ;
    . G4 k# E; b# wBinomial(n, 9) ;
    , m) r4 z  [4 EBinomial(n, 10) ;
    . @$ _+ M& ?( c7 Y7 L8 y0 W' {; f9 JBinomial(n, 11) ;4 m1 X, ]3 _7 v
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    8 T. T& P/ a5 K# X) u4 a; A
    % X4 I: `" I+ o1 `* L) _Fibonacci(n);斐波数Fibonacci(n-1);6 t3 c4 r! K4 y# q3 x/ U
    Fibonacci(n+1);8 X6 z0 d# N) O+ |, Y  W! }
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;
    6 F7 _; r& R) R' g' z) ^( C4 t3 YGeneralizedFibonacciNumber(0, 1, n) ;, H& P8 H$ b, Z2 b+ g# L" [( [
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))/ z5 ?, S3 S/ h5 s* ~
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;1 h; o( S  F7 D. a. M+ Z8 K+ V% z9 M) W
    Catalan(1);
    , W! _* Q$ W- @( q) B" [Catalan(2);
    8 H; \/ S2 k1 UCatalan(3);Catalan(4);1 y" P% J! L; Q* \% c7 r
    Catalan(12);. s6 ^' ~3 }7 Y

    & X5 l& O6 u% ^9 T7 oLucas(n);卢卡斯数) T  ?/ f; r: i3 c* A
    127 {* W$ M+ x" d) |* B7 L
    479001600
    # y% V% e, h' ^7 C( k2 z831600
    5 e  C+ c$ f$ |121 r+ U+ D8 j2 y9 A% {; h' K5 x
    132
    : G: C+ k1 n( H11880& a6 L4 y6 c- H8 Q0 t* z
    479001600! ?, m* w( a* q" T( @  X) {. n& C
    12
    / P, i: \- ~, M5 y66, u" V7 x* I. I6 u- h, ]$ D* M) n
    220: z- S2 w/ e& m0 ?
    220
    ; l) }  ]2 C9 l# v8 J' c) g66
    - A% N. L8 t+ n; }12+ l: c- m5 o+ v+ M- V2 d7 Q9 t
    8316004 f9 @1 z# i8 B4 Z4 C- |0 D
    144
    , k8 h4 S. q. N0 R) w  i% n" E89
    $ `; a6 G- u, Q+ d+ O2332 o4 ~% w& k% A6 c5 A$ }- S% q
    233" ]- p" N6 v; j+ j1 u
    610
    - D+ j1 h. a8 }/ Q0 X, e# y& d# s144
    4 t( f8 D4 |1 ]& W) h- s208012
    4 h6 W" d# j2 ?% H2080120 x# b) `. ?1 M7 g+ k% P8 q$ {
    1* }" B6 _% q  y2 r9 E- P
    2
    ( [. }0 E9 t3 R+ h& ~5
    5 F0 U3 z( J3 _% _# k- M14/ E3 S- b* N% J* q
    208012
    9 ]* B- `. w; o5 N+ c! O322) k% \- {& j, O4 A/ J# ^% N7 L% C

    6 v5 P4 ]. o* }' s; e4 i! ^2 O! x0 s  ?" m9 Q& ~  B9 X
    卡特兰数=(2n)!/((n)!*(n+1)!)): N" `1 H+ y& F1 N0 I% s
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    7 e0 F# _+ d+ ~' u; E: J, K**YYY XYXXYY XYXYXY XXYYXY XXYXYY
    6 N# e6 E9 r5 z$ b- O. n! k9 K将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数:
    * V5 y. Z% b* Q((())) ()(()) ()()() (())() (()())
    0 c0 k2 g9 `% O' \Cn表示有n+1个叶子的二叉树的个数。
    & t$ y+ H- p. X# E4 d
    & b/ p; b/ |% b3 n0 c: K0 aCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    # x2 W& v% ^& m4 ^! HCn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况: ' T. p& o& C6 a; o- _' ?/ j

    8 s& q! f& q7 mCn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    + I4 a* t' |: 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. 6 P0 T: ?  }7 l0 I3 z; ]' y* @5 M
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    5 G% Y+ C4 }+ B* H" k3 R3 b) z9 P
    ! ~. u' y* R) \0 m2 L& x8 a
    . o$ ]: X; t5 o: h
    ; h$ p( L5 B: [' @卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。* k  i* `$ J' ]' V" [$ W3 b& j

    . R8 s  l6 K/ p: O- O+ }但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同
    * m7 z3 Q% W* [9 \, o+ i

    - O3 R5 }) _. P6 s: c, _7 ]# g5 A" O2 v0 r$ \& }3 R8 d0 T6 _

    7 w: \" w$ }6 G8 y8 T, pn:=100;n;, I4 [. V  u4 I8 S# `
    a:=Lucas(n);a;; h# n7 L5 j3 E+ A. S# [% s) X: C
    b:=Fibonacci(n+1)+Fibonacci(n-1);b;) s& m7 C0 E) `4 q3 F' T; z
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    & i1 n7 L0 }" ?) t2 o8 r8 }4 L9 f: |7 g4 n- W& Q, Q  ]
    1009 W2 u( u, J6 M$ q: p% Z
    792070839848372253127
    ! e% F( ^/ M. w! [% A$ @792070839848372253127
    ( L' {+ K+ v% [1 X! U/ C1771124240896309575375. l9 }  G/ J$ P
    1771124240896309575375

    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 编辑
    " f& k% s" P5 [$ F/ p9 F; _' Y, w, N& T5 {  K
    反费波那西数列反费波那西数列的递归公式如下:
      |6 p$ U* S9 U% e9 ^1 n
    * e) n, h2 C# ~1 i/ T0 tGn + 2 = Gn − Gn + 1 ) A) P' Z3 m- T$ \- @' S, ~% r9 s! [0 f
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...4 P+ h4 z' J* E! P) G9 X4 V
    + |" N1 J5 z0 n
    即是F2n + 1 = G2n + 1,F2n = − G2n。
    4 G+ b% M* C# s* T7 N: c: Y3 o4 S0 a- X  L2 d/ a! K: l% N
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    # m) s, D8 q$ @1 e, rStirlingFirst(4,3);
    3 q9 Y' r- ]7 g* h( s5 v- F& f# hStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
    / p6 W2 k* a( C5 VStirlingSecond(4, 3);
    ( e. d' a" D4 A2 d2 ~4 G: o9 y. s23 L: w& q, ]9 B; U
    52
    6 ^7 }2 b. a- _9 h5+ V. S) v2 ^9 e% m, A- Z
    15$ g% l* I' o; d3 }- p
    -6
    3 R4 e! ^5 g4 s% s0 f11$ k! e8 t5 m* ]3 |  L* |& i4 x* h
    -6
    * u" |  A+ l% m) p8 U1
    % `3 e; M& q4 w0 e8 E. }* ?$ t77 w$ |8 p2 {& q! G- T
    68 s3 T' N& \/ S
    & M/ v/ v, f7 B
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:+ D- o- u* A! s' W

    9 u, I9 f9 [! K( W: I+ N{{a}, {b}, {c}} ( ]+ r5 W/ B- k" ]. ]
    {{a}, {b, c}} 1 v: k7 Y; u# v2 n8 v
    {{b}, {a, c}} ; \# G2 b' J( e, ~, `+ N
    {{c}, {a, b}} * o4 s; z+ Q  u1 j" @2 D6 ~
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
    $ v3 B5 k. t, es(n,k)是递降阶乘多项式的系数/ z8 ?& v! e+ w0 F8 k
    有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)) `6 w( z& c" R8 V3 {4 H9 Y% W
    * b1 D/ d+ u) X/ R
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):
    $ d1 Z) }7 ?1 y! u  h6 n  O0 h; z3 @8 F7 }
    {A,B},{C,D} 4 j! a+ }+ F2 N2 P& y. g
    {A,C},{B,D}
    6 s. l2 u( {% R+ z  e* M{A,D},{B,C}
    * R7 S" \! ]% P- ~6 T7 O{A},{B,C,D} . e  X# ~4 r* ~- z! R3 H, }/ t+ y
    {A},{B,D,C} ) H/ u6 {& s8 T2 r0 f6 C
    {B},{A,C,D} 2 [% G# f; z- f( d
    {B},{A,D,C} + h9 h# `3 {8 v' q* `2 ]0 d9 i
    {C},{A,B,D}
    & X3 E7 k) d) S2 r* [4 B& L{C},{A,D,B}
    # j* C4 @- |* F% y* _. ~{D},{A,B,C} , q' _/ G6 ]( Q4 f) K! b& v
    {D},{A,C,B}
    - E' [6 `- \; g% K3 g2 B
    第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    ' p+ D, o, V5 A6 Y给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) & t# W6 l* S" n$ k2 w9 D  P
    S(n,n − 1) = C(n,2) = n(n − 1) / 2 6 ^3 W1 s1 K% ]3 K0 U8 n
    S(n,2) = 2n − 1 − 1 9 C, g9 f7 }, X$ q$ \" a( D

    - U- A# g# Z+ O( I( e" K换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:+ Z6 g/ ]' p( _/ }( d/ z8 Y8 ^

    1 ~( v: D/ x: p) W( }) }{A,B},{C,D} 0 S( }0 `% n) j5 `7 n  \, q! Z
    {A,C},{B,D} ( }/ ~  ?5 l+ r4 J5 ^: S2 d8 O2 s
    {A,D},{B,C} 1 R' w, {# T! ?4 t' t" ], ]' a
    {A},{B,C,D} ) t$ ~& `' K! I
    {B},{A,C,D} , ]3 s1 L; e8 T9 m. U: j6 d8 @6 e/ J
    {C},{A,B,D}
    - r% ^7 `0 ^; l  k$ Y* w7 C0 s{D},{A,B,C} ' \2 J$ p" {. z% 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 编辑
    * T! P4 A0 Q, S, t
    # @# ]  ^* L: jn:=5;r:=3;0 y) T! r( f+ P2 `) x7 J" v
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    * R. y+ K" g. y5 g  V! yBernoulliPolynomial(n) ;伯努利多项式6 ?) X8 G# j- Q0 o
    2 S) a& q9 C" f2 y" S" o/ w8 m
    26
    : m6 U7 ~- A: A0 q$ l! l137/60, _- l) Z$ [: Q2 A# Z2 X$ y# A+ t9 P
    0
    # @$ U! k2 o, Z. u! G0.0000000000000000000000000000004 y* c- o* t* M3 R9 R" A
    $.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 编辑 3 d6 S! V9 z& A+ d

    , T9 J" K+ V2 m$ a4 s6 X; k' D9 T" q( `8 o: ]  W

    $ ~1 M5 I% I! `( j. p8 L& B$ z3 W伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。" m( h' C+ C1 x$ h

    * p  U* C  b& z# x: X- u伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关8 p2 T0 u- I* E8 B! j* e6 c* H9 M3 K
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 * ]* t6 e( N" x! M5 ?2 I+ ^
    8 q- |! d- s$ j6 @2 |
    拆分 。。。。强!
    ' r2 F0 U! e' P) \
    0 C8 D! p/ ^* D5 c9 cNumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;: M" ~, j& x6 R6 N! m5 Q

    . V9 m$ O0 V0 W" Q. u. ^- f# @3 X7$ U3 Y; I$ F7 c. w6 b: B, A  m4 N
    190569292
    . ?' c$ {! ~5 I* N* Z6 d! H% S" C[" j  _1 M9 s: [; f  F/ i, Z
        [ 10 ],
    # ?* {+ A/ ]' D5 D( s% C) F/ O' K# `    [ 9, 1 ],
    6 G! L$ x0 o4 k- B! K4 K    [ 8, 2 ],7 G; w2 @- a, a6 i% k: H. w) y
        [ 8, 1, 1 ],- X, N+ j/ q8 @2 S+ z3 F
        [ 7, 3 ],
    ; N" T- K+ \' ~3 V+ h9 \4 P* u    [ 7, 2, 1 ],& f# X9 ~8 b" P
        [ 7, 1, 1, 1 ],
    " a5 h. n$ L0 Y0 }) x    [ 6, 4 ],# h* s& }. [  L2 x4 d
        [ 6, 3, 1 ],
    ! ?2 r. y6 X- `9 [. u- c    [ 6, 2, 2 ],
    7 u/ v1 p. p* a' O. i    [ 6, 2, 1, 1 ],) p" A+ @  r# E* J% E
        [ 6, 1, 1, 1, 1 ],- a. i# ~: ~$ V6 h- }
        [ 5, 5 ],3 S4 w( |! |6 u* q" B% r) J
        [ 5, 4, 1 ],( B. e. D& J9 w  A0 ~) k
        [ 5, 3, 2 ],7 t3 I4 W. B. [
        [ 5, 3, 1, 1 ],6 Z  \3 u! D9 e6 t! D
        [ 5, 2, 2, 1 ],
    * W8 G8 `4 z* m  A( [* O& D1 N* n8 _    [ 5, 2, 1, 1, 1 ],
    * p! `' x: v0 B! R$ N/ g, R7 T1 f* \    [ 5, 1, 1, 1, 1, 1 ],5 G6 ?4 m: H3 }) n0 E0 [  T) M
        [ 4, 4, 2 ],
    : y# ~, g- D: z    [ 4, 4, 1, 1 ],
    2 {  Z' U) R6 n% O) I" S0 C    [ 4, 3, 3 ],3 F  p" ~6 R" V+ X2 {2 `9 j
        [ 4, 3, 2, 1 ],( M. H# i; m2 _, x, P
        [ 4, 3, 1, 1, 1 ],% b( T7 L5 \& X8 o  i
        [ 4, 2, 2, 2 ],) |* I7 n! b3 [: v5 ]9 b
        [ 4, 2, 2, 1, 1 ],( ?( |. x8 x% s  K2 K
        [ 4, 2, 1, 1, 1, 1 ],
    $ ?  Q! a, I9 l" ?    [ 4, 1, 1, 1, 1, 1, 1 ],
    8 m+ F2 l% @6 q9 b& b& ?# H! K4 K5 h1 {: P    [ 3, 3, 3, 1 ],* h% q) v7 Y! O. o
        [ 3, 3, 2, 2 ]," J/ }9 i4 y+ V. E
        [ 3, 3, 2, 1, 1 ],  q8 L3 x! I" ?9 M
        [ 3, 3, 1, 1, 1, 1 ],
    ! S3 Y$ b+ B" ]4 {2 B" H3 |. [& k    [ 3, 2, 2, 2, 1 ],, ]$ Q* ]* e: F. k
        [ 3, 2, 2, 1, 1, 1 ],
    8 P0 l7 N" D3 q. X/ \    [ 3, 2, 1, 1, 1, 1, 1 ],& G$ Q! L% @) f2 g, d
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],8 d- i1 A& Q: S0 ?6 O7 g4 ]% x
        [ 2, 2, 2, 2, 2 ],& }2 K* n; S+ K6 o2 O
        [ 2, 2, 2, 2, 1, 1 ],
    # [4 P( H7 M/ m, ~" D) W) X    [ 2, 2, 2, 1, 1, 1, 1 ],
    0 E& R4 W" I8 K6 f  |    [ 2, 2, 1, 1, 1, 1, 1, 1 ],
    ( A9 Q+ @9 n" \8 s' q6 C# G    [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],
    3 b7 D' S, `/ l9 m: n    [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    ' p6 b) B: \0 T* d- c" k9 ]& C4 p]
    回复

    使用道具 举报

    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-12 09:21 , Processed in 0.504155 second(s), 92 queries .

    回顶部