H,W,N,D=map(int, input().split()) A=[[-1]*W for _ in range(H)] for i in range(N): h,w=map(int, input().split()) h-=1;w-=1 A[h][w]=i B=[] for i in range(-D,D+1): for j in range(-D,D+1): if abs(i)+abs(j)<=D: B.append((i,j)) par=[i for i in range(N)];rank=[0]*N;size=[1]*N def find(x): if par[x]==x: return x else: par[x]=find(par[x]) return par[x] #同じ集合か判定 def same(x,y): return find(x)==find(y) def union(x,y): x=find(x) y=find(y) if x==y: return if rank[x]>rank[y]: par[y]=x size[x]+=size[y] else: par[x]=y size[y]+=size[x] if rank[x]==rank[y]: rank[y]+=1 for h in range(H): for w in range(W): if A[h][w]==-1: continue p=A[h][w] for x,y in B: nh,nw=h+x,w+y if 0<=nh