QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 5367|回复: 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 编辑
    ! \# W" H9 k5 E! Y! i
    + U) F) \4 L$ r$ xn:=12;n;7 A" [. ^& E  v
    Factorial(n);求阶乘
    8 E+ L3 q) \- K" i1 r% {Factorial(n)/(Factorial(2)*Factorial(2)*Factorial(3)*Factorial(4));
    5 ?; `1 V, M* ]/ ^' _. O  X) ANumberOfPermutations(n, 1);组合数NumberOfPermutations(n, 2);2 A, M& D4 {- p8 V. F6 x2 x5 p
    NumberOfPermutations(n, 4);. D( Y) W1 M3 b/ V) f/ Y
    NumberOfPermutations(n, 11);
      u/ p$ O" ?; k! S$ w* y2 @( P* d/ TBinomial(n, 1) ;二项式系数Binomial(n, 2) ;* |1 W9 J$ l6 A1 ~. [6 Z4 `
    Binomial(n, 3) ;: }- k! f% w; G( \) K4 e2 d; G
    Binomial(n, 9) ;2 P9 `) \: O6 a6 w
    Binomial(n, 10) ;( M/ u6 _) q7 [8 p$ W
    Binomial(n, 11) ;
    : n1 L4 x; N% v9 oMultinomial(n, [1,2,2,3,4]) ;x*y^2*z^2*t^3*k^4系数=12!/1!*2!*2!*3!*4!=831600
    3 x* F# r3 I  O
    ; j0 J) ~8 a/ [/ C' {Fibonacci(n);斐波数Fibonacci(n-1);
    8 d2 Z6 B2 B: [7 ~- WFibonacci(n+1);- h: A5 X2 ^7 N3 f( c/ Q& R
    GeneralizedFibonacciNumber(1, 1, n) ;斐波位数加数GeneralizedFibonacciNumber(2, 3, n) ;
    6 s# M8 t4 y) _0 l3 D4 v+ W* aGeneralizedFibonacciNumber(0, 1, n) ;
    * q9 p2 p" m& Y$ N$ ^Catalan(n);卡特兰数=(2n)!/((n)!*(n+1)!))# R$ j/ [: p9 r" U
    k:=Factorial(24)/Factorial(12);k;m:=k/Factorial(13);m;
    ! l8 o( x+ J3 B, e6 jCatalan(1);( P) u( ]- S" J4 q
    Catalan(2);
    3 M) O) k5 [5 t/ Y; B+ LCatalan(3);Catalan(4);
    1 n$ }7 F- r) C% @: TCatalan(12);
    & |" I5 g8 k9 k  F* K
    ! n0 n8 B& N+ N8 k0 PLucas(n);卢卡斯数
    5 I" E& C$ P; _1 X* ~12
    8 p, {# W% m) d# Q$ I( T479001600
    & ?. s& x! R- c% s% V, U/ L; S' ^8316000 q0 U/ R4 _0 V, R( ^. }6 ~7 a
    12
    ( S2 x' ^; y9 ^132, w6 u' W2 x2 y# l3 h" D
    11880
    # K" J  a! B5 w" u' s9 L/ m6 g479001600
    . C& [+ A1 X7 U; G; l12
    $ I, Z: @6 M6 {3 o. b- G1 }66
    ! E6 b) W' l" d, l220; l$ r# `1 X7 }" F( j9 P/ c6 N
    220
    ! }, n; h: \, ^/ u% E66
    # g3 r1 F" p; g4 T12
    5 }* v7 f' ?+ T, i: K831600
    # K' I) j& ]$ [+ W144
    - \7 K6 R" ?- _: ^89
    ( |& L3 y3 |0 i* \233
    ; M' v$ `4 a* b- I' J! q233
    . R# q( B0 L& `- X9 l610
    ( D9 Y. h  \2 F! X$ c144
    & u1 P8 S. n% J# D6 s5 X208012( p2 p$ N3 K( K
    208012* l" A( m* J1 S
    1
    / H! G1 B/ g6 Q9 e3 I2
    . L/ }! B1 k! ?: D2 C0 O2 ~55 M% g3 }, {9 o+ `! j
    14
    , ]6 X# O2 U2 N' s8 S208012
    % j$ {" t/ _( o322
    7 p6 H/ p' `. c  ?! E2 V
    , {! S) V# ]3 Z& ?" T. G. D* M! i7 G: w; I$ Y
    卡特兰数=(2n)!/((n)!*(n+1)!))% \7 T1 Y5 z1 A' o2 u  {% d/ V
    Cn表示长度2n的dyck word的个数。Dyck word是一个有n个X和n个Y组成的字串,且所有的部分字串皆满足X的个数大于等于Y的个数。以下为长度为6的dyck words:
    + R/ ^( f7 W& ?4 |8 A" p# {**YYY XYXXYY XYXYXY XXYYXY XXYXYY4 d5 _, ~1 d/ K/ u# L
    将上例的X换成左括号,Y换成右括号,Cn表示所有包含n组括号的合法运算式的个数:
    % u% j$ S! W* f((())) ()(()) ()()() (())() (()())
    & P3 Q  {! c  X# A$ Q" t( UCn表示有n+1个叶子的二叉树的个数。 - ^8 F+ G. }& {! Q
    3 o, u. {+ B) A  `
    Cn表示所有不同构的含n个分枝结点的满二叉树的个数。(一个有根二叉树是满的当且仅当每个结点都有两个子树或没有子树。) 6 a8 \8 V) R! j% N4 A2 a+ R) [
    Cn表示通过连结顶点而将n + 2边的凸多边形分成三角形的方法个数。下图中为n = 4的情况:
    ; h# U3 G8 m; k4 I8 c0 E3 p9 W6 ~  J. @
    Cn表示对{1, ..., n}依序进出栈的置换个数。一个置换w是依序进出栈的当S(w) = (1, ..., n), 其中S(w)递归定义如下:令w = unv,其中n为w的最大元素,u和v为更短的数列;再令S(w) = S(u)S(v)n,其中S为所有含一个元素的数列的单位元。
    3 c; K! B9 k* p3 g2 MCn表示集合{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. - q0 u3 q9 ~# g
    Cn表示用n个长方形填充一个高度为n的阶梯状图形的方法个数。下图为 n = 4的情况:

    3 g% M% {  o, D7 }- e7 I
    " k& A' \8 v  y
    + z- ~7 F/ ^1 J& ]5 x- f. L3 D; s3 B' I6 c5 ^% K
    卢卡斯数是一个以数学家爱德华·卢卡斯命名的整数序列,他既研究了这个数列,也研究了有密切关系的斐波那契数(两个数列都是卢卡斯数列)。与斐波那契数一样,每一个卢卡斯数都定义为前两项之和,也就是说,它是一个斐波那契整数序列。两个相邻的卢卡斯数之比收敛于黄金分割比。, W: W. a) Z( a
    + C0 c9 `6 v( i* r$ E: |( m
    但是,最初两个卢卡斯数是L0 = 2和L1 = 1,而不是0和1。所以,卢卡斯数的性质与斐波那契数的性质有些不同

    & {" J0 W$ }8 l! G% S% o% Z' C+ K  ]7 H# v+ `2 g, q2 _2 J1 [7 g) M% D

    0 J1 n: R: @0 D! A5 U; o  Z0 o& a2 P* D* E
    n:=100;n;
    $ S9 d: L+ P$ J1 A( L; x" }a:=Lucas(n);a;
    6 V8 F/ A( _5 V. Eb:=Fibonacci(n+1)+Fibonacci(n-1);b;
    . T. c1 @2 l1 k4 ?; L& nLucas(n+1)+Lucas(n-1);5*Fibonacci(n);: ~8 o; z% j. [
    ) T' i# ?% F) T
    1001 y0 W5 X, ]* `8 q1 x
    7920708398483722531272 |& S( u0 o5 N/ n' m' Q% j. \  e
    7920708398483722531273 r7 k7 u) A. {3 O
    1771124240896309575375' H: m! q8 P' c8 S/ y# \
    1771124240896309575375

    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 编辑
    ; L% Z/ v1 Q1 I, K2 M9 ^/ a9 C2 T; z! Y, I( r- _7 d# I: w
    反费波那西数列反费波那西数列的递归公式如下:
    : k( D' C) {  Y( g3 Y! {4 q: C
      @$ b0 }3 a2 NGn + 2 = Gn − Gn + 1
    ! W  f2 N# S% ?- J1 D如果它以1,-1,之后的数是:1,-1,2,-3,5,-8, ...
    7 ^/ Z% L' P" U- f2 m) {3 F8 n
    4 \' H& M, R0 t4 Z+ x, b3 P即是F2n + 1 = G2n + 1,F2n = − G2n。
    ! {, ^# }2 e, b- n, u$ K1 J, C+ m) A. @. n/ D  ?% t& Y# s
    Bell(2);Bell(5);Bell(3);Bell(4);贝尔数StirlingFirst(4,1);第一Stirling数StirlingFirst(4,2);2 J5 H  {) y5 i! ~' c
    StirlingFirst(4,3);
    , M0 w2 _0 b. LStirlingSecond(4, 1);第二Stirling数StirlingSecond(4, 2);" q# Y* Y; \2 j9 H* b" y
    StirlingSecond(4, 3);
    4 ?  N; G! s! {2' y, E6 z! B# g
    52! W. W, x# `4 A# \8 _. S$ y% ^
    5& N. K& o# ?1 J! G2 O
    15
    7 ?1 S8 L% q7 d9 ], s. ]-60 v( s. a' l6 G* j. S9 y2 v
    11
    , Q; d0 l8 Z5 s* ]! Y-6% |. H/ A; o, U
    1* j; K/ V# n0 F. B9 E+ E* c
    7  O  A' Q6 s# K+ v
    65 L7 x4 R/ p+ A/ w8 k
    4 x8 U! T+ P* d/ m
    Bn是基数为n的集合的划分方法的数目。集合S的一个划分是定义为S的两两不相交的非空子集的族,它们的并是S。例如B3 = 5因为3个元素的集合{a, b, c}有5种不同的划分方法:
    & P7 z; R4 v" F/ S2 U
    5 X" Z1 N# C$ q{{a}, {b}, {c}} % C% D, D' Y/ m7 x, s0 B
    {{a}, {b, c}} % s! k% g. d$ u% S/ \  N  j" ~' Y
    {{b}, {a, c}} 5 E, t1 [* j) y2 W# h/ f. w
    {{c}, {a, b}} 3 Z8 p8 c: q9 T2 K3 v" l/ y' _
    {{''a'', ''b'', ''c''}};
    第一类Stirling数是有正负的,其绝对值是n个元素的项目分作k个环排列的方法数目用小写s
    ; B& B1 B2 a6 K- {! x! N* ys(n,k)是递降阶乘多项式的系数
    / I- M" O) X' b( o! ?6 H有递归关系S(n,k) = S(n +1,k) + S(n ,k-1) -n*s(n.k)) H( V$ p; n. V1 z
    3 [* @5 K: x4 ~# z  S
    换个较生活化的说法,就是有n个人分成k组,每组内再按特定顺序围圈的分组方法的数目。例如s(4,2):! Z" G$ B+ B+ p5 M) j

    ! K8 i2 k3 y0 P! U{A,B},{C,D}
    * c. w3 a# ^! k1 _" W8 ]{A,C},{B,D}
    % `, L9 q3 ?/ _" o+ q{A,D},{B,C}
    $ E2 @3 U8 W+ b9 c7 z- T4 F{A},{B,C,D}
    ) U" {" ^( B, Y* {{A},{B,D,C}
    , E2 M& k4 [8 C' z4 V4 B1 ?{B},{A,C,D} 0 u) n7 ~) ?5 d( {8 o# L9 C
    {B},{A,D,C} $ U; v: Z- c8 B) I* [0 [
    {C},{A,B,D}
    & Z9 h9 t9 o' {* z7 t/ k{C},{A,D,B}
    % v& `2 J0 Q& K1 @: T" J{D},{A,B,C}
    $ s1 K* p$ U8 t" M) M" R& s' Y{D},{A,C,B}

    ; |/ z$ l8 M# I. Z  j第二类Stirling数是n个元素的集定义k个等价类的方法数目。用大写S
    $ k- \1 M5 ~  q; p" c+ q# s+ }8 Y给定S(n,n) = S(n,1) = 1,有递归关系S(n,k) = S(n − 1,k − 1) + kS(n − 1,k)
    & c7 }  D7 c  P: w* sS(n,n − 1) = C(n,2) = n(n − 1) / 2 ) O+ {' f  m* a* a: Y; c
    S(n,2) = 2n − 1 − 1
    - H+ j4 l6 P. y2 y
    4 K- q2 h) @" B: D* d换个较生活化的说法,就是有n个人分成k组的分组方法的数目。例如有甲、乙、丙、丁四人,若所有人分成1组,只有所有人在同一组这个方法,因此S(4,1) = 1;若所有人分成4组,只可以人人独立一组,因此S(4,4) = 1;若分成2组,可以是甲乙一组、丙丁一组,或甲丙一组、乙丁一组,或甲丁一组、乙丙一组,或其中三人同一组另一人独立一组,即是:3 Q( }- P# z* `* B. r) ?+ ?
    % k  p* F$ L8 m/ o% C6 |, i! x
    {A,B},{C,D} + A+ h8 K# Z/ j2 u' ^4 ]
    {A,C},{B,D}
    - T% s# p" i6 R{A,D},{B,C} 2 D9 L) l) R0 v/ R
    {A},{B,C,D}
    4 f( [$ i( }+ Z' Y, g3 x4 _{B},{A,C,D}
    ; t( ?2 K7 I+ F8 i3 @7 Q2 y" Z3 A{C},{A,B,D} / o% V$ W" |8 S/ ^
    {D},{A,B,C}
    . m& S' B( f, s3 f* W% C* 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 编辑
    ( ]1 {3 Q8 b3 A4 }" ]
    ! ^- P' W* B" U! V, p% a. d+ h. xn:=5;r:=3;' T# S5 ~6 U+ q( f: D- Y
    EulerianNumber(n, r) ;欧拉数HarmonicNumber(n) ;调和数列和BernoulliNumber(n) ;伯努利数有时会写成小写bn,以便与贝尔数分别开。BernoulliApproximation(n) ;
    4 V; T$ [8 W. {, fBernoulliPolynomial(n) ;伯努利多项式
    ' @; o7 d% h3 \+ }; {, C1 q5 w! z5 N- d# D5 z2 U" }! Q& Y' S* _$ M+ l
    269 B4 }; |, t. V; O
    137/60) q: Z* W8 c3 y8 [/ d8 o: V
    0
    # x1 N5 l/ v- f# \! r2 N7 O0 E6 w0.000000000000000000000000000000" ]- U% {. L1 y
    $.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 编辑
    9 V4 w2 F( A- p9 M+ l3 [
    + L$ n7 F3 c, f/ A* g. m9 O3 \$ E( @& d/ d4 ?
    " g& L  k7 N) d2 g) k' X
    伯努利数可以用黎曼ζ函数表达为Bn = − nζ(1 − n),也就说明它们本质上是这函数在负整数的值。因此,可推测它们有深刻的算术性质,事实也的确如此,这是库默尔(Kummer)研究费马最后定理时发现的。
    1 Y% G& g. h; A2 b9 ~( W$ g6 O: d0 y8 d& M7 r" _: }
    伯努利数的可整除性是与分圆域的理想类群有关。这关系由库默尔的一道定理和更强的埃尔贝朗-里贝定理(Herbrand-Ribet)描述。而这性质与实二次域的关系由安克尼-阿廷-乔拉猜想(Ankeny-Artin-Chowla)给出。伯努利数还和代数K理论有关  \. g  u  g& ]9 s2 ~1 T
    回复

    使用道具 举报

    lilianjie        

    43

    主题

    4

    听众

    204

    积分

    升级  52%

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

    [LV.4]偶尔看看III

    本帖最后由 lilianjie 于 2012-1-12 20:38 编辑 % ~0 I2 G- M: z1 u
    % f  @6 I- H! Z, ?9 T8 q
    拆分 。。。。强!$ Z0 f7 z/ X( p* m
    " s1 a5 _# L0 {: M) P
    NumberOfPartitions(5);NumberOfPartitions(100)artitions(10) ;  p- N1 Q2 ^$ ?* m

    # |; n* Q, G2 b% D" N6 z; I7
    # J; l% b* @; S  Q190569292  N! {1 v! I8 @8 Z6 N/ l
    [
    0 l9 S; E+ v' q4 K+ A/ w. @    [ 10 ]," o8 G& N% d9 t6 ?; G6 S5 Z
        [ 9, 1 ],! k2 `+ \7 h5 f5 x  f( f2 R
        [ 8, 2 ],% s$ D) q  ?5 a0 Z# L8 r4 I6 y
        [ 8, 1, 1 ],% _' U/ q! i2 h$ r( ~* M* n
        [ 7, 3 ],  A/ c) b5 S' a0 z9 R; q& y1 ?
        [ 7, 2, 1 ],2 a/ |: w: L7 e2 o. b! K
        [ 7, 1, 1, 1 ],
    " x# E5 i+ B# f+ X. I; e8 i; l; G- n    [ 6, 4 ],' L; h& i1 t, u+ u, O& P
        [ 6, 3, 1 ],6 ~8 v( @* ^! p
        [ 6, 2, 2 ],
    - G8 B$ ^9 u+ q$ U: E; G    [ 6, 2, 1, 1 ],
    7 B: w: P$ c3 V* @: @7 c. B2 I    [ 6, 1, 1, 1, 1 ],4 p, b/ {: s: G3 n/ r% a. f2 b
        [ 5, 5 ],% g! f% l5 [4 S2 ]8 m
        [ 5, 4, 1 ],& j. O" q" \5 j9 w; W; P
        [ 5, 3, 2 ],
    3 o( z; c/ g& p- v    [ 5, 3, 1, 1 ],4 y- }6 ]! @" \4 d
        [ 5, 2, 2, 1 ],; z3 w' ]% v8 [4 l6 g& h' v
        [ 5, 2, 1, 1, 1 ],
    : l+ F* y  O" O    [ 5, 1, 1, 1, 1, 1 ],
    ( P3 m3 R5 r' O0 E    [ 4, 4, 2 ],
    - x( n; _* E$ }: M    [ 4, 4, 1, 1 ],
    5 [+ C- \2 N( T+ p) T: H    [ 4, 3, 3 ],
    - [: B# t9 F" C* w. r% O, S    [ 4, 3, 2, 1 ],! J$ O9 R' P0 y, N9 e7 _
        [ 4, 3, 1, 1, 1 ],
    # d- a  c; Q0 c+ I    [ 4, 2, 2, 2 ],4 O& B% e: o& I
        [ 4, 2, 2, 1, 1 ],
    : Z* j4 ^0 C2 c* }. p8 n8 h4 {    [ 4, 2, 1, 1, 1, 1 ],* N7 U& N1 R1 Y# p7 H2 R
        [ 4, 1, 1, 1, 1, 1, 1 ],6 U- Y0 R! b& B# P% V
        [ 3, 3, 3, 1 ],& r8 H; V  S8 T$ D# p2 @; g) n8 R7 G
        [ 3, 3, 2, 2 ],) Q: i& h% ^) y2 \! O; O7 P' H! U
        [ 3, 3, 2, 1, 1 ],
    - s- L" _1 A7 X! S2 l    [ 3, 3, 1, 1, 1, 1 ],5 G, t0 b! M8 d. ]( i9 A+ H
        [ 3, 2, 2, 2, 1 ],0 [( i9 S( ^5 v- j6 |
        [ 3, 2, 2, 1, 1, 1 ],
    5 G% {% q# F& P3 d2 @8 q    [ 3, 2, 1, 1, 1, 1, 1 ],
    / Y6 f: q! f* d* m8 x8 Q( Y5 q    [ 3, 1, 1, 1, 1, 1, 1, 1 ],' c* R# j; |, b' `! z3 |- H
        [ 2, 2, 2, 2, 2 ],
    ! V% i* n2 P& v: S    [ 2, 2, 2, 2, 1, 1 ],
    5 T8 ~8 r' T1 u7 e    [ 2, 2, 2, 1, 1, 1, 1 ],
    9 _, Y& y! D" w9 c! _    [ 2, 2, 1, 1, 1, 1, 1, 1 ],
    % k' @0 Y' h" _8 v  P    [ 2, 1, 1, 1, 1, 1, 1, 1, 1 ],) n# a6 `* G" l# Z1 p: o
        [ 1, 1, 1, 1, 1, 1, 1, 1, 1, 1 ]
    4 T; W" ?0 q3 x' ]& k9 @]
    回复

    使用道具 举报

    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-11 09:13 , Processed in 0.701006 second(s), 92 queries .

    回顶部