OFFSET
1,1
PROG
(PARI) isok(newp, v, n) = {for (k=1, n, if (!issquare(Mod(newp, v[k])), return (0)); ); return (1); }
lista(nn) = {my(v=vector(nn), lastp=2); v[1] = lastp; for (n=2, nn, my(newp = nextprime(lastp+1)); while (! isok(newp, v, n-1), newp = nextprime(newp+1)); v[n] = newp; lastp = newp; ); v; } \\ Michel Marcus, Sep 24 2020
CROSSREFS
KEYWORD
nonn,more
AUTHOR
EXTENSIONS
Name edited by Michel Marcus, Sep 24 2020
STATUS
approved