Contents

2019杭电多校6_1005算法日常[6/100]

题目

题目链接

2019杭电多校6_1005_HDU6638

思路

自己

一开想偏了,想着像以前0,1矩阵那种同权值点一样叠加成为矩形的最大面积来求解,这样子就能让复杂度在O(n^2)的样子

然后并非如此,这里的权值w是一个可正可负的整数…自己隐隐约约地感觉复杂度要达到O(n^2*log(n)),但是没有往下想,虽然想了也不一定会,,,emmmm,继续努力吧

正解

  • 首先将纵坐标离散化到 O(n) 的范围内,方便后续的处理。
  • 将所有点按照横坐标排序,枚举矩形的上边界,然后往后依次加入每个点,这样就确定了矩形的上下边界。
  • 设 v[y] 表示矩形内部纵坐标为 y 的点的权值和,则答案为 v 的最大子段和,用线段树维护带修改的最大子段和即可。
  • 时间复杂度 O(n^2*log(n))

std理解版

今天比较晚了,明天要出一趟远门,所以暂时没有手写,对不起自己啊,等回家一定要好好把这个重新写几遍

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=2010,M=4100;
int Case,n,m,i,j,k,cb,b[N],pos[N];ll pre[M],suf[M],s[M],v[M],ans;
struct E{int x,y,z;}e[N];
inline bool cmp(const E&a,const E&b){return a.x<b.x;}

/*用纵坐标建的线段树,大佬对于线段树的理解以及如同我对1+1的理解一样了!*/
void build(int x,int a,int b){
  pre[x]=suf[x]=s[x]=v[x]=0;
  if(a==b){
    pos[a]=x;
    return;
  }
  int mid=(a+b)>>1;
  build(x<<1,a,mid),build(x<<1|1,mid+1,b);
}
inline void change(int x,int p){
  x=pos[x];
  s[x]+=p;
  if(s[x]>0)pre[x]=suf[x]=v[x]=s[x];else pre[x]=suf[x]=v[x]=0;
  /*上传的操作精辟,orz*/
  for(x>>=1;x;x>>=1){
    /*根的左边的 = max(左子树之前的,左子树+右子树之前的)*/
    pre[x]=max(pre[x<<1],s[x<<1]+pre[x<<1|1]);
    /*根的右边的 = max(右子树右边的,右子树+左子树右边的)*/
    suf[x]=max(suf[x<<1|1],s[x<<1|1]+suf[x<<1]);
    /*s是直接叠加的*/
    s[x]=s[x<<1]+s[x<<1|1];
    /*区间最大值 = max(左子树最大,右子树最大,左子树后面的+左子树前面的)*/
    v[x]=max(max(v[x<<1],v[x<<1|1]),suf[x<<1]+pre[x<<1|1]);
  }
}
int main(){
  scanf("%d",&Case);
  while(Case--){
    scanf("%d",&n);
    for(cb=0,i=1;i<=n;i++){
      scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].z);
      b[++cb]=e[i].y;
    }
    /*离散化纵坐标*/
    sort(b+1,b+cb+1);
    for(m=0,i=1;i<=cb;i++)if(i==1||b[i]!=b[m])b[++m]=b[i];
    /*给横坐标排序*/
    sort(e+1,e+n+1,cmp);
    ans=0;
    /*用离散化后的纵坐标覆盖掉原来的纵坐标*/
    for(i=1;i<=n;i++)e[i].y=lower_bound(b+1,b+m+1,e[i].y)-b;
    /*枚举上边界,x是行号,是上边界*/
    for(i=1;i<=n;i++)if(i==1||e[i].x!=e[i-1].x){
      build(1,1,m);
      /*加入点确定好下边界,这样上下边界都确定好了*/
      for(j=i;j<=n;j=k){
        /*又是逐步插入空树维护区间最大值的操作*/
        for(k=j;k<=n&&e[j].x==e[k].x;k++)change(e[k].y,e[k].z);
        if(ans<v[1])ans=v[1];
      }
    }
    printf("%lld\n",ans);
  }
}

每天一句叨叨

今天是情人节,然而…没得女朋友…不过有队友和我一起大杭电多校还是很开心的啦(还是有点点失落)

虽然

https://cdn.jsdelivr.net/gh/wolfdan666/BlogPic/%E7%AE%97%E6%B3%95/%E7%AE%97%E6%B3%95%E6%97%A5%E5%B8%B8%E7%9A%84%E5%8F%A8%E5%8F%A8/6%E4%B8%83%E5%A4%95%E6%83%85%E4%BA%BA%E8%8A%82/1.jpg

但是

谁终将声震人间,必长久深自缄默

谁终将点燃闪电,必长久如云漂泊

全力以赴打完这一段时光的退役赛一定会是一个大学乃至人生最珍贵的记忆,所以这段时间先不要让自己被一个体内的激素控制,等时机到了,一定会更加美好!加油吧,少年!