题意:

给定一个矩阵,矩阵上有若干点,每个点有正或负的权值,找一个方框框住一些点使得方框中点权值最大。

题解:

离散化横纵坐标,容易将这个问题转化为在矩阵上求最大和子矩阵的问题。

普通的n*n的矩阵的子矩阵最大和正解为$O(n^3)$,枚举上下端点后dp

然而此题是一个稀疏矩阵,n*n矩阵中只有O(n)个点,要求$O(n^2logn)$解法。

正解是枚举上下端点,用线段树维护区间最大和,每枚举到一个下端点,将这个下端点上所有的点的权值更新到线段树上,每次更新logn

由于点的个数是O(n)的,因此每枚举一个上端点,最多更新线段树O(n)次,总时间复杂度$O(n^2logn)$

用线段树维护区间最大和,需要节点上保存如下信息:sum(总和),maxsum(最大子段和),lmax(最大前缀和),rmax(最大后缀和)

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
typedef pair<int, LL>P;
const int M = 2e3 + 5;
const LL mod = 998244353;
const LL lINF = 0x3f3f3f3f3f3f3f3f;
#define ls (rt<<1)
#define rs (rt<<1|1)
LL gcd(LL a, LL b) { return b ? gcd(b, a%b) : a; }
struct node {int l, r;LL sum, lsmx, rsmx, mx;
}tr[M * 4];
LL a[M];
LL pre, ans;//pre前一个搜索区间从右端点开始的最大子段和,ans当前最大值
void pushup(int rt)
{tr[rt].sum = tr[ls].sum + tr[rs].sum;tr[rt].lsmx = max(tr[ls].sum + tr[rs].lsmx,tr[ls].lsmx);tr[rt].rsmx = max(tr[rs].sum + tr[ls].rsmx,tr[rs].rsmx);tr[rt].mx = max(max(tr[ls].mx, tr[rs].mx), tr[ls].rsmx + tr[rs].lsmx);
}
void build(int rt, int l, int r)
{tr[rt].l = l, tr[rt].r = r;if (l == r){tr[rt].sum += a[l];tr[rt].lsmx += a[l];tr[rt].rsmx += a[l];tr[rt].mx += a[l];return;}int mid = (l + r) >> 1;build(ls, l, mid);build(rs, mid + 1, r);pushup(rt);
}
void update(int rt, int l, int r, int pos, LL v)
{if (l == r){tr[rt].sum += v;tr[rt].lsmx += v;tr[rt].rsmx += v;tr[rt].mx += v;return;}int mid = (l + r) >> 1;if (pos <= mid)update(ls, l, mid, pos, v);elseupdate(rs, mid + 1, r, pos, v);pushup(rt);
}
void query(int rt, int ql, int qr)
{if (ql <= tr[rt].l&&tr[rt].r <= qr){ans = max(ans, tr[rt].mx);ans = max(ans, pre + tr[rt].lsmx);pre = max(pre + tr[rt].sum, tr[rt].rsmx);return;}int mid = (tr[rt].l + tr[rt].r) >> 1;if (ql <= mid){query(ls, ql, qr);}if (qr > mid){query(rs, ql, qr);}}
int n;
int q;
LL x[M], y[M], val[M];
LL xid[M], yid[M];
int xsz, ysz;
vector<P>ve[M];
int main()
{int _;scanf("%d", &_);while (_--){scanf("%d", &n);for (int i = 1; i <= n; i++){scanf("%lld%lld%lld", &x[i], &y[i], &val[i]);xid[i] = x[i];yid[i] = y[i];}sort(yid + 1, yid + 1 + n);sort(xid + 1, xid + 1 + n);xsz = unique(xid + 1, xid + 1 + n) - (xid + 1);ysz = unique(yid + 1, yid + 1 + n) - (yid + 1);for (int i = 1; i <= ysz; i++){ve[i].clear();}for (int i = 1; i <= n; i++){int xi, yi;xi = lower_bound(xid + 1, xid + 1 + xsz, x[i]) - xid;yi = lower_bound(yid + 1, yid + 1 + ysz, y[i]) - yid;ve[yi].push_back(make_pair(xi, val[i]));}ans = 0;for (int i = 1; i <= ysz; i++)//下边界
        {memset(tr, 0, sizeof tr);for (int j = i; j >= 1; j--)//上边界
            {for (auto tmp : ve[j]){update(1, 1, xsz, tmp.first, tmp.second);}ans = max(ans, tr[1].mx);}}printf("%lld\n", ans);}
}

