QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5201|回复: 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 编辑 ) \1 L. h$ v% Y* D: u
    / Y# N2 Y( d! l0 k: Y
    n:=12;n;6 e7 {3 L8 Z3 q
    Factorial(n);求阶乘. S# g1 K; _3 a2 m. k/ M
    Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));% }+ f8 W; ?% U, T) j
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    8 h6 X' V) M' a& s4 {& aNumberOfPermutations(n, 4);
    - E% h2 X! g8 lNumberOfPermutations(n, 11);- i3 e1 f' A0 h- r& r% ?3 x2 \
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;+ T! l& K/ |% N' Q
    Binomial(n, 3) ;/ d/ S0 I7 Z" S" u+ M* B
    Binomial(n, 9) ;
    5 {, x  c9 V/ D( {+ p7 EBinomial(n, 10) ;
    % G- _+ S- [- a& a) _5 [Binomial(n, 11) ;+ `0 a8 F. x  P* Y5 x4 M7 l
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    1 q' ~& E% e0 |6 X% ]0 {: q0 s5 i
    Fibonacci(n);斐波数Fibonacci(n-1);
    # s! i1 c- @9 H, k$ G" NFibonacci(n+1);
    ( }' S+ u. }# [+ lGeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;6 O2 P' X- e" t) x8 r0 S
    GeneralizedFibonacciNumber(0, 1, n) ;9 B# ~/ M7 m: E+ @- ~; e! L1 [9 r7 t
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))
      }5 G, B( u# W6 A" K2 {k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    0 T( J/ e# O0 }; r0 jCatalan(1);; L7 J' d6 G$ X7 D  K
    Catalan(2);
    $ f( \1 C5 s9 oCatalan(3);Catalan(4);+ m  ]# T+ r9 C+ P$ e5 {8 ?6 `
    Catalan(12);/ o: m) M% m0 T( O3 z1 t8 g" F

    # l) T( F  i0 V: gLucas(n);卢卡斯数
    8 g* N# M% [+ N; {3 S9 |12% u2 s. L3 |; b. C! h
    479001600
    % u! U- X, N5 _3 w" p8 ]831600: N0 o0 A: Y0 m) a7 C
    12
    " o0 l0 u( ]( d! Z4 E132
    - ]" d$ K( u' k. }11880: w( h, c* H& s; J" q
    4790016007 U# E3 O" O! t7 p- b6 t. i7 }7 u
    12
    . L, H# a: [7 ]) M  i! y66) |3 J; A+ C% Q) k5 J5 S% k
    220
    - X5 J' e8 I  S  |, b' W2204 D2 z2 A# m. _/ y
    66
    5 K2 H& R9 |* E1 r12
    & X8 S- y+ j+ C2 H831600
    ( j$ g. n8 I# w% P, X! l144( l/ Y; J. p: m$ t
    89
    1 h- i8 n" W" C! f: H( Z3 u233( Y: p* @+ d: d" G% S6 M( r% G
    233
    # [# C* J$ `- N/ [9 k; H610
    1 d. Y0 _, z# i. a144% _$ {0 u% G+ @+ P1 Q5 T9 ~4 `$ u4 p
    2080125 y4 ]" a% O+ j0 x8 v7 y# N
    2080120 ]+ D  Y# o' |2 d
    1
    . s8 R( @) }/ z+ l; O% d25 p, v6 P( D) \1 ?, |* F1 a
    5
    & m  S" z0 C: T% s* V( k14
    9 |7 Y4 G- t& @, q4 Z208012, T5 |8 Q+ R1 u* J" z5 t
    322
    ) ]3 d7 K9 x, `# O$ c# G" _( _- V4 j! [

    - \* Y$ T+ q  w7 v卡特兰数=(2n)!/((n)!*(n+1)!))
    ' U  \8 E1 }+ dCn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: 8 }: k5 ^. ^; m( G8 f- l# P( w
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY# Q& w9 J% n/ [8 b; `
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数:
    8 j0 E+ O! m1 ^/ ^- X# I0 i* g((())) ()(()) ()()() (())() (()())
    ' r2 H( g' H# Q9 P* j. R/ KCn表示有n+1个叶子的二叉树的个数。 : ^2 W1 G: W, Q  K& N& B6 M
    8 s1 Y& s/ u2 n! E! ?* |0 Q
    Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) * C( |) |7 J' @; J" @( J7 s
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    * j& I  w9 K* x: X- A0 n, `( v' S1 c/ E" h* b/ I/ i
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    # H4 u, G1 _  l7 u" {) O3 o" _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. ; |3 n5 Z- ]1 _2 B7 }& C4 Q
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:
    ; w+ p$ l* H* {
    * j  m5 d- K* F+ X0 C9 Z

    4 V+ p3 Y, {* D8 p0 O& f' P' f: I8 W- E1 a/ d
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。$ z* u8 G0 t2 e- Z$ |

    ' h+ q7 @, `1 e; ^但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    ' X, V/ P' {$ ?' a6 l* x# i3 m+ R$ n$ a4 Z$ S

    4 Y$ q/ N) y5 o# ]3 C5 _1 f; _1 N$ i1 O
    n:=100;n;0 d. k5 J" ^8 f" r4 l/ a0 O5 Y9 n5 H
    a:=Lucas(n);a;! {$ O4 x8 B: Q& n: h3 s
    b:=Fibonacci(n+1)+Fibonacci(n-1);b;
    / K! l9 f3 ^: o4 J, cLucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    $ l, B& X  ]: m2 d  T: a* u& \
    ; t6 r3 t8 J- H3 P6 A, N: k1003 Z- ~6 B& O* X+ O2 J  s/ Y+ E
    792070839848372253127
    7 |0 M, z! k. F792070839848372253127
    " N1 V6 q$ l" G7 @1771124240896309575375
    2 q+ Z/ y9 L6 R& w; p$ S- n1771124240896309575375

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

    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 编辑
    . |4 T8 ^0 N0 }, \% N( O5 R* r: h) B6 u
    3 F9 e; ]: c. z! x  ^  E9 K反费波那西数列反费波那西数列的递归公式如下:
    , q( t7 a' R: i
    * b" k* \5 X8 M% z7 oGn + 2 = Gn − Gn + 1
    / A5 _# B* E; \' D- o如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...( G( E: k/ w* ]3 b( E5 W/ j: q

    6 \/ K+ X7 z# ?6 `即是F2n + 1 = G2n + 1,F2n = − G2n。
    ! U. j! ^; \8 i1 h  F" m6 l+ Y# M) l& h3 P* L& W9 A5 E3 j# C
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    ! t. W3 o" h6 x, ^StirlingFirst(4,3);
    ) `' S; m9 I6 yStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);5 t% ^* c" q' F' P
    StirlingSecond(4, 3);0 j# V( W5 V' @1 X& M
    23 B9 t6 j* U% v  z
    52
      H% s4 B  M* I( [- N7 |57 k" O8 [" l6 m7 ~
    15
    . m0 t$ H; M" Z% {, o( t3 h-61 X9 k; P  D7 f0 R+ }  g: H* t
    11
    " T( S5 f( O+ A2 ^5 O8 B' q-62 C' o( S6 d) }2 d
    1! a& V6 l7 G% Z# ~3 f
    7+ M6 h" ?  n  R/ n% G
    6
    * K9 h# I0 F1 A9 B( X9 ?6 i) c0 f1 Y0 C9 u- b) g
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:5 K6 `; e  [! F9 |6 q9 n
    ) l) j- f# ^7 j. q
    {{a}, {b}, {c}} ) |4 Z7 ?4 {7 c0 U8 D5 f- S% {
    {{a}, {b, c}}
    - G8 K4 q% k/ D% I0 d; E{{b}, {a, c}}
    ) @% w6 p4 i* e9 Z  d% l1 T{{c}, {a, b}} " R2 E" _" `7 Z" |" I4 E
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s9 u. b9 f2 a) W# m! P; ?% O0 B
    s(n,k)是递降阶乘多项式的系数
    " Y4 x  s  l4 W! I有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)5 g( W( Q: U- B) O

    , d4 r  R: w( R$ j; k9 P1 b% Q1 R, V- q换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):7 d5 Q7 _; A& |0 b, c
    2 I$ h# s3 P4 q( |  Q% ]
    {A,B},{C,D} 7 X0 M+ R5 N- Y1 Y2 T
    {A,C},{B,D}
    - D, X+ T+ s4 o  O4 e{A,D},{B,C} * c6 p2 @. B# j: ~$ a
    {A},{B,C,D}
    4 M0 C' Q2 a8 m8 B$ K{A},{B,D,C} 1 n7 ?/ K- {  E" o) U" m  N; V
    {B},{A,C,D}
    " X6 E# G  A( e" ~( \{B},{A,D,C}
    4 h/ k' C2 M3 e- J$ I) o- I0 k{C},{A,B,D}
    * a; C0 \, A$ i$ r7 B4 n- P( }# p{C},{A,D,B}
    ! s% s3 u: l) |1 v% I{D},{A,B,C}
    ! @1 e0 n: B" K, V2 ?9 y{D},{A,C,B}

    , L( A9 V5 H0 I  d( m6 f第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S& v1 m+ H+ K! e
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k)
    ( I4 Y0 K* A/ NS(n,n − 1) = C(n,2) = n(n − 1) / 2 1 K( }# m& F4 c! d. T
    S(n,2) = 2n − 1 − 1 : n: g; y7 R* ^+ s
    9 }6 y2 Z- L& a
    换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
    3 X2 K' Z. y% D
      ?, Z9 J8 o# n3 x$ H0 @5 {1 B{A,B},{C,D}
    % ]! I, Z  r$ c4 J! ?* r{A,C},{B,D}
      `& |; _  R; o+ \4 P" R% I{A,D},{B,C} 7 F* R0 q, |5 V) a7 [0 U7 r
    {A},{B,C,D} ! M, c+ t) W4 T5 M
    {B},{A,C,D}
    ! |; ^' v- Y6 ~9 f; c8 U{C},{A,B,D}
    9 W* ~7 Y  x* R: v( H{D},{A,B,C} - z/ |+ E1 A4 a3 s- 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 编辑 - z1 Z8 V. _) H

    & m0 o, Q: c' `1 pn:=5;r:=3;
    & s  r) L8 d- b8 ]# NEulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;! G! @8 h! X1 ?! t
    BernoulliPolynomial(n) ;伯努利多项式
    . P, o7 P7 |+ l$ {# o$ Q* x
    ( C4 ?2 ]+ i. H. o) s8 \26
    5 T$ R9 s* e* C1 d+ g3 S$ n137/603 U5 C7 P8 S* F3 m
    0
    0 J! H4 Y" A( [5 k( S3 C, x5 V% c0.000000000000000000000000000000. I$ s- e( F9 F) N
    $.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑
      ?4 N" c/ Q$ \/ H% I* q" X  N
    9 i( p* g/ _0 s- }! f( J0 J6 r# ^5 E4 T. i

    . g- Q0 U4 ^$ [# k+ \7 f: \4 g. ?2 Y伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    . a# U9 `1 U8 H3 y$ u; r3 x: K! w$ ~5 K# ~; C
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
    ' J- D" l* O: W, G+ z" h
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑
    . }0 e3 _& k# j% p! Y0 r0 \& ^; m% g  H
    拆分 。。。。强!1 t) b) e3 R! S) P4 Z, {- Z; S$ E
    5 _" r* x* [# u, m) e: h
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;8 }. G. N2 I+ F8 }' x

    2 O6 p5 v$ D) ]4 Q+ E+ C76 u; J2 I/ B4 N4 r- W9 Z: q2 e
    190569292' D" o& T0 Y: r1 t* }
    [
    9 z: E, o5 D( c# k    [ 10 ],
    ' d( A! Y$ x" [. Q3 Y; J# }. m    [ 9, 1 ],
    3 [# o! h4 L4 N    [ 8, 2 ],
    9 F" y8 }: {* @    [ 8, 1, 1 ],
    + T7 o3 y5 R" Y: I' `( [1 `1 q    [ 7, 3 ],/ F# T1 X* E5 Z" ~
        [ 7, 2, 1 ],+ x3 I$ D/ `3 n' @% V7 @) ~. |
        [ 7, 1, 1, 1 ],9 k1 M- g9 t8 R# m1 x- O7 O
        [ 6, 4 ],9 I- Z* e6 z4 o) q% u* _! N% f
        [ 6, 3, 1 ],$ a  d( @* C8 x5 A1 D3 G: E) V
        [ 6, 2, 2 ],  c/ M: S2 t* p5 t
        [ 6, 2, 1, 1 ],
    . {3 `; m: d# }; r* A6 f    [ 6, 1, 1, 1, 1 ],
    # ^0 i. i8 |: C& X) y+ a: \    [ 5, 5 ],
    " w9 c% y9 @0 s5 H" u    [ 5, 4, 1 ],
    / b0 Z* \) M8 U4 ?" k$ @; [; |/ {7 o    [ 5, 3, 2 ],
    / V/ G9 ^6 a! M9 x    [ 5, 3, 1, 1 ],
    6 Z/ U; N$ H1 L3 h0 `    [ 5, 2, 2, 1 ],
    , S( s8 v6 I4 Z7 d    [ 5, 2, 1, 1, 1 ],
    1 \0 A$ r6 V. I: k' K1 l    [ 5, 1, 1, 1, 1, 1 ],: Z  r* \. B7 O; `: W/ O+ f
        [ 4, 4, 2 ],
    4 F% |% r, k4 |# ~0 V2 ]# a# x9 o    [ 4, 4, 1, 1 ],
    8 P: V" ?; l- Z5 h  V2 p$ W( f4 n8 B    [ 4, 3, 3 ],2 i" c2 A* M) b+ P( X% W
        [ 4, 3, 2, 1 ],
    2 e0 Z( z% o& u4 C6 |& D- R+ H0 o    [ 4, 3, 1, 1, 1 ],! H. a, z0 [, D! A2 p$ G  O8 n  u
        [ 4, 2, 2, 2 ],& F  ~* h' X6 d% \
        [ 4, 2, 2, 1, 1 ],
    ( [( P& U% o( }; ^    [ 4, 2, 1, 1, 1, 1 ],
    8 ^" `$ t( W% p; u. a9 ]- Y/ L# J    [ 4, 1, 1, 1, 1, 1, 1 ],
    * S7 \2 x5 p( s  P* S2 w1 t9 ?2 k8 L    [ 3, 3, 3, 1 ],1 U/ X2 C  P9 b7 f0 `# }
        [ 3, 3, 2, 2 ],4 ^( g$ d% M; U: e1 |: H- @2 _" L
        [ 3, 3, 2, 1, 1 ],( y% ~5 s! W; N
        [ 3, 3, 1, 1, 1, 1 ],( @7 @  G% P+ p6 h3 k- p1 X$ D
        [ 3, 2, 2, 2, 1 ],3 C- c8 q" @0 S# u% o5 ]$ U
        [ 3, 2, 2, 1, 1, 1 ],
    . B5 q% h8 W: u/ r* M' I    [ 3, 2, 1, 1, 1, 1, 1 ],1 p" Z9 n3 c6 w
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    7 x' q2 }2 }* ~3 f$ Y& x4 V& g( f    [ 2, 2, 2, 2, 2 ],  R1 \5 T: f  Z+ z. i+ b. h, q8 m
        [ 2, 2, 2, 2, 1, 1 ],& u+ _% E! `2 x# _5 v
        [ 2, 2, 2, 1, 1, 1, 1 ],
    2 `7 @. H4 X6 e, n2 U    [ 2, 2, 1, 1, 1, 1, 1, 1 ],
    1 ?$ p9 \6 O  M) @( k) C6 W' v    [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],: e3 G  {- W+ k+ o% M/ |
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]- P& q$ I/ m/ {1 {4 C: t
    ]
    回复

    使用道具 举报

    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-20 15:07 , Processed in 0.493935 second(s), 92 queries .

    回顶部