OFFSET
1,1
COMMENTS
Values of k: {2, 2, 3, 2, 2, 3, 3, 2, 3, 3, 5, 3, 3, 6, 3, 3, 3, 3, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 3, 4, 4, 4, 4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, ...}. - Michael De Vlieger, Mar 17 2017
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..150
EXAMPLE
sigma(382) = 576 and d(382) * d(382^2) * d(382^3) = 4 * 9 * 16 = 576;
sigma(9598) = 14400 and d(9598) * d(9598^2) * d(9598^3) * d(9598^4) = 4 * 9 * 16 * 25 = 14400.
MAPLE
with(numtheory): P:=proc(q) local a, k, n; for n from 1 to q do a:=1; k:=0; while a<sigma(n) do k:=k+1; a:=a*tau(n^k); if sigma(n)=a then print(n); break; fi; od; od; end: P(10^5);
MATHEMATICA
Select[Range[2, 40000], Module[{k = 1, d = DivisorSigma[1, #], b}, While[Set[b, Product[DivisorSigma[0, #^j], {j, k}]] < d, k++]; If[b == d, True, False]] &] (* Michael De Vlieger, Mar 17 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Mar 16 2017
STATUS
approved