OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
360360 = 2^3 * 45045 is a term since it is divisible by a cube and A073185(360360) = 733824 > 2*360360.
MATHEMATICA
f[p_, e_] := 1 + p + If[e == 1, 0, p^2]; q[n_] := AnyTrue[(fct = FactorInteger[n])[[;; , 2]], # > 2 &] && Times @@ f @@@ fct > 2*n; Select[Range[2, 5*10^6], q]
PROG
(PARI) is(n) = {my(f = factor(n)); if(n == 1 || vecmax(f[, 2]) < 3, return(0)); prod(i=1, #f~, 1 + f[i, 1] + if(f[i, 2]==1, 0, f[i, 1]^2)) > 2*n};
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 10 2022
STATUS
approved