HDU 1542 POJ 1151 Atlantis(线段树+扫描线)

2024-08-24 21:58

本文主要是介绍HDU 1542 POJ 1151 Atlantis(线段树+扫描线),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目地址:HDOJ地址:HDU 1542  POJ 地址:POJ 1151

第一发扫描线。。费了好大一番功夫。。构思用了半天。。写出来调试成功用了半天。。。真是弱渣。。

所谓扫描线就是从上往下或从下往上扫描,每到一个边,就进行增或删的处理。最后出来的值就是总的面积。对于求面积并的问题,可以参考这篇博客(博客地址),讲的不错。

具体实现过程是用lazy标记此时的边数量,如果大于0,说明这个地方有边,等于0则说明没有边。然后进行更新操作。最后直接使用根节点的值就行。

代码如下:

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <stdlib.h>
#include <math.h>
#include <ctype.h>
#include <queue>
#include <map>
#include <set>
#include <algorithm>using namespace std;
#define lson l, mid, rt<<1
#define rson mid+1, r, rt<<1|1
int lazy[10000], cnt;
double sum[10000], c[10000];
struct node
{double l, r, h;int f;
} edge[1000];
int cmp(node x, node y)
{return x.h<y.h;
}
void add(double l, double r, double h, int f)
{edge[cnt].l=l;edge[cnt].r=r;edge[cnt].h=h;edge[cnt++].f=f;
}
void PushUp(int rt)
{sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void update(int ll, int rr, int x, int l, int r, int rt)
{if(ll<=l&&rr>=r){lazy[rt]+=x;if(lazy[rt]){sum[rt]=c[r+1]-c[l];}else if(l!=r){PushUp(rt);}else if(l==r)sum[rt]=0;return ;}int mid=l+r>>1;if(ll<=mid) update(ll,rr,x,lson);if(rr>mid) update(ll,rr,x,rson);if(lazy[rt]==0){PushUp(rt);}else{sum[rt]=c[r+1]-c[l];}
}
int erfen(double x, int high)
{int low=0, mid;while(low<=high){mid=low+high>>1;if(c[mid]==x)return mid;else if(c[mid]>x){high=mid-1;}elselow=mid+1;}return -1;
}
int main()
{int n, i, j, num=0, k;double x1, x2, y1, y2, ans;while(scanf("%d",&n)!=EOF&&n){num++;memset(sum,0,sizeof(sum));memset(lazy,0,sizeof(lazy));ans=0;cnt=0;k=0;for(i=0; i<n; i++){scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);c[k++]=x1;c[k++]=x2;add(x1,x2,y1,1);add(x1,x2,y2,-1);}sort(edge,edge+2*n,cmp);sort(c,c+k);for(i=0; i<2*n-1; i++){int l=erfen(edge[i].l,2*n);int r=erfen(edge[i].r,2*n);update(l,r-1,edge[i].f,0,2*n-1,1);ans+=sum[1]*(edge[i+1].h-edge[i].h);//printf("%.2lf   --   %.2f   %.2f\n",ans,sum[1],edge[i+1].h-edge[i].h);}printf("Test case #%d\nTotal explored area: %.2lf\n\n",num,ans);}return 0;
}


这篇关于HDU 1542 POJ 1151 Atlantis(线段树+扫描线)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1103714

相关文章

poj3468(线段树成段更新模板题)

题意:包括两个操作:1、将[a.b]上的数字加上v;2、查询区间[a,b]上的和 下面的介绍是下解题思路: 首先介绍  lazy-tag思想:用一个变量记录每一个线段树节点的变化值,当这部分线段的一致性被破坏我们就将这个变化值传递给子区间,大大增加了线段树的效率。 比如现在需要对[a,b]区间值进行加c操作,那么就从根节点[1,n]开始调用update函数进行操作,如果刚好执行到一个子节点,

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1689(线段树成段更新)

两种操作:1、set区间[a,b]上数字为v;2、查询[ 1 , n ]上的sum 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdl

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2093 考试排名(sscanf)

模拟题。 直接从教程里拉解析。 因为表格里的数据格式不统一。有时候有"()",有时候又没有。而它也不会给我们提示。 这种情况下,就只能它它们统一看作字符串来处理了。现在就请出我们的主角sscanf()! sscanf 语法: #include int sscanf( const char *buffer, const char *format, ... ); 函数sscanf()和

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n