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

Author

Michael De Vlieger, Aug 15 2025

Keywords

• nonn

Links

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

Crossrefs

Cf. A386482, A387076.

Programs

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] ]