题目大意:

有一些屋顶,相当于一些线段(不想交)。

问每一条线段能够接到多少水,相对较低的屋顶能够接到高屋顶留下的水(如题图所看到的)。因为y1!=y2,所以保证屋顶是斜的。



解题思路:

扫描线,由于对于同一个x最多有25条线段。所以不须要线段树更新。

在扫描线的过程中构造出线段与线段之间的关系。好在最后计算每一个屋顶能够接多少水。



以下是代码:

#include <set>
#include <map>
#include <queue>
#include <math.h>
#include <vector>
#include <string>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm> #define eps 1e-8
#define pi acos(-1.0)
#define inf 107374182
#define inf64 1152921504606846976
#define lc l,m,tr<<1
#define rc m + 1,r,tr<<1|1
#define iabs(x) ((x) > 0 ? (x) : -(x))
#define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE))
#define clearall(A, X) memset(A, X, sizeof(A))
#define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE))
#define memcopyall(A, X) memcpy(A , X ,sizeof(X))
#define max( x, y ) ( ((x) > (y)) ? (x) : (y) )
#define min( x, y ) ( ((x) < (y)) ? (x) : (y) ) using namespace std; struct node
{
int x1,x2,y1,y2;
int num,next,w,deg;
double k;
} segment[400005]; int n,point[400000<<1],ans[400005],q[400005],a[100]; bool cmp(node a,node b)
{
return a.x1<b.x1;
} int main()
{
int last=0,now,tail=0,k,cnt,u,v,pointcnt,levelcnt;
while(scanf("%d",&n)!=EOF)
{
pointcnt=0;
for(int i=0; i<n; i++)
{
scanf("%d%d%d%d",&segment[i].x1,&segment[i].y1,&segment[i].x2,&segment[i].y2);
segment[i].num=i;
segment[i].next=-1;
segment[i].w=0;
segment[i].deg=0;
segment[i].k=(1.0*segment[i].y2-segment[i].y1)/(segment[i].x2-segment[i].x1);
point[pointcnt++]=segment[i].x1;
point[pointcnt++]=segment[i].x2;
} sort(point,point+pointcnt);
pointcnt=unique(point,point+pointcnt)-point; sort(segment,segment+n,cmp); last=0;
tail=0;
levelcnt=0; for(int i = 0; i < pointcnt ; last=now,i++)
{
now=point[i];
if(levelcnt)
{
segment[a[levelcnt-1]].w+=now-last;
}
while(segment[tail].x1==now&&tail<n)
{
k=-1;
for(int j=0; j<levelcnt; j++)
{
if( (segment[a[j]].y1 + segment[a[j]].k * (now - segment[a[j]].x1))<segment[tail].y1)k=j;
}
for(int j=levelcnt-1; j>k; j--)
{
a[j+1]=a[j];
}
a[k+1]=tail;
levelcnt++;
tail++;
}
for(int j=1; j<levelcnt; j++)
{
if((segment[a[j]].x1==now&&segment[a[j]].y1<segment[a[j]].y2)||(segment[a[j]].x2==now&&segment[a[j]].y1>segment[a[j]].y2))
{
segment[a[j]].next=a[j-1];
segment[a[j-1]].deg++;
}
}
cnt=0;
for(int j=0; j<levelcnt; j++)
{
if(segment[a[j]].x2!=now)a[cnt++]=a[j];
}
levelcnt=cnt;
}
queue <int >q;
for(int i=0; i<n; i++)
{
if(!segment[i].deg)q.push(i);
}
while(!q.empty())
{
u=q.front();
q.pop();
v=segment[u].next;
if(v==-1)continue;
segment[v].deg--;
segment[v].w+=segment[u].w;
if(!segment[v].deg)q.push(v);
}
for(int i=0; i<n; i++)
{
ans[segment[i].num]=segment[i].w;
} for(int i=0; i<n; i++)
{
printf("%d\n",ans[i]);
} }
return 0;
}



最新文章

  1. Kinect开发资源汇总
  2. Python 学习小结
  3. 修改安卓串口蓝牙app问题记录
  4. 1137. Bus Routes(dfs)
  5. OM Price Lists
  6. Bzoj-2820 YY的GCD Mobius反演,分块
  7. BZOJ_1620_[Usaco2008_Nov]_Time_Management_时间管理_(二分+贪心)
  8. uva 10916 Factstone Benchmark(对数函数的活用)
  9. python进度条代码
  10. android中自定义shape
  11. javascript核心概念之——数组
  12. Ajax跨域 CROS处理
  13. Java中空串和null串的区别
  14. 隐写术之steghide的使用
  15. MySQL常用数值函数
  16. NGINX 502错误排查(转)
  17. spring学习(三) ———— spring事务操作
  18. 在Sublime中配置JsFormat
  19. ABP框架系列之二十:(Dependency-Injection-依赖注入)
  20. mysql case when 判断null

热门文章

  1. chrome打开网址但是没有地址栏
  2. 【WaaCaa】一款开源科学作图/数据可视化工具 —— 诞生篇
  3. python fuzzy c-means demo
  4. ROS-导航功能-RVIZ
  5. Core篇——初探Core的认证,授权机制
  6. JS实现掷骰子
  7. 关于数据未渲染完,要获取document高度问题——ajax全局事件
  8. Android 自定义简单控件--星级评价
  9. c#可自定义码表的base64加密解密算法类
  10. LeetCode(11)Container With Most Water