OFFSET
1,1
COMMENTS
Are all the terms of the form 10*k+1?
LINKS
EXAMPLE
MATHEMATICA
pspQ[n_] := PowerMod[2, n - 1, n] == 1; a[1] = 11; a[n_] := a[n] = Module[{p = NextPrime[a[n-1]], r = Product[a[i], {i, 1, n - 1}]}, While[! pspQ[p*r], p = NextPrime[p]]; p]; Array[a, 10]
PROG
(PARI) ispsp(n) = Mod(2, n)^(n-1) == 1;
lista(len) = {my(prd = 1, c = 0, k = 11); while(c < len, while(!ispsp(prd * k), k = nextprime(k+1)); prd *= k; print1(k, ", "); c++); }
(PARI) my(P=List(11), base=2); my(m = vecprod(Vec(P))); my(L = znorder(Mod(base, m))); print1(P[1], ", "); while(1, forstep(p=lift(1/Mod(m, L)), oo, L, if(p > P[#P] && isprime(p) && base % p != 0, if((m*p-1) % znorder(Mod(base, p)) == 0, print1(p, ", "); listput(P, p); L = lcm(L, znorder(Mod(base, p))); m *= p; break)))); \\ Daniel Suteu, Jun 30 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 26 2024
EXTENSIONS
a(20)-a(26) from Daniel Suteu, Jun 30 2024
STATUS
approved