N - Picture
扫描线求周长。
刚学完求面积 然后不会这一题。。
求周长其实和面积差不多POJ 1177 Picture (线段树+离散化+扫描线) 详解这篇博客讲的不错
通常有2中方法。
1.横向扫一遍 求竖线和,然后纵向求横线和
2.只扫一遍 记录每段扫描线中间的缺口数量,没缺一个 横向边要多2个。
懒得写题解了 看看刚才那篇博客+我的代码应该就差不多了。
//#incldue<bits/stdc++.h>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
#include<map>
using namespace std;
vector<double>v;
const int M = 100000+7;
struct node
{
int x,y1,y2;
int k;
}rem[M];
struct nodet
{
int cov,len,lp,rp,num;
}st[M>>2];
bool cmp(node a,node b)
{
return a.x<b.x;
}
int getid(double x)
{
return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}
void pushup(int l,int r,int rt)
{
if(st[rt].cov>0)
{
st[rt].lp=1;
st[rt].rp=1;
st[rt].num=1;
st[rt].len=v[r]-v[l-1];//vector下标从0开始
}
else if(l==r)
{
st[rt].len=st[rt].lp=st[rt].rp=0;
st[rt].num=0;
}
else
{
st[rt].len=st[rt<<1].len+st[rt<<1|1].len;
//如果当前段未完全覆盖 中间缺口数等于2个子树相加 减去子树中点是否连接。
st[rt].num=st[rt<<1].num+st[rt<<1|1].num-st[rt<<1].rp*st[rt<<1|1].lp;
st[rt].lp=st[rt<<1].lp;
st[rt].rp=st[rt<<1|1].rp;
}
}
void update(int l,int r,int rt,int x,int y,int k)
{
if(x<=l&&r<=y)
{
st[rt].cov+=k;
pushup(l,r,rt);
return ;
}
int mid=(l+r)>>1;
if(x<=mid)update(l,mid,rt<<1,x,y,k);
if(y>mid)update(mid+1,r,rt<<1|1,x,y,k);
pushup(l,r,rt);
}
int main()
{
int n,cnt=0,x1,x2,y1,y2;
int ans=0;
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d%d%d%d",&x1,&y1,&x2,&y2),v.push_back(y1),v.push_back(y2),
rem[++cnt].x=x1,rem[cnt].y1=y1,rem[cnt].y2=y2,rem[cnt].k=1,
rem[++cnt].x=x2,rem[cnt].y1=y1,rem[cnt].y2=y2,rem[cnt].k=-1;
sort(v.begin(),v.end()),v.erase(unique(v.begin(),v.end()),v.end());//离散化去重,很方便
sort(rem+1,rem+1+cnt,cmp);
int pre=0;
for(int i=1;i<=cnt;i++)
{
update(1,cnt,1,getid(rem[i].y1),getid(rem[i].y2)-1,rem[i].k);
//cout<<abs(st[1].len-pre)<<" "<< st[1].num*2*(rem[i+1].x-rem[i].x) <<" "<<st[1].num<<endl;
ans+=abs(st[1].len-pre);//横向扫描 求竖线长度和
pre=st[1].len;
if(i<cnt)//注意这里 横线像扫描线求面积一样只加到倒数第二个竖线,而前面的竖线每次都要求差值相加
ans+=st[1].num*2*(rem[i+1].x-rem[i].x);//横线长度和
}
printf("%d\n",ans);
return 0;
}