class unif:
  def __init__(self,n):
    self.pare=[-1]*n
    self.size=[1]*n
  def root(self,x):
    while self.pare[x]!=-1:
      x=self.pare[x]
    return x
  def unite(self,u,v):
    rootu=self.root(u)
    rootv=self.root(v)
    if rootu!=rootv:
      if self.size[rootu]>=self.size[rootv]:
        self.pare[rootv]=rootu
        self.size[rootu]+=self.size[rootv]
      else:
        self.pare[rootu]=rootv
        self.size[rootv]+=self.size[rootu]
  def same(self,s,t):
    return self.root(s)==self.root(t)
H,W,N,D=map(int,input().split())
C=set()
v=[[-1]*(W+1) for i in range(H+1)]
for i in range(N):
  x,y=map(int,input().split())
  v[x][y]=i
Z=unif(N)
for x in range(1,H+1):
  for y in range(1,W+1):
    if v[x][y]==-1:
      continue
    for dx in range(-D,D+1):
      k=D-abs(dx)
      for dy in range(-k,k+1):
        x2,y2=x+dx,y+dy
        if x2<=0 or x2>=H+1 or y2<=0 or y2>W:
          continue
        pos=v[x2][y2]
        if pos>=0:
          Z.unite(v[x][y],pos)
b=set()
for i in range(N):
  p=Z.root(i)
  if Z.size[p]>1:
    b.add(p)
count=len(b)
ans1=10**10
ans2=-10**10
for x in range(1,H+1):
  for y in range(1,W+1):
    if v[x][y]>=0:
      continue
    c=0
    k=set()
    for dx in range(-D,D+1):
      k2=D-abs(dx)
      for dy in range(-k2,k2+1):
        x2,y2=x+dx,y+dy
        if x2<=0 or x2>=H+1 or y2<=0 or y2>W:
          continue
        pos=v[x2][y2]
        if pos>=0:
          l=Z.root(pos)
          if Z.size[l]>1:
            k.add(l)
          else:
            c=1
    if c==0:
      if len(k)==0:
        ans1=min(ans1,count)
        ans2=max(ans2,count)
      else:
        ans1=min(ans1,count+1-len(k))
        ans2=max(ans2,count+1-len(k))
    else:
      ans1=min(ans1,count+1-len(k))
      ans2=max(ans2,count+1-len(k))
print(ans1,ans2)