OFFSET
1,1
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..500
EXAMPLE
Prime factors of 2106 are 2, 3, 3, 3, 3, 13. Then rad(2106) = 2 * 3 * 13 = 78, sopfr(2106) = 2 + 3 + 3 + 3 + 3 + 13 = 27 and 78 * 27 = 2106.
MAPLE
with(numtheory): P:=proc(q) local a, k, n; for n from 1 to q do a:=ifactors(n)[2];
if 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,easy
AUTHOR
Paolo P. Lava, Jan 11 2017
STATUS
approved