USACO 2016 open【Field Reduction】

Description

农夫约翰的N(5<=N<=50000)头牛被定在了平面内的不同的位置。他想用栅栏(平行于x和y轴)围住所有的牛。他想这个栅栏尽可能小(牛在边界上也被视作围住)。

他因为牛奶产量低而感到经费紧张，所以他想卖掉三头牛再围起剩下的牛。请算出栅栏围出的最小面积。(面积有可能为0，例如，最后剩下的牛排成一排或一列。）

Input Format

第一行输入正整数n

剩下2-n+1行，输入每头牛的坐标位置(x,y) ( 1<= x,y<=40000的整数）

Output Format

最小面积

Sample Input

6
1 1
7 8
10 9
8 12
4 100
50 7

Sample Output

12

【题解】

   一题比较水的题目，不知道为何我在考场上脑子抽了233

   只需要给牛按照x，y分别排序，用dfs，每次枚举要卖掉最左或者最右，最上或者最下的牛就可以了，到三头牛都卖掉了就算一下栅栏总面积

   我在考场上打的是先排序，然后算出从左，从右，从上，从下最多能卖掉几条线上的奶牛然后再去dfs每一条线（反正特别麻烦就是了，代码也很长）

   

   先贴一下正解（ 战神）的代码

    

#include <algorithm>
#include <cstring>
#include <cstdio>
const int inf=2000000000;
int b[50010],c[50010],d[50010],n,ans=inf;
struct P{
	int x,y;
}a[50010];
bool cmp1(int x,int y){return a[x].x<a[y].x || (a[x].x==a[y].x && a[x].y<a[y].y);}
bool cmp2(int x,int y){return a[x].y<a[y].y || (a[x].y==a[y].y && a[x].x<a[y].x);}
void work(int x,int y,int z){
	memset(d,0,sizeof d);
	for (int i=1;i<=x;i++) d[b[i]]=1;
	for (int i=1;i<=y;i++) d[b[n-i+1]]=1;
	int m=0;
	for (int i=1;i<=n;i++)
		if (!d[i]) c[++m]=i;
	std::sort(c+1,c+m+1,cmp2);
	for (int i=1;i<=z;i++) d[c[i]]=1;
	for (int i=1;i+x+y+z<=3;i++) d[c[m-i+1]]=1;
	int maxx=0,maxy=0,minx=inf,miny=inf;
	for (int i=1;i<=n;i++)
		if (!d[i]){
			maxx=std::max(maxx,a[i].x);
			minx=std::min(minx,a[i].x);
			maxy=std::max(maxy,a[i].y);
			miny=std::min(miny,a[i].y);
		}
	if (maxx<=minx || maxy<=miny) ans=0;
		else ans=std::min(ans,(maxx-minx)*(maxy-miny));
}
int main(){
	scanf("%d\n",&n);
	for (int i=1;i<=n;i++) scanf("%d%d\n",&a[i].x,&a[i].y);
	if (n<=4){
		printf("0\n");
		return 0;
	}
	for (int i=1;i<=n;i++) b[i]=i;
	std::sort(b+1,b+n+1,cmp1);
	for (int i=0;i<=3;i++)
		for (int j=0;i+j<=3;j++)
			for (int k=0;i+j+k<=3;k++) work(i,j,k);
	printf("%d\n",ans);
}

 

再贴一下我（蒟蒻）的代码（长到爆炸）

 

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <stack>
#include <vector>
#include <queue>
#include <map>
using namespace std;
int i,j,k,l,m,n;
struct info
  {
  	int x,y,u;
  }c[50005];
bool cmp1(int x,int y) {return c[x].x<c[y].x;}
bool cmp2(int x,int y) {return c[x].y<c[y].y;}
int line[15][2],rankx[50005],ranky[50005],num,ans;
void dfs(int k,int z)
  {
  	int i,l,j;
  	 if (k>=3)
  	   {
  	   	 int xmax,xmin,ymax,ymin;
  	   	 for (i=1;i<=4;i++) 
			if (!c[rankx[i]].u) {xmin=c[rankx[i]].x;break;}
  	   	 for (i=n;i>=n-3;i--) if (!c[rankx[i]].u) {xmax=c[rankx[i]].x;break;}
  	   	 for (i=1;i<=4;i++) if (!c[ranky[i]].u) {ymin=c[ranky[i]].y;break;}
  	   	 for (i=n;i>=n-3;i--) if (!c[ranky[i]].u) {ymax=c[ranky[i]].y;break;}
  	     ans=min(ans,(xmax-xmin)*(ymax-ymin));
  	     return;
	   }
	for (j=z;j<=m;j++)
	  {
	  	if (line[j][0]==1) 
	  	  {
			for (i=1;i<=3;i++) if ((!c[rankx[i]].u)&&(c[rankx[i]].x==line[j][1])) k++,c[rankx[i]].u=1;
			                    else if (c[rankx[i]].x==line[j][1]) c[rankx[i]].u++; 
  	   	    for (i=n;i>=n-2;i--) if ((!c[rankx[i]].u)&&(c[rankx[i]].x==line[j][1])) k++,c[rankx[i]].u=1;
  	   	                        else if (c[rankx[i]].x==line[j][1]) c[rankx[i]].u++; 
		  }
		if (line[j][0]==2)
		  {
		  	for (i=1;i<=3;i++) if ((!c[ranky[i]].u)&&(c[ranky[i]].y==line[j][1])) k++,c[ranky[i]].u=1;
		  	                     else if (c[ranky[i]].y==line[j][1]) c[ranky[i]].u++; 
  	   	    for (i=n;i>=n-2;i--) if ((!c[ranky[i]].u)&&(c[ranky[i]].y==line[j][1])) k++,c[ranky[i]].u=1;
  	   	                         else if (c[ranky[i]].y==line[j][1]) c[ranky[i]].u++; 
		  }
		if (k<=3) dfs(k,j+1);
		if (line[j][0]==1) 
	  	  {
			for (i=1;i<=3;i++) if ((c[rankx[i]].u==1)&&(c[rankx[i]].x==line[j][1])) k--,c[rankx[i]].u=0;
			                    else if (c[rankx[i]].x==line[j][1]) c[rankx[i]].u--; 
  	   	    for (i=n;i>=n-2;i--) if ((c[rankx[i]].u==1)&&(c[rankx[i]].x==line[j][1])) k--,c[rankx[i]].u=0;
  	   	                        else if (c[rankx[i]].x==line[j][1]) c[rankx[i]].u--; 
		  }
		if (line[j][0]==2)
		  {
		  	for (i=1;i<=3;i++) if ((c[ranky[i]].u==1)&&(c[ranky[i]].y==line[j][1])) k--,c[ranky[i]].u=0;
		  	                     else if (c[ranky[i]].y==line[j][1]) c[ranky[i]].u--; 
  	   	    for (i=n;i>=n-2;i--) if ((c[ranky[i]].u==1)&&(c[ranky[i]].y==line[j][1])) k--,c[ranky[i]].u=0;
  	   	                         else if (c[ranky[i]].y==line[j][1]) c[ranky[i]].u--; 
		  }
	  }
  } 
int main()
  {
  	scanf("%d",&n);
  	if (n<=3) {printf("0");return 0;}
  	for (i=1;i<=n;i++) scanf("%d%d",&c[i].x,&c[i].y),rankx[i]=ranky[i]=i;
  	sort(rankx+1,rankx+n+1,cmp1);
  	sort(ranky+1,ranky+n+1,cmp2);
  	for (i=1;i<=3;i++) 
  	  	if (c[rankx[i]].x!=c[rankx[i-1]].x) line[++m][1]=c[rankx[i]].x,line[m][0]=1;
  	if (c[rankx[4]].x==line[m][1]) m--;
  	for (i=n;i>=n-2;i--)
  	    if (c[rankx[i]].x!=c[rankx[i+1]].x) line[++m][1]=c[rankx[i]].x,line[m][0]=1;
    if (c[rankx[n-3]].x==line[m][1]) m--;
    for (i=1;i<=3;i++) 
  	  	if (c[ranky[i]].y!=c[ranky[i-1]].y) line[++m][1]=c[ranky[i]].y,line[m][0]=2;
  	if (c[ranky[4]].y==line[m][1]) m--;
  	for (i=n;i>=n-2;i--)
  	    if (c[ranky[i]].y!=c[ranky[i+1]].y) line[++m][1]=c[ranky[i]].y,line[m][0]=2;
    if (c[ranky[n-3]].y==line[m][1]) m--;
    ans=(c[rankx[n]].x-c[rankx[1]].x)*(c[ranky[n]].y-c[ranky[1]].y);
    num=0;
    dfs(0,1);
    printf("%d",ans);
  }