import sys input = sys.stdin.readline class UnionFind(): def __init__(self, n=1): self.par = [i for i in range(n)] self.rank = [0 for _ in range(n)] self.size = [1 for _ in range(n)] def find(self, x): if self.par[x] == x: return x else: self.par[x] = self.find(self.par[x]) return self.par[x] def union(self, x, y): x = self.find(x) y = self.find(y) if x != y: if self.rank[x] < self.rank[y]: x, y = y, x if self.rank[x] == self.rank[y]: self.rank[x] += 1 self.par[y] = x self.size[x] += self.size[y] def is_same(self, x, y): return self.find(x) == self.find(y) def get_size(self, x): x = self.find(x) return self.size[x] H, W, N, D = map(int, input().split()) X, Y = [0] * N, [0] * N U = UnionFind(N) DD = dict() for i in range(N): X[i], Y[i] = map(int, input().split()) X[i], Y[i] = X[i] - 1, Y[i] - 1 DD[(X[i], Y[i])] = i for i in range(N): for a in range(X[i] - D, X[i] + D + 1): for b in range(Y[i] - D, Y[i] + D + 1): if abs(a - X[i]) + abs(b - Y[i]) > D: continue if a < 0 or a > H - 1 or b < 0 or b > W - 1: continue if (a, b) not in DD: continue ind = DD[(a, b)] U.union(i, ind) mi, ma = 0, -10 ** 18 for i in range(H): for j in range(W): if (i, j) in DD: continue S = set() v = 0 flag = 0 for a in range(i - D, i + D + 1): for b in range(j - D, j + D + 1): if abs(a - i) + abs(b - j) > D: continue if a < 0 or a > H - 1 or b < 0 or b > W - 1: continue if (a, b) not in DD: continue ind = U.find(DD[(a, b)]) if ind in S: continue S.add(ind) if U.get_size(ind) >= 2: v -= 1 else: flag = 1 if v < 0: v += 1 mi = min(mi, v + flag) ma = max(ma, v + flag) ans = 0 S = set() for i in range(N): if U.find(i) in S: continue S.add(U.find(i)) if U.get_size(i) >= 2: ans += 1 print(ans + mi, ans + ma)