OFFSET
1,5
EXAMPLE
A003626(8) = 47; 47 + 3*1 - 1 = 49 (not prime), 47 + 3^2 - 1 = 56 (not prime), 47 + 3^3 - 1 = 73, a prime, so a(8) = 3.
MATHEMATICA
s = Select[Prime[Range[400]], Mod[#, 3] == 2 &]; (* A003627 *)
g[n_] := s[[n]];
f[n_] := Select[Range[200], PrimeQ[g[n] + 3^# - 1] &, 1];
Join[{-1}, Flatten[Table[f[n], {n, 2, 200}]]]
CROSSREFS
KEYWORD
sign
AUTHOR
Clark Kimberling, Feb 24 2026
STATUS
approved
