OFFSET
1,2
LINKS
Andrew Howroyd, Table of n, a(n) for n = 1..200
Math.stackexchange.com, Normalizer of the cyclic group in S_n
FORMULA
a(n) = n! * A013939(n). - Andrew Howroyd, Jul 30 2018
PROG
(GAP) List([1..7], n->Sum(Filtered( ConjugacyClassesSubgroups( SymmetricGroup(n)), x->IsPrime( Size( Representative(x))) ), x->Size(x)*Size( Normalizer( SymmetricGroup(n), Representative(x))) )); # Andrew Howroyd, Jul 30 2018
(GAP)
a:=function(n) local total, perm, g, p, k;
total:= 0; g:= SymmetricGroup(n);
for p in Filtered([2..n], IsPrime) do for k in [1..QuoInt(n, p)] do
perm:=PermList(List([0..p*k-1], i->i - (i mod p) + ((i + 1) mod p) + 1));
total:=total + Size(Normalizer(g, perm)) * Factorial(n) / (p^k * (p-1) * Factorial(k) * Factorial(n-k*p));
od; od;
return total;
end; # Andrew Howroyd, Jul 30 2018
(PARI) a(n)={n!*sum(p=2, n, if(isprime(p), n\p))} \\ Andrew Howroyd, Jul 30 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Olivier Gérard, Apr 04 2012
EXTENSIONS
Terms a(8) and beyond from Andrew Howroyd, Jul 30 2018
STATUS
approved