X, Y = map(int, input().split())
n=X*Y
m=X*Y
print(n,m)
for i in range(X):
  u = i*Y+1
  v = ((i+1)*Y+1)%(X*Y)
  print(u,v)
  for j in range(Y-1):
    u = (i*Y)+(j+1)
    v = u+1
    print(u, v)