转载于:https://www.cnblogs.com/isakovsky/p/11323176.html

hdu多校第六场1005 (hdu6638) Snowy Smilel 线段树/区间最大和相关推荐

  1. 吉首大学校赛 K 白山茶与红玫瑰 (线段树区间操作)

    链接:https://ac.nowcoder.com/acm/contest/925/K 来源:牛客网 题目描述 公元2019年6月22日,白山茶王国与红玫瑰王国展开大战,在世外仙境--天空花园处,双 ...

  2. 2018牛客网暑假ACM多校训练赛(第六场)I Team Rocket 线段树

    原文链接https://www.cnblogs.com/zhouzhendong/p/NowCoder-2018-Summer-Round6-I.html 题目传送门 - https://www.no ...

  3. 杭电多校第六场个人补题6 7 9 10 12

    杭电多校第六场个人补题6 7 9 10 12 6 题意 给定一棵有n结点的树,点权为1~n,求对所有结点子树的mex操作和的最大值 思路 其实就是从最底部开始网上找,由于0是唯一的一个,所欲最好给在最 ...

  4. 2018 HDU多校第四场赛后补题

    2018 HDU多校第四场赛后补题 自己学校出的毒瘤场..吃枣药丸 hdu中的题号是6332 - 6343. K. Expression in Memories 题意: 判断一个简化版的算术表达式是否 ...

  5. hdu 5692 Snacks(dfs序+线段树区间更新)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5692 解题思路:这道题是树节点的点权更新,而且涉及到子树,常用的思路是利用dfs序,用线段树来对区间进 ...

  6. hdu 3966(树链剖分+线段树区间更新)

    传送门:Problem 3966 https://www.cnblogs.com/violet-acmer/p/9711441.html 学习资料: [1]线段树区间更新:https://blog.c ...

  7. hdu 1698 Just a Hook 线段树区间更新

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1698 Let us number the consecutive metallic sticks of ...

  8. POJ 2777 ZOJ 1610 HDU 1698 --线段树--区间更新

    直接将这3题 放一起了  今天在做线段树的东西 这3个都是区间更新的 查询方式互相不同 反正都可以放到一起吧 直接先上链接了 touch me touch me touch me 关于涉及到区间的修改 ...

  9. hdu 多校赛 第二场

    slove  3/12 rank  224 补题   6/12 --------------------------------------------------- hdu 6595 http:// ...

  10. HDU 5861 Road 线段树区间更新单点查询

    题目链接: http://acm.split.hdu.edu.cn/showproblem.php?pid=5861 Road Time Limit: 12000/6000 MS (Java/Othe ...

最新文章

  1. 大型网站架构模式之三
  2. 条件概率、全概率公式与贝叶斯公式
  3. kotlin 一个加号引发的血案
  4. nyoj-小明的密钥(362)--数论
  5. LeetCode 最大正方形
  6. pagehelper 不分页的解决方法
  7. 提升普适性,阿里云官方SDK发布支持Go语言SDK
  8. android jni fork()子进程不运行_Android高级面试谈谈Zygote的理解
  9. Tensorflow Estimator之DNNClassifier
  10. 科大讯飞2017年报:营收54亿利润5.9亿,政府补助1.18亿
  11. 转载:用图片搜索图片的几个好网站
  12. macOS中LaTeX报错Package fontspec Error: The font “SimHei“ cannot be found
  13. 单变量微积分笔记—— 积分方法之换元法总结(简单换元和三角换元)
  14. 2022年全球市场木材干燥窑总体规模、主要生产商、主要地区、产品和应用细分研究报告
  15. Ubuntu更换系统字体
  16. oracle pdb监听配置,oracle 12c 监听
  17. ideapad720s在接通电源情况下,关机后自动重启的解决办法
  18. 基础实验——ADC(模数转换)
  19. 中国冶金工业节能减排规划及投资前景预测报告2022年
  20. C# 调用钉钉接口进行发送企业通知消息,适应于网页版

热门文章

  1. 数字签名、电子签名与电子合同
  2. RabbitMQ虚拟主机
  3. linux 薄膜键盘驱动,一份*级薄膜键盘的简历— HHKB 静电容 键盘
  4. Map container is already initialized.
  5. 2015阿里巴巴北京年会——马云“北伐”讲话
  6. PHP图片尺寸调整(全填充,没有黑底)
  7. 360抢票服务器响应,抢票首日 360浏览器成最佳
  8. 修真院教学模式三大阶段之真实项目
  9. 软件测试的系统代码,软件测试工程师管理系统代码
  10. 我一生中最重要的12个人