OFFSET
1,2
PROG
(GAP) Concatenation([1], List([2..7], n->Sum( Filtered( ConjugacyClassesSubgroups( GL(2, Integers mod n)), x->IsCyclic( Representative(x))), Size)));
(PARI)
MatOrder(M)={my(id=matid(#M), k=1, N=M); while(N<>id, k++; N=N*M); k}
a(n)={sum(a=0, n-1, sum(b=0, n-1, sum(c=0, n-1, sum(d=0, n-1, my(M=Mod([a, b; c, d], n)); if(gcd(lift(matdet(M)), n)==1, 1/eulerphi(MatOrder(M)))))))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Howroyd, Jul 06 2018
STATUS
approved