OFFSET
1,2
COMMENTS
The sequence arose while solving puzzle 329 from Carlos Rivera's Prime Puzzles & Problems Connection site.
LINKS
EXAMPLE
a(2)=3 because the second prime being 3, then 3^2 * 5 * 7 * 11 = 3465 and sigma(3465) - 2*3465 = 558, a positive number (i.e., 3465 is abundant), but 3^2 * 5 * 7 = 315 and sigma(315) - 2*315 = -6, a nonpositive number (i.e., 315 is not abundant).
MATHEMATICA
abQ[n_] := DivisorSigma[1, n] > 2n; f[0] = 0; f[n_] := f[n] = Block[{k = f[n - 1]}, p = Fold[Times, Prime[n], Prime[ Range[n, n + k]]]; While[ !abQ[p], k++; p = p*Prime[n + k]]; k]; Table[ f[n], {n, 48}] (* Robert G. Wilson v *)
PROG
(PARI) forprime(p=2, 100, k=0; while(k++, if(sigma(n=p^2*prod(j=1, k, prime(j+primepi(p))))-n>n, print(k); break)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Igor Schein, Sep 13 2005
EXTENSIONS
Edited and extended by Robert G. Wilson v, Sep 15 2005
STATUS
approved