OFFSET
1,6
EXAMPLE
The 2 X 2 idempotent matrices over Z/3Z are {{0, 0}, {0, 0}}, {{0, 0}, {0,1}}, {{0, 0}, {1, 1}}, {{0, 0}, {2, 1}}, {{0, 1}, {0, 1}}, {{0, 2}, {0, 1}}, {{1, 0}, {0, 0}}, {{1, 0}, {0, 1}}, {{1, 0}, {1, 0}}, {{1, 0}, {2, 0}}, {{1, 1}, {0, 0}}, {{1, 2}, {0, 0}}, {{2, 1}, {1, 2}}, {{2, 2}, {2, 2}}. Their sum is {{10, 9}, {9, 10}} == 1*{{1, 0}, {0, 1}} (mod 3) and therefore a(3) = 1.
MATHEMATICA
K[n_] := K[n] = Mod[Sum[If[ Mod[{{a, b}, {c, d}}.{{a, b}, {c, d}} - {{a, b}, {c, d}}, n] == 0{{a, b}, {c, d}}, {{a, b}, {c, d}}, 0], {a, n}, {b, n}, {c, n}, {d, n}], n]; Table[K[n][[1, 1]], {n, 1, 22}]
PROG
(PARI) a(n) = lift(sum(i=0, n-1, sum(j=0, n-1, sum(k=0, n-1, sum(l=0, n-1, m = Mod([i, j; k, l], n); if ((m^2 == m), m[1, 1])))))); \\ Michel Marcus, Apr 04 2016
CROSSREFS
KEYWORD
nonn,more
AUTHOR
José María Grau Ribas, Nov 03 2015
STATUS
approved