A363636 Indices of numbers of the form k^2+1, k >= 0, that can be written as a product of smaller numbers of that same form.
0, 3, 7, 13, 17, 18, 21, 31, 38, 43, 47, 57, 68, 73, 91, 99, 111, 117, 123, 132, 133, 157, 183, 211, 241, 242, 253, 255, 268, 273, 293, 302, 307, 313, 322, 327, 343, 381, 413, 421, 438, 443, 463, 487, 507, 515, 553, 557, 577, 593, 601, 651, 693, 697, 703, 707
Offset: 1
Keywords
Examples
0 is a term because 0^2+1 = 1 equals the empty product. 3 is a term because 3^2+1 = 10 = 2*5 = (1^2+1)*(2^2+1). 38 is a term because 38^2+1 = 1445 = 5*17*17 = (2^2+1)*(4^2+1)^2. (This is the first term that requires more than two factors.)
Links
- David Trimas, Table of n, a(n) for n = 1..2260
Crossrefs
Sequences that list those terms (or their indices or some other key) of a given sequence that are products of smaller terms of the same sequence (in other words, the nonprimitive terms of the multiplicative closure of the sequence):
this sequence (A002522),
Programs
-
Mathematica
g[lst_, p_] := Module[{t, i, j}, Union[Flatten[Table[t = lst[[i]]; t[[j]] = p*t[[j]]; Sort[t], {i, Length[lst]}, {j, Length[lst[[i]]]}], 1], Table[Sort[Append[lst[[i]], p]], {i, Length[lst]}]]]; multPartition[n_] := Module[{i, j, p, e, lst = {{}}}, {p, e} = Transpose[FactorInteger[n]]; Do[lst = g[lst, p[[i]]], {i, Length[p]}, {j, e[[i]]}]; lst]; output = Join[{0}, Flatten[Position[Table[ test = Sqrt[multPartition[n^2 + 1][[2 ;; All]] - 1]; Count[AllTrue[#, IntegerQ] & /@ test, True] > 0 , {n, 707}], True]]] (* David Trimas, Jul 23 2023 *)
Comments