OFFSET
1,1
COMMENTS
If p is the largest prime divisor of a term k, then p^2 divides k. - Max Alekseyev, Feb 03 2024
EXAMPLE
Prime factors of 288 are 2, 2, 2, 2, 2, 3, 3. Then phi(288) = 96, rad(288) = 2 * 3 = 6, sopfr(288) = 2 + 2 + 2 + 2 + 2 + 3 + 3 = 16 and 6 * 16 = 96.
MAPLE
with(numtheory): P:=proc(q) local a, k, n; for n from 1 to q do a:=ifactors(n)[2];
if phi(n)=mul(a[k][1], k=1..nops(a))*add(a[k][1]*a[k][2], k=1..nops(a)) then print(n);
fi; od; end: P(10^9);
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Jan 11 2017
EXTENSIONS
a(9)-a(16) from Max Alekseyev, Feb 03 2024
STATUS
approved