def solve(A, B):
	N = A*A + B*B
	ans_A = []
	ans_B = []
	for i in range(N):
		ans_A.append((A*i % N, i*B % N))
	for i in range(N):
		ans_B.append((A*(i+1) % N, i*B % N))
	return ans_A, ans_B

A, B = map(int, input().split())
N = A*A + B*B

if A < B:
	ans_B, ans_A = solve(B, A)
else:
	ans_A, ans_B = solve(A, B)

for i in range(N):
	print(*ans_A[i])
for i in range(N):
	print(*ans_B[i])