gap> Q:=QuadraticNumberField(-1);;
gap> OQ:=RingOfIntegers(Q);;
gap> I:=QuadraticIdeal(OQ,41+56*Sqrt(-1));
ideal of norm 4817 in O(GaussianRationals)
gap> G:=HAP_CongruenceSubgroupGamma0(I);;
gap> AbelianInvariants(G);
[ 2, 2, 4, 5, 7, 16, 29, 43, 157, 179, 1877, 7741, 22037, 292306033, 
  4078793513671 ]

gap> II:=QuadraticIdeal(OQ,47+61*Sqrt(-1));
ideal of norm 5930 in O(GaussianRationals)
gap> GG:=HAP_CongruenceSubgroupGamma0(II);;
gap> AbelianInvariants(GG);
[ 0, 0, 0, 0, 0, 0, 0, 0, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 
  2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 
  4, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 11, 16, 16, 16, 16, 17, 17, 17, 32, 37, 
  61, 61, 64, 64, 128, 263, 263, 263, 263, 512, 1024, 5099, 5099, 72043, 
  72043 ]

gap> III:=QuadraticIdeal(OQ,49+69*Sqrt(-1));
ideal of norm 7162 in O(GaussianRationals)
gap> GGG:=HAP_CongruenceSubgroupGamma0(III);;
gap> AbelianInvariants(GGG);
[ 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 5, 8, 8, 8, 8, 13, 13, 25, 59, 
  59, 179, 283, 283, 379, 857, 967, 967, 3769, 13537, 25601, 222659, 
  8180323, 8180323, 11450932001, 11450932001 ]
