|
%I
%S 1,3,4,7,13,25,33,48,61,79,91,119,133,157,181,212,230,269,289,331,363,
%T 399,423,483,514,556,596,652,682,754,786,849,897,951,999,1090,1128,
%U 1188,1244,1334,1376,1472,1516,1600,1678,1750,1798,1922,1979,2072,2144,2242
%N a(n)=sigma(n) if sigma(n)>=a(n-1), else a(n)=a(n-1)+sigma(n).
%C Does a(n)=sigma(n) only for n=1,2,3,4?
%o (PARI) { v=vector(60); v[1]=sigma(1); for (i=2,60, if (sigma(i)<v[i-1],v[i]=v[i-1]+sigma(i),v[i]=sigma(i))); v }
%Y Cf. A000203.
%K nonn
%O 1,2
%A _Jon Perry_, Mar 22 2004
|