with(LinearAlgebra):
n := 10;
Cn := Matrix(n,n);
  for i to n do 
    for j to n do 
        Cn[i,j] := x[(i+j-2) mod n + 1 ];
    od;
  od;
Dn := CodeTools[Usage](Determinant(Cn,method=minor)):
printf("n=%d  #Dn=%d\n",n,nops(Dn));
Fn := CodeTools[Usage](factors(Dn)):
printf("#f=%a\n",map(proc(f) nops(f[1]) end,Fn[2]));