(1, 1) 24
(1, 2) 48
(1, 3) 48
(1, 4) 48
(1, 5) 48
(1, 6) 24
(2, 1) 48
(2, 2) 8
(2, 3) 16
(2, 4) 16
(2, 5) 8
(2, 6) 48
(3, 1) 48
(3, 2) 16
(3, 3) 24
(3, 4) 24
(3, 5) 16
(3, 6) 48
(4, 1) 48
(4, 2) 16
(4, 3) 24
(4, 4) 24
(4, 5) 16
(4, 6) 48
(5, 1) 48
(5, 2) 8
(5, 3) 16
(5, 4) 16
(5, 5) 8
(5, 6) 48
(6, 1) 24
(6, 2) 48
(6, 3) 48
(6, 4) 48
(6, 5) 48
(6, 6) 24

def mul(a,b):
	return (a[0]*b[0]-a[1]*b[1], a[0]*b[1]+b[0]*a[1])

def mod(z, p):
	return (z[0] % p, z[1] % p)

def order (z, p, k):
	a = z
	if a[0] == 1:
		if a[1] == 0:
			return 1
	for i in range(2,k):
		a = mul(a, z)
		a = mod(a, p)
		if a[0] == 1:
			if a[1] == 0:
				return i
	else:
		return 0

def pr2(a):
	d = []
	for i in range(1,a):
		for j in range(1,a):
			z = mod((i,j), a)
			y = order(z, a, a**2)
			if y not in d:
				d.append(y)
	d.sort()
	return d