' *********************************************************************************
'
' Somme_{i=1..n} C(i+2,3)^2 = [ 10*n^7 +105*n^6 +427*n^5 +840*n^4 +805*n^3 +315*n^2
' +18*n ]/2520
' ______________________________________________
'
' [(> La Recherche, numéro 319, avril 1999, page 97 <)]
'
' -----------------------------------------------------------------
' Somme_{i=1..n} C(i+2,3)^2 = [ 20*n^7 +210*n^6 +854*n^5 +1680*n^4 +1610*n^3 +630*n^2
' +36*n ]/7!
'
' >>__>>__>> 1^2 +4^2 +10^2 +20^2 +35^2 +56^2 +84^2 ... OEIS A086020 '
' ==================================================================
' Somme_{i=1..n} C(i+3,4)^2 = [ 70*n^9 +1260*n^8 +9420*n^7 +37800*n^6 +87654*n^5
' +117180*n^4 +83720*n^3 +25200*n^2 +576*n ]/9!
'
' >>__>>__>> 1^2 +5^2 +15^2 +35^2 +70^2 +126^2 +210^2 ... OEIS A086023 '
' ==================================================================
' Somme_{i=1..n} C(i+5,6)^2 = [ 924*n^13 +36036*n^12 +623532*n^11 +6306300*n^10
' +41333292*n^9 +183891708*n^8 +565114836*n^7 +1195854660*n^6 +1701884184*n^5
' +1545800256*n^4 +804035232*n^3 +181621440*n^2 +518400*n ]/13!
'
' >>__>>__>> 1^2 +7^2 +28^2 +84^2 +210^2 +462^2 +924^2 ... OEIS A086027 '
' ==================================================================
' Somme_{i=1..n} C(i+1,2)^3 = C(n+2,3)*{70 +420*C(n-1,1) +798*C(n-1,2) +630*C(n-1,3)
' +180*C(n-1,4)}/70
' _________________ {"triangulaires"} -----------------
'
' = [ 15*n^7 +105*n^6 +273*n^5 +315*n^4 +140*n^3 -8*n ]/840
'
' >>__>>__>> 1^3 +3^3 +6^3 +10^3 +15^3 +21^3 +28^3 ... OEIS A085438 '
' ==================================================================
' Somme_{i=1..n} C(i+2,3)^3 = C(n+3,4)*{1 +12*C(n-1,1) +46*C(n-1,2) +84*C(n-1,3) +81*C(n-1,4)
' +40*C(n-1,5) +8*C(n-1,6)} = [ 4*n^10 +60*n^9 +375*n^8 +1260*n^7 +2442*n^6 +2700*n^5
' +1535*n^4 +300*n^3 -36*n^2 ]/8640
'
' >>__>>__>> 1^3 +4^3 +10^3 +20^3 +35^3 +56^3 +84^3 ... OEIS A086021 '
' ==================================================================
' Somme_{i=1..n} C(i+3,4)^3 = [ 385*n^13 +10010*n^12 +114205*n^11 +750750*n^10
' +3138135*n^9 +8678670*n^8 +15996695*n^7 +19269250*n^6 +14370356*n^5 +5885880*n^4
' +960960*n^3 +13824*n ]/69189120
'
' >>__>>__>> 1^3 +5^3 +15^3 +35^3 +70^3 +126^3 +210^3 ... OEIS A086024 '
' ================================================================
' Somme_{i=1..n} C(i+1,2)^5 = C(n+2,3)*{77 +4620*C(n-1,1) +52437*C(n-1,2) +244629*C(n-1,3)
' +596574*C(n-1,4) +831600*C(n-1,5) +669900*C(n-1,6) +291060*C(n-1,7) +52920*C(n-1,8)}/77
'
' Somme_{i=1..n} C(i+1,2)^5 = 3 * [ C(n+2,3) ] * [ 1/3 +(80/4)*C(n-1,1) +(1135/5)*C(n-1,2)
' +(6354/6)*C(n-1,3) +(18078/7)*C(n-1,4) +(28800/8)*C(n-1,5) +(26100/9)*C(n-1,6)
' +(12600/10)*C(n-1,7) +(2520/11)*C(n-1,8) ]
' _________________ {"triangulaires"} -----------------
'
' Somme_{i=1..n} [ ( {(i+1)!}/{[2!]*[(i-1)!]} )^5 ] = [ 113400*n^11 +1247400*n^10 +5544000*n^9
' +12474000*n^8 +14196600*n^7 +6237000*n^6 -831600*n^5 +1108800*n^3 -172800*n ]/11!
'
' >>__>>__>> 1^5 +3^5 +6^5 +10^5 +15^5 +21^5 +28^5 ... OEIS A085440 '
' =========================================================
' a(10^6) = Somme_{i=1..10^6} C(i+1,2)^5 = 1^5 +3^5 +6^5 +..... +499998500001^5
' +499999500000^5 +500000500000^5
' = 2,840,940,341,047,980,110,480,153,634,716,134,539,051,226,551,254,329,004,329,000,000
' = 2.840940341047980*10^63
' _________________ {"triangulaires"} -----------------
'
' >>__>>__>> 1^5 +3^5 +6^5 +10^5 +15^5 +21^5 +28^5 ...
OEIS A085440 '
' =============================================================
' Somme_{i=1..n} C(i+1,2)^7 = [ 681080400*n^15 +10216206000*n^14 +64178730000*n^13
' +214540326000*n^12 +394840882800*n^11 +357567210000*n^10 +101026926000*n^9
' +71513442000*n^8 +223394371200*n^7 -228843014400*n^5 +116237721600*n^3
' -17679513600*n ]/15!
' _________________ {"triangulaires"} -----------------
'
' >>__>>__>> 1^7 +3^7 +6^7 +10^7 +15^7 +21^7 +28^7 ...
OEIS A085442 '
' =========================================================
' Somme_{i=1..n} C(i+1,2)^13
' = [ 49229914688306352000000*n^27 + 1329207696584271504000000*n^26
' + 15666928050406613460480000*n^25 + 103678200333573177312000000*n^24
' + 409954623058114230096000000*n^23 + 950383503057754125360000000*n^22
' + 1304589926102136139065600000*n^21 + 2280920407338609900864000000*n^20
' + 6679361953069057940688000000*n^19 + 1710690305503957425648000000*n^18
' - 37340957475434530857139200000*n^17 + 380153401223101650144000000*n^16
' + 264266613863508763537430400000*n^15 + 17279700055595529552000000*n^14
' - 1405644051017889353691820800000*n^13 + 5554912088752317123293260800000*n^11
' - 15478160145362448144003993600000*n^9 + 28228929004743466359667507200000*n^7
' - 30032019165552710131125780480000*n^5 + 15214403703917992281381273600000*n^3
' - 2312312279803787595350016000000*n ]/27!
'
' a(17) = Somme_{i=1..17} C(i+1,2)^13 = 1^13 +3^13 +6^13 +... +120^13 +136^13 +153^13
' = 31,913,502,112,695,362,871,585,445,449
'
' a(100) = Somme_{i=1..100} C(i+1,2)^13 = 1^13 +3^13 +6^13 +... +4851^13 +4950^13 +5050^13
' = 5,895,612,842,498,185,159,439,515,658,002,406,318,025,269,701,300
'
' a(500) = Somme_{i=1..500} C(i+1,2)^13 = 1^13 +3^13 +6^13 +... +124251^13 +124750^13
' +125250^13
' = 35,547,430,035,304,730,153,508,305,015,206,030,226,163,604,004,226,233,806,473,922,506,500
'
' a(1000) = Somme_{i=1..1000} C(i+1,2)^13 = 1^13 +3^13 +6^13 +... +498501^13 +499500^13
' +500500^13
' = 4,644,641,358,989,254,975,484,323,232,443,760,271,430,898,253,464,331,524,087,016,896,
' 101,432,513,000
' _________________ {"triangulaires"} -----------------
'
' >>__>>__>> 1^13 +3^13 +6^13 +10^13 +15^13 +21^13 +28^13 ...
'
' =========================================================
' Somme_{i=1..n} c(i) =
' Somme_{i=1..n} C(2*i-2,i-1)/i =
' 1/(n-1)! * [ n^(n-2) + C(n,2)*n^C(n-3,1) + { 8*C(n-4,0) +19*C(n-4,1)
' +24*C(n-4,2) +14*C(n-4,3) +3*C(n-4,4) }*n^(n-4) + { 18*C(n-5,0)
' +82*C(n-5,1) +229*C(n-5,2) +323*C(n-5,3) +244*C(n-5,4) +95*C(n-5,5)
' +15*C(n-5,6) }*n^(n-5) + { 24*C(n-6,0) +256*C(n-6,1) +1606*C(n-6,2)
' +4247*C(n-6,3) +6216*C(n-6,4) +5459*C(n-6,5) +2875*C(n-6,6)
' +840*C(n-6,7) +105*C(n-6,8) }*n^(n-6) + { 0*C(n-7,0) +588*C(n-7,1)
' +10172*C(n-7,2) +45156*C(n-7,3) +107047*C(n-7,4) +157749*C(n-7,5)
' +152103*C(n-7,6) +96544*C(n-7,7) +39025*C(n-7,8) +9135*C(n-7,9)
' +945*C(n-7,10) }*n^(n-7) + { 720*C(n-8,0) +5244*C(n-8,1)
' +78584*C(n-8,2) +470406*C(n-8,3) +1584010*C(n-8,4)
' +3412933*C(n-8,5) +5001194*C(n-8,6) +5124412*C(n-8,7)
' +3689819*C(n-8,8) +1836107*C(n-8,9) +602910*C(n-8,10)
' +117810*C(n-8,11) +10395*C(n-8,12) }*n^(n-8) + ... + C(n-3,0)*(n-1)! ]
'
' Sommes partielles des catalans
'
' >>__>>__>> 1 +1 +2 +5 +14 +42 +132 +429 +1430 +4862 ... OEIS A014137 '
' =========================================================
' ! n =
' n + C(n-2,1) + 3*C(n-3,1) + C(n-2,2) + 9*C(n-4,1) + 8*C(n-3,2)
' + 33*C(n-5,1) + 46*C(n-4,2) + 8*C(n-3,3) + 153*C(n-6,1) + 272*C(n-5,2)
' + 101*C(n-4,3) + 3*C(n-3,4) + 873*C(n-7,1) + 1796*C(n-6,2)
' + 975*C(n-5,3) + 114*C(n-4,4) + 5913*C(n-8,1) + 13424*C(n-7,2)
' + 9175*C(n-6,3) + 1935*C(n-5,4) + 65*C(n-4,5) + 46233*C(n-9,1)
' + 112928*C(n-8,2) + 90255*C(n-7,3) + 26795*C(n-6,4) + 2289*C(n-5,5)
' + 15*C(n-4,6) + 409113*C(n-10,1) + 1058864*C(n-9,2) + 949323*C(n-8,3)
' + 353507*C(n-7,4) + 49474*C(n-6,5) + 1615*C(n-5,6) + 4037913*C(n-11,1)
' + 10961744*C(n-10,2) + 10744131*C(n-9,3) + 4709551*C(n-8,4)
' + 902164*C(n-7,5) + 60080*C(n-6,6) + 630*C(n-5,7) + 43954713*C(n-12,1)
' + 124247984*C(n-11,2) + 130864827*C(n-10,3) + 64951499*C(n-9,4)
' + 15506754*C(n-8,5) + 1582455*C(n-7,6) + 48104*C(n-6,7) + 105*C(n-5,8)
' + 522956313*C(n-13,1) + 1530857264*C(n-12,2) + 1711679259*C(n-11,3)
' + 937614803*C(n-10,4) + 263894810*C(n-9,5) + 35980245*C(n-8,6)
' + 1953272*C(n-7,7) + 24535*C(n-6,8) + 6749977113*C(n-14,1)
' + 20375613104*C(n-13,2) + 23968818075*C(n-12,3) + 14233877531*C(n-11,4)
' + 4556231886*C(n-10,5) + 767814355*C(n-9,6) + 60761721*C(n-8,7)
' + 1700860*C(n-7,8) + 7245*C(n-6,9) + 93928268313*C(n-15,1)
' + 291387218864*C(n-14,2) + 358154343579*C(n-13,3)
' + 227595259787*C(n-12,4) + 80837373206*C(n-11,5) + 16034074699*C(n-10,6)
' + 1659663347*C(n-9,7) + 75834563*C(n-8,8) + 1025927*C(n-7,9)
' + 945*C(n-6,10) + 1401602636313*C(n-16,1) + 4456193638064*C(n-15,2)
' + 5692256078619*C(n-14,3) + 3832864320155*C(n-13,4)
' + 1484006742566*C(n-12,5) + 335238037899*C(n-11,6)
' + 42502032383*C(n-10,7) + 2722410516*C(n-9,8) + 70102976*C(n-8,9)
' + 408870*C(n-7,10) + 22324392524313*C(n-17,1) + 72575127295664*C(n-16,2)
' + 95928546784539*C(n-15,3) + 67923228538139*C(n-14,4)
' + 28290972852710*C(n-13,5) + 7110886730459*C(n-12,6)
' + 1057591686983*C(n-11,7) + 86965201932*C(n-10,8) + 3425012224*C(n-9,9)
' + 47489645*C(n-8,10) + 97020*C(n-7,11) + 378011820620313*C(n-18,1)
' + 1254135138098864*C(n-17,2) + 1709251167332379*C(n-16,3)
' + 1264954489761947*C(n-15,4) + 561056789094902*C(n-14,5)
' + 154233447541179*C(n-13,6) + 26113493816703*C(n-12,7)
' + 2606173335276*C(n-11,8) + 139414081273*C(n-10,9) + 3311069421*C(n-9,10)
' + 22947155*C(n-8,11) + 10395*C(n-7,12) + 6780385526348313*C(n-19,1)
' + 22919748340428466*C(n-18,2) + 32114052359081499*C(n-17,3)
' + 24718973715872411*C(n-16,4) + 11582165901190167*C(n-15,5)
' + 3437026107969291*C(n-14,6) + 648154316989263*C(n-13,7)
' + 75563588842868*C(n-12,8) + 5104268853259*C(n-11,9)
' + 176457238662*C(n-10,10) + 2443281775*C(n-9,11) + 7498260*C(n-8,12)
' + 128425485935180313*C(n-20,1) + .....
'
' Sommes partielles des factorielles
'
' >>__>>__>> 1 +1 +2 +6 +24 +120 +720 +5040 +40320 +362880 ... OEIS A003422 '
' =========================================================
' Premier(n_impair) = A000040(n_impair) =
' A115298(n-1) + Somme_{i=1..(n-3)/2} [ A115297(i,n-2) - A115297((n-3)/2,n-2)
' - 2*A115297(i,n-1) + A115297(1,n-1) ] + A000040(n-2) - A000040(n-1) + 2
'
' Premier(n_pair) = A000040(n_pair) =
' A115298(n-1) + Somme_{i=1..(n-4)/2} [ A115297(i,n-2) - 2*A115297(i,n-1) ]
' + A000040(n-2) - A000040(n-1) + 2
'
' >>__>>__>> 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, ... OEIS A000040 '
' =========================================================
' 1, 1, 2, 7, 8, 3, 6, 38, 93, 111, 65, 15, 24, 226, 874, 1821, 2224, 1600, 630,
' 105, 120, 1524, 8200, 24860, 47185, 58465, 47474, 24430, 7245, 945,
' 720, 11628, 81080, 326712, 852690, 1522375, 1905168, 1676325,
' 1018682, 407925, 97020, 10395, 5040, 99504, 859068, 4356044, 14604590, 34397790,
' 58808449, 74133703, 69077049, 47080775, 22850135, 7487865, 1486485, 135135,
' 40320, 945936, 9794808, 60241948, 248388056, 731834110, 1598901626,
' 2646575953, 3354909248, 3263579776, 2420334620, 1345662010, 543543000,
' 150720570, 25675650, 2027025, 362880, 9902880, 120120696, 872663304,
' 4292337076, 15266260344, 40842369036, 84242791416, 135989069049,
' 173146169241, 174179263764, 137842990956, 84865264770, 39833463670,
' 13776803040, 3310807500, 493918425, 34459425, 3628800, 113286240,
' 1580814432, 13296262728, 76281141320, 319203963200, 1015089654600,
' 2519208133344, 4964854771986, 7855656271695, 10037861043000, 10373047112175,
' 8643848682612, 5763278078558, 3032428599200, 1231344264150, 372324642690,
' 78946542675, 10475665200, 654729075, 39916800, 1406609280, 22257138816,
' 213361382256, 1403169369360, 6775648692560, 25055902129720, 72957415507592,
' 170471577504822, 323766787021590, 503968974045405, 645868669875495,
' 682319501142201, 592990687600311, 421592441600330, 242793694613470,
' 111519426298280, 39897964606500, 10716125364945, 2032933777875,
' 242904486825, 13749310575, 479001600, 18844755840, 334185525504,
' 3605269060368, 26806966110144, 147122560810720, 622040823172560,
' 2085073890864888, 5652686401346776, 12569866592166378, 23152570388738010,
' 35551359181917645, 45674567831475504, 49149200941792848, 44233504001611834,
' 33157172615782215, 20550945564949800, 10414427391738370, 4243935735429390,
' 1356888542884875, 327808007326800, 56262178872900, 6113860102350,
' 316234143225, 6227020800, 271011605760, 5334101735040, 64090364217984,
' 532765816242192, 3282792660428112, 15657681457207504, 59526437968631920,
' 184149883848744528, 470536008487194364, 1003841580495932464,
' 1801841116485683988, 2735095330182696585, 3521355039279088089,
' 3848865796297683342, 3567679519438103602, 2795786253119533219,
' 1841999197773089355, 1011785970162587620, 457735276322243920,
' 167611471723270575, 48437013885234975, 10632705460527150, 1666132289271450,
' 166022925193125, 7905853580625,
'
' >>__>>__>> 1, 1, 2, 7, 8, 3, 6, 38, 93, 111, 65, 15, 24, ... OEIS A094216 '
' =========================================================
' a(n,k) =
' s(k,k-n+1) =
' (-1)^(n-1) * S1(k,k-n+1) =
' (-1)^(n-1) * ( k!/{(k-n)!*(n-1)!*2^(n-1)} ) * [ {1}*k^(n-2)
' - {C(n-3,0)/3 +2*C(n-3,1)/3 +C(n-3,2)/3}*k^(n-3)
' + {C(n-4,1)/3 +C(n-4,2) +C(n-4,3) +C(n-4,4)/3}*k^(n-4)
' - {-2*C(n-5,0)/15 -8*C(n-5,1)/15 -11*C(n-5,2)/45 +76*C(n-5,3)/45 +16*C(n-5,4)/5
' +20*C(n-5,5)/9 +5*C(n-5,6)/9}*k^(n-5)
' + {-2*C(n-6,1)/3 -10*C(n-6,2)/3 -145*C(n-6,3)/27 -5*C(n-6,4)/27 +260*C(n-6,5)/27
' +332*C(n-6,6)/27 +175*C(n-6,7)/27 +35*C(n-6,8)/27}*k^(n-6)
' - {16*C(n-7,0)/63 +32*C(n-7,1)/21 +44*C(n-7,2)/63 -856*C(n-7,3)/63 -2455*C(n-7,4)/63
' -2354*C(n-7,5)/63 +751*C(n-7,6)/63 +532*C(n-7,7)/9 +497*C(n-7,8)/9 +70*C(n-7,9)/3
' +35*C(n-7,10)/9}*k^(n-7)
' + {404*C(n-8,1)/135 +2828*C(n-8,2)/135 +2128*C(n-8,3)/45 -112*C(n-8,4)/27
' -623*C(n-8,5)/3 -51541*C(n-8,6)/135 -...}*k^(n-8)
' - {-16*C(n-9,0)/15 -128*C(n-9,1)/15 -124*C(n-9,2)/45 +7072*C(n-9,3)/45
' +5390*C(n-9,4)/9 +34832*C(n-9,5)/45 +...}*k^(n-9)
' + {-208*C(n-10,1)/9 -208*C(n-10,2) -16172*C(n-10,3)/27 +156*C(n-10,4) +...}*k^(n-10)
' - {256*C(n-11,0)/33 +2560*C(n-11,1)/33 +488*C(n-11,2)/165 -83504*C(n-11,3)/33 -...}*k^(n-11)
' + {255968*C(n-12,1)/945 +2815648*C(n-12,2)/945 +...}*k^(n-12)
' - {-353792*C(n-13,0)/4095 -1415168*C(n-13,1)/1365 -...}*k^(n-13) + ..... ]
'
' a(n,k) =
' s(k,k-n+1) =
' (-1)^(n-1) * S1(k,k-n+1) =
' (-1)^(n-1) * ( k!/{(k-n)!*2^(n-1)} ) * [ { 1/(n-1)! }*k^(n-2)
' - { (1/6)*(1/(n-3)!) }*k^(n-3)
' + { (1/72)*(1/(n-5)!) }*k^(n-4)
' - { (1/6480)*(5/(n-7)!-36/(n-5)!) }*k^(n-5)
' + { (1/155520)*(5/(n-9)!-144/(n-7)!) }*k^(n-6)
' - { (1/6531840)*(7/(n-11)!-504/(n-9)!+2304/(n-7)!) }*k^(n-7)
' + { (1/1175731200)*(35/(n-13)!-5040/(n-11)!+87264/(n-9)!) }*k^(n-8)
' - { (1/7054387200)*(5/(n-15)!-1260/(n-13)!+52704/(n-11)!-186624/(n-9)!) }*k^(n-9)
' + { (1/338610585600)*(5/(n-17)!-2016/(n-15)!+164736/(n-13)!-2156544/(n-11)!) }*k^(n-10)
' - { (1/1005673439232000)*(275/(n-19)!-166320/(n-17)!+23379840/(n-15)!
' -726713856/(n-13)!+2149908480/(n-11)!) }*k^(n-11)
' + { (1/84476568895488000)*(385/(n-21)!-332640/(n-19)!+73846080/(n-17)!
' -4425974784/(n-15)!+47769772032/(n-13)!) }*k^(n-12) - ..... ]
'
' a(n,k) =
' s(k,k-n+1) =
' (-1)^(n-1) * S1(k,k-n+1) =
' (-1)^(n-1) * ( k!/{(k-n)!*2^(n-1)} )
' * [ Somme_{i=1..n-1} (-1)^(i-1)
' * { ( 1 / { 6^(i-1)*(i-1)!*(n-2*i+1)! } )
' - ( { (i-3) } / { 5*6^(i-2)*(i-3)!*(n-2*i+3)! } )
' + ( { (i-5)*(21*i-46) } / { 1050*6^(i-3)*(i-5)!*(n-2*i+5)! } )
' - ( { (i-7)*(i-4)*(7*i-11) } / { 5250*6^(i-4)*(i-7)!*(n-2*i+7)! } )
' + ( { (i-9)*6*(5040+2959*C(i-10,1)+...) } / { 1819125*6^(i-5)*(i-9)!*(n-2*i+9)! } )
' - ..... } * k^(n-i-1) ]
'
' Triangle des nombres de Stirling de 1 ère espèce
' ( classés par rangées )
'
' >>__>>__>> 1, 1, -1, 1, -3, 2, 1, -6, 11, -6, 1, -10, 35, -50, 24, ... OEIS A008276 '
' s(n,k) =
' (-1)^(n-k) * S1(n,k) =
' (-1)^(n-k) * ( n!/{(k-1)!*(n-k)!*2^(n-k)} ) * [ {1}*n^(n-k-1)
' - { (n-k)!/(6*(n-k-2)!) }*n^(n-k-2)
' + { (n-k)!/(72*(n-k-4)!) }*n^(n-k-3)
' - { (1/6480)*(5*(n-k)!/(n-k-6)!-36*(n-k)!/(n-k-4)!) }*n^(n-k-4)
' + { (1/155520)*(5*(n-k)!/(n-k-8)!-144*(n-k)!/(n-k-6)!) }*n^(n-k-5)
' - { (1/6531840)*(7*(n-k)!/(n-k-10)!-504*(n-k)!/(n-k-8)!+2304*(n-k)!/(n-k-6)!) }*n^(n-k-6)
' + { (1/1175731200)*(35*(n-k)!/(n-k-12)!-5040*(n-k)!/(n-k-10)!+87264*(n-k)!/(n-k-8)!) }*n^(n-k-7)
' - {-16*C(n-k-8,0)/15 -128*C(n-k-8,1)/15 -124*C(n-k-8,2)/45 +7072*C(n-k-8,3)/45
' +5390*C(n-k-8,4)/9 +34832*C(n-k-8,5)/45 +...}*n^(n-k-8)
' + {-208*C(n-k-9,1)/9 -208*C(n-k-9,2) -16172*C(n-k-9,3)/27
' +156*C(n-k-9,4) +...}*n^(n-k-9)
' - {256*C(n-k-10,0)/33 +2560*C(n-k-10,1)/33 +488*C(n-k-10,2)/165
' -83504*C(n-k-10,3)/33 -...}*n^(n-k-10)
' + {255968*C(n-k-11,1)/945 +2815648*C(n-k-11,2)/945 +...}*n^(n-k-11)
' - {-353792*C(n-k-12,0)/4095 -1415168*C(n-k-12,1)/1365 -...}*n^(n-k-12) + ..... ]
'
' s(n,k) =
' (-1)^(n-k) * S1(n,k) =
' (-1)^(n-k) * ( n!/{(k-1)!*2^(n-k)} ) * [ { 1/(n-k)! }*n^(n-k-1)
' - { (1/6)*(1/(n-k-2)!) }*n^(n-k-2)
' + { (1/72)*(1/(n-k-4)!) }*n^(n-k-3)
' - { (1/6480)*(5/(n-k-6)!-36/(n-k-4)!) }*n^(n-k-4)
' + { (1/155520)*(5/(n-k-8)!-144/(n-k-6)!) }*n^(n-k-5)
' - { (1/6531840)*(7/(n-k-10)!-504/(n-k-8)!+2304/(n-k-6)!) }*n^(n-k-6)
' + { (1/1175731200)*(35/(n-k-12)!-5040/(n-k-10)!+87264/(n-k-8)!) }*n^(n-k-7)
' - { (1/7054387200)*(5/(n-k-14)!-1260/(n-k-12)!+52704/(n-k-10)!-186624/(n-k-8)!) }*n^(n-k-8)
' + { (1/338610585600)*(5/(n-k-16)!-2016/(n-k-14)!+164736/(n-k-12)!-2156544/(n-k-10)!) }*n^(n-k-9)
' - { (1/1005673439232000)*(275/(n-k-18)!-166320/(n-k-16)!+23379840/(n-k-14)!
' -726713856/(n-k-12)!+2149908480/(n-k-10)!) }*n^(n-k-10)
' + { (1/84476568895488000)*(385/(n-k-20)!-332640/(n-k-18)!+73846080/(n-k-16)!
' -4425974784/(n-k-14)!+47769772032/(n-k-12)!) }*n^(n-k-11) - ..... ]
'
' s(n,k) =
' (-1)^(n-k) * S1(n,k) =
' (-1)^(n-k) * ( n!/{(k-1)!*2^(n-k)} )
' * [ Somme_{i=1..n-k} (-1)^(i-1)
' * { ( 1 / { 6^(i-1)*(i-1)!*(n-k-2*i+2)! } )
' - ( { (i-3) } / { 5*6^(i-2)*(i-3)!*(n-k-2*i+4)! } )
' + ( { (i-5)*(21*i-46) } / { 1050*6^(i-3)*(i-5)!*(n-k-2*i+6)! } )
' - ( { (i-7)*(i-4)*(7*i-11) } / { 5250*6^(i-4)*(i-7)!*(n-k-2*i+8)! } )
' + ( { (i-9)*6*(5040+2959*C(i-10,1)+...) } / { 1819125*6^(i-5)*(i-9)!*(n-k-2*i+10)! } )
' - ..... } * n^(n-k-i) ]
'
' ((( Cette "formule explicite" {d'échelon ou de rang 1 selon
Louis Comtet (2.3); et aussi pas
' seulement récurrente comme d'aucunes} couvre plus de 72 % des coefficients présents dans
' le tableau de Francis L. Miksa (1956) tiré du Formulaire ou Sommier d'Abramowitz et Stegun
' [ Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables;
' Table 24.3; pp. 833-834 ] ! )))
'
' {{{ «« . . . mais on ne connaît pas d'expression donnant s(n,k) en fonction de n et k . . . »»
' [ " Stirling " - Dictionnaire des Mathématiques par Alain Bouvier, Michel George et François Le Lionnais;
' Presses universitaires de France (P.U.F.); page 711 (1983, 2 ème édition), pages 800-801
' (1993, 4 ème édition) ] }}}
'
' =========================================================
' a(n,k) =
' S(k,k-n+1) =
' S2(k,k-n+1) =
' ( k!/{(k-n)!*(n-1)!*2^(n-1)} ) * [ {1}*k^(n-2)
' - {5*C(n-3,0)/3 +10*C(n-3,1)/3 +5*C(n-3,2)/3}*k^(n-3)
' + {6*C(n-4,0) +79*C(n-4,1)/3 +43*C(n-4,2) +31*C(n-4,3) +25*C(n-4,4)/3}*k^(n-4)
' - {502*C(n-5,0)/15 +3508*C(n-5,1)/15 +30161*C(n-5,2)/45 +45524*C(n-5,3)/45
' +12752*C(n-5,4)/15 +3400*C(n-5,5)/9 +625*C(n-5,6)/9}*k^(n-5)
' + {760*C(n-6,0)/3 +7390*C(n-6,1)/3 +30800*C(n-6,2)/3 +649625*C(n-6,3)/27
' +939175*C(n-6,4)/27 +...}*k^(n-6)
' - {152696*C(n-7,0)/63 +213060*C(n-7,1)/7 +10586924*C(n-7,2)/63
' +34009180*C(n-7,3)/63 +...}*k^(n-7)
' + {84112*C(n-8,0)/3 +58548164*C(n-8,1)/135 +344857628*C(n-8,2)/135 +...}*k^(n-8)
' - {17120272*C(n-9,0)/45 +315705056*C(n-9,1)/45 +...}*k^(n-9)
' + {29621376*C(n-10,0)/5 +...}*k^(n-10) - ..... ]
'
' Triangle des nombres de
Stirling de 2 ème espèce
' ( classés par rangées )
'
' >>__>>__>> 1, 1, 1, 1, 3, 1, 1, 6, 7, 1, 1, 10, 25, 15, 1, ... OEIS A008278 '
' S(n,k) =
' S2(n,k) =
' ( n!/{(k-1)!*(n-k)!*2^(n-k)} ) * [ {1}*n^(n-k-1)
' - {5*C(n-k-2,0)/3 +10*C(n-k-2,1)/3 +5*C(n-k-2,2)/3}*n^(n-k-2)
' + {6*C(n-k-3,0) +79*C(n-k-3,1)/3 +43*C(n-k-3,2) +31*C(n-k-3,3)
' +25*C(n-k-3,4)/3}*n^(n-k-3)
' - {502*C(n-k-4,0)/15 +3508*C(n-k-4,1)/15 +30161*C(n-k-4,2)/45
' +45524*C(n-k-4,3)/45 +12752*C(n-k-4,4)/15 +3400*C(n-k-4,5)/9
' +625*C(n-k-4,6)/9}*n^(n-k-4)
' + {760*C(n-k-5,0)/3 +7390*C(n-k-5,1)/3 +30800*C(n-k-5,2)/3
' +649625*C(n-k-5,3)/27 +939175*C(n-k-5,4)/27 +...}*n^(n-k-5)
' - {152696*C(n-k-6,0)/63 +213060*C(n-k-6,1)/7 +10586924*C(n-k-6,2)/63
' +34009180*C(n-k-6,3)/63 +...}*n^(n-k-6)
' + {84112*C(n-k-7,0)/3 +58548164*C(n-k-7,1)/135
' +344857628*C(n-k-7,2)/135 +...}*n^(n-k-7)
' - {17120272*C(n-k-8,0)/45 +315705056*C(n-k-8,1)/45 +...}*n^(n-k-8)
' + {29621376*C(n-k-9,0)/5 +...}*n^(n-k-9) - ..... ]
'
' S(n,k) =
' S2(n,k) =
' 1/(k-1)! * Somme_{i=1..[2*k+1+(-1)^(k-1)]/4}
' [ C(k-1,2*i-2) * (k-2*i+2)^(n-1) - C(k-1,2*i-1) * (k-2*i+1)^(n-1) ]
'
' ***************************************************************************