QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4214|回复: 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 编辑
    ; S5 }" T0 K. V; A- Y2 e5 s  w4 M
    n:=12;n;% H& h5 ^7 J( {* K  e
    Factorial(n);求阶乘
    % c- V7 P& e; V  w0 r6 w3 oFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));% H) @% R1 l* Z: N# R0 E
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    3 c, b+ ?8 X3 Z; @2 B) q' ^NumberOfPermutations(n, 4);
    , m! T" \" }' `- f; f# BNumberOfPermutations(n, 11);
    & P; K! g. J& e8 M8 j4 U; A( oBinomial(n, 1) ;二项式系数Binomial(n, 2) ;
    " H/ `8 h' d) W' E" i+ zBinomial(n, 3) ;9 E2 o1 D; N6 I% K
    Binomial(n, 9) ;) r' t) L8 z, c0 _4 J& E
    Binomial(n, 10) ;% U. p1 @/ Z$ L9 }, D# ^- B
    Binomial(n, 11) ;8 q" g0 v6 u& Z- {) x9 u+ A
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600, c' g5 O6 ?  B; z

    ! J& F5 U* q0 l* N: n# eFibonacci(n);斐波数Fibonacci(n-1);
    9 f4 `- ~; N0 X. u. k% WFibonacci(n+1);7 ?$ [5 u8 y: q+ H/ P2 a' f7 I
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;- j1 C5 Z5 ?, U. l6 v
    GeneralizedFibonacciNumber(0, 1, n) ;9 c. x& F: u) _- _" w
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))& }( ^9 k& I0 P7 T
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    ! B2 y* D* X# c- FCatalan(1);
    4 v2 [4 T1 M; h6 sCatalan(2);- g8 S* U+ g4 c$ Z; r# f5 p8 D- c
    Catalan(3);Catalan(4);
    ; Z/ D3 E% E; Z, K# ACatalan(12);
    / G1 e# N; W6 {2 b# Y7 L6 [
    5 ~. H3 L" m3 {7 A2 XLucas(n);卢卡斯数
    7 I' ^6 j1 R+ A; N4 ]12
    6 k! I% ^0 U( E3 R0 p479001600
    + u) {/ m8 j  N8316001 K0 q- M" O9 Q) [
    12/ Y9 [" b; H" Y& C% O& J4 b
    132* G2 l0 R9 X! v( S& C
    11880- g( M; V( L; ^9 f4 y
    479001600
    5 K) i4 H  N: t( A4 d128 @9 A4 A3 L" j% x9 T  D0 e$ C
    66
    9 U4 `: w* J2 q; `' c* U220
    $ A6 ?+ G% k; V4 X: B2203 O. A4 i# g9 \% K: i
    66- s0 F( a  u0 ~. P
    12/ I0 Z. m( v8 ^5 _  W" x
    831600- K2 Y9 i0 w) N( _. k
    144
    5 p% J3 ~; L' G; k: p89# j$ {2 A6 f% q& f
    233
    * E- M% z4 h" N233
    2 b1 [4 c4 c8 z+ J3 ^. ^610
    ) R; u( o6 u( N/ K144$ s! |) ^$ ~5 [4 ?1 f! y# B
    208012) u+ P( ~7 r: N. S( t+ {+ n5 r- G
    208012
    & P8 w+ `4 W4 o' o19 I" m+ m$ ^4 t2 x
    2
    " W" I" i, d8 S& ]. e% K4 \+ W58 x# u8 h* m( f3 f7 j. ^5 R0 h
    14
    # B6 s9 V! J- A. a$ w208012" Q- O" R( @6 Z. Y4 ?" d
    322
    $ W& I" x! {& }9 ^2 {
    % L. [: o# A- W% U/ Q2 _2 k, E- u' A5 [& P% k
    卡特兰数=(2n)!/((n)!*(n+1)!))
    2 W! O  `! X" K1 F+ }) _  A- x# \Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: 4 s+ x4 D& `5 p! }' t
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY$ ?& \" J& t; M3 I4 x
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: , E  F$ K, x/ `7 o; l
    ((())) ()(()) ()()() (())() (()())9 z6 S+ x7 q. n6 w) y; ]1 Y9 s
    Cn表示有n+1个叶子的二叉树的个数。
    0 p8 v2 o1 o$ K8 F; r) q! M+ g( f8 I; x& O7 e
    Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) 5 L& ?8 X1 a" U! d5 J# ]( L: @
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
      L$ Y1 b! J% M6 H; q# e4 K! x; h8 P! ]/ m8 D' b& y: w: b4 @
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    ' L/ ?% ~, m- A4 BCn表示集合{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. 7 \$ k9 b1 i2 I# Y- ~: U
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:
    * g% k4 k% D* C

    . z5 H0 d$ Z# W* @! g. w* B9 P- D6 P) b1 d
    & S' D. I, G0 m6 z7 n
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。
    # U% B4 H$ G. u. i2 ]
    , {" \8 `: ^- G. f; |但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同
    * f% a4 {2 J; f/ _  C, j" w0 ~  a

    4 q- f/ Y# I7 F* ?! c  N2 J$ Z$ Z. F0 X1 l8 [" `0 K
    3 _  Z; \2 ^& O) E+ o$ l
    n:=100;n;
    2 @1 O6 C( `: a3 b7 F$ ta:=Lucas(n);a;
    - V$ O2 n8 v% |9 y' ~, J) n( V$ jb:=Fibonacci(n+1)+Fibonacci(n-1);b;
    " t/ k- |* C' ^' p; p: p! ]" vLucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    0 q9 w! l" z, m! Q, ^, N* v2 d$ v) U: |7 I& y
    100
    . Y: p3 O$ k- k0 m) `; y  k/ X792070839848372253127
    " ^2 p6 N/ f1 u$ x+ K6 M792070839848372253127/ M, H+ h3 H7 ?# R: T$ C
    17711242408963095753751 L& }; Z, P5 g# |* S. }  N9 t
    1771124240896309575375

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

    12.JPG

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

    74

    主题

    6

    听众

    3283

    积分

    升级  42.77%

  • 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 编辑
    9 [; L) \$ _9 r1 i9 d$ R7 m5 b- Q+ {( y
    % C( p4 u( X% B0 h反费波那西数列反费波那西数列的递归公式如下:
    1 G  U  d; R- ~3 a4 f1 p8 Y! z" z
    Gn + 2 = Gn − Gn + 1
    ( @2 d! o* t  C9 H* a! ~  ?如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...: v# M& u+ w) H6 j
    / y9 h) Q# ^0 i
    即是F2n + 1 = G2n + 1,F2n = − G2n。
      m) {' w7 Z& ^, A. z7 o; }, {! {
    : @5 U' b5 q) r7 R# Q8 u! I5 xBell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    5 U! f" J4 N$ ~! z% HStirlingFirst(4,3);
    3 x) r9 O8 f) RStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
    & {0 [" J" P4 T5 y9 h$ `8 O* z# lStirlingSecond(4, 3);
    3 {$ m' v3 ?' e0 w6 s: |) x2 R2
    8 H! [) y2 u/ u  ]/ L520 L7 g  Q# n6 {) h4 s
    5
    ) l  V# f1 u6 |) j0 E3 {" k1 p9 I15
    9 Z" m) c2 I8 ]* u-6
      R( t. S2 n& ]3 W& }0 w! M# Z; T119 I  |( `( L% B' `
    -69 f) S- Y; x* Q! w
    1
    4 K$ j7 Z7 ]1 C3 ~- b  d7
    + N' y9 P5 l- G) z3 o; ~6
    9 ~4 u  z+ O6 v4 D, d# t& A. x2 N$ W0 U) V! l9 {* Z( r/ R! a
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    : _% k2 x! z9 }8 Q8 w
    / F$ {$ I" P' f) g{{a}, {b}, {c}} 8 V" W* w! |3 S" f& F+ a
    {{a}, {b, c}}
    # D9 m2 `! v- d! Z# Z0 k{{b}, {a, c}}
    1 j) x( b/ [7 Z: @3 t* c{{c}, {a, b}}
    / p# E; I7 H7 |! ?) w( u{{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s& d: F* C" i3 Y$ v/ N0 I$ @
    s(n,k)是递降阶乘多项式的系数8 Y5 o* F" s1 @$ ~
    有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)1 I) p1 ?2 Z. l8 X* f3 H: u$ `
    ) ~5 ?2 e+ H  T3 e5 b+ U( `5 k
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):% m2 N8 k0 ?( ]5 e5 a* _5 Y

    4 K) R# `# h* i5 g' l! ?2 t{A,B},{C,D}   R( J# u8 m- K3 w. ~2 [
    {A,C},{B,D}
    % f$ z, B  e6 {9 G( q. ~( C1 P- E{A,D},{B,C} - ]$ W( N; g3 s) v% {
    {A},{B,C,D}
    . G1 _" E1 V" F: s: X$ A& t- j- S4 s{A},{B,D,C}
    $ ?+ ^! @9 v" a4 c6 Q{B},{A,C,D} 1 ]9 S7 [/ H/ i' p/ o' o
    {B},{A,D,C} " Z' y6 L/ E9 `
    {C},{A,B,D} 7 d8 ^2 S1 T" @# h8 x+ r
    {C},{A,D,B} 3 q! u7 k* n' e* ]$ E) ~' l: G  q
    {D},{A,B,C} ' l7 m8 ]9 h# v" ~
    {D},{A,C,B}
    ; n5 g) Z  H$ _. b1 j7 r, D
    第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    ' Y" L$ `' y3 W3 T; i7 ~3 t给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k)
    6 k* a+ l9 q6 @( wS(n,n − 1) = C(n,2) = n(n − 1) / 2 0 V. f( B& S, i" F# j0 ~
    S(n,2) = 2n − 1 − 1
    & {4 y& J2 u8 `" ^$ C3 `' a8 r, N7 w; A3 N
    换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:6 A* E! V. O+ @7 O

    % j- Y7 ]( z8 {. G% C0 m{A,B},{C,D}
    : E2 j0 U! r7 H) |( z/ k{A,C},{B,D} ; L' \0 Y, S" T& z0 W/ J/ b# y' N* v
    {A,D},{B,C}
    : ?5 c) J6 T% z5 E6 _9 ^# y- x{A},{B,C,D}
    5 Q- @) Q' \. Y+ @# z1 D! |{B},{A,C,D}
    2 y1 C5 i9 @/ G5 z4 \{C},{A,B,D} 5 D: @1 \" A/ k, W
    {D},{A,B,C}
    # G# `! K- L# J: T' M- r因此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 编辑 2 f- C  I" s0 N% a
    + W$ @8 ~$ D6 P' L1 P( O
    n:=5;r:=3;7 G. N4 \9 E7 Z+ ]
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    7 W3 I6 n8 ]2 W* p' n$ xBernoulliPolynomial(n) ;伯努利多项式
      E8 S' b5 b! |/ n3 k. j& c
    " Z4 O5 v8 }! s; w" F+ V. t0 g* l26
    7 S$ F: \/ g" D4 p8 }! Z137/60$ r- `; J; Q/ k2 W
    0; ~# O. X7 L+ i5 _$ I
    0.000000000000000000000000000000  l1 R2 r, E# p+ y. i
    $.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑 1 v, L4 I, V) h
    8 C8 o5 v* L" B* j
    9 _, i3 W) _. H7 q8 B! a

    5 E5 E; o# c! e4 q$ U# x伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    & Q5 P( e1 i3 f* e  F) k8 U2 ~) U* C3 v
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关3 }- p' u; u- _& f& q) G
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 ! x3 C% g# ~. f! {
    : \! N- c% f% X, p6 I4 p
    拆分 。。。。强!8 r' T& d, o. O% X  Q2 ?$ B
    " v+ H/ ]& \1 Q4 w% T8 T
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;
    $ F. H  e/ T  i2 A1 y
    7 _# p( g) p* s3 F3 V70 Y1 A! b3 V$ @
    190569292
    0 O/ P; [& i$ U8 A6 \# Z[
    8 f3 t6 o& h$ b) F. d    [ 10 ],
    " c# a$ l/ }) l5 q4 h    [ 9, 1 ],
    ; k* D& O, J5 @: G* ^    [ 8, 2 ],
    8 T' w' c$ E. z; U5 \& [1 d+ q& i( C) E    [ 8, 1, 1 ],
    1 B+ v( w! l4 {5 e1 I    [ 7, 3 ],
    % F/ Z3 Q: Y- _! A% w    [ 7, 2, 1 ],- S5 D6 a2 M  |
        [ 7, 1, 1, 1 ],
    9 E! o, q5 E. c6 q) \    [ 6, 4 ],; U. B" |7 @8 K2 p4 q3 \& l) X, t
        [ 6, 3, 1 ],
    1 H! J. S6 K4 I8 k; u9 \    [ 6, 2, 2 ],4 ]5 d# Y4 z0 P; ~& S
        [ 6, 2, 1, 1 ],( a8 O5 X9 y+ A$ ?8 Z' @% X
        [ 6, 1, 1, 1, 1 ],& ^; J, `  f. M: r# X; B
        [ 5, 5 ],
    * |) r7 f3 r1 f% h6 z    [ 5, 4, 1 ],
    8 R6 R2 u4 r! h8 p7 Y    [ 5, 3, 2 ],/ y' R7 u- ~& \! t% v
        [ 5, 3, 1, 1 ],0 l2 W3 S- H+ g2 v( F
        [ 5, 2, 2, 1 ],# g& S3 ~! o0 u# ?6 B/ H
        [ 5, 2, 1, 1, 1 ],+ g  U9 a$ y, A9 r
        [ 5, 1, 1, 1, 1, 1 ],9 D, U0 u; L* k: P2 P* ]! f
        [ 4, 4, 2 ],
    3 a+ A5 m6 X% j% N; L3 M# o% S    [ 4, 4, 1, 1 ],( ~. ^; c; w2 S# i
        [ 4, 3, 3 ],( t5 I$ F" u' `9 C& N/ H
        [ 4, 3, 2, 1 ],
    $ u+ f4 x" E, n7 T/ |: R    [ 4, 3, 1, 1, 1 ],
    5 Y& W& G9 ^  X  Z* g. D4 t5 M    [ 4, 2, 2, 2 ],
    # C, B$ j; c% }9 w) ~    [ 4, 2, 2, 1, 1 ],' u7 U# w1 M- u0 q8 x& N! v
        [ 4, 2, 1, 1, 1, 1 ],% C2 c3 Y$ Y8 d# ]8 h+ o- ]
        [ 4, 1, 1, 1, 1, 1, 1 ],
    & T/ M2 Y( O$ a    [ 3, 3, 3, 1 ],
    4 _6 d7 C4 A% A: s* P: l. v    [ 3, 3, 2, 2 ],: g- K8 |9 {: {0 h: m9 S( T8 E
        [ 3, 3, 2, 1, 1 ],
    , Y5 {- ^; f$ U. o' t! z/ q    [ 3, 3, 1, 1, 1, 1 ],* F/ y+ P5 Z* [+ a, t. T% l- h
        [ 3, 2, 2, 2, 1 ],
    ! [# A2 h" r4 S7 x+ m0 }    [ 3, 2, 2, 1, 1, 1 ],
    ( p3 B$ ?0 X0 L' t8 I4 L    [ 3, 2, 1, 1, 1, 1, 1 ],
    ' Z, z  y* n% K8 g    [ 3, 1, 1, 1, 1, 1, 1, 1 ],7 m$ P  }6 c( U( c
        [ 2, 2, 2, 2, 2 ],+ O6 ?1 m- _8 B4 \* r
        [ 2, 2, 2, 2, 1, 1 ],
    & a' p. i3 v6 X    [ 2, 2, 2, 1, 1, 1, 1 ],
    4 V2 l8 q" `+ u    [ 2, 2, 1, 1, 1, 1, 1, 1 ],: t  |( i6 Z7 [8 i4 ]" q
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],
    0 T4 l; T4 G# G6 s    [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    . x$ W( T) G+ S9 V1 a]
    回复

    使用道具 举报

    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, 2024-6-20 01:13 , Processed in 0.535960 second(s), 89 queries .

    回顶部