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

print( Z , Z )

for i in range(X) :
    P = i * Y
    print( P + 1 , ( P + Y ) % Z + 1 )
    for j in range(Y-1) :
        print( P + j + 1 , P + j + 1 + 1 )