from collections import defaultdict,deque

class UnionFind():
    def __init__(self, n):
        self.n = n
        self.parents = [-1] * n

    def find(self, x):
        if self.parents[x] < 0:
            return x
        else:
            self.parents[x] = self.find(self.parents[x])
            return self.parents[x]

    def union(self, x, y):
        x = self.find(x)
        y = self.find(y)

        if x == y:
            return

        if self.parents[x] > self.parents[y]:
            x, y = y, x

        self.parents[x] += self.parents[y]
        self.parents[y] = x

    def size(self, x):
        return -self.parents[self.find(x)]

    def same(self, x, y):
        return self.find(x) == self.find(y)

    def members(self, x):
        root = self.find(x)
        return [i for i in range(self.n) if self.find(i) == root]

    def roots(self):
        return [i for i, x in enumerate(self.parents) if x < 0]

    def group_count(self):
        return len(self.roots())

    def all_group_members(self):
        group_members = defaultdict(list)
        for member in range(self.n):
            group_members[self.find(member)].append(member)
        return group_members

H,W,N,D=map(int,input().split())
XY=[tuple(map(int,input().split())) for _ in range(N)]

mp=[[N for _ in range(W+1)] for _ in range(H+1)]
for (i,(x,y)) in enumerate(XY):mp[x][y]=i

nxt=[]
for i in range(-5,6):
    for j in range(-5,6):
        if i==0 and j==0:continue
        if abs(i)+abs(j)<=D:nxt.append((i,j))

uf=UnionFind(N)
for (i,j) in XY:
    for (di,dj) in nxt:
        ni=i+di
        nj=j+dj
        if 1<=ni<=H and 1<=nj<=W and mp[ni][nj]!=N:
            uf.union(mp[i][j],mp[ni][nj])

cand=set()
star=set()
sign=set()
star_cand=[set() for _ in range(N)]
sign_cand=[set() for _ in range(N)]
for (i,j) in XY:
    root=uf.find(mp[i][j])
    size=uf.size(root)
    if size==1:star.add(root)
    else:sign.add(root)

    for (di,dj) in nxt:
        ni=i+di
        nj=j+dj
        if 1<=ni<=H and 1<=nj<=W and mp[ni][nj]==N:
            cand.add((ni,nj))
            if size==1:star_cand[root].add((ni,nj))
            else:sign_cand[root].add((ni,nj))

star_num=len(star)
sign_num=len(sign)

star_cand_cnt=[[0 for _ in range(W+1)] for _ in range(H+1)]
sign_cand_cnt=[[0 for _ in range(W+1)] for _ in range(H+1)]
for v in star_cand:
    for (x,y) in v:star_cand_cnt[x][y]+=1
for v in sign_cand:
    for (x,y) in v:sign_cand_cnt[x][y]+=1

inc=0
dec=0
for (x,y) in cand:
    if star_cand_cnt[x][y]>0 and sign_cand_cnt[x][y]==0:inc+=1
    elif star_cand_cnt[x][y]==0 and sign_cand_cnt[x][y]>0:dec=max(dec,sign_cand_cnt[x][y])

mx=sign_num
mn=sign_num
if inc>0:mx+=1
if dec>0:mn-=dec-1

print(mn,mx)