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

V = [[] for _ in range(X*Y)]

for i in range(X):
    V[i].append((i+1)%X)

for i in range(Y-1):
    for j in range(X):
        V[(i*X)+j].append(((i+1)*X)+j)

print(X*Y,Y*X)
for i,v in enumerate(V):
    for j,e in enumerate(v):
        print(i+1,e+1)