A135525 Row sums of terms > 1 in A137152.

0, 2, 5, 7, 12, 19, 23, 29, 40, 53, 61, 78, 97, 120, 140, 158, 187, 218, 234, 271, 312, 355, 402, 444, 497, 556, 617, 649, 716, 787, 860, 939, 993, 1076, 1165, 1262, 1363, 1466, 1573, 1682, 1795, 1905, 2005, 2132, 2196, 2327, 2464, 2603, 2752, 2903, 3060, 3223

Offset

1,2

Links

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

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: A135525 := proc(n) local a, k, i ; a := 0 ; k :=A137152(n) ; for i in k do if i <> 1 then a := a+i; fi; od; a ; end: for n from 1 to 80 do printf("%d, ", A135525(n)) ; od: # R. J. Mathar, Dec 17 2008

Crossrefs

Cf. A137152.

Sequence in context: A173088 A005895 A238661 * A319142 A350497 A117538

Adjacent sequences: A135522 A135523 A135524 * A135526 A135527 A135528

Keyword

nonn

Author

Mats Granvik, Feb 19 2008

Extensions

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

Status

approved