OFFSET
1,1
COMMENTS
Smallest k such that sopf(k) = n*q where q is a square.
LINKS
Michel Lagneau, Table of n, a(n) for n = 1..1000
EXAMPLE
a(55) = 2631 because 2631 = 3*877 and 3 + 877 = 880 = 55*16 where 16 is a square.
MAPLE
with (numtheory):
sopf:= proc(n) option remember;
add(i, i=factorset(n))
end:
a:= proc(n) local k, p;
for k from 2 while irem(sopf(k), n, 'p')>0 or
sqrt(p)<>floor(sqrt(p)) or p=1 do od; k
end:
seq (a(n), n=1..100);
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jun 10 2012
STATUS
approved