OFFSET
1,1
REFERENCES
R. Osborn, Tables of All Primitive Roots of Odd Primes Less Than 1000, Univ. Texas Press, 1961.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
A. E. Western and J. C. P. Miller, Tables of Indices and Primitive Roots. Royal Society Mathematical Tables, Vol. 9, Cambridge Univ. Press, 1968, p. XLIV.
LINKS
Michel Marcus, Table of n, a(n) for n = 1..38 (using McGown and Sorenson).
Stephen D. Cohen, Tomás Oliveira e Silva, and Tim Trudgian, On Grosswald's conjecture on primitive roots, arXiv:1503.04519 [math.NT], 2015.
R. K. Guy and N. J. A. Sloane, Correspondence, 1988.
Kevin J. McGown and Jonathan P. Sorenson, Computation of the least primitive root, arXiv:2206.14193 [math.NT], 2022.
Tomás Oliveira e Silva, Least prime primitive root of prime numbers
A. E. Western and J. C. P. Miller, Tables of Indices and Primitive Roots, Royal Society Mathematical Tables, Vol. 9, Cambridge Univ. Press, 1968 [Annotated scans of selected pages]
MATHEMATICA
s = {2}; rm = 1; Do[p = Prime[k]; r = PrimitiveRoot[p]; If[r > rm, Print[p]; AppendTo[s, p]; rm = r], {k, 10^6}]; s (* Jean-François Alcover, Apr 05 2011 *)
DeleteDuplicates[Table[{p, PrimitiveRoot[p, 1]}, {p, Prime[Range[61100]]}], GreaterEqual[ #1[[2]], #2[[2]]]&][[All, 1]] (* The program generates the first 15 terms of the sequence. *) (* Harvey P. Dale, Aug 22 2022 *)
PROG
(Python)
from sympy import isprime, primitive_root
from itertools import count, islice
def f(n): return 0 if not isprime(n) or (r:=primitive_root(n))==None else r
def agen(r=0): yield from ((m, r:=f(m))[0] for m in count(1) if f(m) > r)
print(list(islice(agen(), 15))) # Michael S. Branicky, Feb 13 2023
CROSSREFS
KEYWORD
nonn,easy,nice
AUTHOR
EXTENSIONS
More terms from Scott Lindhurst (ScottL(AT)alumni.princeton.edu)
STATUS
approved