poj1177——Picture

线段树+离散化+扫描：

离散化的概念可以参考99年国家集训队陈宏的论文。

重点理解两个公式：[1]   ans+= st[ 1 ]. segnum * 2 *( yline[ i+1 ]. x- yline[ i]. x);

                           [2]    ans+= abs ( st[ 1 ]. sum- lsum);

首先，读入数据，存入yline[]数组中，并将y坐标存入INDEX[]数组中，其中，yline[]的结构应该为：

struct line { int st,ed;//记录竖段的上下坐标 int x;//记录其x坐标 bool inout;//记录左边还是右边 }yline[max];  ；

接下来，对yline[]数组以x为关键字进行排序，INDEX[]数组同样得排序，并对INDEX[]数组进行去重操作，得到剩下的不同的y值为cnt个；

现在，到了建树阶段：这里建树，是以[0,cnt-1]为起始点建树的，也就是说st[1].l=0,st[1].r=cnt-1;因为这样建的树的话，每个节点的l,r值是和INDEX[]数组里面的y值相对应的，而且，用的是超元线段，这样，可以节省很多空间和操作的时间，当然，建完树，应该顺便初始化的，st[max]的结构是：

struct stnode { int l,r;//当前节点的左右区间范围 int len;//当前节点的区间长度 int segnum;//当前节点里不相交的段数 int cover;//当前节点区间被覆盖的次数 int sum;//当前节点被覆盖的总长度 bool lcover,rcover;//当前节点的左右端点是否被覆盖 };  

现在，从0开始，对yline[]数组进行扫描，扫描到左边，进行Insert操作，如果是右边，则进行Delete操作，并且，每次都对数进行更新，即getlen ( index)和 getsegnum ( index)两个操作 ；

对于ans+= abs ( st[ 1 ]. sum- lsum)，这个等式，应该不难理解，每次加上|线段树中所覆盖的总长度-上次还剩的长度|即更新的长度；

而 ans+= st[ 1 ]. segnum * 2 *( yline[ i+1 ]. x- yline[ i]. x)中， 难理解的是 st[ 1 ]. segnum 。对于这个的值，就是该区间中，有多少的不相交的段数。而求这个值的过程中，就用的了cover、st[i].lcover、st[i].rcover，而且，对轴进行了离散化。自己理解下，不是很好明白，但你必须弄懂！

现在贴上代码：

#include<cstdio> #include<iostream> #include<string> #include<algorithm> using namespace std; #define max 10005 int n,cnt=0; int INDEX[max]; struct stnode { int l,r; int len; int segnum; int cover; int sum; bool lcover,rcover; stnode() { l=r=0; len=segnum=cover=sum=0; lcover=rcover=false; } }; stnode st[max*4]; struct line { int st,ed; int x; bool inout; bool operator <(line o)const { return x<o.x; } }; line yline[max]; void build(int l,int r,int index) { st[index].l=l; st[index].r=r; st[index].len =INDEX[r]-INDEX[l]; st[index].cover =0; st[index].segnum =0; st[index].sum =0; st[index].lcover =st[index].rcover =false; if(r-l>1) { int mid=(l+r)>>1; build(l,mid,index*2); build(mid,r,index*2+1); } } int getindex(int i) { return lower_bound(INDEX,INDEX+cnt,i)-INDEX;//lower_bound函数返回一个元素在容器中的迭代器，数组可以看成特殊的容器 //，所以这里返回的迭代器就是指针 } void getlen(int index) { if(st[index].cover >0) st[index].sum =st[index].len ; else if(st[index].r -st[index].l>1) st[index].sum=st[index*2].sum +st[index*2+1].sum ; else st[index].sum =0; } void getsegnum(int index) { if(st[index].cover>0) { st[index].lcover=st[index].rcover=true; st[index].segnum =1; } else if(st[index].r-st[index].l>1) { st[index].lcover=st[index*2].lcover; st[index].rcover=st[index*2+1].rcover; st[index].segnum=st[index*2].segnum+st[index*2+1].segnum-st[index*2].rcover*st[index*2+1].lcover; } else { st[index].lcover=st[index].rcover=false; st[index].segnum=0; } } void Insert(int l,int r,int index) { if(st[index].l==l&&st[index].r ==r) st[index].cover ++; else { int mid=(st[index].l +st[index].r )>>1; if(l>=mid) Insert(l,r,index*2+1); else if(r<=mid) Insert(l,r,index*2); else { Insert(l,mid,index*2); Insert(mid,r,index*2+1); } } getlen(index); getsegnum(index); } void Delete(int l,int r,int index) { if(st[index].l==l&&st[index].r==r) st[index].cover --; else { int mid=(st[index].l +st[index].r )>>1; if(r<=mid) Delete(l,r,index*2); else if(l>=mid) Delete(l,r,index*2+1); else { Delete(l,mid,index*2); Delete(mid,r,index*2+1); } } getlen(index); getsegnum(index); } int main() { while(scanf("%d",&n)!=EOF) { int x1,x2,y2,y1; int i,j,k; for(i=0;i<n;i++) { scanf("%d%d%d%d",&x1,&y1,&x2,&y2); yline[i*2].x=x1;yline[i*2].st =y1;yline[i*2].ed =y2; yline[i*2].inout =true; yline[i*2+1].x =x2;yline[i*2+1].st =y1;yline[i*2+1].ed =y2; yline[i*2+1].inout =false; INDEX[i*2]=y1; INDEX[i*2+1]=y2; } sort(INDEX,INDEX+n*2); sort(yline,yline+n*2); cnt=0; for(i=1;i<n*2;i++) { if(INDEX[i]!=INDEX[i-1]) INDEX[cnt++]=INDEX[i-1]; } INDEX[cnt++]=INDEX[2*n-1]; build(0,cnt-1,1); int ans=0; int lsum=0; for(i=0;i<2*n-1;i++) { if(yline[i].inout) Insert(getindex(yline[i].st),getindex(yline[i].ed),1); else Delete(getindex(yline[i].st),getindex(yline[i].ed),1); ans+=st[1].segnum *2*(yline[i+1].x-yline[i].x); ans+=abs(st[1].sum-lsum); lsum=st[1].sum ; } Delete(getindex(yline[2*n-1].st ),getindex(yline[2*n-1].ed ),1); ans+=abs(st[1].sum -lsum); printf("%d/n",ans); } return 0; }