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

edge = []

cycle = list(range(1, X+1))
for i in range(X):
    edge.append((cycle[i], cycle[(i+1)%X]))
pre = list(range(1, X+1))
for i in range(X*(Y-1)):
    edge.append((pre[i%X], pre[i%X]+X))
    pre[i%X] += X

print(X*Y, X*Y)
for u, v in edge:
    print(u, v)