OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..200
PROG
(PARI) us(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d));
for(n=1, 10^7, if(us(n)==eulerphi(n)*omega(n), print(n)))
(PARI) us(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d))
{ n=0; for (m=1, 10^9, if(us(m)==eulerphi(m)*omega(m), write("b063795.txt", n++, " ", m); if (n==200, break)) ) } \\ Harry J. Smith, Aug 31 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Aug 18 2001
STATUS
approved