#define rep(i, n) for (int i = 0; i < (int)(n); i++) #define ALL(v) v.begin(), v.end() typedef long long ll; #include using namespace std; class Dis{ public: vector rank,p,siz; Dis(int s){ rank.resize(s,0); p.resize(s,0); siz.resize(s,1); rep(i,s) makeSet(i); } void makeSet(int x){ p[x]=x; rank[x]=0; } bool same(int x,int y){ return root(x)==root(y); } void unite(int x,int y){ if(same(x,y)) return; link(root(x),root(y)); } void link(int x,int y){ if(rank[x]>rank[y]){ p[y]=x; siz[x]+=siz[y]; } else{ p[x]=y; siz[y]+=siz[x]; if(rank[x]==rank[y]) rank[y]++; } } int root(int x){ if(x != p[x]) p[x]=root(p[x]); return p[x]; } int size(int x){ return siz[root(x)]; } }; int main(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); int h,w,n,d; cin>>h>>w>>n>>d; vector X(n),Y(n); rep(i,n) cin>>X[i]>>Y[i]; map,set> m; map,int> m1; rep(k,n){ m1[{X[k],Y[k]}]=k; } Dis ds=Dis(n); rep(k,n){ for(int i=-d;i<=d;i++){ for(int j=-d-abs(i);j<=d+abs(i);j++){ if(m1.count({X[k]+i,Y[k]+j})) ds.unite(k,m1[{X[k]+i,Y[k]+j}]); m[{X[k]+i,Y[k]+j}].insert(k); } } } int ma=0; rep(i,n){ if(ds.root(i)==i && ds.size(i)>1) ma++; } int tmp=1; for(auto [a,b]:m){ if(m1.count(a)) continue; set s; for(auto x:b){ s.insert(ds.root(x)); } tmp=max(tmp,(int)s.size()); } if(ma>0) cout<