题意:给n个矩形,求他们的总周长
跟面积并类似,唯一不同的是每个节点还需要保存它所在的区间有多少个不连续的区间,则宽的个数就为该个数的二倍。具体做法是用两个数组分别保存区间的左端点和右端点是否被覆盖。如果当前区间整个被覆盖,那么左右端点都被覆盖,宽的个数就为2.如果该区间没有被覆盖,那么宽的个数就为左子结点宽的个数+右子节点宽的个数-左右子区间是否连起来,如果连起来了,那么宽的个数应该减2,画个图很容易理解。判断左右子区间是否连起来就看左子区间的右端点和右子区间的左端点是否都被覆盖,如果都被覆盖就说明左右子区间被连起来了。
代码:
#include <cstdlib>
#include <cctype>
#include <cstring>
#include <cstdio>
#include <cmath>
#include<climits>
#include <algorithm>
#include <vector>
#include <string>
#include <iostream>
#include <sstream>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <fstream>
#include <numeric>
#include <iomanip>
#include <bitset>
#include <list>
#include <stdexcept>
#include <functional>
#include <utility>
#include <ctime>
using namespace std;
#define PB push_back
#define MP make_pair
#define REP(i,x,n) for(int i=x;i<(n);++i)
#define FOR(i,l,h) for(int i=(l);i<=(h);++i)
#define FORD(i,h,l) for(int i=(h);i>=(l);--i)
#define SZ(X) ((int)(X).size())
#define ALL(X) (X).begin(), (X).end()
#define RI(X) scanf("%d", &(X))
#define RII(X, Y) scanf("%d%d", &(X), &(Y))
#define RIII(X, Y, Z) scanf("%d%d%d", &(X), &(Y), &(Z))
#define DRI(X) int (X); scanf("%d", &X)
#define DRII(X, Y) int X, Y; scanf("%d%d", &X, &Y)
#define DRIII(X, Y, Z) int X, Y, Z; scanf("%d%d%d", &X, &Y, &Z)
#define OI(X) printf("%d",X);
#define RS(X) scanf("%s", (X))
#define MS0(X) memset((X), 0, sizeof((X)))
#define MS1(X) memset((X), -1, sizeof((X)))
#define LEN(X) strlen(X)
#define F first
#define S second
#define Swap(a, b) (a ^= b, b ^= a, a ^= b)
#define Dpoint strcut node{int x,y}
#define cmpd int cmp(const int &a,const int &b){return a>b;}
/*#ifdef HOME
freopen("in.txt","r",stdin);
#endif*/
const int MOD = 1e9+7;
typedef vector<int> VI;
typedef vector<string> VS;
typedef vector<double> VD;
typedef long long LL;
typedef pair<int,int> PII;
//#define HOME
int Scan()
{
int res = 0, ch, flag = 0;
if((ch = getchar()) == '-') //判断正负
flag = 1;
else if(ch >= '0' && ch <= '9') //得到完整的数
res = ch - '0';
while((ch = getchar()) >= '0' && ch <= '9' )
res = res * 10 + ch - '0';
return flag ? -res : res;
}
/*----------------PLEASE-----DO-----NOT-----HACK-----ME--------------------*/
#define maxn 20050
struct seg
{
int l, r, h, s;
seg(){}
seg(int _l,int _r,int _h,int _s):l(_l),r(_r),h(_h),s(_s){}
bool operator<(const seg&cmp)const
{
if(cmp.h==h)
return s>cmp.s;
else
return h<cmp.h;
}
};
seg ss[maxn<<2];
int lp[maxn<<2];
int rp[maxn<<2];
int cnt[maxn<<2];
int sum[maxn<<2];
int numseg[maxn<<2];
void build(int l,int r,int rt)
{
lp[rt]=rp[rt]=cnt[rt]=sum[rt]=numseg[rt]=0;
int m=(l+r)>>1;
if(l==r)
return;
build(l,m,rt<<1);
build(m+1,r,rt<<1|1);
}
void pushup(int rt,int l,int r)
{
if(cnt[rt])
{sum[rt]=r-l+1;
lp[rt]=1;
rp[rt]=1;
numseg[rt]=2;
}
else
if(l==r)
{
sum[rt]=lp[rt]=rp[rt]=numseg[rt]=0;
}
else
{
lp[rt]=lp[rt<<1];
rp[rt]=rp[rt<<1|1];
sum[rt]=sum[rt<<1]+sum[rt<<1|1];
numseg[rt]=numseg[rt<<1]+numseg[rt<<1|1];
if(lp[rt<<1|1]&&rp[rt<<1])
numseg[rt]-=2;
}
}
void update(int L,int R,int l,int r,int rt,int c)
{
if(L<=l&&r<=R)
{
cnt[rt]+=c;
pushup(rt,l,r);
return;
}
int m=(l+r)>>1;
if(L<=m)
update(L,R,l,m,rt<<1,c);
if(R>m)
update(L,R,m+1,r,rt<<1|1,c);
pushup(rt,l,r);
}
int main()
{
int n;
while(RI(n)!=EOF)
{int Min=10001;
int Max=-10001;
int m=0;
for(int i=0;i<n;i++)
{
int a,b,c,d;
scanf("%d%d%d%d",&a,&b,&c,&d);
ss[m++]=seg(a,c,b,1);
ss[m++]=seg(a,c,d,-1);
if(a<Min)
Min=a;
if(c>Max)
Max=c;
}
build(-10000,10000,1);
sort(ss,ss+m);
int ans=0;
int last=0;
REP(i,0,m)
{
update(ss[i].l,ss[i].r-1,-10000,10000,1,ss[i].s);
ans+=numseg[1]*(ss[i+1].h-ss[i].h);
ans+=abs(sum[1]-last);
last=sum[1];
}
printf("%d\n",ans);
}
return 0;
}