|
%I #10 Nov 27 2017 00:38:05
%S 1,5,25,35,175,385,1225,1925,9625,13475,48125,55825,148225,279125,
%T 390775,1395625,1730575,4298525,8652875,12114025,43264375,60570125,
%U 133254275,302850625,423990875,642043325,2119954375
%N Indices of records in A192013: Sum_{d|n} Kronecker(-6, d).
%H Charles R Greathouse IV, <a href="/A279541/b279541.txt">Table of n, a(n) for n = 1..400</a>
%o (PARI) smooth(P:vec, lim)=my(v=List([1]), nxt=vector(#P, i, 1), indx, t); while(1, t=vecmin(vector(#P, i, v[nxt[i]]*P[i]), &indx); if(t>lim, break); if(t>v[#v], listput(v, t)); nxt[indx]++); Vec(v);
%o ok(p)=p=p%24; p==1 || p==5 || p==7 || p==11;
%o A192013(n)=sumdivmult(n, d, kronecker(-6, d));
%o list(lim)=my(u=List([1]),v=[5],r,t); forprime(p=7,, if(ok(p), if(v[#v]*p>lim, break); v=concat(v, v[#v]*p))); v=smooth(v,lim); for(i=2,#v, t=A192013(v[i]); if(t>r, r=t; listput(u,v[i]))); Vec(u)
%Y Cf. A192013, A279542, A109017.
%K nonn
%O 1,2
%A _Charles R Greathouse IV_, Dec 14 2016
|
