OFFSET
1,2
COMMENTS
A harmonic number is a harmonic seed if it does not have a proper unitary divisor (a divisor d > 1 such that gcd(d,n/d)=1), which is harmonic. - minor correction by Antti Karttunen, Sep 19 2025
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)}
(PARI)
is_A001599(n) = (1==denominator(n*numdiv(n)/sigma(n)));
is_A035527(n) = if(!is_A001599(n), 0, fordiv(n, d, if(d>1 && d<n && 1==gcd(d, n/d) && is_A001599(d), return(0))); (1)); \\ Antti Karttunen, Sep 19 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Offset corrected by Donovan Johnson, Aug 21 2012
STATUS
approved
