24小时热门版块排行榜    

查看: 614  |  回复: 0

mcfroo

铁虫 (初入文坛)

[交流] 新概念新方法新成果,顺便解决 1,2类stirling 数计算问题

简单到一目了然, 新在于前所未见, 实用性显而易见

一, 数形
M个不同正整数(K1,K2...Km),Ki 从小到大排序,相邻数间共M-1个间隔,用A代表连续,B代表空隙(称为洞),记做M-1个字符的串: AABB...代表一个类型,记为PX
1,2,4),(1,2,8),(2,3,6)是PX=AB型,(1,2,5,6),(2,3,6,7)是ABA型,(1,3,5),(1,3,6),(2,4,6)是BB型,(2,3443865)是B型,(2,3)是PX=A型

MIN(PX)为PX型中的最小积,如MIN(AA)=1*2*3,MIN(AB)=1*2*4
IDX(PX)=2+(A的数量)+2*(B的数量),如IDX(AA)=4,IDX(AB)=5
SUM(N,PX)为1到N-1中所有PX型项积的和,如 SUM(6,AB)=1*2*4+1*2*5+2*3*5,SUM(7,AB)=SUM(6,AB)+(1*2+2*3+3*4)*6

--->SUM(N,PX)=Min(PX)*C(N,IDX(PX))

for example:
SUM(6,AA)=1*2*3+2*3*4+3*4*5=1*2*3*C(6,4)=90
SUM(6,AB)=1*2*4+1*2*5+2*3*5=1*2*4*C(6,5)=48
SUM(6,BA)=1*3*4+1*4*5+2*4*5=1*3*4*C(6,5)=72
SUM(7,BB)=1*3*5+1*3*6+1*4*6+2*4*6=1*3*5*C(7,6)=105
SUM(8,BB)=SUM(7,BB)+1*(3+4+5)*7+2*(4+5)*7+3*5*7=1*3*5*C(8,6)=420
SUM(8,BAB)=1*3*4*6+1*(3*4+4*5)*7+2*4*5*7=576=1*3*4*6*C(8,7)
SUM(9,BAB)=SUM(8,BAB)+1*(3*4+4*5+5*6)*8+2*(4*5+5*6)*8+3*5*6*8=2592=Min*C(9,7)

二, 1类stirling 数 s(N,N-K)的 直接算法
定义 F1(N,M)=s(N,N-M), s(N,M)为第1类Stirling数,  F1 就是 1到N-1中 所有M个数乘积之和
从定义即知 F1(N,M)=∑MIN(PX)*C(N,IDX(PX)),求和遍历PM=M的PX

for example:
(A+B)^3=AAA+(AAB+ABA+BAA)+(ABB+BAB+BBA)+BBB
-->F1(N,4)=1*2*3*4*C(N,5)+(1*2*3*5+1*2*4*5+1*3*4*5)*C(N,6)+(1*2*4*6+1*3*4*6+1*3*5*6)*C(N,7)+1*3*5*7*C(N,8)

三, 2类stirling 数 S2(N,N-K)的 简单递推算法
用组合计算天昏地暗, 递推算法类似杨辉3角
S2(N,N-8)=C(N,9)+501*C(N,10)+22935*C(N,11)+302995*C(N,12)+1636635*C(N,13)+4099095*C(N,14)+4729725*C(N,15)+2027025*C(N,16)
其中 302995=56980*(12-8)+ 6825*11,4729725=135135*(15-8)+270270*14,2027025=135135*15
S2(N,N-9),S2(N,N-10)...也能在几分钟算出

四,新的同余关等式

M相同且B相同的{PX},且B>0,IDX(PX)=P,P>3,则∑MIN(PX)≡0 MOD P*(P-1)
例如:
ABB,BAB,BBAM=4,Idx=7--->1*2*4*6+1*3*4*6+1*3*5*6=5*6*7≡0MOD7*6
AAAB,AABA,ABAA,BAAAM=5,Idx=7-->1*2*3*4*6+1*2*3*5*6+1*2*4*5*6+1*3*4*5*6=2*11*7*6≡0MOD7*6

详见:https://kdocs.cn/l/srECYr2HM?f=501  或附件
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 数形,Stirling数的一般计算公式.pdf
  • 2019-12-29 05:40:15, 1.34 M

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mcfroo 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见