线段树扫描线~~哥的非主流写法。。。hdu1542

最新推荐文章于 2021-01-19 22:05:01 发布

最新推荐文章于 2021-01-19 22:05:01 发布 · 86 阅读

本文介绍了一种利用集合数据结构简化区间更新与查询算法的方法。通过使用 set 和自定义结构体，解决了区间更新与查询的问题，并实现了高效的面积求解算法。文中详细展示了算法的实现过程，包括数据结构的设计、区间更新的处理方式等。

用了个set，所以不用写蛋疼的二分了，哈哈哈~~多方便啊！！

一开始output 没搞空行，PE了两次才ac。。T T

#include<iostream>

#include<vector>
#include<string>
#include<algorithm>
#include<cstdio>
#include<set>
#include<map>
#include<iomanip>
using namespace std;
const int maxn=2222;
struct pp
{
double l;
double r;
double h;
int d;
pp(){}
pp (double a,double b,double c,int dd) : l(a) , r(b) , h(c) , d(dd) {};
bool operator < (const pp &cmp) const
{
return h < cmp.h;
}
}px;
struct dd
{
double x;
int b;
dd(){}
dd(double a)
{
x=a;
b=0;
}
dd(double a,int bb)
{
x=a;
b=bb;
}
bool operator < (const dd &cmp) const
{
return x < cmp.x;
}
bool operator == (const dd &cmp) const
{
return x==cmp.x;
}
}dx;
double sum[maxn<<2];
int add[maxn<<2];
int lend[maxn<<2];
int rend[maxn<<2];
vector<double>vx;
vector<double>v;
set<dd> sx;
set<dd>::iterator si;
vector<pp>vp;
double a,b,c,d,t,total,dps;
int n,temp,tt,lstart,rstart;
int find(double x)
{
si=sx.find(x);
return si->b;
}
void push_up(int root)
{
sum[root]=sum[root<<1]+sum[root<<1|1];
return ;
}
void build(int l,int r,int root=1)
{
lend[root]=l;
rend[root]=r;
int mid=(r+l)>>1;
if(r-l==1)
{
sum[root]=v[r]-v[l];
return ;
}
else
{
build(l,mid,root<<1);
build(mid,r,root<<1|1);
push_up(root);
return ;
}
}
void update(int l,int r,int d,int root=1)
{
if(l<=lend[root] && r>=rend[root])
{
add[root]+=d;
return ;
}
else
{
if(l<rend[root<<1] && r>lend[root<<1])
{
update(l,r,d,root<<1);
}
if(r>lend[root<<1|1] && l<rend[root<<1|1])
{
update(l,r,d,root<<1|1);
}
return ;
}
}
void query(int root=1)
{
if(add[root]>0)
{
dps+=sum[root];
return ;
}
else
{
if(rend[root<<1]-lend[root<<1]>=1)
{
query(root<<1);
}
if(rend[root<<1|1]-lend[root<<1|1]>=1)
{
query(root<<1|1);
}
return ;
}
}
int main()
{
tt=1;
while(cin>>n && n)
{
memset(sum,0,sizeof(sum));
memset(add,0,sizeof(add));
total=0.0;
dps=0.0;
sx.clear();
vx.clear();
vp.clear();
v.clear();
v.push_back(0.0);
for(int i=1;i<=n;i++)
{
scanf("%lf%lf%lf%lf",&a,&b,&c,&d);
vx.push_back(a);
vx.push_back(c);
px=pp(a,c,b,1);
vp.push_back(px);
px=pp(a,c,d,-1);
vp.push_back(px);
}
sort(vp.begin(),vp.end());
sort(vx.begin(),vx.end());
temp=1;
t=-0.1;
for(int i=0;i<vx.size();i++)
{
if(vx[i]!=t)
{
dx=dd(vx[i],temp);
temp++;
t=vx[i];
v.push_back(vx[i]);
sx.insert(dx);
}
else
{
continue;
}
}
build(1,v.size()-1);
t=vp[0].h;
for(int i=0;i<vp.size();)
{

if(vp[i].h==t)
{
update(find(vp[i].l),find(vp[i].r),vp[i].d);
i++;
}
else
{
dps=0;
query();
total+=dps*(vp[i].h - t);
t=vp[i].h;
}
}
cout<<"Test case #"<<tt++<<endl;
cout<<"Total explored area: ";
//<<setiosflags(ios::fixed)<<setprecision(2)<<total<<endl<<endl;
printf("%.2lf\n\n",total);
}
return 0;
}
最后。。不知道为什么时间不是0ms。。╮(╯▽╰)╭~~