A287883 Partial sums of A037276.
1, 3, 6, 28, 33, 56, 63, 285, 318, 343, 354, 577, 590, 617, 652, 2874, 2891, 3124, 3143, 3368, 3405, 3616, 3639, 5862, 5917, 6130, 6463, 6690, 6719, 6954, 6985, 29207, 29518, 29735, 29792, 32025, 32062, 32281, 32594, 34819, 34860, 35097, 35140, 37351, 37686, 37909, 37956, 60179, 60256, 60511
Offset: 1
Links
- N. J. A. Sloane, Table of n, a(n) for n = 1..20000
Crossrefs
Cf. A037276.
Programs
-
Mathematica
co[n_, k_] := Nest[Flatten[IntegerDigits[{#, n}]] &, n, k - 1]; A037276 = Table[FromDigits[Flatten[IntegerDigits[co @@@ FactorInteger[n]]]], {n, 50}]; Table[Sum[A037276[[k]], {k, 1, n}], {n, 1, 25}] (* G. C. Greubel, Jun 23 2017 *)
-
Python
from sympy import factorint def a037276(n): f=factorint(n) l=sorted([i for i in f]) return 1 if n==1 else int("".join(str(i)*f[i] for i in l)) l=[0, 1] for n in range(2, 101): l.append(l[n - 1] + a037276(n)) print(l[1:]) # Indranil Ghosh, Jun 23 2017