poj 1177 Picture 線段樹+離散化+線掃描

慚愧啊，還是不怎麼懂，主要是lcover，和rcover，看了很多解題報告都沒怎麼寫清楚，看陳宏的論文，感覺自己理解到了，不過寫不出代碼，悲劇啊。看來是我的iq有問題~

下面這個人寫的還可以：

http://www.cppblog.com/abilitytao/archive/2010/07/21/120927.html

/* * File: main.cpp * Author: Mi * * Created on 2011年3月24日, 下午7:02 */ #include <cstdlib> #include <stdio.h> #include <stdio.h> #include <algorithm> #define N 10005 using namespace std; /* * */ struct tree { int l,r; int len; int seg; int cover; int sum; bool lcover,rcover; }T[N*3]; struct line { int st,ed; int x; bool inout; bool operator<(line o) const { return x<o.x; } }Y[N]; int Index[N],cnt,n; void build(int l,int r,int root) { T[root].l=l,T[root].r=r; T[root].len=Index[r]-Index[l]; T[root].seg=0; T[root].lcover=T[root].rcover=false; T[root].cover=0; T[root].sum=0; if(r>l+1) { int mid=(r+l)>>1; build(l,mid,root<<1); build(mid,r,root<<1|1); } } void update(int root) { if(T[root].cover>0) { T[root].sum=T[root].len; T[root].lcover=T[root].rcover=true; T[root].seg=1; } else if(T[root].r-T[root].l>1) { T[root].sum=T[root<<1].sum+T[root<<1|1].sum; T[root].lcover=T[root<<1].lcover; T[root].rcover=T[root<<1|1].rcover; T[root].seg=T[root<<1].seg+T[root<<1|1].seg-T[root<<1].rcover*T[root<<1|1].lcover; } else { T[root].sum=0; T[root].lcover=T[root].rcover=false; T[root].seg=0; } } void Insert(int l,int r,int root) { if(l==T[root].l&&r==T[root].r) T[root].cover++; else { int mid=(T[root].l+T[root].r)>>1; if(r<=mid) Insert(l,r,root<<1); else if(l>=mid) Insert(l,r,root<<1|1); else { Insert(l,mid,root<<1); Insert(mid,r,root<<1|1); } } update(root); } void Delete(int l,int r,int root) { if(l==T[root].l&&r==T[root].r) T[root].cover--; else { int mid=(T[root].l+T[root].r)>>1; if(r<=mid) Delete(l,r,root<<1); else if(l>=mid) Delete(l,r,root<<1|1); else { Delete(l,mid,root<<1); Delete(mid,r,root<<1|1); } } update(root); } int GetIndex(int x) { int t=lower_bound(Index,Index + cnt,x)-Index; return t; } int main(int argc, char** argv) { while(scanf("%d",&n)==1) { cnt=0; int x1,y1,x2,y2; for(int i=0;i<n;i++) { scanf("%d%d%d%d",&x1,&y1,&x2,&y2); Y[i<<1].st=Y[i<<1|1].st=y1; Y[i<<1].ed=Y[i<<1|1].ed=y2; Y[i<<1].x=x1,Y[i<<1|1].x=x2; Y[i<<1].inout=true; Y[i<<1|1].inout=false; Index[i<<1]=y1; Index[i<<1|1]=y2; } sort(Y,Y+2*n); sort(Index,Index+2*n); for(int i=1;i<2*n;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,back=0; for(int i=0;i<2*n-1;i++) { if(Y[i].inout) Insert(GetIndex(Y[i].st),GetIndex(Y[i].ed),1); else Delete(GetIndex(Y[i].st),GetIndex(Y[i].ed),1); ans+=T[1].seg*(Y[i+1].x-Y[i].x)*2; ans+=abs(T[1].sum-back); back=T[1].sum; } Delete(GetIndex(Y[2*n-1].st),GetIndex(Y[2*n-1].ed),1); ans+=abs(T[1].sum-back); printf("%d/n",ans); } return 0; }