x,y=map(int,input().split())
n=x*y
ans=[]
nj=x
for i in range(x):
  ans+=(i,(i+1)%x),
  ni=i
  for j in range(y-1):
    ans+=(ni,nj),
    ni,nj=nj,nj+1
print(n,n)
for i,j in ans:
  print(i+1,j+1)