login
A397165
Apply PrimePi to A397164.
2
33, 120, 171, 201, 203, 210, 213, 230, 238, 264, 286, 288, 294, 322, 354, 373, 395, 479, 486, 506, 542, 553, 563, 583, 609, 628, 646, 659, 688, 700, 719, 724, 730, 744, 749, 770, 814, 862, 883, 963, 975, 1009, 1020, 1056, 1103, 1115, 1119, 1182, 1247, 1260, 1273, 1324, 1358
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
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 *)
KEYWORD
nonn,new
AUTHOR
N. J. A. Sloane, Jun 17 2026
STATUS
approved