OFFSET
1,1
COMMENTS
a(20) > A000960(10^7). - Donovan Johnson, Aug 23 2012
PROG
(PARI) A181488(Nmax) = { my(S=[]); for(n=1, Nmax, my(F=n); for(i=1, n-1, F=(n-i)*ceil(F/(n-i)+1)); S=setunion(S, Set(F)); isprime(F) & next; fordiv(F, d, setsearch(S, d)||next(2)); print1(F", "))}
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Oct 29 2010
EXTENSIONS
a(13)-a(16) from Donovan Johnson, Oct 30 2010
a(17)-a(18) from Donovan Johnson, Nov 11 2010
a(19) from Donovan Johnson, Aug 23 2012
STATUS
approved