%I#22 2018年9月18日11:31:36
%S 132264396899110545359642555301559844929990747755078,
%电话:89599104167264942837971624739526039296960706379993620,
%U 1000999899111111111102222222033333333330444444444405555555066666666660777777708888888099999999990
%N Osiris或Digit re-assembly numbers:等于自身数字子样本排列总和的数字。
%C该序列与A241754不同,因为该序列只使用一次置换。
%C排列的长度k相同,允许使用前导零。
%C序列中的k是:2、2、2,3、3、4、4、5、5、5,5、6、7、10、10、10,10、10,10,10、6,7、7,8、7,9、9。
%H Giovanni Resta,n表,n=1..33的a(n)(术语<10^16)
%e 10545=014+015+041+045+051+054+055+055+104+105+105+140+145+150+154+155+401+405+410+415+450+451+450+450+510+514+515+540+545+550+551+540+550+550+5501+551+5401+554。
%o(Python)
%o导入itertools
%o定义getData(a,b):
%o dig=(itertools.permutations(str(a),b))
%o表示挖掘中的d:
%o产量d
%o表示w在范围(2.6)内:
%o kk=整数(w*'1')
%o对于范围内的i(kk,10**(w+3),kk):
%o m=[]
%o get=获取数据(i,w)
%o为True时:
%o尝试:
%o n=下一个(get)
%o ee=int(“”.join((n)))
%o如果ee不以m为单位:
%o m.附录(ee)
%o除StopIteration外:
%o如果总和(m)==i,长度(m)>1:
%o m.分拣()
%o打印(总和(m),长度(m)、m、i)
%o中断
%Y参见A047726、A179239、A241754、A241899。
%K基,nonn
%O 1,1号机组
%2018年9月16日发送邮件
%E a(12)-a(26)摘自_Giovanni Resta_,2018年9月16日
|