OFFSET
1,1
COMMENTS
No other terms < 2700000000. - Jud McCranie, Nov 04 2001
PROG
(PARI) us(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d));
for(n=1, 10^8, if(us(n)==n+eulerphi(n), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Jason Earls, Aug 18 2001
EXTENSIONS
More terms from Jud McCranie, Nov 04 2001
Offset corrected by Mohammed Yaseen, Jul 17 2023
STATUS
approved