OFFSET
1,2
COMMENTS
Positions k where R(k) = A290732(k)/k, achieves a new minimum.
LINKS
Hugo Pfoertner, Table of n, a(n) for n = 1..182
Andreas Enge, William Hart, Fredrik Johansson, Short addition sequences for theta functions, arXiv:1608.06810 [math.NT], 2016-2018. See Table 6.
MATHEMATICA
a[n_] := Product[{p, e} = pe; If[p <= 3, p^e, (p^e - p^(e-1))/2 + (p^(e-1) - p^(Mod[e+1, 2]))/(2*(p+1)) + 1], {pe, FactorInteger[n]}];
r = 2; Reap[For[j=1, j <= 10^4, j = j+1, t = a[j]/j; If[t<r, r = t; Sow[j] ]]][[2, 1]] (* Jean-François Alcover, Oct 02 2018, after Hugo Pfoertner *)
PROG
(PARI) a290732(n)={my(f=factor(n)); prod(k=1, #f~, my([p, e]=f[k, ]); if(p<=3, p^e, (p^e-p^(e-1))/2+(p^(e-1)-p^((e+1)%2))/(2*(p+1))+1))}
my(r=2); for(j=1, 10001, my(t=a290732(j)/j); if(t<r, r=t; print1(j, ", "))) \\ Hugo Pfoertner, Aug 26 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Aug 10 2017
EXTENSIONS
a(1) corrected by Hugo Pfoertner, Aug 26 2018
STATUS
approved