OFFSET
1,1
COMMENTS
Because 257 is a Fermat prime, these numbers are all primitive roots (mod 257). Complement of A136803.
LINKS
Nathaniel Johnston, Table of n, a(n) for n = 1..128 (full sequence)
FORMULA
a(n) + a(129-n) = 257.
MAPLE
A136804 := {$(0..256)}: for n from 0 to 256 do A136804 := A136804 minus {n^2 mod 257}: od: op(sort(convert(A136804, list))); # Nathaniel Johnston, Jun 23 2011
MATHEMATICA
p=257; Select[Range[0, p-1], JacobiSymbol[ #, p]==-1&]
PROG
(PARI) for (n=0, 256, if (! issquare(Mod(n, 257)), print1(n, ", "))) \\ Michel Marcus, Mar 12 2017
(PARI) A136804=setminus(S=[0..256], Set([k^2 | k <- S]%257)); \\ M. F. Hasler, Nov 15 2017
CROSSREFS
KEYWORD
fini,full,easy,nonn
AUTHOR
T. D. Noe, Jan 22 2008
STATUS
approved