A387077 Indices of prime terms in A386482.

2, 5, 11, 17, 24, 31, 44, 52, 57, 73, 76, 115, 126, 144, 189, 207, 236, 287, 310, 320, 368, 453, 479, 652, 667, 674, 678, 684, 809, 821, 832, 837, 996, 1016, 1034, 1088, 1206, 1289, 1425, 1497, 1532, 2020, 2026, 2053, 2079, 2425, 2442, 2445, 2522, 2542, 2578, 2637

Offset

1,1

Links

Michael De Vlieger, Table of n, a(n) for n = 1..4310

Mathematica

Block[{c, j, k, m, p, r, nn},
  nn = 3000; c[_] := False; m[_] := 1; j = 2; c[1] = c[2] = True; r = 1;
  {1}~Join~Monitor[Most@ Reap[Do[
    If[PrimePowerQ[j],
      Set[{p, k, m}, {#1, #1^(#2 - 1), #1^(#2 - 1)}] & @@
        FactorInteger[j][[1]]; While[And[c[k*p], k != 0], k--];
        If[k == 0, k = m; While[c[k*p], k++]]; k *= p,
      k = j - 1; While[And[Or[c[k], CoprimeQ[j, k]], k != 1], k--];
        If[k == 1, k += j; While[Or[c[k], CoprimeQ[j, k] ], k++] ] ];
    If[PrimeQ[k], Sow[n]];
    Set[{c[k], j}, {True, k}], {n, 3, nn}] ][[-1, 1]], n] ]

Crossrefs

Cf. A386482, A387076.

Sequence in context: A217303 A053033 A136244 * A115057 A228344 A157421

Adjacent sequences: A387074 A387075 A387076 * A387078 A387079 A387080

Keyword

nonn

Author

Michael De Vlieger, Aug 15 2025

Status

approved