login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A211537
Smallest number k such that the sum of prime factors of k (counted with multiplicity) equals n times a nontrivial integer power.
1
4, 15, 35, 39, 51, 95, 115, 87, 155, 111, 123, 215, 235, 159, 371, 183, 302, 335, 219, 471, 395, 415, 267, 623, 291, 303, 482, 327, 339, 791, 554, 1255, 635, 655, 411, 695, 662, 447, 698, 471, 734, 815, 835, 519, 1211, 543, 842, 1895, 579, 591, 914, 2167, 1263
OFFSET
1,1
COMMENTS
Smallest k such that sopfr(k) = n * m^q where m, q >= 2.
a(n) = A211144(n) except for n = 55, 63, 73, ... Example: a(55) = 1964 = 2^2*491 but A211144(55) = 2631 = 3*877.
LINKS
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
Sequence in context: A124150 A054556 A113693 * A213420 A209345 A323452
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jun 27 2012
STATUS
approved