OFFSET
1,1
COMMENTS
Essentially (except the first term) the same as A007684, where the product is only required to be non-deficient, i.e., possibly a perfect number. This happens for the first term, but can't happen later any more. - M. F. Hasler, Jul 30 2016
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
MATHEMATICA
a[n_] := Module[{p = Prime[n]}, k = n; r = 1 + 1/p; While[r <= 2, p = NextPrime[p]; r *= 1 + 1/p; k++]; k]; Array[a, 46] (* Amiram Eldar, Jun 29 2019 *)
PROG
(PARI) a(n) = {p = prime(n); sig = p+1; prd = p; np = n; while (sig <= 2*prd, p = nextprime(p+1); sig *= p+1; prd *= p; np++; ); return (np); } \\ Michel Marcus, Mar 10 2013
(PARI) a=1; i=0; for(n=1, 99, until(2<a*=1+1/prime(i++), ); print1(i", "); a/=1+1/prime(n)) \\ M. F. Hasler, Jul 30 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Don Reble, Nov 10 2005
STATUS
approved