首页 > 代码库 > BZOJ 4581: [Usaco2016 Open]Field Reduction
BZOJ 4581: [Usaco2016 Open]Field Reduction
Description
有 \(n\) 个点,删掉三个点后,求最小能围住的面积.
Sol
搜索.
找出 左边/右边/上边/下边 的几个点枚举就可以了.
我找了 12 个点,统计一下坐标的个数,然后找到最大最小坐标.
PS:md.Cena 好像对 STL 不太友好啊,用 map 居然T了???exm???
昨天用 map 离散化就T的飞起.md..mama我再也不用 map 离散化了...以后用lower_bound()好了...
一开始写了个DFS来枚举...发现不是很好...又删掉变成三维循环枚举...代码奇丑无比...
Code
/************************************************************** Problem: 4581 User: BeiYu Language: C++ Result: Accepted Time:440 ms Memory:5168 kb****************************************************************/ #include <cstdio>#include <map>#include <utility>#include <algorithm>#include <iostream>using namespace std; typedef long long LL;const int N = 50050;const int M = 15;int C = 12; struct P{ int x,y,id; };bool operator < (const P &a,const P &b){ return a.x==b.x?a.y<b.y:a.x<b.x; }bool operator == (const P &a,const P &b){ return (a.x==b.x) && (a.y==b.y); } int n;LL ans;int b[N],xx[N],yy[N],cx[N],cy[N];P a[N],p[N];P mxx[M],mxy[M],mix[M],miy[M];P del[M]; map< int,int > mpx,mpy; inline int in(int x=0){ scanf("%d",&x);return x; }/*LL calc(){// cout<<"--------------"<<endl;// for(int i=1;i<=3;i++) cout<<del[i].x<<" "<<del[i].y<<endl; int minx=1,maxx=n,miny=1,maxy=n; for(int i=3;i>=1;i--){ if(del[i].x == xx[minx]) minx++; if(del[i].x == xx[maxx]) maxx--; if(del[i].y == yy[miny]) miny++; if(del[i].y == yy[maxy]) maxy--; } for(int i=3;i>=1;i--){ if(del[i].x == xx[minx]) minx++; if(del[i].x == xx[maxx]) maxx--; if(del[i].y == yy[miny]) miny++; if(del[i].y == yy[maxy]) maxy--; } for(int i=3;i>=1;i--){ if(del[i].x == xx[minx]) minx++; if(del[i].x == xx[maxx]) maxx--; if(del[i].y == yy[miny]) miny++; if(del[i].y == yy[maxy]) maxy--; }// cout<<xx[maxx]<<" "<<xx[minx]<<" "<<yy[maxy]<<" "<<yy[miny]<<endl;// cout<<1LL*(xx[maxx]-xx[minx])*(yy[maxy]-yy[miny])<<endl; return 1LL*(xx[maxx]-xx[minx])*(yy[maxy]-yy[miny]);}void DFS(int w,int d){ if(w==4 || d==0){ if(!d) ans=min(ans,calc());return; } int c=0; DFS(w+1,d); int tmp[M]; switch(w){ case 0: for(int i=1;i<=C;i++){ if(c>d) break; if(!b[mix[i].id]) tmp[++c]=mix[i].id,del[d-c+1]=mix[i],b[mix[i].id]=1,DFS(w+1,d-c); }break; case 1: for(int i=1;i<=C;i++){ if(c>d) break; if(!b[mxx[i].id]) tmp[++c]=mxx[i].id,del[d-c+1]=mxx[i],b[mxx[i].id]=1,DFS(w+1,d-c); }break; case 2: for(int i=1;i<=C;i++){ if(c>d) break; if(!b[miy[i].id]) tmp[++c]=miy[i].id,del[d-c+1]=miy[i],b[miy[i].id]=1,DFS(w+1,d-c); }break; case 3: for(int i=1;i<=C;i++){ if(c>d) break; if(!b[mxy[i].id]) tmp[++c]=mxy[i].id,del[d-c+1]=mxy[i],b[mxy[i].id]=1,DFS(w+1,d-c); }break; } for(int i=1;i<=c;i++) b[tmp[i]]=0; return;}*/int main(){ n=in();C=min(C,n),ans=1e18; for(int i=1;i<=n;i++) a[i].x=in(),a[i].y=in(),a[i].id=i; for(int i=1;i<=n;i++) xx[i]=a[i].x,yy[i]=a[i].y; sort(xx+1,xx+n+1),sort(yy+1,yy+n+1); cx[0]=unique(xx+1,xx+n+1)-(xx+1),cy[0]=unique(yy+1,yy+n+1)-(yy+1); for(int i=1;i<=cx[0];i++) mpx[xx[i]]=i; for(int i=1;i<=cy[0];i++) mpy[yy[i]]=i; for(int i=1;i<=n;i++) cx[mpx[a[i].x]]++,cy[mpy[a[i].y]]++; int cnt=0; sort(a+1,a+n+1); for(int i=1;i<=C;i++) p[++cnt]=a[i]; for(int i=1;i<=n;i++) a[i].x=-a[i].x; sort(a+1,a+n+1); for(int i=1;i<=C;i++) p[++cnt]=a[i],p[cnt].x=-p[cnt].x; for(int i=1;i<=n;i++) a[i].x=-a[i].x,swap(a[i].x,a[i].y); sort(a+1,a+n+1); for(int i=1;i<=C;i++) p[++cnt]=a[i],swap(p[cnt].x,p[cnt].y); for(int i=1;i<=n;i++) a[i].x=-a[i].x; sort(a+1,a+n+1); for(int i=1;i<=C;i++) p[++cnt]=a[i],p[cnt].x=-p[cnt].x,swap(p[cnt].x,p[cnt].y); sort(p+1,p+cnt+1);// for(int i=1;i<=cnt;i++) cout<<p[i].x<<" "<<p[i].y<<endl; cnt=unique(p+1,p+cnt+1)-(p+1); // for(int i=1;i<=cnt;i++) cout<<p[i].x<<" "<<p[i].y<<endl; for(int i=1;i<=cnt;i++) for(int j=i+1;j<=cnt;j++) for(int k=j+1;k<=cnt;k++){ cx[mpx[p[i].x]]--,cx[mpx[p[j].x]]--,cx[mpx[p[k].x]]--; cy[mpy[p[i].y]]--,cy[mpy[p[j].y]]--,cy[mpy[p[k].y]]--; int minx=1,maxx=cx[0],miny=1,maxy=cy[0]; while(!cx[minx]) minx++; while(!cx[maxx]) maxx--; while(!cy[miny]) miny++; while(!cy[maxy]) maxy--; ans=min(ans,1LL*(xx[maxx]-xx[minx])*(yy[maxy]-yy[miny])); cx[mpx[p[i].x]]++,cx[mpx[p[j].x]]++,cx[mpx[p[k].x]]++; cy[mpy[p[i].y]]++,cy[mpy[p[j].y]]++,cy[mpy[p[k].y]]++; } // for(int i=1;i<=n;i++) cout<<xx[i]<<" ";cout<<endl;// for(int i=1;i<=n;i++) cout<<yy[i]<<" ";cout<<endl; // sort(a+1,a+n+1);// for(int i=1;i<=C;i++) mix[i]=a[i];// // for(int i=1;i<=n;i++) a[i].x=-a[i].x;// sort(a+1,a+n+1);// for(int i=1;i<=C;i++) mxx[i]=a[i],mxx[i].x=-mxx[i].x;// // for(int i=1;i<=n;i++) a[i].x=-a[i].x,swap(a[i].x,a[i].y);// sort(a+1,a+n+1);// for(int i=1;i<=C;i++) miy[i]=a[i],swap(miy[i].x,miy[i].y);// // for(int i=1;i<=n;i++) a[i].x=-a[i].x;// sort(a+1,a+n+1);// for(int i=1;i<=C;i++) mxy[i]=a[i],mxy[i].x=-mxy[i].x,swap(mxy[i].x,mxy[i].y); // cout<<"--------------"<<endl;// for(int i=1;i<=C;i++) cout<<mix[i].x<<" "<<mix[i].y<<endl;// cout<<"--------------"<<endl;// for(int i=1;i<=C;i++) cout<<mxx[i].x<<" "<<mxx[i].y<<endl;// cout<<"--------------"<<endl;// for(int i=1;i<=C;i++) cout<<miy[i].x<<" "<<miy[i].y<<endl;// cout<<"--------------"<<endl;// for(int i=1;i<=C;i++) cout<<mxy[i].x<<" "<<mxy[i].y<<endl;// cout<<"--------------"<<endl; // DFS(0,3); cout<<ans<<endl; return 0;}
BZOJ 4581: [Usaco2016 Open]Field Reduction
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。