OFFSET
1,2
MATHEMATICA
z = 1000; s[n_] := Total[Prime[1 + Range[n]]];
t = Table[s[n], {n, 1, z}] (* A071148 *)
a[n_] := Select[Range[z], IntegerQ[t[[#]]/n] &, 1];
Flatten[Table[a[n], {n, 1, z/4}]]
PROG
(PARI) a(n) = my(k=2); while ((vecsum(primes(k))-2) % n, k++); k-1; \\ Michel Marcus, Jan 12 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jan 05 2026
STATUS
approved
