\\
subcyclopclgp(22220, 101)
subcyclopclgp(44440, 101)
subcyclopclgp(44440, 101, 1)
subcyclopclgp(22220, 11)
subcyclopclgp(22220, 11, 1)

\\
T = polcompositum(x^2-36322, polcyclo(5), 2)
bnfinit(T).cyc
subcyclopclgp(T, 5)
subcyclopclgp(T, 5, 1+4)
default(parisizemax, "4G")
subcyclopclgp(T, 5, 1)

Q = bnfinit(y);
[,bnr1,H1] = rnfconductor(Q, x^2-36322);
[,bnr2,H2] = rnfconductor(Q, polcyclo(5));
[bnr,H] = bnrcompositum([bnr1,H1], [bnr2,H2]);
subcyclopclgp([bnr,H], 5)

bnr1 = bnrinit(Q, 8); H1 = Mat(2);
bnr2 = bnrinit(Q, [53^2, [1]]); H2 = Mat(53);
[bnr,H] = bnrcompositum([bnr1, H1], [bnr2, H2]);
subcyclopclgp([bnr,H], 107)

\\ disaster !!!
T = polcompositum(x^2-2,polsubcyclo(53^2,53), 2)
subcyclopclgp(T, 107)

p = 7860079; G = znstar(p, 1);
subcyclohminus([G, Mat(13122)], 3)

subcycloiwasawa(x^2+1501391, 3)