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的更多相关文章
- BZOJ 4576: [Usaco2016 Open]262144
Description 一个序列,每次可以将两个相同的数合成一个数,价值+1,求最后最大价值 \(n \leqslant 262144\) Sol DP. 这道题是 BZOJ 4580: [Usaco ...
- 【BZOJ 4581】【Usaco2016 Open】Field Reduction
http://www.lydsy.com/JudgeOnline/problem.php?id=4581 考虑\(O(n^3)\)暴力. 实际上枚举最靠边的三个点就可以了,最多有12个点. 还是暴力= ...
- BZOJ 4742: [Usaco2016 Dec]Team Building
4742: [Usaco2016 Dec]Team Building Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 21 Solved: 16[Su ...
- bzoj 4506: [Usaco2016 Jan]Fort Moo
4506: [Usaco2016 Jan]Fort Moo Description Bessie is building a fort with her friend Elsie. Like any ...
- bzoj 4412: [Usaco2016 Feb]Circular Barn
4412: [Usaco2016 Feb]Circular Barn Description 有一个N个点的环,相邻两个点距离是1.点顺时针标号为1..N.每一个点有ci头牛,保证∑ci=N.每头牛都 ...
- Bzoj 4582 [Usaco2016 Open] Diamond Collector 题解
4582: [Usaco2016 Open]Diamond Collector Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 204 Solved: ...
- BZOJ 4582: [Usaco2016 Open]Diamond Collector
Descrirption 给你一个长度为 \(n\) 的序列,求将它分成两个序列后最多个数,每个序列最大值最小值不能超过 \(k\) Sol 二分+DP. 排一下序,找出以这个点结尾和开始的位置. 这 ...
- BZOJ 4580: [Usaco2016 Open]248
Description 一个序列,每次可以把相邻的两个数合为一个,价值+1,求最后的最大价值. Sol 区间DP. \(f[i][j]\) 表示 \(i-j\) 中合成一个数字为多少,转移就是枚举断点 ...
- BZOJ 4579: [Usaco2016 Open]Closing the Farm
Description 依次删去一个点和它的边,问当前图是否连通. Sol 并查集. 倒着做就可以了. 每次将一个点及其的边加入,如果当前集合个数大于 1,那么就不连通. Code /******** ...
随机推荐
- 嵌入式Linux驱动学习之路(十三)按键驱动-异步通知
之前的按键方式: 查询: 极度占用CPU资源 中断: 在读的时候产生休眠,在没有信号的时候永远不会返回. poll机制: 在中断的基础上加上超时时间. 异步通知就是通过信号来传送. 首先在应用程序中有 ...
- Visual Studio 2013支持Xamarin的解决方案
转自博客园[遗忘的代码] Xamarin的Visual Studio插件目前还不支持VS 2013,所以需要在安装Xamarin的VS插件时把2010和2012全选上 (由于我的电脑里只剩2013,而 ...
- jquery的add()用法总结
<!DOCTYPE html> <html> <head lang="en"> <meta charset="utf-8&quo ...
- zookeeper安装
http://blog.itpub.net/27099995/viewspace-1394831/ http://blog.csdn.net/huwei2003/article/details/491 ...
- 3942: [Usaco2015 Feb]Censoring [KMP]
3942: [Usaco2015 Feb]Censoring Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 375 Solved: 206[Subm ...
- javaScript中的页面传值
function getURIParam(name) { var search = window.location.search; search = search.substring(1); if ( ...
- mvc SelectList selected失效的解决方法
ViewBag.QuestionnaireType = new SelectList(questionCollectionTypeList, "CodeID", "Cod ...
- .NET Web API 新手遇到的那些“坑”
以前一直做web项目,用的mvc.最近公司需要上线APP,准备用webAPI做接口,用习惯了MVC,API还是有好多细节的地方感觉很“坑”... 第一个坑,webAPI的默认路由规则不一样,通过命名来 ...
- 可变长子网掩码 VLSM
2014-05-05 17:26:08 标签:IP地址 路由器 技术部 路由表 长子 VLSM的介绍: VLSM(VLSM(Variable Length Subnetwork Mask,可变长子网掩 ...
- 利用Microsoft.Practices.Unity的拦截技术,实现.NET中的AOP
1.记住这个单词的意思:Interception(拦截) 2.首先说一下原理和背景 原理:所谓的AOP就是面向切面编程,这里不多说,百度搜索. 目的:个人认为是为了解耦,部分代码跟业务代码分离,业务代 ...