QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4926|回复: 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, L" A; w0 n; c7 H
    ( ~+ P7 V+ Z6 @$ k3 B0 u  h
    n:=12;n;
    % c2 u' I; a1 T- XFactorial(n);求阶乘
    " \1 q7 S9 }' M$ R7 v5 Y5 z( `Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
    , K' {9 K: `: HNumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);
    % u) r7 [# _- y" F5 U) yNumberOfPermutations(n, 4);
    ) l, x  l9 |* |/ QNumberOfPermutations(n, 11);
    & K& }  W) ~9 p& K3 T8 z! d+ NBinomial(n, 1) ;二项式系数Binomial(n, 2) ;
    8 Z% |9 I0 W) m$ D$ uBinomial(n, 3) ;
    " a$ ^. d8 f( J+ q& uBinomial(n, 9) ;( z; W7 x* r0 Z1 U. h9 O& u
    Binomial(n, 10) ;
    4 G0 X5 ?7 Q3 ^; g5 tBinomial(n, 11) ;
    ; b) H5 V6 F. p0 C9 y2 AMultinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    , D; Z7 ~0 S$ d% V+ i. z+ I4 z& D1 Z
    Fibonacci(n);斐波数Fibonacci(n-1);
    , f  D" P# U) lFibonacci(n+1);
    ; m- C- k1 [" |GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;. \4 L/ ^( D. s! k9 z/ c! X& O' w  r) R
    GeneralizedFibonacciNumber(0, 1, n) ;- p2 R1 x. S3 I9 L
    Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))4 w1 J6 R8 {$ M0 y$ N, \  t
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    ! e* v) ]/ \3 _9 gCatalan(1);# f8 w7 Y. Q: t
    Catalan(2);
    ) s  T2 X% r, ~* RCatalan(3);Catalan(4);
    - ~: J4 l6 V7 T& }$ iCatalan(12);: K) g# c* W( G9 R/ c" X1 c
    ; D2 h6 ~3 p& \9 _/ A
    Lucas(n);卢卡斯数
    6 |" h9 d  f+ K: i( g12
    4 Y' p% G1 A. q% |3 g) D. G479001600$ d, w1 U6 ^; l9 ?8 G! p' Z( l* k
    8316003 D6 |5 r. s# `( g: U$ u' O
    12
    % T# ^' b" \% h9 m' @# X6 }( w1324 [: X6 u7 _0 p! F, J
    11880
    ( T! m: J! h1 M479001600
    2 T3 G: `  v! i  O3 `# L12
    : k. e7 y/ q/ T6 s8 ?1 M5 ?66
    . \- d$ D$ n8 h0 m9 d" h  L  M220
    5 H2 L( N: c0 o) _220. F# V- Y+ X" H1 O- E. [* V5 `
    66
    1 S0 h$ ?. F% r12
    0 `  F0 _$ `2 C: R831600
    % P& [8 G% I8 j6 {144
    8 X* V' ~# y5 M4 i0 T89
    9 [) G1 ?. Z6 y1 n, D0 v233
    & B3 v0 C3 e9 L+ C4 w  G0 m233
    & O. }! r; Q, k, A- q. S610* \' u+ v4 W% H+ K9 S1 W" |: J9 V
    144# [1 x4 o: h7 t2 O
    208012
    2 {/ ]7 h8 ~) ^- Y$ j208012
    6 R5 F9 x- m) `1  ?. p! c4 C2 c; q* D# i- [5 n
    21 E4 f1 }. s2 Z- W* E; L4 S. T
    5
    % M  g, ?3 _( o# W5 k14' p4 u$ ~0 ]; t
    208012
    1 w) G1 o# M3 Q322
    ' S% K1 L/ e0 s/ I$ H0 q# d  C' b5 I, H) e" y! ]

    & [4 a+ [* f* _卡特兰数=(2n)!/((n)!*(n+1)!))
    8 o7 E! q- S+ N* K8 s( s% @: zCn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words: 3 w& o' l' t2 B
    **YYY XYXXYY XYXYXY XXYYXY XXYXYY
    7 T& ]( K) O% S8 O* ?' ]: a将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数:
    9 K; O8 G2 K" \" R((())) ()(()) ()()() (())() (()())/ X# l0 S! g; W/ H5 H0 `
    Cn表示有n+1个叶子的二叉树的个数。 : [1 B5 p2 C" g$ t8 l7 P

    ( a( i: Y& P; V$ F$ G/ C9 h% i2 DCn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) 6 J4 n, ?, a4 b
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    1 t& f) I9 q, g, @3 \3 D
    # }- ]1 E$ N& ^: D9 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为所有含一个元素的数列的单位元。 2 L1 U  u. Z* s! I! K7 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. ; S9 |  t3 H" @' N6 I% c
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    ) [. |% @7 F1 ^- _( F2 y  B0 O& A. e* w7 Q+ X! \! x8 e3 z

    7 a- v1 ^$ Y( G& Z* n! h0 ]. W: x' `4 _
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。( D. i/ C; G1 g- `  z' B
    & t% n, P# K3 t5 M. `1 O
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同
    5 u7 |" A) B$ k

    ) [; P7 [6 V7 Q- {8 H0 D8 w# d( Z' h6 k" a9 P" x7 ^
    & [9 f3 H7 ~. v' ?# p+ G
    n:=100;n;
    0 M) B* Q* E) [( `4 `5 Ra:=Lucas(n);a;
    ( O1 J" ]2 j7 F3 s6 D* {" Tb:=Fibonacci(n+1)+Fibonacci(n-1);b;6 c6 O5 x% S# y
    Lucas(n+1)+Lucas(n-1);5*Fibonacci(n);
    ( E# U% H, P  Z  j# h2 y1 U& y; ]: P9 v
    1004 J" K8 L9 n! z
    7920708398483722531276 {* ?+ G$ U2 q$ R. f7 h/ H, A
    792070839848372253127( a6 D  E5 ]2 @$ E* j/ `, N) W3 p% j
    1771124240896309575375
    6 ?2 ~& e9 T$ d" m* r0 r# `1771124240896309575375

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

    12.JPG

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

    74

    主题

    6

    听众

    3292

    积分

    升级  43.07%

  • 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 W8 K& r1 H0 Y& s1 t' W- O- C9 r
    . x! ]3 t: e% W9 _+ h& K% b反费波那西数列反费波那西数列的递归公式如下:# v0 C% Z0 V' A( U% c

    3 @( E7 }+ b9 B6 V8 T6 J5 n( R7 H* |Gn + 2 = Gn − Gn + 1
      X0 `8 r; K+ e如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...
    - M3 \: k8 B% [) c) @  ]3 l* L8 d/ _* l; q; Y
    即是F2n + 1 = G2n + 1,F2n = − G2n。- z! n( F$ q& h

      d" L% h) G; MBell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);
    ; j1 p3 I4 N5 g( X! H4 AStirlingFirst(4,3);
    ; g/ t( F+ {1 m2 mStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);
    " z( w( i3 l# g4 P8 j2 J- c% jStirlingSecond(4, 3);1 b: Q$ p: G7 n
    2
    0 O2 b1 k' c8 w52( g. o5 z4 ]; B( g% A
    5% J/ s" g% m+ @! b( }" g
    15
    ! c* I  S# w$ Y$ x# ^+ R; d-69 G: L) H) v- @! n3 I
    11
    0 N/ s6 f5 c0 @6 `% o-6
    ( d4 j' P( Z' ~17 g5 o2 O* C7 Z2 ?) ]7 F+ B2 n: @' m1 @
    7
    , H( K" k# U5 s" y' |+ q  d6/ e; U! p' X/ b- k. ]4 d

    5 x/ a1 U0 u) J: Y# hBn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    # @7 _) m4 Z5 Y  m
    # ~7 v8 \" I1 M# ^7 D( k+ v. ]{{a}, {b}, {c}} 7 ~& @) z' v% s( T; h  ^
    {{a}, {b, c}}
    ) \% b8 ^. [4 V% U$ f( L; T4 l{{b}, {a, c}} ( m- m: u" k/ X$ \' ~; l+ n
    {{c}, {a, b}} 2 y: Q0 C- \9 _$ U+ a" k! l# n
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
    2 ^+ R# H4 s% @( {, n! @. Rs(n,k)是递降阶乘多项式的系数
    9 l0 l" B: z; Y; F3 X有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k), y' X" m  a6 v
    , P/ |& p# L6 m/ E% F" L: {
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):  n& }5 ^9 v, x/ I4 }7 ]

      u7 I  S6 i, j5 r' d{A,B},{C,D} 6 L4 j! ]& P9 ]
    {A,C},{B,D} ( i0 F$ Y( [' E8 W" w7 j  R3 T  D
    {A,D},{B,C} 7 |- f& O/ }0 N# ]2 ~% ?! T
    {A},{B,C,D}
    5 H5 T# {' w$ }) T{A},{B,D,C}
    5 s! J9 B/ c% L! e{B},{A,C,D} / N, h% i  ?" ?* U7 b4 |3 \8 D
    {B},{A,D,C}   P5 x5 K/ _0 M
    {C},{A,B,D}
    6 S' v% |( e3 o{C},{A,D,B}
    " [0 G, V# i4 p7 Q( i$ t# @{D},{A,B,C} : e) V% s* E9 n  j9 H) }& S
    {D},{A,C,B}

    * E8 ^- X4 d+ d& D. s第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    9 h9 o) r& L; y; k* c给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k) $ T. L. m3 G' p0 p5 U
    S(n,n − 1) = C(n,2) = n(n − 1) / 2 ) B8 s- f# F# B
    S(n,2) = 2n − 1 − 1 ' r8 t0 @$ z  D. w* v

    & w0 z) r' {7 b* o2 Z, i" L- [换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:
    # T% [9 C+ R' X# G& H- B* x( C* R
    4 ^; T0 d, e7 @# `! U{A,B},{C,D} 6 e; ~3 @% W6 s0 m
    {A,C},{B,D}
    8 a/ L& Y- L; P9 j- B{A,D},{B,C}
    3 M' s# F2 I$ a, `! m* k) A% f{A},{B,C,D}
    ! t2 _1 c. G. y. X2 ]{B},{A,C,D}
    $ j# K8 @3 a+ }5 j+ u  }{C},{A,B,D}
    " A- F( @% F" ]. t- X{D},{A,B,C}
    % w+ t# g, B: V* @1 w% V! x$ t( a因此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 编辑 ' s3 x8 k0 V- V% ]- J. u

      P9 k1 H9 g, Y: Gn:=5;r:=3;; A( x$ F4 x& o/ S6 a% B) G& S
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    $ B. I2 A, I. j; u$ Y6 J# iBernoulliPolynomial(n) ;伯努利多项式
    - z- I; Q+ l/ E: @* T* k, S, S* P! p  E" i+ f+ ^9 p; P- E
    26' ]7 k' y+ R$ _; L" a+ M7 N
    137/60
    % A& t) V9 K1 o9 @4 q$ d! q0
    ' Y$ Y" Y" I& x2 j4 @0.000000000000000000000000000000  L* F! K# ?3 S% h& A7 d
    $.1^5 - 5/2*$.1^4 + 5/3*$.1^3 - 1/6*$.1

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

    22.JPG

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

    33.JPG

    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 19:56 编辑
    5 L( G% Y9 d! Q- t" B. B( C7 x4 a: i3 X; {, X
    , R5 }& L1 ]2 b; J5 l
    % ]$ o* u3 r3 P# n* m# A
    伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    6 y; Q  }2 S) j2 S$ _' F/ H2 j! v$ x' S3 ], `
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关
    ) e& U9 O2 g) J) n
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 / b8 Y4 c( l# x4 T! H% K3 W

    & k1 g1 w. I; V/ _1 a+ j2 }0 T拆分 。。。。强!4 X) T0 Y1 w! o0 m7 h

    ! s: O! V) S2 M# {" G; [NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;
    $ N4 b6 W; z; E+ F3 A; d/ O, J8 l$ ?/ o5 A, A/ L; z( d" U  K
    7; _$ X7 F2 I- C) O; Q3 o2 K
    1905692921 g" ~* Y) `8 A, M4 O3 u
    [
    9 l; I7 s2 r+ |0 ~( b2 ?* f! L6 y- u    [ 10 ],1 S5 {1 s/ V" {  V
        [ 9, 1 ],
      X$ s3 m3 c- S    [ 8, 2 ],
    ; }/ b1 n# E; h  @: B, r3 I$ T- M. U; W    [ 8, 1, 1 ],. d& c) s1 Z: Z$ D3 P, ~) ^
        [ 7, 3 ],
    9 b+ t! A4 j" J    [ 7, 2, 1 ],' d5 d0 J  r) Q3 E! i+ }8 M* b1 [
        [ 7, 1, 1, 1 ],
    9 a0 \5 v+ I6 x1 u3 ?- E% s    [ 6, 4 ],! a6 c& h6 ]6 \) l7 O$ [
        [ 6, 3, 1 ],  m6 x. K5 c  [' y" b$ [( u
        [ 6, 2, 2 ],
    ) [- H5 f$ Y- {6 [# Y    [ 6, 2, 1, 1 ],% a* v4 ]' k  ~1 Q$ L
        [ 6, 1, 1, 1, 1 ],% Z7 I$ r7 J" @( Z' H
        [ 5, 5 ],$ _1 q' i5 E4 f* W
        [ 5, 4, 1 ],
    1 l5 H4 H3 i4 @' U  @- i& k1 y    [ 5, 3, 2 ],2 Z% [  q! u0 C; {
        [ 5, 3, 1, 1 ],
    5 H* I7 o* N3 G/ e* h/ J( e4 l  o    [ 5, 2, 2, 1 ],1 z- L3 f4 w7 [; @
        [ 5, 2, 1, 1, 1 ],3 q2 {, b, S$ |+ Z5 s
        [ 5, 1, 1, 1, 1, 1 ],: D7 w9 b9 D0 ~
        [ 4, 4, 2 ],
    / R$ f% k# a( u# W    [ 4, 4, 1, 1 ],, Q( E# F+ S3 z7 m1 Q
        [ 4, 3, 3 ],
    + g, w$ b9 M+ n6 i) P2 N" }6 Z- o    [ 4, 3, 2, 1 ],
    ' A+ D) q! z5 B    [ 4, 3, 1, 1, 1 ],3 l) Z8 K8 @- u1 q; h0 d
        [ 4, 2, 2, 2 ],
    8 u% X3 y, n" @' b- ?( P+ G    [ 4, 2, 2, 1, 1 ],
    3 l2 k9 n) z1 E; c' ~5 S3 Y# L    [ 4, 2, 1, 1, 1, 1 ],
    " j& L+ \% q- C: C! B    [ 4, 1, 1, 1, 1, 1, 1 ],% T3 o  c0 v+ S$ }% b; t
        [ 3, 3, 3, 1 ],+ G; t8 I8 x" Z& W7 G
        [ 3, 3, 2, 2 ],
    1 I8 _# V" |0 x" D$ [0 L    [ 3, 3, 2, 1, 1 ],
    5 p/ I! _7 U" y! W9 `* ?: \    [ 3, 3, 1, 1, 1, 1 ],+ p* f4 O2 t8 P4 w% u) P
        [ 3, 2, 2, 2, 1 ],7 J  K/ T* h6 M3 W" _
        [ 3, 2, 2, 1, 1, 1 ],
    ( L4 p/ ~" l/ Q    [ 3, 2, 1, 1, 1, 1, 1 ],6 a3 ^4 i, D. g* h# ?
        [ 3, 1, 1, 1, 1, 1, 1, 1 ],
    ; [+ E" s/ C9 t# o" V6 E! O    [ 2, 2, 2, 2, 2 ],
    , V5 ^& D/ Q! m    [ 2, 2, 2, 2, 1, 1 ],
    0 J; J. ?5 T6 a    [ 2, 2, 2, 1, 1, 1, 1 ],
    1 K* s) Q0 S' f3 @4 z" j3 O. k    [ 2, 2, 1, 1, 1, 1, 1, 1 ],( W: O2 i# u2 \
        [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],/ N& P% w7 B0 t1 T0 Z+ {, P1 Q# g1 y
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]2 M: {6 m8 G" \0 j. {
    ]
    回复

    使用道具 举报

    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, 2025-9-15 01:48 , Processed in 0.704408 second(s), 91 queries .

    回顶部