X, Y = map(int, input().split())

edges = []

for i in range(X):
    edges.append((i, (i + 1) % X))

cnt = 0

for i in range(X, X * Y):
    edges.append((cnt // (Y - 1), i))
    cnt += 1

N, M = X * Y, len(edges)
print(N, M)

for u, v in edges:
    print(u + 1, v + 1)