2019杭电暑假多校训练 第六场 Snowy Smile HDU - 6638
很多题解都是简单带过,所以打算自己写一篇,顺便也加深自己理解
前置知识:线段树、线段树维护最大字段和、二维坐标离散化
题解:
1.很容易想到我们需要枚举所有子矩阵来得到一个最大子矩阵,所以我们的任务是 “枚举所有子矩阵”,
二维前缀和暴力枚举达到O(n^4), DP结合前缀和枚举也需要O(n^3),然而我们的时间需要更少.
2.可以看到坐标范围为 -10^9~10^,但是点 n<=2000个,所以我们需要先将点 离散化。
3.将点 按y轴高度排序,枚举矩阵的上下界,这将达到O(n^2)了。
4.最重点的一步,将矩阵内的点 加入线段树维护。下面解释下这一步。
首先假设我们枚举的这一个矩阵的 上界为 up ,下界为 down ,目前的矩阵的宽度就已经知道是 up-down。
所以现在我们剩下的任务就是“枚举矩阵宽度” 来达到 “枚举宽度为up-down的所有子矩阵,找出宽度为up-down的最大子矩阵”。
我们宽度已知,所有要枚举也就是长度,这样我们可以把 “矩阵压缩称为一条线”。
这时候线段树的功能就能解决这个问题了。
用线段树来维护最大字段和,其实维护的是 “宽度为up-down”的最大矩阵和。
这样我们的时间复杂度就可以达到O(n^2 logn)了。
记得每次改变下界的时候要初始化线段树,关于这个初始化代码中还有一个小技巧,差不多减少了1.5s左右的时间。
#include <bits/stdc++.h> using namespace std;
typedef double dou;
typedef long long ll;
typedef pair<int, int> pii;
typedef map<int, int> mii; #define pai acos(-1.0)
#define M 4005
#define inf 0x3f3f3f3f
#define mod 1000000007
#define IN inline
#define left k<<1
#define right k<<1|1
#define lson L, mid, left
#define rson mid + 1, R, right
#define W(a) while(a)
#define lowbit(a) a&(-a)
#define ms(a,b) memset(a,b,sizeof(a))
#define Abs(a) (a ^ (a >> 31)) - (a >> 31)
#define false_stdio ios::sync_with_stdio(false),cin.tie(0),cout.tie(0) int T, n;
ll vx[M], vy[M];
ll xlen, ylen, pos, ans;
struct Data {
int x, y, val;
bool operator <(Data& t) {
return y < t.y;
}
}node[M];
struct Data_t {
ll sum;
ll Lmax, Rmax, Max;
}tree[M << ]; IN void Updata(int L, int R, int k, int id,int add) {
if (L == R) {
tree[k].sum += (ll)add;
tree[k].Lmax = tree[k].Rmax = tree[k].Max = tree[k].sum;
return;
}
int mid = L + R >> ;
if (id <= mid)Updata(lson, id, add);
else Updata(rson, id, add); //维护最大字段和
tree[k].sum = tree[left].sum + tree[right].sum;
tree[k].Lmax = max(tree[left].Lmax, tree[left].sum + tree[right].Lmax);
tree[k].Rmax = max(tree[right].Rmax, tree[right].sum + tree[left].Rmax);
tree[k].Max = max(max(tree[left].Max, tree[right].Max), tree[left].Rmax + tree[right].Lmax);
} int main() {
false_stdio;
cin >> T;
W(T--) {
cin >> n;
for (int i = ; i <= n; i++) {
cin >> node[i].x >> node[i].y >> node[i].val;
vx[i] = node[i].x, vy[i] = node[i].y;
} //二维坐标离散化
sort(vx + , vx + n + );
sort(vy + , vy + n + );
xlen = unique(vx + , vx + n + ) - vx - ;
ylen = unique(vy + , vy + n + ) - vy - ;
for (int i = ; i <= n; i++) {
node[i].x = lower_bound(vx + , vx + xlen + , node[i].x) - vx;
node[i].y = lower_bound(vy + , vy + ylen + , node[i].y) - vy;
}
sort(node + , node + n + ); ans = ; //首先枚举下界
for (int dw = ; dw <= ylen; dw++) {
pos = ;
memset(tree, , (xlen * + ) * sizeof(Data_t));//初始化线段树,离散化完有多少个点就初始化多大 W(node[pos].y < dw && pos <= n)pos++;//直接跳过小于下界的点 for (int up = dw; up <= ylen; up++) {//枚举上界
W(pos <= n && node[pos].y <= up) {//将上界与下届之间的点加入线段树中
Updata(, xlen, , node[pos].x, node[pos].val);
pos++;
}
ans = max(ans, tree[].Max);
}
}
cout << ans << endl;
}
return ;
}
2019杭电暑假多校训练 第六场 Snowy Smile HDU - 6638的更多相关文章
- 2019牛客网暑假多校训练第四场 K —number
链接:https://ac.nowcoder.com/acm/contest/884/K来源:牛客网 题目描述 300iq loves numbers who are multiple of 300. ...
- HDU 4639 hehe 杭电2013多校联赛第四场1008题
解题报告:题目的意思是输入一个字符串,并规定,里面的“hehe”可以用"wqnmlgb"来代替,也可以不代替,问输入的这个字符串在经过相关的代替之后可以有多少种不同的形态.先打一个 ...
- 2019牛客暑假多校赛(第二场) F和H(单调栈)
F-Partition problem https://ac.nowcoder.com/acm/contest/882/F 题意:输入一个数n,代表总共有2n个人,然后每个人对所有人有个贡献值,然后问 ...
- 2019牛客暑期多校训练营(第六场)C - Palindrome Mouse (回文自动机)
https://ac.nowcoder.com/acm/contest/886/C 题意: 给出一个串A , 集合S里面为A串的回文字串 , 现在在集合S里面找出多少对(a,b),b为a的字串 分析: ...
- 2019牛客暑期多校训练营(第六场)J Upgrading Technology
传送门 题意: 就是给你n个技能,每个技能最高升到m级,每升一级就是耗费Cij钱,这个Cij可能是负的,如果所有技能都升到或者说超过j等级,就会获得Dj钱,这个Dj也有可能是负值,让你求你最多得到多少 ...
- 山东省ACM多校联盟省赛个人训练第六场 poj 3335 D Rotating Scoreboard
山东省ACM多校联盟省赛个人训练第六场 D Rotating Scoreboard https://vjudge.net/problem/POJ-3335 时间限制:C/C++ 1秒,其他语言2秒 空 ...
- HDU6578 2019HDU多校训练赛第一场 1001 (dp)
HDU6578 2019HDU多校训练赛第一场 1001 (dp) 传送门:http://acm.hdu.edu.cn/showproblem.php?pid=6578 题意: 你有n个空需要去填,有 ...
- HDU6579 2019HDU多校训练赛第一场1002 (线性基)
HDU6579 2019HDU多校训练赛第一场1002 (线性基) 传送门:http://acm.hdu.edu.cn/showproblem.php?pid=6579 题意: 两种操作 1.在序列末 ...
- HDU6621 K-th Closest Distance HDU2019多校训练第四场 1008(主席树+二分)
HDU6621 K-th Closest Distance HDU2019多校训练第四场 1008(主席树+二分) 传送门:http://acm.hdu.edu.cn/showproblem.php? ...
随机推荐
- 7 —— node —— 响应图片
const http = require('http'); const fs = require('fs'); const server = http.createServer(); server ...
- JPA#复杂查询#自定义查询
编写自定义SQL基于下面信息:1. SpringData JPA 在为Repository接口生成实现的时候,会查找是否有 "接口名称"+"Impl"的类,如果 ...
- 105-PHP使用var_dump查看类的类型
<?php class ren{ //定义人类 } class mao{ //定义猫类 } $ren1=new ren(); //实例化人类 $mao1=new mao(); //实例化猫类 $ ...
- ActiveMQ消息队列和SignalR之日志实时监控及警报小实例
主要技术: log4net-生成日志. ActiveMQ-生成日志的时候发送消息,并实时监控日志. SignalR-将ActiveMQ监控的日志实时显示到浏览器上,而不用刷新浏览器. 小实例介绍: 左 ...
- jqGrid 添加 合计行 footDate
jQuery(table_id).jqGrid({ url : url,//组件创建完成之后请求数据的url datatype : "json",//请求数据返回的类型.可选jso ...
- ISO处理jq事件
jq事件在IOS上,如果标签本身不具备某些功能,而我们通过document或者body绑定上去的话,有些事件是不起作用的,同时在IOS上jq事件也存在延迟. 事件不起作用 这里点击事件为例,在IOS中 ...
- mysql视图使用方法
1.为什么要使用视图 对于复杂的查询,往往是有多个数据表进行关联查询而得到,而这种语句往往比较复杂,也可能非常频繁的使用.比如: select 字段一,字段二.字段三, from 数据表1 join ...
- Java算法练习——整数反转
题目链接 题目描述 给出一个 32 位的有符号整数,你需要将这个整数中每位上的数字进行反转. 备注 注意: 假设我们的环境只能存储得下 32 位的有符号整数,则其数值范围为 $[−2^{31}, 2^ ...
- CSS font-family 各字体一览表
windows常见内置中文字体字体中文名 字体英文名宋体 SimSun(浏览器默认) 黑体 SimHei 微软雅黑 ...
- ssh and scp从远程服务器下载文件
scp -r root@172.16.252.32:/home/files /home/files 下载目录 -r root是用户172.16.252.32是ip:/home/files 是你要 ...