/*
考虑统计关键点小于 k 个的子矩形数量。
枚举上边界，然后从小到大扫下边界，可以维护 c[i] 表示第 i 列的关键点数。
就是对左右边界做双指针，答案即 sum(l) r-l，其中 r 为 [l...r] 的 c 之和不超过 k 的最大 r。
对每个 l 维护 r[l]。c=0 的部分形成一些连续段 [L,R]，每段内的 r 等于 R+1 的 r。
维护每个 c>0 的 l 的 r，可以用链表找前后下一个 c>0 的位置。O(N(R+C)k)。
*/
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int inf=1e9;
int N,M,C,k;
vector<int> rimp[3100];
int c[3100],nxt[3100],pre[3100],rp[3100];
int res;
void twoPointers(int ll,int lr)
{
	for(int l=ll,r=ll,cnt=0;l<=lr;cnt-=c[l],l=nxt[l])
	{
		// cout<<"l="<<l<<" cnt="<<cnt<<endl;
		while(cnt+c[r]<=k)
			cnt+=c[r],r=nxt[r];
		// cout<<"l="<<l<<" r <- "<<r-1<<endl;
		// cout<<"cnt="<<cnt<<" c[r]="<<c[r]<<endl;
		res+=(l-pre[l]-(l==M+1))*((r-1)-rp[l]);
		rp[l]=r-1;
		if(nxt[l]==l)
			break;
	}
}
void build(int U)
{
	// cout<<"Build U="<<U<<endl;
	memset(c,0,sizeof(c));
	memset(nxt,0,sizeof(nxt));
	memset(pre,0,sizeof(pre));
	memset(rp,0,sizeof(rp));

	for(int i=U;i<=N;i++)
	{
		for(int p:rimp[i])
			c[p]+=1;
	}
	c[0]=c[M+1]=inf;
	pre[0]=0,nxt[M+1]=M+1;
	for(int i=1,p=0;i<=M+1;i++)
	{
		if(c[i])
			tie(pre[i],p)=make_pair(p,i);
	}
	for(int i=M,p=M+1;i>=0;i--)
	{
		if(c[i])
			tie(nxt[i],p)=make_pair(p,i);
	}

	// cout<<"c   :";
	// for(int i=0;i<=M+1;i++)
	// 	cout<<" "<<c[i];
	// cout<<endl;
	// cout<<"pre :";
	// for(int i=0;i<=M+1;i++)
	// 	cout<<" "<<pre[i];
	// cout<<endl;
	// cout<<"nxt :";
	// for(int i=0;i<=M+1;i++)
	// 	cout<<" "<<nxt[i];
	// cout<<endl;

	res=0;
	twoPointers(0,M+1);
	res-=(M-1)*M/2;

	// cout<<"rp  :";
	// for(int i=0;i<=M+1;i++)
	// 	cout<<" "<<rp[i];
	// cout<<endl;
}
void decreaseC(int p)
{
	// cout<<"DecC "<<p<<endl;
	if((--c[p])==0)
	{
		res+=(p-pre[p])*(rp[nxt[p]]-rp[p]);
		tie(nxt[pre[p]],pre[nxt[p]])=make_pair(nxt[p],pre[p]);
		tie(pre[p],nxt[p],p)=make_tuple(0,0,pre[p]);
	}
	int lr=p,ll=p;
	for(int i=1;i<=k;i++)
		ll=pre[ll];
	twoPointers(ll,lr);
}
signed main()
{
	cin>>N>>M>>C>>k;
	k--;
	for(int i=1;i<=C;i++)
	{
		int x,y;
		cin>>x>>y;
		rimp[x].push_back(y);
	}
	int ans=0;
	for(int U=1;U<=N;U++)
	{
		build(U);
		for(int D=N;D>=U;D--)
		{
			// cout<<"D="<<D<<" : res="<<res<<endl;
			ans+=res;
			for(int p:rimp[D])
				decreaseC(p);

			// cout<<"c   :";
			// for(int i=0;i<=M+1;i++)
			// 	cout<<" "<<c[i];
			// cout<<endl;
			// cout<<"pre :";
			// for(int i=0;i<=M+1;i++)
			// 	cout<<" "<<pre[i];
			// cout<<endl;
			// cout<<"nxt :";
			// for(int i=0;i<=M+1;i++)
			// 	cout<<" "<<nxt[i];
			// cout<<endl;
			// cout<<"rp  :";
			// for(int i=0;i<=M+1;i++)
			// 	cout<<" "<<rp[i];
			// cout<<endl;
		}
	}
	// cout<<ans<<endl;
	cout<<(N*(N+1)/2)*(M*(M+1)/2)-ans<<'\n';
}
/*
2 2 1 1
1 2

3 2 3 2
1 1
3 1
2 2
*/