p = 5
m = 1
F = GF(p)
Rx.<x> = PolynomialRing(F)
f = x
C = superelliptic(f, 1)
AS1 = as_cover(C, [C.x^3], prec = 200)
print(AS1.genus())