OFFSET
1,2
COMMENTS
a(n) = n if n prime.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
EXAMPLE
a(18) = 4 because k^18 == 0, 1, 9, 10 (mod 18) => 4 distinct residues.
From R. J. Mathar, Aug 27 2013: (Start)
The triangle of k^n (mod n) starts in row n=1 with columns k>=0 as:
0;
0, 1;
0, 1, 2;
0, 1, 0, 1;
0, 1, 2, 3, 4;
0, 1, 4, 3, 4, 1;
0, 1, 2, 3, 4, 5, 6;
0, 1, 0, 1, 0, 1, 0, 1;
0, 1, 8, 0, 1, 8, 0, 1, 8;
0, 1, 4, 9, 6, 5, 6, 9, 4, 1;
0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10;
0, 1, 4, 9, 4, 1, 0, 1, 4, 9, 4, 1;
0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12;
Its row sums are 0, 1, 3, 2, 10, 13, 21, 4, 27, 45, 55, 38, 78, 77, 105, 8,.... (End)
MAPLE
a:= n-> nops({seq(k&^n mod n, k=0..n-1)}):
seq(a(n), n=1..100);
MATHEMATICA
Table[Length[Union[PowerMod[Range[0, n - 1], n, n]]], {n, 100}] (* T. D. Noe, Sep 21 2011 *)
PROG
(PARI) a(n)=if(isprime(n), n, #Set(vector(n, i, lift(Mod(i-1, n)^n)))) \\ Charles R Greathouse IV, Jul 31 2016
(Python)
def A195637(n): return len({pow(x, n, n) for x in range(n)}) # Chai Wah Wu, Aug 22 2023
CROSSREFS
KEYWORD
nonn,nice
AUTHOR
Michel Lagneau, Sep 21 2011
STATUS
approved