#include<bits/stdc++.h>
using namespace std;

//* ATCODER
#include<atcoder/all>
using namespace atcoder;
typedef modint998244353 mint;
//*/

/* BOOST MULTIPRECISION
#include<boost/multiprecision/cpp_int.hpp>
using namespace boost::multiprecision;
//*/

typedef long long ll;

#define rep(i, s, n) for (int i = (int)(s); i < (int)(n); i++)
#define rrep(i, s, n) for (int i = (int)(n)-1; i >= (int)(s); i--)

template <typename T> bool chmin(T &a, const T &b) {
	if (a <= b) return false;
	a = b;
	return true;
}

template <typename T> bool chmax(T &a, const T &b) {
	if (a >= b) return false;
	a = b;
	return true;
}

template <typename T> T max(vector<T> &a){
	assert(!a.empty());
	T ret = a[0];
	for (int i=0; i<(int)a.size(); i++) chmax(ret, a[i]);
	return ret;
}

template <typename T> T min(vector<T> &a){
	assert(!a.empty());
	T ret = a[0];
	for (int i=0; i<(int)a.size(); i++) chmin(ret, a[i]);
	return ret;
}

template <typename T> T sum(vector<T> &a){
	T ret = 0;
	for (int i=0; i<(int)a.size(); i++) ret += a[i];
	return ret;
}

struct UndoUnionFind{
	int n;
	vector<int> par;
	vector<pair<int,int>> history;

	explicit UndoUnionFind(int v): n(v), par(v,-1) {}

	void merge(int a, int b){
		int x = find(a);
		int y = find(b);
		history.push_back(pair(x, par[x]));
		history.push_back(pair(y, par[y]));
		if (x == y) return;
		if (-par[x] < -par[y]) swap(x, y);
		par[x] += par[y];
		par[y] = x;
		return;
	}

	int find(int x){
		if (par[x] < 0) return x;
		return find(par[x]);
	}

	void undo(){
		par[history.back().first] = history.back().second;
		history.pop_back();
		par[history.back().first] = history.back().second;
		history.pop_back();
		return;
	};
};

int main(){
	int h, w; cin >> h >> w;
	int n, d; cin >> n >> d;

	vector<int> x(n), y(n);
	rep(i,0,n) cin >> x[i] >> y[i];

	vector taio(h+1, vector<int>(w+1, -1));

	auto exist = [&](int i, int j) -> int {
		if (! (0 <= i && i <= h)) return -1;
		if (! (0 <= j && j <= w)) return -1;	
		return taio[i][j];
	};

	rep(i,0,n){
		taio[x[i]][y[i]] = i;
	}

	UndoUnionFind uf(n+1);
	int val = 0;
	rep(i,0,n){
		rep(dx,-d,d+1){
			rep(dy,-d,d+1){
				if (abs(dx)+abs(dy)>d) continue;
				if (exist(x[i]+dx,y[i]+dy)==-1){
					continue;
				}
				if (uf.find(exist(x[i]+dx,y[i]+dy)) != uf.find(i)){
					if (-uf.par[uf.find(exist(x[i]+dx,y[i]+dy))] >= 2){
						val--;
					}
					if (-uf.par[uf.find(i)] >= 2){
						val--;
					}
					uf.merge(exist(x[i]+dx,y[i]+dy), i);
					if (-uf.par[uf.find(i)] >= 2){
						val++;
					}
				}
			}
		}
	}

	int mn = 1e9;
	int mx = -1;

	rep(i,1,h+1){
		rep(j,1,w+1){
			if (exist(i,j) >= 0){
				continue;
			}
			int cnt = 0;
			int orval = val;
			rep(dx,-d,d+1){
				rep(dy,-d,d+1){
					if (abs(dx)+abs(dy) > d) continue;
					if (exist(i+dx,j+dy) == -1){
						continue;
					}
					if (uf.find(exist(i+dx, j+dy)) != uf.find(n)){
						if (-uf.par[uf.find(exist(i+dx, j+dy))] >= 2){
							val--;
						}
						if (-uf.par[uf.find(n)] >= 2){
							val--;
						}
						uf.merge(exist(i+dx, j+dy), n);
						
						if (-uf.par[uf.find(n)] >= 2){
							val++;
						}
						cnt++;
					}
				}
			}


			chmax(mx,val);
			chmin(mn,val);

			rep(num,0,cnt){
				uf.undo();
			}
			val = orval;
		}
	}
	cout << mn << ' ' << mx << '\n';
}