洛谷P1034 矩形覆盖

题目链接

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
int N,K,x[55],y[55],vis[55];
int ans=250000,disx[55][55],disy[55][55];
struct node
{
	int xa,ya,xb,yb;
}a[10];
bool judge(int i,int k)
{
	for(int j=1;j<=k;j++)
	{
		if(i!=j)
		{
			if(a[i].xa<=a[j].xa && a[i].xb>=a[j].xa && a[i].ya<=a[j].ya && a[i].yb>=a[j].ya) return 0;
			if(a[i].xa<=a[j].xb && a[i].xb>=a[j].xb && a[i].ya<=a[j].ya && a[i].yb>=a[j].ya) return 0;
			if(a[i].xa<=a[j].xa && a[i].xb>=a[j].xa && a[i].ya<=a[j].yb && a[i].yb>=a[j].yb) return 0;
			if(a[i].xa<=a[j].xb && a[i].xb>=a[j].xb && a[i].ya<=a[j].yb && a[i].yb>=a[j].yb) return 0;
		}
	}
	return 1;
}
void dfs(int n,int k,int sum)//?
{
	if(sum>ans || k>K) return;//剪枝
	if(n==N+1)
	{
		ans=min(ans,sum); return;
	}
	for(int i=1;i<=k;i++)//点被已有矩形包含 
	{
		if(a[i].xa<=x[n] && x[n]<=a[i].xb && a[i].ya<=y[n] && y[n]<=a[i].yb) 
		{
			dfs(n+1,k,sum);
			return;//?
		}
	} 
	if(k<K)//增加一个新的矩形 
	{
		a[k+1].xa=a[k+1].xb=x[n]; a[k+1].ya=a[k+1].yb=y[n];
		dfs(n+1,k+1,sum); 
	}
	for(int i=1;i<=k;i++)//扩展矩形 
	{
		node pre=a[i];
		if(a[i].xa>x[n]) 	  a[i].xa=x[n];
		else if(a[i].xb<x[n]) a[i].xb=x[n];
		if(a[i].ya>y[n]) 	  a[i].ya=y[n];
		else if(a[i].yb<y[n]) a[i].yb=y[n];
		if(judge(i,k))  
		{
			dfs(n+1,k,sum-(pre.xb-pre.xa)*(pre.yb-pre.ya)+(a[i].xb-a[i].xa)*(a[i].yb-a[i].ya));
		}
		a[i]=pre;//尝试结束，回溯 
	 } 
	/*蒟蒻的瓷代码 
	if(k==K+1&&n==N+1)
	{
		ans=min(ans,sum);
		return;
	}
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)
	{
		disx[i][j]=abs(x[i]-x[j]);
		disy[i][j]=abs(y[i]-y[j]);
		if(!vis[i]&&!vis[j])
		{
			vis[i]=1; vis[j]=1;
			dfs(n+2,k+1,sum+disx[i][j]*disy[i][j]);
			vis[i]=0; vis[j]=0;
		}
	}
	*/
}
int main()
{
	cin>>N>>K;
	for(int i=1;i<=N;i++)
	{
		cin>>x[i]>>y[i];
	}
	dfs(1,0,0);//? 
	printf("%d\n",ans);
}```