#include<bits/stdc++.h>
using namespace std;
int a[20][20],ABS[20][20],l[20],f[20][20],hABS[20],tot=1;
int n,m,r,c;
void init(){
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			ABS[i][j]=0;
			for(int k=1;k<=c;k++){
				ABS[i][j]+=abs(a[i][l[k]]-a[j][l[k]]);
			}
		}
	}
	for(int i=1;i<=n;i++){
		hABS[i]=0;
		for(int j=2;j<=c;j++){
			hABS[i]+=abs(a[i][l[j-1]]-a[i][l[j]]);
		}
	}
	return;
}
int cmin=2e9;
void DP(){
	for(int i=0;i<=19;i++){
		for(int j=0;j<=19;j++) f[i][j]=1e9;
	}
	for(int i=0;i<=19;i++){
		f[i][0]=0;
	}
	for(int i=1;i<=n;i++){
		for(int j=1;j<=min(i,r);j++){
			if(j==1) f[i][j]=hABS[i];
			else if(i==j){
				f[i][j]=f[i-1][j-1]+hABS[i]+ABS[i-1][i];
			}
			else for(int k=j-1;k<i;k++){
				f[i][j]=min(f[i][j],f[k][j-1]+hABS[i]+ABS[k][i]);
			}
			if(j==r)
			cmin=min(f[n][r],cmin);
		}
	}
}
void dfs(int x){
	if(x>n){
		init();
		DP();
		return;
	}
	if(c-tot+1==n-x+1){
		l[tot++]=x;
		dfs(x+1);
		l[tot--]=0;
		return;
	}
	dfs(x+1);
	if(tot<=c){
		l[tot++]=x;
		dfs(x+1);
		l[tot--]=0;
		return;
	}
	return;
}
int main(){
	scanf("%d%d%d%d",&n,&m,&r,&c);
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++) scanf("%d",&a[i][j]);
	}
	dfs(1);
	printf("%d",cmin);
	return 0;
}