Description

见上帝动了恻隐之心,天后也想显示一下慈悲之怀,随即从口袋中取出一块魔术方巾,让身边的美神维纳斯拿到后堂的屏风上去试试,屏风是正方形的,高和宽方向上各划有m条鱼屏风的边平行的直线,平行直线间的距离为1厘米。这2m条直线共有m*m个交点,在某些交点上镶嵌着宝石。如果魔术方巾的边与屏风的边平行且魔术方巾触碰到屏风上镶嵌着的宝石,就将与这些宝石等值的金银送给人们。维纳斯想让魔术方巾触碰到的宝石的价值最多,可要在短短的1秒钟之内解决问题,也感到力不从心,你能帮帮她吗?

Input

输入文件gem.in的第一行有三个正整数m,n,k,数与数之间用一个空格分隔。其中m为屏风在高和宽方向上被划分出的直线数。魔术方巾为正方形的,它的边长为k厘米。N为屏风上宝石的个数。
接下来的n行,每行三个正整数,依次表示宝石所在直线的行号、列号、宝石的价值,数与数之间用一个空格分隔。

Output

输出文件gem.out只有一个正整数,为魔术方巾触碰到的宝石的最大价值总数。

Sample Input

10    4    4
1 1 9
2 3 5
6 2 12
4 5 6

Sample Output

23

Hint

【输入样例2】
10 4 3
1 1 9
2 3 5
6 2 12
4 5 6
【输出样例2】
18

【数据限制】
30%的数据,1≤m≤500,1≤n≤10000,1≤k≤100;
60%的数据,1≤m≤3000,1≤n≤10000,1≤k≤1000;
100%的数据,1≤m≤50000,1≤n≤50000,1≤k≤10000;

题目大意就是在一个正方形内给定一些点,用一个小正方形去覆盖,由于每个点的分数不同,因此肯定有所谓最大得分,算出最大得分就是我们的任务

我们放正方形,有一个点放在正方形的一角上是最优的,这里我们采用左下角

把输入的点按照x离散化,再按照y进行排序,考虑从下到上扫描线。线怎么来呢?我们把每个点置于小矩形的左下角,每个点就拥有了一个对应的范围,对于每一个范围我们把上边界所在直线进行扫描。

想象各点范围相交的地方肯定是这些点的分数都可以取到的。

在向上扫描的过程中,我们相当于是将一个长为m的矩形(离散化之后就是cnt),宽为k的矩形不断想上移动,然后去其中长为k的得分最大的一段。具体实现就是在扫描线的线段树上维护区间最大值就可以了。

下面附上代码:

#include<algorithm>
#include<cstdio>
using namespace std; const int maxn=5e4+;
int m,n,k,tot,cnt;
int left[maxn<<],right[maxn<<];
struct STAR
{
int id,x,y,l;
bool fl;
}a[maxn<<],b[maxn<<];
struct node
{
int l,r,m,lazy;
}t[maxn<<];
inline int read()
{
char ch=getchar();
int s=,f=;
while (!(ch>=''&&ch<='')) {if (ch=='-') f=-;ch=getchar();}
while (ch>=''&&ch<='') {s=(s<<)+(s<<)+ch-'';ch=getchar();}
return s*f;
}
bool cmp(STAR a,STAR b) {return a.x<b.x;}
bool CMP(STAR a,STAR b) {if (a.y==b.y) return a.l>b.l;else return a.y<b.y;}
void build(int u,int l,int r)
{
t[u].l=l;
t[u].r=r;
t[u].lazy=t[u].m=;
if (l==r) return;
int mid=l+r>>;
build(u<<,l,mid);
build(u<<|,mid+,r);
}
void down(int u)
{
if (t[u].lazy)
{
int k=t[u].lazy;
t[u].lazy=;
t[u<<].lazy+=k;
t[u<<|].lazy+=k;
t[u<<].m+=k;
t[u<<|].m+=k;
}
}
void add(int u,int l,int r,int k)
{
if (t[u].l>=l&&t[u].r<=r) {
t[u].m+=k;
t[u].lazy+=k;
down(u);
return;
}
down(u);
int mid=t[u].l+t[u].r>>;
if (l<=mid) add(u<<,l,r,k);
if (r>mid) add(u<<|,l,r,k);
t[u].m=max(t[u<<].m,t[u<<|].m);
return;
}
int query(int u,int l,int r)
{
if (t[u].l>=l&&t[u].r<=r) return t[u].m;
down(u);
int ans=,mid=t[u].l+t[u].r>>;
if (l<=mid) ans=query(u<<,l,r);
if (r>mid) ans=max(ans,query(u<<|,l,r));
return ans;
}
int main()
{
m=read();n=read();k=read();
for (int i=;i<=n;i++)//a数组用来离散化
{
a[i].x=read();a[i].y=read();a[i].l=read();
a[i].id=i;
a[i].fl=;
b[i].y=a[i].y;//扫描线左端点
b[i].id=i;
b[i].l=a[i].l;
}
tot=n;
for (int i=;i<=n;i++)
{
a[++tot].id=i;
a[tot].x=a[i].x+k;
a[tot].fl=;
}
sort(a+,a++tot,cmp);//按x排序
for (int i=;i<=tot;i++)
{
int tt=;
if (a[i].x==a[i-].x&&i!=) tt=cnt;
else tt=++cnt;
if (!a[i].fl) left[a[i].id]=tt;
else right[a[i].id]=tt;
}
tot=n;
for (int i=;i<=n;i++)
{
b[++tot].id=b[i].id;//扫描线右端点
b[tot].l=-b[i].l;
b[tot].y=b[i].y+k;
}
sort(b+,b++tot,CMP);//按y排序
build(,,cnt);
int ans=;
for (int i=;i<=tot;i++)
{
if (b[i].l>) {
int sum=query(,left[b[i].id],right[b[i].id]);//只有在b[i]的范围内才能加上b[i].l
if (sum+b[i].l>ans) ans=b[i].l+sum;
}
add(,left[b[i].id],right[b[i].id],b[i].l);
}
printf("%d",ans);
return ;
}

