|
(MAGMA) /* For n a prime integer */ [NthPrime(n)-1 +Factorial(2*NthPrime(n)-1) div Factorial(NthPrime(n)): n in [0..10]]; // Vincenzo Librandi, Aug 01 2015
(MAGMA) /* For all n */
nbofdblecos := function(G, H, K);
CG := Classes(G); nCG := #CG; oG := #G; CH := Classes(H); nCH := #CH; oH := #H; CK := Classes(K); nCK := #CK; oK := #K;
resH := []; for mu in [1..nCG] do Gmurep := CG[mu][3]; Hmupositions := {j: j in [1..nCH] | CycleStructure(CH[j][3]) eq CycleStructure(Gmurep)};
Hmugoodpositions := {j : j in Hmupositions | IsConjugate(G, CH[j][3], Gmurep) eq true}; bide := 0; for j in Hmugoodpositions do bide := bide + CH[j][2]; end for; Append(~resH, bide); end for;
resK := []; for mu in [1..nCG] do Gmurep := CG[mu][3]; Kmupositions := {j: j in [1..nCK] | CycleStructure(CK[j][3]) eq CycleStructure(Gmurep)};
Kmugoodpositions := {j : j in Kmupositions | IsConjugate(G, CK[j][3], Gmurep) eq true}; bide := 0; for j in Kmugoodpositions do bide := bide + CK[j][2]; end for; Append(~resK, bide); end for;
ndcl := 0; tot := 0; for mu in [1..nCG] do tot := tot + resH[mu]* resK[mu]/CG[mu][2]; end for; ndcl:= tot * oG/(oH * oK); return ndcl;
end function;
OOfull := function(n); G:=Sym(2*n); genH:={}; for j in [1..(n-1)] do v := G!(1, 2*j+1)(2, 2*j+2); Include(~genH, v) ; end for;
H := PermutationGroup< 2*n |genH>;
beta:=G!Append([2..2*n], 1); Cbeta:=Centralizer(G, beta);
return nbofdblecos(G, H, Cbeta); end function;
[OOfull(n) : n in [1..10]];
// Robert Coquereaux, Aug 01 2015
| 