OFFSET
1,1
COMMENTS
If m <> 0 (mod 24) then Sum_{(a+b*i)^m: 0<=a<m, 0<=b<m} == 0 (mod m).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
MATHEMATICA
aa[n_] := aa[n] = Mod[Sum[PowerMod[a + b *I, n, n], {a, n}, {b, n}], n]; Table[aa[24*n], {n, 1, 10}]
PROG
(PARI) a(n)=my(N=24*n, a); lift(sum(A=0, N-1, a=Mod(A, N); sum(b=0, N-1, (a+b*I)^N))) \\ Charles R Greathouse IV, Nov 05 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Oct 16 2013
STATUS
approved