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 10005int 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);elseDelete(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;}