OFFSET
1,1
PROG
(PARI) s = 3; S = Set([1, 2]); i = 3; while (1, d = divisors(s); j = 1; while (setsearch(S, d[j]), j++); n = d[j]; if (n == s, print1(s, ", ")); s += n; S = setunion(S, Set(n)); i++); \\ David Wasserman, Jan 09 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Nov 04 2005
EXTENSIONS
More terms from David Wasserman, Jan 09 2009
STATUS
approved