QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5202|回复: 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 编辑
    , X& P7 u( ^5 r& u' b; F  W* i
    7 c" X5 x, K3 y& \) L) on:=12;n;2 m# b: f2 P# B; |1 u
    Factorial(n);求阶乘
    2 t, P8 p, T6 a9 r& cFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));, K% g/ w/ h( ?- K9 l! d2 t
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    " p. e0 c( @" L% a1 aNumberOfPermutations(n, 4);1 z; T- I1 E/ m: p: K% J6 b1 T8 }
    NumberOfPermutations(n, 11);
    ' i* W; j9 J) L' d  m: {' S: z( PBinomial(n, 1) ;二项式系数Binomial(n, 2) ;
    3 ~; g* z4 a- x4 Z1 V$ c( TBinomial(n, 3) ;
    + z3 r* \' {4 p$ }) p! |4 N. t( h0 VBinomial(n, 9) ;) H7 e' H5 G6 S  c7 W7 P
    Binomial(n, 10) ;' Z2 E* Q0 h9 h7 s, S0 y6 y. u* ~
    Binomial(n, 11) ;6 t$ T  r- G3 @
    Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    5 O# ]) f8 D% V& L4 J% B
    + e( d! M$ N7 b. I+ LFibonacci(n);斐波数Fibonacci(n-1);3 c/ A1 j9 _" ]% S. T! r; W) G
    Fibonacci(n+1);$ N+ _6 \% w; x8 r- L  R. v
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;# y4 P* R+ G$ M
    GeneralizedFibonacciNumber(0, 1, n) ;
    ! W% Q& p( O3 t, H; K! U4 V+ z- i* eCatalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))- @$ r- Z' n" n
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    ) W7 J( O3 o6 u3 J. n7 s/ Z% _4 eCatalan(1);! f" r( {4 I1 M( F
    Catalan(2);: b! D' r& [; X# v/ ~! k
    Catalan(3);Catalan(4);/ K6 n3 F, s6 n( e" `7 X. e
    Catalan(12);' P1 t# q7 \6 L0 Y- a7 U7 g
    # j7 r: w, H5 x# c
    Lucas(n);卢卡斯数
    ( b& c$ M  a; U# n123 J# r) @: ?6 q
    479001600( v2 l1 c- s7 w4 H  V' b8 h
    831600
    : }9 H# i: Y/ J# {& c12
    * w; a5 G$ b1 [, v$ M3 |! G132
    9 s9 M3 a, k7 s" X0 X8 y11880
    % M, M0 |8 S. E2 u479001600
    6 F+ \0 K$ i% X12
    % O: [7 O4 e/ w6 j- Z& y: O66
    " b! Y3 X. L/ |+ j& Z5 Y7 s# A7 u220& k3 q8 \! y* v' e! C7 ^
    220
    2 D# X+ z/ H! ~8 }- E66
    1 x4 l* M& @9 J% N9 u7 F/ D: y) n12: t* b. [+ F% I* b
    831600* Z$ N3 e" k7 e; L$ U5 e: }
    1447 [; F+ `2 F9 L6 y: M
    89
    # e' F6 U/ _( k! l3 _+ J1 M& i$ z7 B% C2331 V; G6 @* ]: u3 C+ z
    233
    . Y# z% ?- o7 ^- |- Z610
    ( ~0 J  A! v- w9 w6 D# e$ c/ d, g144
    + v/ p1 M( f$ U  |3 y208012
    ; K, P- }% z' p2 u. R7 G$ P# C$ M2080128 F0 q: S, N8 F& q& O2 p
    1
    4 b% Y4 s  i) X& V2
    % _1 f& N2 r9 S9 K9 Z5 O  A5, Y# O+ A9 Z: X% L1 O# V0 g$ E
    146 p: }' d/ T0 e, p
    208012
    * a9 M1 F$ v& b( ]. z$ O322
    + l3 E! `! t: s( x1 X
    * H) ^% e/ f$ q) z9 p2 h  S7 o/ |$ Q0 `- \6 |3 k+ j2 U
    卡特兰数=(2n)!/((n)!*(n+1)!))
    0 }* ^" ?5 P* W. P' _* @- a' BCn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: & K5 e: f0 g# i6 \# u5 i
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY
    ( {; O% Q5 m! n2 |将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: : ?% c# x( s% I: h- F2 p
    ((())) ()(()) ()()() (())() (()())
    ( ^; H: v0 @$ X3 @& r* eCn表示有n+1个叶子的二叉树的个数。 / `2 r0 L7 C( _( h

    8 K5 K. F% }* u7 u4 I' QCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) : L4 e9 R8 l4 p  s0 z0 a* C
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况: / u4 r% E7 d" E4 \, T  i( x7 k
    7 O$ a% m4 p& w7 Z0 ]  |
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。   h. j4 q4 ?2 `: h
    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. ! U. }* `: x$ Z5 W4 `% i; Y
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    # w9 Y8 F7 Z- t/ @+ b" t. T# m5 g0 C' s9 V% O
    + _) k# z# Z3 k, C/ T

    1 B3 w: N2 l$ ?4 [0 }2 Z卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。3 }5 c4 J; G9 P5 K

      P! |8 T+ B; c) ?但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    & J, V; n. j. p& G) P6 n; g# T
    ! n, u2 a. b( v1 U+ |
    " v1 S/ q8 w8 q5 f) @8 p" Q2 I, S- E% h% P
    n:=100;n;
    ; |; _0 }- b# x1 r- g( D$ @9 ^9 Wa:=Lucas(n);a;) n7 }7 E6 v- W; w
    b:=Fibonacci(n+1)+Fibonacci(n-1);b;
    9 ]0 n, L4 X  i, e9 S5 H0 YLucas(n+1)+Lucas(n-1);5*Fibonacci(n);& ^8 N. m# a* }7 I
    / g4 E$ g2 a# G! e" p# ~9 ~6 a3 S
    100
    3 e+ {/ z; H: K/ t4 w: H) z3 K4 S792070839848372253127
    * E% z5 l" p; s5 O792070839848372253127! L! U+ ]* W9 g
    1771124240896309575375( q" r# E% ?# N  N8 G+ k
    1771124240896309575375

    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 编辑 # \, O6 d2 D1 H9 B- H, q

    , B; ]$ x1 E1 j+ }2 V) w1 c1 j反费波那西数列反费波那西数列的递归公式如下:6 j  X) G( |" i7 [6 \2 T2 G
    ' K) }, l; L4 v, G5 B. J4 E
    Gn + 2 = Gn − Gn + 1   G4 ^9 C) A! O, G) R+ P( i8 V& v
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...7 _5 i) P1 @1 A

    ! w; j9 S" r0 O0 n+ O' z即是F2n + 1 = G2n + 1,F2n = − G2n。
    , R0 J) R; g3 A. u
    3 Y- B& b+ q- Q) hBell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    ! {# o3 |8 S2 DStirlingFirst(4,3);
    : P" D! ?2 h. i( ~StirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
    * X7 C0 j  q& T: Q. K4 HStirlingSecond(4, 3);( a5 a5 M6 {0 W* b2 W( w
    2
    # F6 K6 y1 q) ^. c! ]52
    1 @, b5 Q- A( Y( O8 X5
    3 J$ W$ F$ {. O; X/ }, B3 c& S15
    + H; o4 f3 W$ k9 v+ W-6
    * {, H& M6 b9 j) \7 _- u: d4 j11$ b( R) x8 r2 q7 p7 P% G; U  B1 P
    -60 \. w4 [: B" f: A  W' U
    1+ _/ o$ q( q# k
    7
    $ k/ r4 T! Q1 i6+ z) e/ X( r6 p& B& f, v

    3 T2 ~+ E) e" o8 mBn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    & E. P' [3 ?$ \0 l& L& i. f  p* n1 A- ^1 b" |) {% Z
    {{a}, {b}, {c}}
    8 C0 W7 \$ |; |1 J; G  ?{{a}, {b, c}}
    ) g6 J  |1 x6 U# o" p{{b}, {a, c}}
    ! g  I- f: H' [1 q) p{{c}, {a, b}} % T: g% x; ^& u) a
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
      g$ g; d# C6 p5 Cs(n,k)是递降阶乘多项式的系数
    1 s. }3 @! r+ h1 h( {& K, n& U有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)
    / i% `. ]8 f0 _6 A7 n" K+ j4 ~1 h* G! `% X& p+ b$ u: X
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):
    % g: p& r4 w  D5 H1 M: {1 J( k. ?3 C% K  Z
    {A,B},{C,D}
    ) |0 R& p- n/ n; b7 K{A,C},{B,D}
    1 n& g" [. q/ r% `3 O! X{A,D},{B,C} 7 C% h/ e  m9 _5 ^* o' }" a9 ~4 ?
    {A},{B,C,D} " G9 W# h5 ~0 u$ J
    {A},{B,D,C}
    , [+ E0 n: y& F, J+ T1 S0 c6 K- X9 V{B},{A,C,D} 7 u' y. z# ]% r
    {B},{A,D,C} 7 H3 h: g/ h" y- X3 f, G- Z( d
    {C},{A,B,D} : \3 \3 O1 N3 k) O9 O2 m  r0 V+ [( a0 U
    {C},{A,D,B}
    3 F. o8 P/ B$ }' N, S& u{D},{A,B,C}
    , t4 N* g( E5 i9 `% w{D},{A,C,B}

    6 [5 x9 Z7 ?2 d7 f/ q第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S/ N8 l! Y: u1 E8 ?  U3 T
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) 4 f% b8 c* |( O2 q1 ^5 h
    S(n,n − 1) = C(n,2) = n(n − 1) / 2
    " n$ r/ g7 o. r# ]! Z/ L# uS(n,2) = 2n − 1 − 1 0 ?! l8 [  y3 j) a0 H4 f

    & r( u4 c8 D- ?: w4 i换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
    ) l5 ~, Z4 r, F/ R8 p* m6 E3 R$ W9 P; @" D' i1 H0 N( }7 t" ^
    {A,B},{C,D} ) f+ a% i2 l3 c
    {A,C},{B,D} * G; q7 r  e. }0 Y2 m3 i, Z
    {A,D},{B,C} ( m3 t0 D' F! k
    {A},{B,C,D} 4 O4 P9 g! m1 D- j$ l
    {B},{A,C,D} - Y" J, i" n- E$ p: |* l
    {C},{A,B,D}
    ; w: N1 U( q( r3 f7 g. T' L, T{D},{A,B,C}
    # c8 m/ h3 X2 Q5 w因此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 Z# E8 B" G" S5 p/ m: i# g& B) p* z, X" @
    n:=5;r:=3;
    2 i1 ?1 {" e1 Y0 s# c" Z& `EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;( o. M1 _7 N  Q* R
    BernoulliPolynomial(n) ;伯努利多项式
    : a# j4 E% @+ W* p
    $ \) W0 C. B- ?$ f( [26* j( }; s) X" x4 \
    137/60
    ! B6 l* P6 w" \6 a2 Z, F05 L3 p" b* a# ?  ?! C' }+ v& ~
    0.000000000000000000000000000000
    , G2 K- k8 h: X" q1 \: s$.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 编辑
    . r4 g$ a  ~$ n7 G; D" T0 P5 g8 G/ d6 S) P) R
    $ o  O: W+ ^/ o& X9 S" Z6 ?

      m  c! ]; U# _  A/ t1 D$ T0 M/ r) v- O伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    1 d8 R9 z# `/ n! l1 e
    7 o$ a9 |- v6 X* l伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
    6 U# K6 b+ w# J0 x- e
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑
    " x4 g" p3 V8 D7 ~# ^# n  i* f  B  G* R, R5 P+ H
    拆分 。。。。强!, n# D+ J& f, a$ y
    : T0 z9 e# U9 P" z3 g
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;
    ' F9 C' m: }' A' o7 Y# U3 Z: @  _# e) q. H* @
    7
    . G9 [! [% g% m1 J& N1905692928 W1 q, W4 @' e& A- `* |' _
    [- F5 a1 R% J6 {
        [ 10 ],
    " S) w$ f, w+ j    [ 9, 1 ],6 ?$ h  j; c2 h
        [ 8, 2 ]," j( j" m7 e5 y) k1 k/ u
        [ 8, 1, 1 ]," Q% ?! w4 F8 _; d! x+ E% w( h
        [ 7, 3 ],
    4 ?# t! t/ g+ T5 K    [ 7, 2, 1 ],
    * Y/ a! j% h7 s4 Y4 J1 I    [ 7, 1, 1, 1 ],
    3 M1 k, Q4 o+ U- E  D- |3 b    [ 6, 4 ],+ F% o) g6 M4 R* K6 s& s
        [ 6, 3, 1 ],1 _' Y: B7 m- D7 ~* m: I
        [ 6, 2, 2 ],) Z& K0 }; Y! a! d' f0 c/ z
        [ 6, 2, 1, 1 ]," l5 U# U* g( {' J9 w4 y
        [ 6, 1, 1, 1, 1 ],- i( {" t* L; l( s3 F5 x- V% m
        [ 5, 5 ],
    . p. V/ E# z7 t' |    [ 5, 4, 1 ],
    3 U/ F: _- t; _9 M  U    [ 5, 3, 2 ],( k% h, Z. B, `; i& Q, i/ o: `
        [ 5, 3, 1, 1 ],/ ?- Y4 b! k8 j! e1 Y' u6 K
        [ 5, 2, 2, 1 ],
    / n4 Q7 `0 n* q/ ^5 d" t    [ 5, 2, 1, 1, 1 ],
    0 q, L5 {- b" j& V1 H% k# n    [ 5, 1, 1, 1, 1, 1 ],
    $ @7 x$ [) Q$ V8 _4 K  S    [ 4, 4, 2 ],) i! c: \, J1 m% e
        [ 4, 4, 1, 1 ],4 ~. h& E7 a! z5 \1 ]
        [ 4, 3, 3 ],/ `: Q3 u3 y6 R$ [" g9 r
        [ 4, 3, 2, 1 ],
    7 ^0 c9 \' `, i9 D" a( m    [ 4, 3, 1, 1, 1 ],
    1 K/ _- ~2 C1 ]) I4 J# g    [ 4, 2, 2, 2 ],
    ( T$ A+ W+ X6 N    [ 4, 2, 2, 1, 1 ],
    5 C. Y3 e: o) R( L' v7 V    [ 4, 2, 1, 1, 1, 1 ],
    9 S, L1 T# g# W7 D9 b    [ 4, 1, 1, 1, 1, 1, 1 ],
    ' l, O" Y: K% B! z9 @    [ 3, 3, 3, 1 ],$ f, f, v7 @' |* v: I: ^! i
        [ 3, 3, 2, 2 ],: t% d5 a, s0 w
        [ 3, 3, 2, 1, 1 ],
    * Y* q( P& d0 Z/ D6 R    [ 3, 3, 1, 1, 1, 1 ],
    3 b6 V* W# H7 Y0 N+ P' y    [ 3, 2, 2, 2, 1 ],5 {% G, t& ?! x" L
        [ 3, 2, 2, 1, 1, 1 ],) }- r; K2 K2 T" R3 Y' P; U
        [ 3, 2, 1, 1, 1, 1, 1 ],0 H/ C7 u- k1 E* }5 c6 {" ], @/ {
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    + ~( Z: Y/ r  O" d+ L" P" ?# w' V    [ 2, 2, 2, 2, 2 ],
    : R' b$ Z: G, X( J& ]    [ 2, 2, 2, 2, 1, 1 ],
    $ _4 ?4 g7 E' X    [ 2, 2, 2, 1, 1, 1, 1 ],
    1 b' V( m0 i: v$ D/ t; `; F    [ 2, 2, 1, 1, 1, 1, 1, 1 ],; l5 ~+ y* l; {' }; n6 D; F
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],( r* M) F0 q5 E6 E" A1 w: o, _
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]; ~' N: a6 S7 f
    ]
    回复

    使用道具 举报

    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:23 , Processed in 0.515829 second(s), 90 queries .

    回顶部