jzoj1495 宝石 解题报告[扫描线]的更多相关文章

  1. ACM-ICPC 2017 Asia HongKong 解题报告

    ACM-ICPC 2017 Asia HongKong 解题报告 任意门:https://nanti.jisuanke.com/?kw=ACM-ICPC%202017%20Asia%20HongKon ...

  2. CH Round #56 - 国庆节欢乐赛解题报告

    最近CH上的比赛很多,在此会全部写出解题报告,与大家交流一下解题方法与技巧. T1 魔幻森林 描述 Cortana来到了一片魔幻森林,这片森林可以被视作一个N*M的矩阵,矩阵中的每个位置上都长着一棵树 ...

  3. 二模13day1解题报告

    二模13day1解题报告 T1.发射站(station) N个发射站,每个发射站有高度hi,发射信号强度vi,每个发射站的信号只会被左和右第一个比他高的收到.现在求收到信号最强的发射站. 我用了时间复 ...

  4. BZOJ 1051 最受欢迎的牛 解题报告

    题目直接摆在这里! 1051: [HAOI2006]受欢迎的牛 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 4438  Solved: 2353[S ...

  5. 习题:codevs 2822 爱在心中 解题报告

    这次的解题报告是有关tarjan算法的一道思维量比较大的题目(真的是原创文章,希望管理员不要再把文章移出首页). 这道题蒟蒻以前做过,但是今天由于要复习tarjan算法,于是就看到codevs分类强联 ...

  6. 习题:codevs 1035 火车停留解题报告

    本蒟蒻又来写解题报告了.这次的题目是codevs 1035 火车停留. 题目大意就是给m个火车的到达时间.停留时间和车载货物的价值,车站有n个车道,而火车停留一次车站就会从车载货物价值中获得1%的利润 ...

  7. 习题: codevs 2492 上帝造题的七分钟2 解题报告

    这道题是受到大犇MagHSK的启发我才得以想出来的,蒟蒻觉得自己的代码跟MagHSK大犇的代码完全比不上,所以这里蒟蒻就套用了MagHSK大犇的代码(大家可以关注下我的博客,友情链接就是大犇MagHS ...

  8. 习题:codevs 1519 过路费 解题报告

    今天拿了这道题目练练手,感觉自己代码能力又增强了不少: 我的思路跟别人可能不一样. 首先我们很容易就能看出,我们需要的边就是最小生成树算法kruskal算法求出来的边,其余的边都可以删掉,于是就有了这 ...

  9. NOIP2016提高组解题报告

    NOIP2016提高组解题报告 更正:NOIP day1 T2天天爱跑步 解题思路见代码. NOIP2016代码整合

随机推荐

  1. python多线程编程代码

    参考了这篇文章,写的挺好的. http://blog.csdn.net/abcjennifer/article/details/49474897 import time import threadin ...

  2. 走进windows编程的世界-----入门篇

    1   Windows编程基础 1.1Win32应用程序基本类型 1)  控制台程序 不须要完好的windows窗体,能够使用DOS窗体方式显示 2)  Win32窗体程序 包括窗体的程序,能够通过窗 ...

  3. Android与server通信的方法之中的一个(json)效率不高安全性不好

    http通信.以xml或者json为载体.相互通信数据. Android对于http的网络通信,提供了标准的java接口--httpURLConnection接口.以及apache的接口--httpc ...

  4. VMware虚拟机的CentOS7安装Nginx后本机用CentOS的IP地址无法访问

    因为CentOS7的默认防火墙改成了Firewall,不再使用iptables为默认防火墙了 所以需要使用以下命令添加80端口 firewall-cmd --zone=public --add-por ...

  5. SPOJ 694/705 后缀数组

    思路: 论文题*n Σn-i-ht[i]+1 就是结果 O(n)搞定~ //By SiriusRen #include <cstdio> #include <cstring> ...

  6. BZOJ 2190 欧拉函数

    思路: 递推出来欧拉函数 搞个前缀和 sum[n-1]*2+3就是答案 假设仪仗队是从零开始的 视线能看见的地方就是gcd(x,y)=1的地方 倒过来一样 刨掉(1,1) 就是ans*2+1 再加一下 ...

  7. 在CentOS下安装tomcat并配置环境变量(改默认端口8080为8081)

    不多说,直接上干货! 第一步:下载tomcat压缩包 http://archive.apache.org/dist/tomcat/tomcat-7/v7.0.73/bin/ 第二步:上传tomcat压 ...

  8. 使用JS方法使页面滚动到指定元素+优化+API介绍(动画)

    前言 当页面最上部有顶部菜单是,使用锚点跳转的方法很容易挡住想要呈现的内容(如下图技能两个字被挡住了一半),为避免出现这样的问题,故滚动到指定元素使用用JS的方法来实现. 目录 使用的API简介 初版 ...

  9. 微信小程序------开发测试

    一.注册小程序 注:微信小程序注册的邮箱不能被其他微信公众平台注册,未被微信开放平台注册,未被给人微信号绑定的微信号. 二.注册完小程序后,下载开发者工具 开发者工具的使用: 1.打开开发者工具:用已 ...

  10. ActiveMQ学习笔记(18)----Message高级特性(二)

    1. Blob Message 有些时候,我们需要传递Blob(Binary Large Objects)消息,在5.14之前,(5.12和5.13需要在jetty.xml中手动开启)可以按照如下的方 ...