OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
a(55) = 1964 = 2^2*491, since the sum of the prime divisors counted with multiplicity is 491+4 = 495 = 55*3^2.
MAPLE
sopfr:= proc(n) option remember;
add(i[1]*i[2], i=ifactors(n)[2])
end:
a:= proc(n) local k, q;
for k while irem(sopfr(k), n, 'q')>0 or
igcd (map(i->i[2], ifactors(q)[2])[])<2 do od; k
end:
seq (a(n), n=1..100);
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jun 27 2012
STATUS
approved