def H1(i,j,p): if i==0: return 1 if j==0: return -1 if i==j: return 1 return kronecker_symbol(i-j,p)
def H2(i,j,p): if i==0 and j==0: return 0 if i==0 or j==0: return 1 if j==0: return -1 if i==j: return 0 return kronecker_symbol(i-j,p)