求减枝方法

#include<bits/stdc++.h>
using namespace std;
int n,m,r,c,x[1000],y[1000],cnt1,cnt2,ans=99999999;
int p[100][100],vis1[1000],vis2[1000];
void dfsc(int o)
{
	if(o>c)
	{
		return ;
	}
	if(o==c)
	{
		int	sum=0;
		int k[100][100];
		for(int i=1;i<=cnt1;i++)
		{
			for(int j=1;j<=cnt2;j++)
			{
				k[i][j]=p[x[i]][y[j]];
			}
		}
		for(int i=1;i<=cnt2;i++)
		{
			for(int j=1;j<=cnt1-1;j++)
			{
				sum+=abs(k[j][i]-k[j+1][i]);
			}
		}
		for(int i=1;i<=cnt1;i++)
		{
			for(int j=1;j<=cnt2-1;j++)
			{
				sum+=abs(k[i][j]-k[i][j+1]);
			}
		}		
		ans=min(ans,sum);
	}
	for(int i=1;i<=m;i++)
	{
		if(vis2[i]==0)
		{
			if(i<y[cnt2])
			{
				continue;
			}
			y[++cnt2]=i;
			vis2[i]=1;
			dfsc(o+1);
			y[cnt2]=0;
			cnt2--;
			vis2[i]=0;
		}
	}
	return ;
}
void dfsr(int o)
{
	if(o>r)
	{
		return ;
	}
	if(o==r)
	{
		dfsc(0);
	}
	for(int i=1;i<=n;i++)
	{
		if(vis1[i]==0)
		{
			if(i<x[cnt1])
			{
				continue;
			}
			x[++cnt1]=i;
			vis1[i]=1;
			dfsr(o+1);
			x[cnt1]=0;
			cnt1--;
			vis1[i]=0;
		}
	}
}
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",&p[i][j]);
		}
	}
	dfsr(0);
	cout<<ans;
	return 0;
}
/*
5 5 2 3
9 3 3 3 9
9 4 8 7 4
1 7 4 6 6
6 8 5 6 9
7 4 5 6 1
*/