%I #10 Nov 05 2020 06:40:40
%S 33,3,155,3,77,5,51,3,77,3,35,5,50,3,187,6,21,5,39,3,145,33,39,5,69,
%T 39,91,3,33,7,15,12,221,3,28,7,21,3,55,3,33,5,91,66,209,69,35,5,50,3,
%U 115,39,141,5,51,6,145,85,15,7,21,93,95,3,57,5,51,3,65,15,35,7,69,55,287,6
%N Least k such that k, k+n and k+2n have the same prime signature.
%C Third row of A113456.
%H Alois P. Heinz, <a href="/A113458/b113458.txt">Table of n, a(n) for n = 1..20000</a>
%e a(4) = 3 because 3, 7 and 11 have the same prime signature.
%p s:= n-> sort(map(i-> i[2], ifactors(n)[2])):
%p a:= proc(n) option remember; local k; for k
%p while s(k)<>s(k+n) or s(k)<>s(k+2*n) do od; k
%p end:
%p seq(a(n), n=1..100); # _Alois P. Heinz_, Feb 28 2018
%t s[n_] := FactorInteger[n][[All, 2]] // Sort;
%t a[n_] := Module[{k}, For[k = 2, True, k++, If[s[k] == s[k+n] == s[k+2n], Return[k]]]];
%t Array[a, 100] (* _Jean-François Alcover_, Nov 05 2020 *)
%Y Cf. A052214, A113456, A113467.
%K easy,nonn
%O 1,1
%A _David Wasserman_, Jan 08 2006
|