QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5193|回复: 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 编辑
    - q) z! ?4 b" b# r. Z  b, G  k1 k1 N8 e8 f) m! l
    n:=12;n;- |  V7 j% ^' l  {  C9 J$ O
    Factorial(n);求阶乘" J4 l4 j* {/ Q$ G3 b6 s; n6 [
    Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
    / @1 r, a( m9 o( R, jNumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    ! u/ e0 u' G: ?: C# FNumberOfPermutations(n, 4);" d" Z9 n1 e7 Y+ i, Q. u  @, f
    NumberOfPermutations(n, 11);$ `1 W9 ~0 f) |* \" m
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;* [" S7 m% p( ]- p6 a* d# V
    Binomial(n, 3) ;/ @! c; |$ o! F# y! Y% R4 C
    Binomial(n, 9) ;# s$ j5 t" B" ^, \
    Binomial(n, 10) ;
    . \5 r  r) q" t  u/ }0 m( QBinomial(n, 11) ;* b* \( O' m1 o' Y& n+ ~! _* m
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=8316005 Q% s# ~3 Y: J- z0 m" C
    ! J- d0 }  B' `8 h
    Fibonacci(n);斐波数Fibonacci(n-1);
    0 O' |. I" ~  C: Z- qFibonacci(n+1);- }) f) E9 C! k/ \* Q
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;! w) M4 s# U+ b  b6 Z; r/ e# e0 n. Y
    GeneralizedFibonacciNumber(0, 1, n) ;
    . V' G$ f1 R' KCatalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))5 A' \9 ]- K7 w$ t  n
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;! [3 l. p7 j4 J8 b) X! |; K  E1 e
    Catalan(1);! C; i6 b8 u* ~6 [: d& y
    Catalan(2);9 r( q1 E7 q- M7 V& s& Z
    Catalan(3);Catalan(4);. G% C+ N: B) c" J5 P
    Catalan(12);
    5 }7 g3 m" Y8 ]: G+ W  A& G5 p, P5 v& e+ \
    Lucas(n);卢卡斯数4 a0 @' H" P5 X; X! e
    12
      `6 o7 p+ j, D: E0 V2 _7 u3 O3 I479001600
    $ \1 @& W5 T7 E5 u9 m831600
    ( L+ x6 O  L2 X) `) G4 M/ e12
    . z1 w( u6 W5 {2 m: x" s132
    % P0 a* ?: |$ n7 \8 ?118805 \* A$ y: I+ f3 V7 @4 [9 l
    479001600
    ' H1 U' C. s# A( g- x. `, ]$ a0 ?12; P0 k  b, s7 z! i
    66* C5 }! ^  K' w+ x- m
    220
    / ], u8 d. ^4 O4 t; j2201 L. b4 k' x. q( d
    66
    $ }7 n& {$ _& n; ]. ?5 q& D12
    : {1 V, G; w* X8 `& V831600" S0 {/ {* E1 l  `4 q4 W  @
    144
    7 G, v8 y4 v9 R1 {  X89" h& m& I) K3 n' b
    233+ E) s0 t/ k+ D
    233
    2 e: a; |- R/ v4 b610
    / Q  ], M3 t7 N; v* [  v" q144
    % h6 F3 Y) ~) ~* v7 B" N# Q1 D9 K* H1 P208012
    % f) Q' [* m2 o  g! n/ A208012
    / ~& {$ H) o* o/ L, w1) S5 R* T4 {+ @  @
    2' ]) Y7 ]4 `8 C
    5
    1 ^7 f) N. P* Y8 n6 k( l14
    ; W3 j' x2 g1 h) R# u208012
    3 S( ^5 T/ U; e; k: l322% b" s$ |/ b. T3 V/ E( |% K+ X/ J
    3 F8 X5 J: h, N- L/ N' o
    - K7 `0 X, E  ]1 L2 u1 x
    卡特兰数=(2n)!/((n)!*(n+1)!))
    5 z) ~' k- y3 U0 MCn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    : b) C' C& T; Q$ M6 E3 n1 [/ V**YYY XYXXYY XYXYXY XXYYXY XXYXYY
    , `7 h4 U7 }5 e  `% E" ]- F) I将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: 6 n6 \$ S5 ?  r
    ((())) ()(()) ()()() (())() (()())
      [; l0 b! r- f( _; {Cn表示有n+1个叶子的二叉树的个数。
    1 w, M& J( [) c, _1 i1 u3 `% d- b# ]" e# h6 b5 d7 m
    Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    ) ]9 A' N' S+ }+ P% [' {0 |2 i6 \Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况: 1 m9 U8 g% _( B5 _+ M
    9 L& I+ f) v0 C4 r( p
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    ; R- Q1 k; K: }3 N/ [% ]Cn表示集合{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.
    9 d+ o9 N0 L; C$ W; _6 {Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    5 ~. {: _7 C8 i- E9 Y% n# W% l5 r: p; n0 i* P5 d
    0 K" _% z6 C8 V0 l( N$ M6 v3 U
    * j( V- \# ^% b) k  f5 Q
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。
    " t1 p: M1 S& t& T
    6 H  l, {% i$ _6 A但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    ( `) Q; G$ A, C5 }8 R" ]: d1 _  q8 w2 W

    - w. {3 e7 A& @. P+ D2 N, H, r" n! }1 A( g' L" |
    n:=100;n;3 H3 P. `4 o5 s: n  D2 [  {
    a:=Lucas(n);a;' c& ^' Y# j3 l# `' q, E* j/ r
    b:=Fibonacci(n+1)+Fibonacci(n-1);b;
    ) Q; m' R- R* _  @Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);, `2 e- c3 `& l2 F3 q3 x
    , E* p0 |  \8 ~" w
    1004 U8 I$ y$ Z& a( Y. Q6 |. j
    792070839848372253127
    / X% e9 [* j6 z; B. s2 _' W3 w9 J792070839848372253127
    , q, ]+ Z1 M# M1771124240896309575375
    6 p4 f: H6 K+ u; R% q1771124240896309575375

    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 编辑
    2 S* z! V7 p# G  ?7 S  K' l" g# K8 O2 v; Y; H5 V
    反费波那西数列反费波那西数列的递归公式如下:
    3 a3 F5 D# }  m* o, L( _2 q7 N* j% V# O/ J$ b) J
    Gn + 2 = Gn − Gn + 1
    3 J& Q9 v3 z7 \5 F& l如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, .../ H6 t* }6 ?$ f, k+ W

    , K) C/ Z9 ?3 z8 v( V: C2 M, o0 t即是F2n + 1 = G2n + 1,F2n = − G2n。9 X+ H* ]( a) `, S( Z+ }& T1 H. y
    # @6 w* l: R/ x; w. e
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);" ^: A4 A( q) v& ]' A, L4 q
    StirlingFirst(4,3);
    8 c/ A. k4 Z9 q6 ^3 c# nStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);  h2 z, I# F4 C. ?) J4 }1 Z$ M/ ?- A; N
    StirlingSecond(4, 3);
    ! z6 y6 ~5 k! |3 H3 p2$ ?6 @! R1 L& w  f
    52
    # S- y8 L) f9 m, M0 e7 |5
    1 P; [8 U) j9 {/ G  O15+ t! z" Y' a' k) y1 G
    -6
    7 R, G3 H" h6 u% x& M) N3 J11- i2 N% G$ A9 p6 Q  Z3 f
    -6; d# w) i' K) f# O
    1$ q* W9 V  L% m4 G  @6 b% {
    7
      ]. R& W4 q; C7 j) x7 v+ B3 E6 r8 A6
    ) F& {' G6 \$ i5 `2 _' l
    & D+ `6 a( V' p6 K/ U. N. _, f8 IBn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:# |) |) R# }& f& u& M
    , [) G1 M2 i. I+ \; W
    {{a}, {b}, {c}} 3 Y- W( V4 u, d2 q5 r
    {{a}, {b, c}} 2 i8 i5 X) ?- B2 e
    {{b}, {a, c}} ( R; O4 f9 x+ ]. v' C
    {{c}, {a, b}} % l' Y: I' }" Y' ~7 t# D; o) v
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s8 ^- w2 l) W& T" i2 u
    s(n,k)是递降阶乘多项式的系数0 p4 P4 |7 W9 b) B5 k1 ?* ^9 _
    有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)7 r% y- i9 b4 `0 L
    ) @: K. n& E& k" g. K& X9 b
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):+ Y. r8 E1 c% t9 Q
    7 s# K1 I) L5 h6 j
    {A,B},{C,D} ) O7 ]2 N4 ~5 d+ v* s# ~$ A
    {A,C},{B,D}
    . p$ J7 z, O, K, V7 F0 |" A{A,D},{B,C}
    ! g8 q  L: i& Q) c* C3 g+ [{A},{B,C,D}
    + R! {: w9 A* z  O{A},{B,D,C} $ _6 [2 \. H' G6 O
    {B},{A,C,D}
    : F" A" V- c" \0 u{B},{A,D,C}
    ; o' u4 Y' q. p) }{C},{A,B,D} ! m1 e( j$ m" g) o/ F0 Z/ R% {
    {C},{A,D,B}
    5 U" Z7 T. Q- x$ M  V! ?9 D) @{D},{A,B,C}
    * \4 A( I' n) K8 c% G. L{D},{A,C,B}
    - v& L4 g( y3 P; H' z
    第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    4 K. s$ v' A" [+ H给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k)
    4 G$ h% Q( z& Z) rS(n,n − 1) = C(n,2) = n(n − 1) / 2 8 f: Z' K* z6 p1 _& H$ o" J) v+ o
    S(n,2) = 2n − 1 − 1
    ) s( U+ z; f. G
      A5 W' r0 |2 \6 w0 E2 v换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:  Z( O* |( X! W. [7 J' H- m

    ( ?4 S6 L7 z$ q3 ]2 F{A,B},{C,D} , L( E6 ]8 q- u& @8 U  r/ r
    {A,C},{B,D}
    / o, e- S; V) @: \+ _# |{A,D},{B,C}
    9 ~: ~9 |- a- B2 ~+ ^. g/ e{A},{B,C,D} 2 Q( d8 C' L" |7 i, T* b. d& x
    {B},{A,C,D} ' s! `" O7 {/ B+ M
    {C},{A,B,D} ' A5 n2 i- f; Q8 J( ?
    {D},{A,B,C}
    : R7 }9 [% F; \7 s2 q2 f因此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 编辑   ]0 t( W/ U5 f! ^
    ) l- A* e( H) ]
    n:=5;r:=3;
    9 b( P3 M+ f, F) U1 S: lEulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    1 W! c% P( W! w  RBernoulliPolynomial(n) ;伯努利多项式
    0 Q4 f) R, h7 i* p$ X8 j
    ; }. N  z8 N& V4 ]- |26
    & R8 r* G  X# u/ t& H137/60
    & |: V( u. U% d8 v( n$ ~& R+ g0& Z/ p0 ^8 v  \
    0.000000000000000000000000000000
    # i4 \& d% ]7 P" i8 P$.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 编辑
    5 V( G) x! j4 r! ~, Z
      ?8 g- X% v; p6 c/ v* J) r4 X( j5 A, a' u- p, N  X5 Z

    % K/ X# Y0 v2 }3 j5 ~伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    & t0 ^5 L+ Z; I; d, V( e$ t7 Y5 `' w
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
    & w& ]: V( P, f0 ?9 w
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 ; w: S; Z" W0 T; z

    ( [6 D/ D+ V* }拆分 。。。。强!( l( E# z/ g& S9 ^: k$ {- P8 d" L2 t
    7 v7 [2 k: X+ d, J
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;4 a7 D% t: _; |# s4 i1 A+ S# }5 X* \

    5 C3 y* e; A0 _3 L4 {7( W9 \( h2 @1 @
    190569292
    5 U+ l. d) Q& u" s) Q[
    : I1 ]% L# r( V0 n4 W2 M    [ 10 ],
    ) B6 p- M7 s# j7 T$ l! f    [ 9, 1 ],( F& ?9 D: l3 }6 l, X0 M5 F# w& G
        [ 8, 2 ],
    1 F" a/ w  G( n0 o8 T% n9 N    [ 8, 1, 1 ],- w& J) K& m. X# \: z, ^+ ^) E# l
        [ 7, 3 ],5 b. x* |5 K! I* B
        [ 7, 2, 1 ],
    1 B0 v! T4 F% Q& m8 b; e; Y    [ 7, 1, 1, 1 ],
    9 t  e" x2 R$ |+ U1 T  C    [ 6, 4 ],
      ~' j$ t5 a) J* ]    [ 6, 3, 1 ],! }5 d; d9 @: `. _. M% T
        [ 6, 2, 2 ],
    + ?/ e& L; _1 Y, h3 l6 w    [ 6, 2, 1, 1 ],
    - E( h" |% f7 i: Y& ^    [ 6, 1, 1, 1, 1 ],1 D& h' x- P& |$ Q& e" w
        [ 5, 5 ],+ |* X4 v' F0 U' F  u# M6 v0 ^9 z
        [ 5, 4, 1 ],* j- I/ p+ |4 P1 Y0 m( e8 i4 @0 r
        [ 5, 3, 2 ],
    * W, M) P# G/ S3 \    [ 5, 3, 1, 1 ],
    - [. j' H/ t2 p, K% v" K+ A    [ 5, 2, 2, 1 ],) M6 _9 R* i: i# B9 z
        [ 5, 2, 1, 1, 1 ],
    3 V: A) c/ R8 c2 |0 s9 k    [ 5, 1, 1, 1, 1, 1 ],
    ) V6 z+ g& s& e' A1 Q! `& y    [ 4, 4, 2 ],) g& U0 k% X5 M. D, I: j# u# K- `
        [ 4, 4, 1, 1 ],/ V3 F( @5 C; A: j$ @
        [ 4, 3, 3 ],2 M# t- q& Z% |" U; I6 m
        [ 4, 3, 2, 1 ],
    ! J$ f$ P  O% G- u+ I4 Q% E1 j    [ 4, 3, 1, 1, 1 ],2 _& f- m3 U" {/ Y% P: q
        [ 4, 2, 2, 2 ],
    / v6 j% P2 ?' X0 E0 H    [ 4, 2, 2, 1, 1 ],
    2 g4 Z( W# i0 n    [ 4, 2, 1, 1, 1, 1 ],( Z, u2 h( T/ m! c; z  l. j* f) `
        [ 4, 1, 1, 1, 1, 1, 1 ],
    + R# Y% n% ]2 A7 G    [ 3, 3, 3, 1 ],( X3 A! f4 `6 o3 a
        [ 3, 3, 2, 2 ],
    ' t: C* C7 E, c( z; I" H& e    [ 3, 3, 2, 1, 1 ],
    ) K2 R2 `$ D/ d    [ 3, 3, 1, 1, 1, 1 ],
    + y; q2 b, f  G1 T2 U7 x    [ 3, 2, 2, 2, 1 ],
    : b0 O+ Z/ g7 u. D6 w. y    [ 3, 2, 2, 1, 1, 1 ],* K% D6 H- V  T
        [ 3, 2, 1, 1, 1, 1, 1 ],
    4 e6 K3 {4 s# R) d    [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    . [& t7 B+ [* V0 v% d$ B* V    [ 2, 2, 2, 2, 2 ],
    ! X+ Q; V0 C' y3 o7 j; [2 h5 }    [ 2, 2, 2, 2, 1, 1 ],
    ! S) ^& I! k0 O3 S  I- l& y) {    [ 2, 2, 2, 1, 1, 1, 1 ],
    8 A# m& T+ m8 `' e6 m0 I$ M% i8 }    [ 2, 2, 1, 1, 1, 1, 1, 1 ],
    + F# A9 b1 V+ \! g+ h    [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],
    # |( u) Y" E1 Y  J" G5 L9 L6 [    [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    # w# {% S9 E' c0 ]8 Z]
    回复

    使用道具 举报

    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-13 05:11 , Processed in 0.445771 second(s), 91 queries .

    回顶部