import sys

sys.setrecursionlimit(2 * 10**3)

if sys.implementation.name == "pypy":
    import pypyjit

    pypyjit.set_param("max_unroll_recursion=-1")


from collections import defaultdict


class UnionFind:
    def __init__(self, N):
        self.N = N
        self.parent_or_size = [-1] * N
        self.group_count = N

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

    def unite(self, x, y):
        root_x = self.root(x)
        root_y = self.root(y)
        if root_x == root_y:
            return
        if -self.parent_or_size[root_x] < -self.parent_or_size[root_y]:
            root_x, root_y = root_y, root_x
        self.parent_or_size[root_x] += self.parent_or_size[root_y]
        self.parent_or_size[root_y] = root_x
        self.group_count -= 1

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

    def size(self, x):
        return -self.parent_or_size[self.root(x)]

    def groups(self):
        ret = defaultdict(list)
        for i in range(self.N):
            ret[self.root(i)].append(i)
        return [group for group in ret.values()]


input = lambda: sys.stdin.readline().rstrip()

SHIFT = 20


def conv(x, y):
    return x << SHIFT | y


H, W, N, D = map(int, input().split())

star = {}
XY = []
for i in range(N):
    x, y = map(int, input().split())
    XY.append((x, y))
    star[conv(x, y)] = i

uf = UnionFind(N)

cand = []
for i, (x, y) in enumerate(XY):
    for xd in range(-D, D + 1):
        for yd in range(-D, D + 1):
            d = abs(xd) + abs(yd)
            if d > D:
                continue
            xn = x + xd
            yn = y + yd
            cand.append((xn, yn))
            posn = conv(xn, yn)
            if posn in star:
                uf.unite(i, star[posn])

stella = uf.groups()
is_seiza = [0] * N

count = 0
for block in stella:
    if len(block) > 1:
        count += 1
        for i in block:
            is_seiza[i] = 1

mi = count
ma = count

for x, y in set(cand):
    seen = set()
    for xd in range(-D, D + 1):
        for yd in range(-D, D + 1):
            d = abs(xd) + abs(yd)
            if d > D:
                continue
            xn = x + xd
            yn = y + yd
            posn = conv(xn, yn)
            if posn not in star:
                continue
            seen.add(uf.root(star[posn]))
    c = count + 1
    for s in seen:
        if is_seiza[s]:
            c -= 1
    mi = min(mi, c)
    ma = max(ma, c)
print(mi, ma)