%I #15 Feb 12 2021 04:06:26
%S 4,4,121,4,4,4,4,4,121,4,4,4,4,4,121,4,4,4,4,4,121,4,4,169,4,4,4,121,
%T 4,4,4,4,4,121,4,4,4,4,4,121,4,4,4,4,4,121,4,4,4,4,4,121,4,4,4,4,4,
%U 121,4,4,4,4,4,121,4,4,4,4,4,121,4,169,4,4,4,4,121
%N The smallest square referenced in A249025 (Numbers k such that 3^k - 1 is not squarefree).
%H Amiram Eldar, <a href="/A283453/b283453.txt">Table of n, a(n) for n = 1..407</a>
%F a(n) = A283454(n)^2. - _Amiram Eldar_, Feb 12 2021
%e A249025(3)=5, 3^5-1 = 242 = 2*11*11. 242 is not squarefree the square being 11*11 = 121.
%t psq[n_] := If[(f = Select[FactorInteger[n], Last[#] > 1 &]) == {}, 1, f[[1, 1]]^2]; psq /@ Select[3^Range[100] - 1, !SquareFreeQ[#] &] (* _Amiram Eldar_, Feb 12 2021 *)
%o (PARI) lista(nn) = {for (n=1, nn, if (!issquarefree(k = 3^n-1), f = factor(k/core(k)); vsq = select(x->((x%2) == 0), f[,2], 1); print1(f[vsq[1], 1]^2, ", ");););} \\ _Michel Marcus_, Mar 11 2017
%Y Cf. A249025, A283454.
%K nonn
%O 1,1
%A _Robert Price_, Mar 07 2017
%E More terms from _Michel Marcus_, Mar 11 2017
|