X,Y=map(int,input().split())
print(X*Y,X*Y)
mx=[0]*X
for i in range(X*(Y-1)):
    print(i+1,i+1+X)
    mx[i%X]=max(mx[i%X],i+1+X)
for i in range(X):
    print(mx[i],mx[(i+1)%X])