OFFSET
1,2
COMMENTS
A harmonic number is a harmonic seed if it does not have a proper unitary divisor (a divisor d such that gcd(d,n/d)=1) which is harmonic.
LINKS
Donovan Johnson, Table of n, a(n) for n = 1..138 (terms < 10^14)
G. L. Cohen and R. M. Sorli, Harmonic seeds, Fibonacci Quart., 36 No. 5 (1998) 386-390 MR1657575 (99j:11002)
PROG
(PARI) harm(n)=if(n<1, 0, 0==(sigma(n, 0)*n)%sigma(n, 1))
{a(n)=local(t); if(n<0, 0, n=a(n-1); until(if(t=1, if(harm(n), fordiv(n, d, t*=((d==n)||(d==1)||(gcd(d, n/d)!=1)||!harm(d))); t, 0)), n++); n)}
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Offset corrected by Donovan Johnson, Aug 21 2012
STATUS
approved