x,y = map(int, input().split())

n = x + x*(y-2) + x
m = n

print(n,m)

for i in range(x):
    a = i + 1
    b = (i+1) % x + 1
    print(a,b)

id = x+1

for i in range(x):
    root = i+1

    li = [root]
    for j in range(y-1):
        li.append(id)
        id += 1
    
    for j in range(len(li)-1):
        print(li[j], li[j+1])