|
%I #6 Sep 18 2015 04:04:10
%S 1,3,6,9,14,21,26,33,44,57,66,83,102,125,146,165,194,225,242,279,320,
%T 363,410,453,506,565,626,659,726,797,870,949,1004,1087,1176,1273,1374,
%U 1477,1584,1693,1806,1917,2018,2145,2210,2341,2478,2617,2766,2917,3074
%N Row sums of A137152.
%p 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
%Y Cf. A137152.
%K nonn
%O 1,2
%A _Mats Granvik_, Feb 19 2008
%E More terms from _R. J. Mathar_, Dec 17 2008
|
