QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5134|回复: 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 编辑
    ) S* v( M! f' Q% v  |. u4 ^0 t7 M7 [9 a, C6 T2 I; T
    n:=12;n;
    $ ^- b/ Z$ n, A: _Factorial(n);求阶乘
    0 e4 O$ U; |" v6 kFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));- @* `1 f  |5 N9 E" M6 J' s) F. z) L
    NumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    $ I/ v( E( B2 p" T* p( N  aNumberOfPermutations(n, 4);, U( o8 A. [( A$ x  u/ `2 o
    NumberOfPermutations(n, 11);
    / c& H3 c; c8 p& u" X+ OBinomial(n, 1) ;二项式系数Binomial(n, 2) ;
    9 \+ q: b' ^6 K. XBinomial(n, 3) ;4 ?, D& H+ R7 K# |2 K& Q6 j
    Binomial(n, 9) ;0 `# a. p: G& n$ r
    Binomial(n, 10) ;
    5 d: A6 Z' b  [5 F! Y" ~- g$ GBinomial(n, 11) ;
    1 t5 k4 c" s7 m, K- |Multinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=8316009 o1 g  [+ Q: O" o/ m

    3 U) W% F% }4 G' `- O( }Fibonacci(n);斐波数Fibonacci(n-1);; i8 H" d' {, v( b' e  _# R
    Fibonacci(n+1);" X6 Z( `+ q3 ?# n# s5 z$ L
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;
    0 _# p& r" M* b8 y8 N( D, |GeneralizedFibonacciNumber(0, 1, n) ;
    " {& D7 z' S7 T: SCatalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))
    9 N6 o3 J% X1 O% Y5 Y/ Q& f  L+ P; Ok:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;9 y& e. i1 X+ o6 x
    Catalan(1);0 c1 M" V. \- Q" g
    Catalan(2);' W% P( b" r4 y5 j5 `" i/ _
    Catalan(3);Catalan(4);
    5 m7 m6 o- v& {( t, L0 Y+ ]Catalan(12);3 ^( @* U) o9 S* @: a; @

    7 w6 x0 O. J1 o+ ?; VLucas(n);卢卡斯数3 r* V" m8 [  F- T( T0 S, U8 o
    12
    " n1 A1 T, v2 z0 L5 ~  B( X479001600- h+ `+ f" W" N9 ?2 G2 H
    831600  W6 U: G1 a* z8 H& i% V
    12& X4 {2 D0 h8 X5 i% |3 n
    132
      {4 E/ y! L9 p3 j11880
    , O& X# f3 x, R1 b479001600
    0 ?+ z: F' P, [7 n! h12
    ) U5 _" x/ E. r7 z4 g! x$ j661 B* g) C+ X6 N  ]% W1 A9 E
    2208 }3 C/ ^  l4 s; G/ k: g
    220" Z& M0 [; o/ @0 }  L
    661 w5 f7 z6 Z' `: o6 k" Q
    12
    5 x" {+ a+ t6 [0 c( d8 U831600/ ]0 d. c( X7 Z$ ~2 X; V% F
    144
    , b- a% P5 J; Y: W, H1 |$ ?' E89. Q1 o8 b% o  T6 ]7 L8 Z  w
    2333 q" e5 H# p2 q, f1 d7 |
    233  v$ }4 W8 B: }
    610: H) x1 O) {; U+ ]# u
    144
    . u3 Y2 J8 V: E7 @/ b8 a; t/ x+ v208012
    4 b$ _2 h. Q( L208012
    # v6 P" v5 L+ b9 B8 ]8 q1
    ! Z0 C# y6 N2 j% C4 n# [2
    % S6 e* o& T$ d5" d% B2 t  k! l' t
    14
    , _9 }6 C* `# |6 o+ W% d208012- E5 {5 Q" Q+ a" {) H
    322
    8 L; \' N  d8 q# p& D
    $ |+ v/ x$ c* _5 I$ L5 h+ C
      ~5 j  d' `( `: F# w% }5 V/ Y卡特兰数=(2n)!/((n)!*(n+1)!))
    : T6 y. V  J. rCn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    ! h6 z/ q% M0 Q& Q0 W9 ]+ f: F+ o( `**YYY XYXXYY XYXYXY XXYYXY XXYXYY; D0 X) O# B, Q) g
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: - i8 E) d7 W2 f
    ((())) ()(()) ()()() (())() (()())' r# f9 t- [$ i* T* ~
    Cn表示有n+1个叶子的二叉树的个数。
    1 i4 V  i0 C5 `1 @4 Y/ G# I
    8 [4 _! ~5 z: D3 uCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    ; N4 V- P5 o. ?& }8 m7 x0 `! |+ `Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    1 s5 l! p- b9 |3 k8 p( e% D
    ' p' I% [3 f( g( E! U0 [# HCn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    / o$ b4 a3 K( T; ]  qCn表示集合{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 d6 [" D- @5 V0 p1 p+ A; a
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    4 e2 j7 @, o, r* j
    : t8 ~2 ^: D8 I8 ]
    0 k$ X2 D' T; S' C. c$ F: ?& g- z) w  S( o4 H
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。
    ) `' H% p6 P# e& o' e+ r. e  c6 G8 E! n
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    * Q! v* A6 c7 V4 P/ x3 U" {2 C" s) l2 Q) z/ f  }
    , ~( ?$ ?7 h% |) b& I

    " ]1 W+ l+ v# j+ D) zn:=100;n;
    6 e+ W4 I0 Y- E6 m4 Ba:=Lucas(n);a;
    7 U8 c4 \0 w7 x& s4 c6 K/ Ub:=Fibonacci(n+1)+Fibonacci(n-1);b;. Y- Y' C  d& u2 Z4 X, h
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    & |" u9 M# j: f$ T+ a. _9 j1 X1 F; t1 u+ [2 N# G
    100+ ^4 b: f' u, m. |2 Y
    792070839848372253127/ r2 p2 j5 Z' Q" Y- g1 W8 f
    792070839848372253127! M2 N# x9 P5 y( g0 H; Z# C
    1771124240896309575375; i( h% ]/ X! u, ~" o; }) l" |
    1771124240896309575375

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

    12.JPG

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

    74

    主题

    6

    听众

    3299

    积分

    升级  43.3%

  • 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 编辑
    * x) O: X$ f: T% T5 K/ {3 r8 K1 V; R& T; J- t6 Q( h
    反费波那西数列反费波那西数列的递归公式如下:
    0 g( w9 l. j, R' {7 I! w8 d! M$ a  R% V% C
    Gn + 2 = Gn − Gn + 1 - v4 \1 n* n% G; H3 |, K
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...
    8 h' z1 J, d. d" j2 Z
    0 {' s+ ?$ c# R: [即是F2n + 1 = G2n + 1,F2n = − G2n。6 K7 `& K7 c* W8 l) ~# }# ]
    : g9 X9 f% ~) y' z; @+ g
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    % p8 ~% ~" }, F6 b8 h5 A5 `StirlingFirst(4,3);
    % }3 @0 z1 v3 v3 A5 ~2 `3 WStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
      F$ u* n7 w0 q) W' n8 bStirlingSecond(4, 3);
    " d( V4 m- R& N2
    3 d( V* T# U. d# E! J$ c( Y52* C/ d1 p7 x$ T; u! v, y3 n/ I* z
    5
    ! V+ j  j" t8 O15
    1 |% k4 x  T" L* y4 j# O-6# e8 F' `7 V) M
    11% Q$ S3 g0 c. A' r0 t
    -6# ?' `2 S1 P$ c9 o
    1
    1 u. o/ S, [' R! i* n- s- x7
    ) u* i3 j& u2 K$ r# a66 C/ N  i4 x& V/ q

    * q+ l5 \+ Q( ^6 IBn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    ) f$ n: S5 a1 Q' T# m/ _5 A/ C& H3 y
    {{a}, {b}, {c}}
    9 x/ e: @! l$ j: p& s& ^+ o. T{{a}, {b, c}}
    : b. [1 B2 H0 J8 {/ F6 C{{b}, {a, c}}
    3 P! M$ A; J* |3 I' W) F5 [8 S+ m{{c}, {a, b}} ( _; X/ A) P( {
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
    4 Z8 b( ]& r9 L, c; {  A* {" ss(n,k)是递降阶乘多项式的系数
      ?1 ~, U- ^# d, j; j! k# `有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)
    2 E/ i% e% X9 m$ }* N
    # {, c5 Q1 J' ]+ d2 B9 u3 q换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):* Z! e& B" B. v! W& Z

    * K# v2 v8 u) {{A,B},{C,D} 8 i% }& c+ O+ D1 M
    {A,C},{B,D} 2 P2 r5 l2 C! \/ g0 ?7 r
    {A,D},{B,C} ) n. l, E- d3 l! i, U5 C! E9 C
    {A},{B,C,D} 5 u& _9 G4 k# q, x/ y/ N7 A) q
    {A},{B,D,C} ' B. c' m4 |5 i1 N) {& x' v
    {B},{A,C,D} ( Z0 S, R% F$ a  |9 W
    {B},{A,D,C}
    9 z. w- L8 K$ g' `3 S' J9 p1 ~3 a8 W{C},{A,B,D} # b! A  Q$ _6 a/ L2 e+ ^6 r
    {C},{A,D,B}
    4 Z, ?  G" M  K' U- l; m$ e{D},{A,B,C}
    5 d* d- R( M+ g! h6 v* ^5 x{D},{A,C,B}

    $ \9 P! n# z6 a& p" J' |* `1 i  _第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S5 X# Q5 G7 c( m. q5 U) A
    给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k)
    1 c( M/ n+ u/ U+ H; SS(n,n − 1) = C(n,2) = n(n − 1) / 2 6 p( f+ q, }1 Z
    S(n,2) = 2n − 1 − 1 3 y/ `% E+ E+ }" E
    ( |+ O1 z/ @/ E: V* o
    换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
    7 d8 M$ i6 z8 z1 ?0 Z$ v
    + n5 f( i: r2 r' B6 N- M- E6 P5 q( z{A,B},{C,D} 5 ^" x" B& A$ J; R: A" I3 ^
    {A,C},{B,D} 6 `; i+ M* z; \: ^) ]( |% [
    {A,D},{B,C} ' N# |, I+ g; {# s/ H
    {A},{B,C,D}
    ' J: B. H5 T+ `{B},{A,C,D}
    : Q1 V( `$ I5 B" j{C},{A,B,D}
    7 s% ~) C6 z4 S; _: _6 m/ f& J. `{D},{A,B,C}
    ( \$ F# j3 U7 J/ d因此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 编辑 & J' A. X+ E- f2 F8 v- F( B8 J

    # ~0 t& l4 X8 ~1 `4 Dn:=5;r:=3;
    . S2 V2 n! i2 {% S1 {EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    & Y! S& E2 L2 D5 Z6 ^$ [% g1 [BernoulliPolynomial(n) ;伯努利多项式' k3 i' L& f7 h+ {
    % q1 C/ t) I) a9 t
    26
    % @  k3 t: Q" E2 }+ `3 M: }. p137/60
    + c% o* \7 k. b5 I) t+ I0
    ! f0 d3 U, g" _& u6 c$ O9 {: u& y$ z1 D5 ~0.000000000000000000000000000000
    6 ?& f5 O1 B' o( D9 `4 i) i* w& b$.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑
    , _5 B6 q; X4 n. o
    5 e6 e: E6 {3 E
    " m! M& D8 H! }, q& S" v) p2 _
    4 D! H! j1 m/ `2 J$ W伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    8 e+ |, Z0 N+ P$ M6 z- C( \" A; h0 U+ p4 V) h5 y( R
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
    ) M$ y9 s0 I0 F0 B# u
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑
    / n; N6 ?2 P. O. t& ~' V
    % h% |; g; y' M3 Q* K! E/ n+ B3 ]拆分 。。。。强!9 f9 q% P& Q& ]4 \1 e; S
    : k( g; J0 l; s0 ^% `8 _
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;2 a0 E% S4 ~6 o& l  v; M
    ! y$ ~. e! I( V- a  K2 g
    7' h( H- |; x! F
    190569292
    - G4 e7 t8 H4 Y. u[- N' d& R4 e2 a% _( }
        [ 10 ],
    6 ~' n8 m8 G9 V: Q0 [0 f    [ 9, 1 ],
    5 ?9 X4 o% O/ g" A3 D# j    [ 8, 2 ],
    # h8 x. q5 H& [5 O) j, n$ j    [ 8, 1, 1 ],
    $ C) n2 S( W6 _" p8 Q    [ 7, 3 ],
    , M- r# d( ^# q3 K" A: l) C    [ 7, 2, 1 ],
    ) S8 T  ?- h* K    [ 7, 1, 1, 1 ],) J2 ^) E/ T' d" ~( j
        [ 6, 4 ],
    ! _" e! X* `! F4 F4 w& \2 y    [ 6, 3, 1 ],( X* o# I' u& m* B5 l7 ]( q: `3 w  i
        [ 6, 2, 2 ],
    + w0 \( n& d& Q8 C1 }0 B    [ 6, 2, 1, 1 ],& i# B- d# Y( x# Q7 e+ }
        [ 6, 1, 1, 1, 1 ],
    ' n2 w0 x/ c$ n    [ 5, 5 ],9 o$ |0 }- A) W6 a9 e# |
        [ 5, 4, 1 ],, b- ?5 n+ Q1 p1 j2 Q
        [ 5, 3, 2 ]," T; O) Q2 }6 Y! k6 D: u1 c7 s
        [ 5, 3, 1, 1 ],
    4 @8 `/ o- d, t% }    [ 5, 2, 2, 1 ],/ F6 w! G+ x8 C4 S2 I. f3 l  [% [
        [ 5, 2, 1, 1, 1 ],
    ' s! j( H7 M) c8 v    [ 5, 1, 1, 1, 1, 1 ],, O9 B% Y8 t, s1 S) T# l4 x" g
        [ 4, 4, 2 ],: M& H( W/ D* _
        [ 4, 4, 1, 1 ],7 Y& v! d$ v( V8 w; g* s
        [ 4, 3, 3 ],- g/ _) y7 d2 O- T3 I2 }& V. z
        [ 4, 3, 2, 1 ],
    ! [3 c+ S7 V1 D) v0 v, C* @8 o    [ 4, 3, 1, 1, 1 ],5 G/ _6 |9 |8 Y' G, h* R& ]
        [ 4, 2, 2, 2 ],, k  b8 Q0 z% _7 }- \, X6 v
        [ 4, 2, 2, 1, 1 ],
    , z! Q, k: {/ y, F7 v5 y' R    [ 4, 2, 1, 1, 1, 1 ]," R" x4 x5 {( {" L# x" B
        [ 4, 1, 1, 1, 1, 1, 1 ],
    & l' n$ Y' ]5 V6 ]: j3 @    [ 3, 3, 3, 1 ],7 ?& v; n0 s  N
        [ 3, 3, 2, 2 ],- X$ G/ M2 c  x% i7 K% I; R
        [ 3, 3, 2, 1, 1 ],0 e! ^7 O9 c8 Q: o
        [ 3, 3, 1, 1, 1, 1 ],
    2 u: n+ i& q% d0 J2 w. t    [ 3, 2, 2, 2, 1 ],
    : ^: {: T2 {! v$ [& w    [ 3, 2, 2, 1, 1, 1 ],5 P$ U+ g" U  Y1 }) {" u
        [ 3, 2, 1, 1, 1, 1, 1 ],  u6 Z3 G. E/ L! s' F4 |' G+ s1 }
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    9 B7 y+ @% H$ ]! {/ [    [ 2, 2, 2, 2, 2 ],
    $ J% r7 z  u+ V/ d' H# G    [ 2, 2, 2, 2, 1, 1 ],3 A6 W; q1 }- M5 G: ^
        [ 2, 2, 2, 1, 1, 1, 1 ],
    7 {7 g2 `2 I" ~* O, r% C    [ 2, 2, 1, 1, 1, 1, 1, 1 ],2 _1 e- d. D3 m7 F$ K- N
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],
      v; ]7 J; e5 a, t) W* p    [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    6 b2 j( V9 S2 t: X! Q6 y$ i]
    回复

    使用道具 举报

    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-4-10 10:42 , Processed in 0.389888 second(s), 92 queries .

    回顶部