QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5195|回复: 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 编辑 % v! K& B' j  T5 O- [, h
    1 o6 p6 I; v& Z' y
    n:=12;n;
    3 J+ |2 |, u& K  e. d0 t! P+ ^1 n; M( ~Factorial(n);求阶乘
    % ?# d0 |0 Q$ J4 RFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));7 n; f" D% V. s3 t
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    ! l9 ~) D. h" zNumberOfPermutations(n, 4);0 ~9 |3 B7 e* `7 F% ~
    NumberOfPermutations(n, 11);
    9 v0 o$ K. ?5 b) OBinomial(n, 1) ;二项式系数Binomial(n, 2) ;7 U, G) a, h9 ?4 P
    Binomial(n, 3) ;4 `+ n* L+ m) |5 \- o
    Binomial(n, 9) ;) J' _" s2 Q( H: {
    Binomial(n, 10) ;
    4 z: ~; Y% I, k. u0 u( {Binomial(n, 11) ;3 z0 O- c3 |2 e8 q3 D
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    / P: Z4 y  G/ v' x6 ^6 P2 Q; m% V3 G6 v  D
    Fibonacci(n);斐波数Fibonacci(n-1);
    / D" f0 |2 I- q( w1 M% bFibonacci(n+1);$ F- g$ q9 y3 D7 X. x" i
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;
    * q  B& {, g# o6 l+ s# z3 e. r/ I; iGeneralizedFibonacciNumber(0, 1, n) ;
    1 n5 m2 V6 S' ~; S* N& }' jCatalan(n);卡特兰数=(2n)!/((n)!*(n+1)!)); Q) e4 u! L$ b  B! K- L* }! C& u
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;& n, @, L, E  k; }+ i
    Catalan(1);
    " L3 I2 b; V; e0 DCatalan(2);8 N& o+ R$ n7 {, {
    Catalan(3);Catalan(4);, T8 j7 v" F; y6 D, z2 ~- [: {
    Catalan(12);8 I/ r$ a* ?' g. G3 N- C
    + C% j) F1 s" f7 C$ y7 t4 [, p
    Lucas(n);卢卡斯数
    ! e4 X, k; p4 U" V& e12* N3 L% g( i# v( v
    4790016002 I. _5 o+ ?" R3 J' j
    831600; b$ I! z  l9 h0 e. ~1 ~' i4 L
    12
    7 N9 ~5 [# d$ e" G9 J+ b$ N; o% Y132
    ! y. q# Y- P' b, N) [11880
    " s8 }& [$ F( u! P4790016004 s: ~) Z( u% T5 R& ?9 z; l
    12
    + I3 w6 e' J& f$ ?' _66
    . P6 J7 p2 f' f; q. ?* b2202 z+ T% Z4 z& T
    2205 @/ a. q4 u1 [
    66
    / c  o. D7 U2 M8 S12/ u: N+ P) a" ]
    831600( p4 ]* F8 d! M/ w( b- s/ W
    1448 [7 P. U; z8 Z+ }0 H( ~
    89
    ( Z! ?" U/ U0 j  y) s0 v# [233+ @! S5 d" g: R
    233
    ) s4 I) N  x! R( V4 @610& S5 A" u+ K2 X9 G; w
    144- e" S% h2 E* r. a8 g6 N4 ]
    208012
    7 J9 Q6 C1 }1 e/ _6 X( h$ ~0 p8 L208012
    ; \8 }# Q8 u& L8 T' _, W3 m1" I5 S  A$ s# \$ i* D
    2
    + r8 a0 R: l3 k' X55 f9 ^. s# N: W
    14) ~1 K- `7 s; |# o
    208012
    / B! B- f, A, |3 r0 r/ }6 O322
    1 v1 `" X+ O+ E. D/ ^5 J0 o6 J) U! z( l' l' q

    3 x0 W! b. E+ u% x8 O/ O8 ^* ~  ~卡特兰数=(2n)!/((n)!*(n+1)!))
    # M4 ~2 @9 G4 W: e, k; c. vCn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:   D; P/ b+ a$ a5 J- A  Z
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY2 a4 {% y! I& K1 R
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: 5 X# L9 c1 p6 ~& |' r# b
    ((())) ()(()) ()()() (())() (()()), i# D0 f0 U5 K# H! f5 T
    Cn表示有n+1个叶子的二叉树的个数。 $ G0 ~) q0 t' {* M

    8 F; j6 x1 d' w, K$ l  _7 sCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    * n: w! g1 `+ u. mCn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    ! O+ B2 _: l- Z4 L/ K( @* l
    , P  O3 u+ @+ R  X& @: n0 M; P/ B: LCn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
      Q& K$ b. f- s! FCn表示集合{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. / ^& R. H- o: j& _
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:
    9 I4 p; T0 _% b1 i: W
    + n" U0 u% q. e' C
    6 \# \& y6 N* P
    * J- z7 ]9 p, `: Z* O1 [4 t/ w
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。
    0 F# I/ u5 u0 Y, }
    $ N& W+ t6 c9 h- M' ?但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同
    / P: `! m9 w, {2 R

    . k  B) ?) P' E; M6 b" A" O
    : x0 u+ S" ^3 E# K9 I
    " j+ S$ Z, s9 y9 x2 `0 {4 Nn:=100;n;: j  X& Q1 p9 ?
    a:=Lucas(n);a;
    9 X; q: P0 J! _$ T; D. T3 Zb:=Fibonacci(n+1)+Fibonacci(n-1);b;" q1 H1 T9 p( W/ M& K& z, f( @& E
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);! [0 M! f; c: i+ V

    ! d5 G8 M& y# u0 f) j, n1005 q; C+ X: G1 o4 g
    7920708398483722531275 j3 W" R  f2 n8 U- V
    792070839848372253127# `" c; m# w3 V
    1771124240896309575375/ o7 B* U; Q9 X2 B3 F
    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 编辑
    8 A+ [3 Y; n* P, o9 o9 r1 W
    + ]: d* I& G: |' [0 s反费波那西数列反费波那西数列的递归公式如下:4 H4 j" H( |" M8 X1 }: o
    3 a6 g) f) n  J' r& d6 \- y- u
    Gn + 2 = Gn − Gn + 1 ( `, v% k' v  W) F8 j" F
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...: a2 T2 D9 {' r) C, N& E
    # p9 s) }7 Q7 @0 F; k: `: ^
    即是F2n + 1 = G2n + 1,F2n = − G2n。0 e2 C, i$ V* q8 x" w& F2 _. p7 M
    4 E6 g& j$ N1 y2 A/ b, f0 x, U
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    , H+ i" d7 \9 z4 T$ `StirlingFirst(4,3);
    5 Q7 L: f7 g& ~9 ?8 f" AStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
    " M0 }9 g- N* b% S9 BStirlingSecond(4, 3);1 p# e% i8 T3 n2 g
    22 t! |4 R4 Y  G9 p' s4 l
    52& y2 d/ f6 h& K' A7 c8 o2 [
    5
    7 W. \" @7 R9 g15
    / Y+ Z* l) t+ \* A: G-6
    + M* _- H" |4 R( y& b- V, i. t" q11
    9 _* w9 o. Q6 {. }* w2 }8 @. z-6
    / J3 y& k9 D4 c* Q2 X1 j* P1
    4 M9 R1 i$ X+ X2 ^7# W& t4 z. D( P! Z3 `' u
    6+ D0 Y% R# W* x  O! j
    ' F$ g  @2 L7 i, o5 n# t$ w
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    ) q, ?; A! f8 H& c3 m
    / f5 V3 T  t# ^2 f$ x3 t" I. b{{a}, {b}, {c}} $ D; a' H% |9 C3 |* N7 J, r
    {{a}, {b, c}} 1 j; ~& Z* D+ y  S! j+ ?
    {{b}, {a, c}} 0 s! J% t, l5 {2 Z7 O" F
    {{c}, {a, b}}
    5 O3 j1 \% q4 R9 Q; M2 I* y{{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
    3 |; u$ ], m8 f. D. ts(n,k)是递降阶乘多项式的系数
    . M5 ]+ R% @1 H1 j有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)7 e& m2 v5 Q/ q# ~7 A9 w1 s

    " c0 n+ d: w9 O) F换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):
    , F% S9 M! h( \) W: \& ^. A% }
    0 @. @; p7 k1 ^{A,B},{C,D}
      V  O$ r9 v2 j9 q{A,C},{B,D} 5 J+ d* C# y4 H6 y, w
    {A,D},{B,C} % `& ?. j* U, |( f, W
    {A},{B,C,D} # T9 ]' c: V6 o
    {A},{B,D,C}
    ( P) }8 i* x4 i4 x" X! t{B},{A,C,D}
    ( o2 g* j' U; I& I. b% `{B},{A,D,C}
    ! @0 v$ Z7 ~2 z% i{C},{A,B,D}
    + i& \( @+ B. _8 O+ ]; z$ ^{C},{A,D,B} 2 i0 y  A: M$ u) d, b* b  t' _$ L
    {D},{A,B,C}
      w! E3 G3 b" q7 Z# b% S4 ~{D},{A,C,B}

    : N" E9 I0 i6 a4 s+ _4 o第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S: s% G5 ]$ J9 e/ T! B* r
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) + o" K! v" x+ s# {" }
    S(n,n − 1) = C(n,2) = n(n − 1) / 2
    % a& I0 _/ }/ q2 ?# D3 ]S(n,2) = 2n − 1 − 1 4 {2 l3 Q1 G# [. }6 L

    9 w: u( p0 Q* h) B; S+ v换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
      y5 |. Q! I. w' Q% t1 w
    9 c0 g( x5 b. K' q$ v{A,B},{C,D} ! v: X3 r. v, V, N! P
    {A,C},{B,D}
    2 g$ N9 @3 P7 D" T% Z# F* O{A,D},{B,C}
    3 ~! e7 Y& _! H{A},{B,C,D}
    ' N" m7 q! @# w, E{B},{A,C,D}
    , r# k- \, G% O1 Y{C},{A,B,D} $ j: V: N8 f, H* P6 I: b
    {D},{A,B,C} 6 G7 r. r# ]! w, H
    因此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 编辑   a. }# N! Y/ ]) i' {- W
    , N! j# L) r: t, F6 r
    n:=5;r:=3;
      K/ @( y) ?' m) |; GEulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;3 s3 |' `8 ~; @1 [4 w) G* _7 `! ^
    BernoulliPolynomial(n) ;伯努利多项式( A* u1 k( I" I
    ( {0 r% V  H$ M6 p. i2 m# O1 Q3 M0 T7 P
    26
    : z) ^& e  U; {! W137/608 g! K# A3 x1 j* X( ?7 N% N0 l
    0
    / f! `7 c1 B- Z1 C) |4 }2 Y7 o0.0000000000000000000000000000002 ^* K5 j& `; l
    $.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 编辑
    & t! v+ b' m! }7 W  Z  l% a! ]& B5 q( j
    - \+ L' i& F- n& C& O' V9 r

    ! e2 P9 q% b: e- Q) ?伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    5 H8 [$ P  [$ k1 i8 O
    6 Q3 x4 h* k: T9 t8 o$ L伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
    & Q; W( q% _' H; k. \* E
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 ) ?1 e/ `- n$ y" i
    $ k% _7 \" }# V2 u4 H9 N: O0 O) A- z
    拆分 。。。。强!
    2 @3 n: \) K; ]1 k1 Y/ \* N. f, M
    $ @% \+ G: `2 M' k0 A, hNumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;
    4 J  ]9 V) w1 m) n; Y0 r1 ^$ t
    4 x" v4 t& L" i/ s4 ?2 C7; M& L) _( `- i0 m4 X$ A1 v
    190569292( Y# ]" y" F  [7 U
    [
    2 c- I9 O& m2 }. g) u    [ 10 ],
    ! e2 }/ T. A: U* w4 f9 }( t9 s2 v. ]    [ 9, 1 ],
    9 u( g$ J# ^/ [; I+ x    [ 8, 2 ],. w1 X/ K9 y* _
        [ 8, 1, 1 ]," @2 C$ \0 }0 o0 {
        [ 7, 3 ],
    : F1 t4 Q2 I) m! h5 n% N6 V    [ 7, 2, 1 ],
    * I& s3 p0 Z. J" \; C1 c    [ 7, 1, 1, 1 ],
    ) [& O" D# E/ E( L) `9 r    [ 6, 4 ],
      j: `5 |, e4 M, F5 h; {/ q    [ 6, 3, 1 ],( u) U/ L9 A2 t; L0 X6 k# X; ?
        [ 6, 2, 2 ],4 @7 i2 C* w1 }# `/ l& d0 m* S
        [ 6, 2, 1, 1 ],* {2 V* i6 z0 q! x) X" \  Z
        [ 6, 1, 1, 1, 1 ],
    ! P4 h4 E+ W- l4 T" x1 r, L    [ 5, 5 ],# L( R$ i. {# O* F' w. g7 i
        [ 5, 4, 1 ],  u4 Q3 X- o6 T* H) J( r
        [ 5, 3, 2 ],
    * [6 D- X' w" G5 G! L- s7 P5 y1 T    [ 5, 3, 1, 1 ],
    ; c' A  q% B1 a, D( S! g. ~. y$ Y* ^    [ 5, 2, 2, 1 ],/ }# u( N3 L' \8 \' C
        [ 5, 2, 1, 1, 1 ],; ?6 m# \% s7 z# @# A
        [ 5, 1, 1, 1, 1, 1 ],
    5 U& y  a+ Z1 ]5 {+ r. a- e    [ 4, 4, 2 ],
    : Y5 ]6 c% m9 A: j& ^) j    [ 4, 4, 1, 1 ],0 r4 K3 ^8 Y$ X
        [ 4, 3, 3 ],
    # O# }* f; s4 {  Z, v: P    [ 4, 3, 2, 1 ],
    / ]8 Q9 E9 p8 h2 j7 `% U    [ 4, 3, 1, 1, 1 ],( v$ P5 t5 T) i9 W8 _" f
        [ 4, 2, 2, 2 ],
    3 J+ O' `# S; A3 U$ [    [ 4, 2, 2, 1, 1 ],
    0 @4 t/ m8 k2 x  L    [ 4, 2, 1, 1, 1, 1 ],
    3 P9 u# r/ U  J) h; h$ d    [ 4, 1, 1, 1, 1, 1, 1 ],) m( l% Z1 i( U8 K  q" d, N/ i/ v* u
        [ 3, 3, 3, 1 ],2 m# S/ n9 L  n2 d2 P8 ~
        [ 3, 3, 2, 2 ],3 l7 y% A- o* U0 ~/ |) F
        [ 3, 3, 2, 1, 1 ],/ `- f! E. R0 k" ~& ?
        [ 3, 3, 1, 1, 1, 1 ],
    : j8 {! _: ~2 A7 G: I7 s1 r" D    [ 3, 2, 2, 2, 1 ],4 \  l9 i6 |) T$ V8 T3 Y
        [ 3, 2, 2, 1, 1, 1 ],
    6 v# m9 G; a8 c0 n/ t% m" g0 c( b    [ 3, 2, 1, 1, 1, 1, 1 ],, F) O1 h' N9 k2 W& l
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],* Z* i7 y0 X$ a9 B' A9 j8 P
        [ 2, 2, 2, 2, 2 ],
    1 y7 J5 u5 n) T; L, I0 n    [ 2, 2, 2, 2, 1, 1 ],
    ' o' r/ v/ l4 J/ V) M$ y8 w" i    [ 2, 2, 2, 1, 1, 1, 1 ],
    ! B0 D" B6 l0 m9 l  |; D( s- G    [ 2, 2, 1, 1, 1, 1, 1, 1 ],, p  ?' w$ C$ O2 ?
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],
    4 h5 p$ u- o; a) X7 ^    [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    6 p# r8 C! j/ g5 ~( A2 n# N5 o]
    回复

    使用道具 举报

    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 19:46 , Processed in 0.508325 second(s), 92 queries .

    回顶部