A135524 Row sums of A137152.

1, 3, 6, 9, 14, 21, 26, 33, 44, 57, 66, 83, 102, 125, 146, 165, 194, 225, 242, 279, 320, 363, 410, 453, 506, 565, 626, 659, 726, 797, 870, 949, 1004, 1087, 1176, 1273, 1374, 1477, 1584, 1693, 1806, 1917, 2018, 2145, 2210, 2341, 2478, 2617, 2766, 2917, 3074

Offset

1,2

Links

Table of n, a(n) for n=1..51.

Maple

A000961 := proc(n) option remember ; local a; if n = 1 then 1; else for a from procname(n-1)+1 do if nops( ifactors(a)[2] ) = 1 then RETURN(a); fi; od: fi; end: A025473 := proc(n) option remember ; if n <= 2 then n; else ifactors( A000961(n))[2] ; op(1, op(1, %)) ; fi; end: A137152 := proc(n) option remember ; local a, m, i; if n = 1 then RETURN([1]) ; else a := procname(n-1) ; m := A025473(n) ; for i from 1 to nops(a) do if gcd(op(i, a), m) <> 1 then m := m*op(i, a) ; a := subsop(i=1, a) ; fi; od; a := [op(a), m] ; fi; RETURN(a) ; end: A135524 := proc(n) add(k, k=A137152(n)) ; end: for n from 1 to 80 do printf("%d, ", A135524(n)) ; od: # R. J. Mathar, Dec 17 2008

Crossrefs

Cf. A137152.

Sequence in context: A082643 A310172 A227849 * A083853 A310173 A310174

Adjacent sequences: A135521 A135522 A135523 * A135525 A135526 A135527

Keyword

nonn

Author

Mats Granvik, Feb 19 2008

Extensions

More terms from R. J. Mathar, Dec 17 2008

Status

approved