hdu 4052 Adding New Machine
2012-09-16 00:42
453 查看
题意:给你W*H大小的矩形,其中有N个地区不能使用(给出了这个地区的两个顶点的坐标即(x1,y1)和(x2,y2)),问能下多少个1*M的矩形。
分别统计1*M的矩形横着放能放多少个,竖着放能放多少个。可以转化成统计哪些位置可以放这矩形的第一个方格。比如说,当我们统计横着能放多少个的时候,问题转化成找出有多少个这样的方格(i,j),使得(i,j+1),(i,j+2)……(i,j+m-1)都是可以被使用的。易得,(x1,y1-m+1)到(x2,y2)这个区域内的方格都是不满足要求的。把不满足条件的方格抠掉,剩下就是满足条件的方格了,问题就变成用扫描线求矩形的面积了。要注意M=1的情况,并且边界上也有一些方格不满足条件,也要抠掉。
分别统计1*M的矩形横着放能放多少个,竖着放能放多少个。可以转化成统计哪些位置可以放这矩形的第一个方格。比如说,当我们统计横着能放多少个的时候,问题转化成找出有多少个这样的方格(i,j),使得(i,j+1),(i,j+2)……(i,j+m-1)都是可以被使用的。易得,(x1,y1-m+1)到(x2,y2)这个区域内的方格都是不满足要求的。把不满足条件的方格抠掉,剩下就是满足条件的方格了,问题就变成用扫描线求矩形的面积了。要注意M=1的情况,并且边界上也有一些方格不满足条件,也要抠掉。
#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <vector> #include <map> using namespace std; #define LL(x) (x<<1) #define RR(x) (x<<1|1) typedef long long LL; const int N=100005; struct Line { int x,y1,y2,valu; Line(){} Line(int a,int b,int c,int d){x=a;y1=b;y2=c;valu=d;} bool operator<(const Line &b)const { return x<b.x; } }; struct node { int lft,rht,sum,valu; int mid(){return lft+(rht-lft)/2;} }; map<int,int> imap; vector<int> y; vector<Line> line; int data [4]; struct Segtree { node tree[N*4]; void relax(int ind) { if(tree[ind].valu>0) tree[ind].sum=y[tree[ind].rht]-y[tree[ind].lft]; else { if(tree[ind].lft+1==tree[ind].rht) tree[ind].sum=0; else tree[ind].sum=tree[LL(ind)].sum+tree[RR(ind)].sum; } } void build(int lft,int rht,int ind) { tree[ind].lft=lft; tree[ind].rht=rht; tree[ind].sum=0; tree[ind].valu=0; if(lft+1!=rht) { int mid=tree[ind].mid(); build(lft,mid,LL(ind)); build(mid,rht,RR(ind)); } } void updata(int be,int end,int ind,int valu) { int lft=tree[ind].lft,rht=tree[ind].rht; if(be<=lft&&rht<=end) { tree[ind].valu+=valu; relax(ind); } else { int mid=tree[ind].mid(); if(be<mid) updata(be,end,LL(ind),valu); if(end>mid) updata(be,end,RR(ind),valu); relax(ind); } } }seg; LL solve(int n,int w,int h,int m) { y.clear(); line.clear(); imap.clear(); y.push_back(1); y.push_back(h); line.push_back(Line(max(1,w-m),1,h,1)); line.push_back(Line(w,1,h,-1)); for(int i=0;i<n;i++) { int x1=max(1,data[i][0]-m),y1=data[i][1]; int x2=data[i][2],y2=data[i][3]; line.push_back(Line(x1,y1,y2,1)); line.push_back(Line(x2,y1,y2,-1)); y.push_back(y1);y.push_back(y2); } sort(y.begin(),y.end()); y.erase(unique(y.begin(),y.end()),y.end()); for(int i=0;i<(int)y.size();i++) imap.insert(make_pair(y[i],i)); sort(line.begin(),line.end()); LL res=0; seg.build(0,(int)y.size(),1); for(int i=0;i<(int)line.size();i++) { if(i!=0) res+=(LL)seg.tree[1].sum*(line[i].x-line[i-1].x); seg.updata(imap[line[i].y1],imap[line[i].y2],1,line[i].valu); } return res; } int main() { int w,h,n,m; while(scanf("%d%d%d%d",&w,&h,&n,&m)!=EOF) { for(int i=0;i<n;i++) for(int j=0;j<4;j++) { scanf("%d",&data[i][j]); if(j==2||j==3) data[i][j]++; } LL res1=(LL)w*h-solve(n,w+1,h+1,m-1); for(int i=0;i<n;i++) { swap(data[i][0],data[i][1]); swap(data[i][2],data[i][3]); } LL res2=(LL)w*h-solve(n,h+1,w+1,m-1); if(m!=1) printf("%I64d\n",res1+res2); else printf("%I64d\n",res1); } return 0; }
相关文章推荐
- HDU 4052 Adding New Machine (线段树+离散化)
- Hdu 4052 Adding New Machine(给你W*H大小的矩形,其中有N个地区不能使用(给出了这个地区的两个顶点的坐标即(x1,y1)和(x2,y2)),问能下多少个1*M的矩形)
- HDU 4052 Adding New Machine(矩形面积并)
- HDU 4052 Adding New Machine(矩形面积并)
- HDU 4052 Adding New Machine(矩形面积并)
- hdu 4052 Adding New Machine 扫描线求矩形面积并
- hdu 4052 Adding New Machine,set
- HDU 4052 Adding New Machine (线段树)
- hdu 4052 Adding New Machine(线段树+扫描线)
- hdu-4052/ LA 5694-Adding New Machine(线段树矩形面积并)
- [线段树]HDU 4052 Adding New Machine
- zoj 3540 Adding New Machine
- zoj 3540 Adding New Machine(map+离散化)
- 【HDU4052】【ZOJ3540】Adding New Machine 线段树+扫描线
- Adding a new disk to a VMWare Virtual Machine in Linux
- UVALive 5694 Adding New Machine
- UVA 1492 - Adding New Machine(线段树)
- uva 1492 - Adding New Machine(线段树)
- 在VMWare 中增加磁盘(Adding a new disk to a VMware Virtual Machine in Linux)
- zoj 3540 Adding New Machine