x, y = map(int,input().split())
n = x*y
m = n
ans = []
for i in range(x):
    for j in range(1,y):
        ans.append([i*y+j,i*y+j+1])

for i in range(1,x+1):
    if y*(i+1) > n:
        b = y*(i+1)%n
    else:
        b = y*(i+1)
    ans.append([y*i,b])

print(n,m)
for i in range(n):
    print(*ans[i])