|
%I #18 Nov 21 2020 16:16:33
%S 3,7,6,21,19,14,12,21,22,27,43,33,63,28,24,66,67,30,98,57,44,129,367,
%T 42,199,63,85,84,463,54,48,93,86,201,76,66,219,111,99,120,163,60,1285,
%U 129,88,274,751,105,156,199,134,198,211,102,327,84,147,346,1765
%N Smallest number k such that the sum of divisors of k equals n times a nontrivial integer power.
%C Smallest k such that sigma(k) = n * m^q where m, q >= 2.
%H Alois P. Heinz, <a href="/A213931/b213931.txt">Table of n, a(n) for n = 1..5000</a>
%e a(34) = 201 because sigma(201) = 272 = 34*2^3.
%p with(numtheory):
%p a:= proc(n) local k, q;
%p for k while irem(sigma(k), n, 'q')>0 or
%p igcd(map(i->i[2], ifactors(q)[2])[])<2 do od; k
%p end:
%p seq(a(n), n=1..100); # _Alois P. Heinz_, Jun 26 2012
%t a[n_] := Module[{k, q, r}, For[k = 1, {q, r} = QuotientRemainder[ DivisorSigma[1, k], n]; r>0 || GCD @@ FactorInteger[q][[All, 2]]<2, k++]; k];
%t Array[a, 100] (* _Jean-François Alcover_, Nov 21 2020, after _Alois P. Heinz_ *)
%o (PARI) a(n)=my(k);while(sigma(k++)%n || !ispower(sigma(k)/n), ); k \\ _Charles R Greathouse IV_, Jun 26 2012
%Y Cf. A000203, A001597, A213401, A213386, A213420.
%K nonn
%O 1,1
%A _Michel Lagneau_, Jun 25 2012
|
