import sys
input = lambda: sys.stdin.readline().rstrip()
# sys.setrecursionlimit(10**7)
# sys.set_int_max_str_digits(10**6)
# import pypyjit
# pypyjit.set_param('max_unroll_recursion=-1')
def mp():return map(int,input().split())
def lmp():return list(map(int,input().split()))
def lm1(LIST): return list(map(lambda x:x-1, LIST))
def mps(A):return [tuple(map(int, input().split())) for _ in range(A)]
def stoi(LIST):return list(map(int,LIST))
def itos(LIST):return list(map(str,LIST))
def atoi(LIST): return [ord(i)-ord("a") for i in LIST]
def Atoi(LIST): return [ord(i)-ord("A") for i in LIST]
def LT(LIST,N): return LIST[bisect.bisect_left(LIST,N)-1]
def LE(LIST,N): return LIST[bisect.bisect_right(LIST,N)-1]
def GT(LIST,N): return LIST[bisect.bisect_right(LIST,N)]
def GE(LIST,N): return LIST[bisect.bisect_left(LIST,N)]
def bitA(X,A):return X & 1<<A == 1<<A
def gtoi(x,y,h,w):return x*w+y
import math
import bisect
import heapq
import time
import random as rd
import itertools
from copy import copy as cc
from copy import deepcopy as dc
from itertools import accumulate, product
from collections import Counter, defaultdict, deque
# from atcoder.dsu import DSU
# from atcoder.fenwicktree import FenwickTree
# from atcoder.segtree import SegTree  # (op, ide_ele, LIST)
# from atcoder.lazysegtree import LazySegTree  # (op, ide_ele, mapping, composition, _id, lst)
def ceil(U,V):return (U+V-1)//V
def modf1(N,MOD):return (N-1)%MOD+1
def pmat(list):
    for i in list:print(*i)
m4 = [[1,0],[0,1],[-1,0],[0,-1]]
m8 = [[-1,-1],[-1,0],[-1,1],[0,-1],[0,1],[1,-1],[1,0],[1,1]]
inf = (1<<63)-1
mod = 998244353

# # n,m = mp()
# # ans = 0
# # for i in range(n):
# #     u,v = map(str,input().split())
# #     v = int(v)
# #     if v >= 1200 and u[:4].count("o") != 4:
# #         ans += 1
# print(ans)

# ok = 0
# ng = 10**9+1
# while (ok-ng) != 1:
#     m = (ok+ng)//2
#     print(m)
#     ans = int(input())
#     if ans == 1:
#         exit()
#     elif ans == 0:
#         ng = m

x,y = mp()
print(x*y,x*y)
for i in range(x-1):
    print(i+1,i+2)
print(x,1)
for i in range(1,x+1):
    for j in range(1,y):
        print(i+(j-1)*x,j*x+i)