|
Search: id:A135524
|
|
| |
|
| 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
(list; graph; listen)
|
|
|
OFFSET
|
1,2
|
|
|
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: [From R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Dec 17 2008]
|
|
CROSSREFS
|
Cf. A137152.
Sequence in context: A134031 A130473 A082643 this_sequence A083853 A058597 A146562
Adjacent sequences: A135521 A135522 A135523 this_sequence A135525 A135526 A135527
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Mats Granvik (mgranvik(AT)abo.fi), Feb 19 2008
|
|
EXTENSIONS
|
More terms from R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Dec 17 2008
|
|
|
Search completed in 0.002 seconds
|
