OFFSET
1,1
COMMENTS
These are the indices of the primes in A280864 that begin a block of length 4 that both begins and ends with a prime.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..2856
MATHEMATICA
Block[{c, j, k, m, q, r, s, t, w, nn, rad},
nn = 24000; c[_] := False; q = 0; r = s = 1; t = {};
rad[n_] := rad[n] = Times @@ FactorInteger[n][[All, 1]];
Monitor[Reap[Do[w = GCD[r, s]; m = r/w; k = m;
While[Or[c[k], ! CoprimeQ[w, k]], k += m];
c[k] = True; s = r; r = rad[k];
If[CoprimeQ[s, r],
If[Length[t] == 4,
If[AllTrue[{First[t], Last[t]}, PrimeQ], Sow[PrimePi@ First[t] ] ] ];
q++; t = {k}, AppendTo[t, k] ], {i, nn}] ][[-1, 1]], i] ] (* Michael De Vlieger, Jun 17 2026 *)
CROSSREFS
KEYWORD
nonn,new
AUTHOR
N. J. A. Sloane, Jun 17 2026
STATUS
approved
