ZOJ3273 POJ3832 HDU3265 Posters

线段树求矩形并。

由于数据范围不大，而都是整数，所以不需要离散化。算法思想是从左往右扫描，node.length保存的是整条y轴可用的长度；LINE中的flag用于记录矩形的左边和右边，左边为1，进入node后就代表LINE的在y轴的这段可以用，右边为-1，进入node后表示LINE的在y轴的这段不能用了。这个要画图才能理解。

/*******************************************************************************
 # Author : Neo Fung
 # Email : [email protected]
 # Last modified: 2012-01-18 20:11
 # Filename: HDU1689 Just a Hook.cpp
 # Description : 线段树求矩形并
 ******************************************************************************/
#ifdef _MSC_VER
#define DEBUG
#endif

#include <fstream>
#include <stdio.h>
#include <iostream>
#include <string.h>
#include <string>
#include <limits.h>
#include <algorithm>
#include <math.h>
#include <numeric>
#include <functional>
#include <ctype.h>
#define L(x) (x<<1)
#define R(x) (x<<1|1)
#define MAX 200010
using namespace std;

struct NODE
{
	int l,r,length,cover;
}node[MAX*4];

struct LINE
{
	int x,y1,y2,flag;
}line[MAX*2];

int y_coord[MAX*2];

void add_line(const int &x1,const int &x2,const int &y1,const int &y2,int &cnt)
{
	line[cnt].x=x1;
	line[cnt].y1=y1;
	line[cnt].y2=y2;
	line[cnt++].flag=1;
	line[cnt].x=x2;
	line[cnt].y1=y1;
	line[cnt].y2=y2;
	line[cnt++].flag=-1;
}

bool inline cmp(const LINE &lhs,const LINE &rhs)
{
	return lhs.x<rhs.x;
}

void init()
{
	memset(node,0,sizeof(node));
}

void build(const int &t,const int &l,const int &r)
{
	node[t].l=l;
	node[t].r=r;
	if(l==r-1) 
		return;
	int mid=(l+r)>>1;
	build(L(t),l,mid);
	build(R(t),mid,r);
}

void len(const int &t)
{
	if(node[t].cover>0)
		node[t].length = y_coord[node[t].r]-y_coord[node[t].l];
	else if(node[t].l==node[t].r-1)
		node[t].length=0;
	else
		node[t].length = node[R(t)].length+node[L(t)].length;
}

void update(const int &t,const LINE &lne)
{
	if(y_coord[node[t].l]>=lne.y1 && y_coord[node[t].r]<=lne.y2)
	{
		node[t].cover+=lne.flag;
		len(t);
		return;
	}
	if(node[t].l==node[t].r-1) return;
	int mid = (node[t].l+node[t].r)>>1;
	if(lne.y1<=y_coord[mid])
		update(L(t),lne);
	if(lne.y2>y_coord[mid])
		update(R(t),lne);
	len(t);
}

__int64 solve(const int &n,const int &cnt)
{
	init();
	build(1,0,cnt-1);
	update(1,line[0]);
	__int64 sum=0ll;
	for(int i=1;i<n;++i)
	{
		sum+=(line[i].x-line[i-1].x)*1ll*node[1].length;
		update(1,line[i]);
	}
	return sum;
}

int main(void)
{
#ifdef DEBUG  
  freopen("../stdin.txt","r",stdin);
  freopen("../stdout.txt","w",stdout); 
#endif  

  int n;
  int x1,y1,x2,y2,x3,y3,x4,y4;

  while(scanf("%d",&n) && n)
  {
	  int cnt=0,y_cnt=0;
	  for(int i=0;i<n;++i)
	  {
		  scanf("%d%d%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&x3,&y3,&x4,&y4);
		  add_line(x1,x3,y1,y2,cnt);
		  add_line(x3,x4,y1,y3,cnt);
		  add_line(x3,x4,y4,y2,cnt);
		  add_line(x4,x2,y1,y2,cnt);
		  y_coord[y_cnt++]=y1;
		  y_coord[y_cnt++]=y2;
		  y_coord[y_cnt++]=y3;
		  y_coord[y_cnt++]=y4;
	  }
	  sort(line,line+cnt,cmp);
	  sort(y_coord,y_coord+y_cnt);
	  y_cnt=unique(y_coord,y_coord+y_cnt)-y_coord;
	  printf("%I64d\n",solve(cnt,y_cnt));
  }

  return 0;
}