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

for i,j in P:
    print(i+1,j+1)