QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5375|回复: 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 编辑
    . G/ T# w( u6 U# f+ ]$ G, h' V9 ]% f9 U5 D
    $ [! h0 u1 W1 }- U% f# x8 Rn:=12;n;$ r4 N8 L( |& ]7 w/ t" H+ d
    Factorial(n);求阶乘
    7 @8 m8 |% }' x8 x2 lFactorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
    ! X0 t0 B1 ~4 P1 a. jNumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);  _4 A/ j8 f0 L. p+ B
    NumberOfPermutations(n, 4);
    7 G' L) C2 Q4 c* SNumberOfPermutations(n, 11);4 Z! j1 D- d3 `& I: Y
    Binomial(n, 1) ;二项式系数Binomial(n, 2) ;
    * k3 {) q7 ?, O% P" F2 d% r. aBinomial(n, 3) ;
    % u6 c- ]$ E- G+ A. O! |Binomial(n, 9) ;
    2 k! B; Z6 w) l$ B. ]Binomial(n, 10) ;
    ! ^( G+ \& o6 V: n, ?Binomial(n, 11) ;
    9 x6 h# V0 t' j+ v3 LMultinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600$ r7 K7 y: h1 r& F
    6 p* y, n9 }" G
    Fibonacci(n);斐波数Fibonacci(n-1);; Q: O! x$ T& c$ _; W; ?
    Fibonacci(n+1);: Y; V( Y# e: c6 v
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;
    ' w/ q( \: V* n1 v, yGeneralizedFibonacciNumber(0, 1, n) ;: x+ m# m2 k+ A1 S3 {
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!)). W) j+ }/ L  j# \1 o
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;" J& o* q5 K$ ]" E+ R( K
    Catalan(1);
    ! {6 Q, g" ^, S/ v4 b& E0 iCatalan(2);4 m0 C9 A5 o- B- t; Z+ a5 G
    Catalan(3);Catalan(4);1 ^: g4 y5 e5 ~7 B1 {$ f; ~' B
    Catalan(12);
    / h7 w$ T6 z: Q" o# i/ L  p# Y6 b/ }  h
    Lucas(n);卢卡斯数
    1 U7 s/ _1 }; p- K0 R$ ~2 l) O12
    8 P. P5 d& i: M( ~% O9 T479001600- f5 f7 ~) b8 i* A# s! Q/ b
    831600
    # i, G" `. D/ f! D12
    9 U1 {: O0 l3 r. F/ y6 v  g132
    1 E  U) i( v/ z0 u; c! ^; C11880
    . E! r- z1 y3 t8 C4790016004 K- Z% C: O, ~6 N: Q) Z4 g
    126 i- L2 {8 W& Q( ?+ v1 n
    66
    5 O* n% b4 q: E220
    2 ?+ u- k8 E3 K/ r' [220
    " x, x' C* u- C66: P* l) \* @& {3 [0 {
    12
    9 E  \$ o* W2 }/ V1 w831600" _. R) ^: X2 u& Z
    144
    9 L4 i( a# I4 C897 C  g+ R: S" Y
    233
    * T! e0 `) Z3 H4 U, L3 N233/ C, M4 y9 |" X$ i5 E0 g& d, A
    610
    # O/ N6 D  I8 M0 F% }8 m144" q* v; V$ O  |# q1 @  {9 e: I
    208012
    - y# }4 Y5 ]9 W* {' F7 Z( U9 `208012
    / c- C" F% U( d; f6 j/ O1 J) @1: o; o4 Q- v0 {
    24 d, Q: {) `- y, d8 o  w
    5
    1 c9 @4 \! ]9 ?, T) X0 z14
    ; B- E" X5 U6 }# ?208012
    ) ?9 @5 a4 F# X# K- {; N& s3227 W4 [6 L; f: u* ?
    3 ^4 M2 n1 E' O( R- Y* K/ q

    , B' L" p' `) g3 C% C* r8 z1 }卡特兰数=(2n)!/((n)!*(n+1)!))- P) x+ U6 c+ f' c
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: ' y' o# M, G! s; W
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY
    ! j: t" h! e" E1 q: i  x将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数: " K; [6 B4 B2 L! p9 O; {
    ((())) ()(()) ()()() (())() (()())# N  e/ n8 |$ {) N; b  j0 ?7 O
    Cn表示有n+1个叶子的二叉树的个数。
    ! k! f# F! J  Y& ?  O
      y# O. P' o" {. D3 }' a% s; y5 H$ CCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。)
    " `, E3 y) V$ K0 H* xCn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    , {; b5 e4 }0 w& ]8 \
    ! A* ^6 n! S! {2 T) V0 aCn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。 , V7 }6 l2 J- W5 V2 M1 f  q
    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.
    ; v# b' Z0 Q* iCn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:
    " |/ R+ X; i' Y

    , X; @$ s, x  G9 v8 L" c4 p5 V# j+ b9 r  d5 R& A! [4 Y# a2 C6 A) p" V
    ( t( M) l7 v$ a& Y  Y) {! q
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。. h6 s* D; q! A5 ?: e( P; }
    # _  b' o5 ]0 n
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同
    3 r) ], S3 J7 H5 N) ~  Q4 C3 D
    ( L8 D  o( @6 |+ C  T
    1 k: |5 |, Q: I
    6 n( f( o) E1 e  ~3 x1 ^
    n:=100;n;" Y( g( y$ |/ e( `5 _4 b' f
    a:=Lucas(n);a;6 `. j1 y" v: Y2 y0 e; ]+ p9 Y4 R
    b:=Fibonacci(n+1)+Fibonacci(n-1);b;
    ' t; X) e3 \  G5 hLucas(n+1)+Lucas(n-1);5*Fibonacci(n);( G7 ~# f1 Z% a/ c7 M. U$ M

    7 ?0 Q6 x" M; w) c100
    + Q* E/ b, ~* \  |2 L7920708398483722531271 ]! g. K# O8 i6 r$ f
    792070839848372253127$ F" D0 ~+ [! w1 I
    1771124240896309575375
    : I3 T4 M' V% b1771124240896309575375

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

    12.JPG

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

    74

    主题

    6

    听众

    3304

    积分

    升级  43.47%

  • 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 编辑 * R8 ~% A/ f- O! i: M% x( k
    ; x" l/ L  V7 X- c
    反费波那西数列反费波那西数列的递归公式如下:
    # I  M# y" M3 d2 B0 Z8 j) d; ~7 i: o& C! g$ X, \. N5 X, j2 W) n$ G
    Gn + 2 = Gn − Gn + 1 3 l$ Q7 a* F4 K8 U) _' D
    如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...$ c2 Y* x  l/ ~  p

    - m5 X6 F/ l' [. L* t3 \即是F2n + 1 = G2n + 1,F2n = − G2n。
    & X2 u: ?3 u+ }+ [
    9 l; \3 C6 _$ D; {! S" ~Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);# Z7 g$ n- E( I  f9 z9 b# V
    StirlingFirst(4,3);! r- p3 @- B& E9 l
    StirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);5 d! p9 n8 z) |; h; b" ?+ z. g+ W; D8 h
    StirlingSecond(4, 3);
    ! I( ?: C" Z! z& B: ~; _% Z+ U2
    ' R* K$ H0 w3 O52
    " R8 U$ c# E6 D+ [& W5% K% D* `$ v  D+ C* W$ G
    15
    3 f0 S% J7 y  e" E* I# J-6
    0 D, m- `' d' S+ b1 A! A11" o# D6 W) z" m' a
    -65 y) {. F3 m2 y, t
    11 f7 S, w$ d. a8 x9 Q) J
    7
    # T% |4 S: w7 c# B( U! m' F& ?- V6
    $ v/ b; r: [. E6 h1 j5 j  O2 T$ @' w$ e# e( A0 B
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    ' }! Q! a0 p) z/ y' q  F1 Q( ?. K' r$ y
    {{a}, {b}, {c}} ) F( _# }! Y% V7 }
    {{a}, {b, c}}
    8 A' V* V/ d6 I( X{{b}, {a, c}}
      O2 l6 V5 v7 J0 F( j1 F{{c}, {a, b}} ! i+ c  C0 v" p( |4 i$ `# d
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s, c3 B4 W: D+ X1 ~9 {: N
    s(n,k)是递降阶乘多项式的系数
    " H* }2 o( R! c0 \6 g有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)+ @0 B: H) I3 Q' Y+ G* v: Z

    6 P. {2 r' O; e/ a换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):
      r+ L- D* B' d- S+ N# I" {# A+ {8 k; f7 T$ m2 S
    {A,B},{C,D} - n8 [" }9 w+ F9 D
    {A,C},{B,D} 4 I2 e8 d( K% ]5 x# Z( m
    {A,D},{B,C} 4 W9 g2 q9 [# L3 Q2 c
    {A},{B,C,D} * X+ K2 I$ `. h3 C3 z" ?1 \
    {A},{B,D,C}
    + C) }1 {1 D) n9 Z{B},{A,C,D} - q* B3 `2 [' D' k
    {B},{A,D,C}
    3 r. v6 v8 u+ f+ |6 }{C},{A,B,D} ; R9 B' ^  I" R; ?+ u3 a, S# a
    {C},{A,D,B}
    # d) m: O: t* {& J5 T' L, n; J{D},{A,B,C} # g' I/ ~+ K; q; N4 M. ~  m
    {D},{A,C,B}

    : q2 {& J( T" S( G  n第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    3 p* U' ~! \3 T( X  b给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) 6 @* S% u7 S) \1 @. Y- a  ~; R: m
    S(n,n − 1) = C(n,2) = n(n − 1) / 2
    * t) x3 w: m/ ~6 k; cS(n,2) = 2n − 1 − 1
    % A: L* l) j; a* `. j( K6 Y% S$ |. g& `- Q- @
    换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:, V0 C  c' @' w4 I
    . y- O6 k$ }; S6 T3 Z+ g
    {A,B},{C,D} ! I! ]: g, \8 F; }" m+ q5 c
    {A,C},{B,D} . r6 l1 t4 B/ V! O5 j/ n
    {A,D},{B,C}
    6 F3 j% }. N) J{A},{B,C,D} 2 \/ c2 b8 o" _+ P+ ^( B1 P/ e
    {B},{A,C,D}
    9 }+ D- g9 j5 l8 k9 @7 T{C},{A,B,D} , L* y/ `7 j. f; D! N/ y8 E6 c
    {D},{A,B,C} ( H: T/ A# `! J) G  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 编辑
    + _& Q; l/ }6 C4 E; i! ?( r( u. V8 a0 K$ p% c  C) T1 w
    n:=5;r:=3;
    ; f/ [; F3 F5 Y0 }EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;3 F0 Y2 E9 p$ ~
    BernoulliPolynomial(n) ;伯努利多项式
    4 Q, a* d& ^7 h( ?% o
    & K2 o! g& p# x26' G7 ~5 I2 k: Y' a0 f  o# _
    137/60. k/ n8 x9 Z; {& k0 y
    0
    , R+ x; K9 s4 s  w) d0.000000000000000000000000000000' j! \3 ?) v* D7 L  A8 z, {
    $.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑 0 N8 u+ W9 p- y2 w% F+ A1 {) x( N, e$ x

    % w' m- c6 K( `; D  k2 a/ L1 J9 ]9 d3 d
    " b7 {- U6 X9 q7 \% K4 l& N
    伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    / s0 T2 L- U* ?4 x0 j
    - k6 ^& n, U& }0 N% H) M伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关8 ^+ W+ s! N: [3 g$ o
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 ) g6 W3 c2 o( X1 T, T
    8 P8 A3 J+ ^! n: a& y& _
    拆分 。。。。强!
    + o2 T( ]3 ^0 U! ^& f, j  C, {8 C
    $ F$ Q4 A! s" A2 W  P8 s6 PNumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;
    " \* J' \2 J+ `/ c1 {! t2 H- m  u8 K7 x! @
    79 r4 b6 u! T7 q; N
    190569292/ }7 I3 I- ^' I. k- ?
    [2 X: F7 |; c; T1 A
        [ 10 ],4 k1 J0 x7 K1 d5 \& i' e
        [ 9, 1 ],
    2 ^! c4 u1 q& T4 O    [ 8, 2 ]," T# |: ^; K' U  g# T; c6 Z. j
        [ 8, 1, 1 ],
    # _8 I! g5 [6 _    [ 7, 3 ],+ X) c+ \: A5 n) j# R% \
        [ 7, 2, 1 ],) y  \8 K0 C6 L
        [ 7, 1, 1, 1 ],
    3 S% K. D+ |- ~* |    [ 6, 4 ],/ _0 K' U5 j7 O9 J- _, X
        [ 6, 3, 1 ],
    & g6 t8 N* E. o( o$ j! h4 x. V0 U    [ 6, 2, 2 ],' s" M5 A/ }; g" u4 m
        [ 6, 2, 1, 1 ],
    & a/ r$ U* G- |+ l    [ 6, 1, 1, 1, 1 ],
    ! f9 c* |) [# h0 _/ p( W    [ 5, 5 ],8 q6 e$ f3 r- v  W+ v: L7 |) Z
        [ 5, 4, 1 ],
    & F  r% B! d) s/ R    [ 5, 3, 2 ],
    ; g0 i9 s1 T  \    [ 5, 3, 1, 1 ],
    ( m' ]4 @6 N4 [1 H9 u" |8 }    [ 5, 2, 2, 1 ],
    % Q0 u. t% t1 S: ~7 v" {" d    [ 5, 2, 1, 1, 1 ],
    - G! T6 I- l2 g3 y0 ~; O* V    [ 5, 1, 1, 1, 1, 1 ],6 j' O# K0 f% \1 Y' |- d
        [ 4, 4, 2 ],; g  A0 Z+ O3 X
        [ 4, 4, 1, 1 ],+ p8 C' f# ]7 W$ i& T) }8 Q
        [ 4, 3, 3 ],
    4 o1 L2 j9 R9 X. q; @- z    [ 4, 3, 2, 1 ],
    . [; E* X. }8 `8 M9 p1 S/ a    [ 4, 3, 1, 1, 1 ],
    ( I( E8 q6 [8 F+ K    [ 4, 2, 2, 2 ],
    2 F( V+ ]2 L- a7 M! @$ L+ s    [ 4, 2, 2, 1, 1 ],1 Z$ ^" m$ B3 r( A: U& L  m! @8 b; X
        [ 4, 2, 1, 1, 1, 1 ],
    , R% T& F5 B5 r+ e5 B2 r    [ 4, 1, 1, 1, 1, 1, 1 ],
    9 t% X/ l1 W$ T4 m  {    [ 3, 3, 3, 1 ],0 A  S* A8 ?5 b( _
        [ 3, 3, 2, 2 ],3 T+ R4 @' Z. F1 j) u
        [ 3, 3, 2, 1, 1 ],7 x2 l( W; ]! C# X, O+ C2 e
        [ 3, 3, 1, 1, 1, 1 ],5 F* k. k8 g+ N: L6 m2 {9 Q% G
        [ 3, 2, 2, 2, 1 ],
    : K3 e, T( _! b" o+ `* b    [ 3, 2, 2, 1, 1, 1 ],
    ) G, r- r% M5 j  [/ N    [ 3, 2, 1, 1, 1, 1, 1 ],) R% x+ `: w/ [1 E) T
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],( _6 I% c; u  d( f
        [ 2, 2, 2, 2, 2 ],- z9 ^" R/ L* P$ e2 @8 S
        [ 2, 2, 2, 2, 1, 1 ],* ?+ d; x6 k' i" r1 j
        [ 2, 2, 2, 1, 1, 1, 1 ],
    * s4 t- H/ L8 z% A& y. N    [ 2, 2, 1, 1, 1, 1, 1, 1 ],
    0 [) M8 v8 u' a) \) o$ _    [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],' U. G/ N8 F: a- ^# e& k( U
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    & Z! e4 t! ?+ V' M% t/ Y; ?]
    回复

    使用道具 举报

    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-10-15 06:10 , Processed in 0.394011 second(s), 91 queries .

    回顶部