Codeforces Round #620 (Div. 2) F2. Animal Observation (hard version) dp + 线段树

传送门

文章目录

  • 题意:
  • 思路:

题意:

在这里插入图片描述
比如下面这个图:
在这里插入图片描述

思路:

对于这个题,比较容易就能考虑到dpdpdp,设f[i][j]f[i][j]f[i][j]为到了第iii行,覆盖了[j,j+k−1][j,j+k-1][j,j+k1]范围时候的最大值。
先不考虑重叠的情况,那么转移方程就是f[i][j]=max(f[i−1][t]+∑x=tt+k−1a[i][x])+∑x=jj+k−1a[i][x]f[i][j]=max(f[i-1][t]+\sum _{x=t}^{t+k-1}a[i][x])+\sum _{x=j}^{j+k-1}a[i][x]f[i][j]=max(f[i1][t]+x=tt+k1a[i][x])+x=jj+k1a[i][x]
考虑重叠的情况,将区间[j,j+k−1][j,j+k-1][j,j+k1]表示为[l,r][l,r][l,r]
对于easyeasyeasy版本,因为k≤20k\le 20k20,所以重叠的部分不会多于202020,所以我们维护一个前后缀的最大值,先取max(pmx[l−1],smx[r+1])max(pmx[l-1],smx[r+1])max(pmx[l1],smx[r+1]),让后对于[l−k+1,r][l-k+1,r][lk+1,r]的部分,我们删掉其重叠的部份,让后转移即可,复杂度O(nmk)O(nmk)O(nmk)
对于hardhardhard版本,由于k≤mk\le mkm,所以上面的复杂度不可行,考虑优化。对于区间滑动问题,我们考虑线段树或者单调队列,由于线段树比较好理解,说一下线段树的吧。
我们设g[i][j]=f[i][j]+∑x=jj+k−1a[i][j]g[i][j]=f[i][j]+\sum _{x=j}^{j+k-1}a[i][j]g[i][j]=f[i][j]+x=jj+k1a[i][j]
考虑a[x][y]a[x][y]a[x][y]都在哪里被重复计算了,可以发现在g[x−1][y−k+1],...,g[x−1][y]g[x-1][y-k+1],...,g[x-1][y]g[x1][yk+1],...,g[x1][y]都被重复计算了,我们只需要减去他们,就可以取maxmaxmax加到f[x][y]f[x][y]f[x][y]上了。这个显然可以用线段树来实现,但是复杂度还是O(nmk)O(nmk)O(nmk)的,但是我们可以发现,当从jjjj+1j+1j+1的时候,jjj不再重叠,j+1+k−1j+1+k-1j+1+k1产生重叠的部分,分别将这两个地方处理一下就好了。

// Problem: F2. Animal Observation (hard version)
// Contest: Codeforces - Codeforces Round #620 (Div. 2)
// URL: https://codeforces.com/contest/1304/problem/F2
// Memory Limit: 512 MB
// Time Limit: 3000 ms
// 
// Powered by CP Editor (https://cpeditor.org)//#pragma GCC optimize("Ofast,no-stack-protector,unroll-loops,fast-math")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4.1,sse4.2,avx,avx2,popcnt,tune=native")
//#pragma GCC optimize(2)
#include<cstdio>
#include<iostream>
#include<string>
#include<cstring>
#include<map>
#include<cmath>
#include<cctype>
#include<vector>
#include<set>
#include<queue>
#include<algorithm>
#include<sstream>
#include<ctime>
#include<cstdlib>
#define X first
#define Y second
#define L (u<<1)
#define R (u<<1|1)
#define pb push_back
#define mk make_pair
#define Mid (tr[u].l+tr[u].r>>1)
#define Len(u) (tr[u].r-tr[u].l+1)
#define random(a,b) ((a)+rand()%((b)-(a)+1))
#define db puts("---")
using namespace std;//void rd_cre() { freopen("d://dp//data.txt","w",stdout); srand(time(NULL)); }
//void rd_ac() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//AC.txt","w",stdout); }
//void rd_wa() { freopen("d://dp//data.txt","r",stdin); freopen("d://dp//WA.txt","w",stdout); }typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;const int N=55,M=20010,mod=1e9+7,INF=0x3f3f3f3f;
const double eps=1e-6;int n,m,k;
int a[N][M],s[N][M],f[N][M];
struct Node
{int l,r;int val,lazy;
}tr[M<<2];void pushup(int u)
{tr[u].val=max(tr[L].val,tr[R].val);
}void pushdown(int u)
{int lazy=tr[u].lazy; tr[u].lazy=0;tr[L].val+=lazy; tr[L].lazy+=lazy;tr[R].val+=lazy; tr[R].lazy+=lazy;
}void build(int u,int l,int r,int pos)
{tr[u]={l,r};if(l==r){tr[u].val=f[pos][l]+s[pos+1][l+k-1]-s[pos+1][l-1];tr[u].lazy=0;return;}build(L,l,Mid,pos); build(R,Mid+1,r,pos);pushup(u);
}void modify(int u,int l,int r,int val)
{if(tr[u].l>=l&&tr[u].r<=r){tr[u].val+=val;tr[u].lazy+=val;return;}pushdown(u);if(l<=Mid) modify(L,l,r,val);if(r>Mid) modify(R,l,r,val);pushup(u);
}void add(int val,int pos)
{if(pos>=1) modify(1,max(1,pos-k+1),min(pos,m-k+1),-a[val][pos]);
}void del(int val,int pos)
{if(pos>=1) modify(1,max(1,pos-k+1),min(pos,m-k+1),a[val][pos]);
}int main()
{
//	ios::sync_with_stdio(false);
//	cin.tie(0);int ans=0;scanf("%d%d%d",&n,&m,&k);for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&a[i][j]),s[i][j]=s[i][j-1]+a[i][j];for(int i=1;i<=n;i++){if(i>1) for(int j=1;j<k;j++) add(i,j);for(int j=1;j<=m-k+1;j++){int l=j,r=j+k-1;f[i][j]=s[i][r]-s[i][l-1];if(i!=1) del(i,l-1),add(i,r),f[i][j]+=tr[1].val;}build(1,1,m-k+1,i);}for(int i=1;i<=m;i++) ans=max(ans,f[n][i]);printf("%d\n",ans);return 0;
}
/**/

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315285.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

设计模式之总体介绍

1. 背景与介绍设计模式是经过反复使用、经过分类的代码总结。设计模式的目的是提高代码可重用性和可靠性&#xff0c;并使代码条理清晰、易于理解、易于维护。设计模式描述了在各种情况下&#xff0c;要选择什么样的方案来解决问题。设计模式通常以类和对象来描述其中的关系和相…

回文自动机:从入门到只会打板

写在前面 如果你会SAMSAMSAM&#xff0c;相信回文自动机不会难懂。 如果你不会&#xff0c;你可以参考我的上一篇文章。 至少回文自动机是治愈系的吧。 作用 回文自动机&#xff0c;也叫回文树&#xff0c;简称PAMPAMPAM实际上它既不是自动机也不是树 处理回文串的有力工…

2021牛客第一场 I. Increasing Subsequence-前缀和优化dp

https://ac.nowcoder.com/acm/contest/11166/I 思路&#xff1a;dp[i][j] 是表示上上步走在i点&#xff0c;上一步走在j点的期望。首先我们很容易想到n^3的做法&#xff0c;那我们必须考虑去优化一维的时间复杂度。我们可以考虑使用前缀和优化dp转移。 我们枚举i点&#xff0c…

Codeforces Round #620 (Div. 2) E. 1-Trees and Queries 思维 + LCA

传送门 文章目录题意思路&#xff1a;题意 思路&#xff1a; 照例&#xff0c;先考虑不加边怎么做。由于可以经过重复的边或点&#xff0c;设aaa与bbb之间长度为lenlenlen&#xff0c;那么需要len<klen<klen<k并且还需要(k−len)mod20(k-len) \bmod 20(k−len)mod20&…

.NET Core 微服务之Polly熔断策略

紧接着上一篇说&#xff0c;咱们继续介绍Polly这个类库熔断策略&#xff08;Circuit-breaker&#xff09;如果调用某个目标服务出现过多超时、异常等情况&#xff0c;可以采取一定时间内熔断该服务的调用&#xff0c;熔断期间的请求将不再继续调用目标服务&#xff0c;而是直接…

牛客第二场 G.League of Legends-单调队列优化dp

https://ac.nowcoder.com/acm/contest/11253/G 上面出题人给的题解&#xff1a; 思路基本差不多&#xff0c;这里主要说一下合并小区间的dp&#xff0c; dp[i][j]代表前i个分成j组最大的时间max 我们首先将区间排好序&#xff0c;如果满足a[k]>b[i] ,则有 j都是由j-1转…

.NET中扩展方法和Enumerable(System.Linq)

LINQ是我最喜欢的功能之一&#xff0c;程序中到处是data.Where(xx>5).Select(x)等等的代码&#xff0c;她使代码看起来更好&#xff0c;更容易编写&#xff0c;使用起来也超级方便&#xff0c;foreach使循环更加容易&#xff0c;而不用for int..&#xff0c;linq用起来那么爽…

Planes, Trains, but not Automobiles-求最小路径覆盖的起点终点

https://vjudge.net/problem/Kattis-planestrainsbutnotautomobiles 题意&#xff1a;给一个有向图&#xff0c;火车可以由任意一个起点开始&#xff0c;每一个点只能经过一次&#xff0c;在坐火车的时候你可以选择坐飞机到另外一个点&#xff0c;求坐飞机的最小次数&#xff…

译 | .NET Core 3.0 对诊断的改进

原文&#xff1a;Sourabh Shirhatti翻译&#xff1a;Edi Wang在 .NET Core 3.0 中&#xff0c;我们将引入一套工具&#xff0c;这些工具利用 .NET 运行时中的新功能&#xff0c;使诊断和解决性能问题变得更加容易。这些运行时功能可帮助您回答一些常见的诊断问题&#xff1a;我…

【AT987】高橋君【组合数】【莫队】

传送门 题意&#xff1a;TTT组询问NNN个相同物品选不超过KKK个的方案数&#xff0c;T,N≤1e5T,N \leq 1e5T,N≤1e5 设f(x,y)∑i0yCxif(x,y)\sum_{i0}^{y}C_x^if(x,y)∑i0y​Cxi​即所求 直接求并没有很好的性质 但我们发现:f(x,y)2f(x−1,y)−Cx−1yf(x,y)2f(x-1,y)-C_{x-1…

2021牛客第一场H.Hash Function—FFT求差值的卷

https://ac.nowcoder.com/acm/contest/11166/H 官方题解。 比赛时&#xff0c;我们都是用暴力写的&#xff0c;数据太弱了&#xff0c;今天突然想起来&#xff0c;用fft写了一下。 主要使用fft求差值的卷&#xff0c;乘法我们在多项式的系数上做&#xff0c;加减法 我们在多项…

设计模式之策略者模式

策略者模式简介策略者模式定义一个算法接口&#xff0c;并由其实现类去实现&#xff0c;使得每一个算法都得到封装&#xff0c;并让他们可以相互替换。这是一种行为型模式。策略者模式降低了算法行为和环境角色的耦合度&#xff0c;使得算法可以独立发生变化。策略者模式在现实…

Codeforces Round #622 (Div. 2) D. Happy New Year 状压dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; n≤1e5,m≤1e9,k≤8.n\le 1e5,m\le 1e9,k\le 8.n≤1e5,m≤1e9,k≤8. 思路&#xff1a; 注意到题目中保证了每个孩子至多收到kkk个&#xff0c;且k≤8k\le 8k≤8&#xff0c;注意到这是题目保证的&#xff0…

2021牛客第四场-E Tree Xor-线段树区间异或

https://ac.nowcoder.com/acm/contest/11255/E 题目大意&#xff1a;有n个点&#xff0c;每一点有一个区间[l,r]范围的权值可能&#xff0c;有n-1条边&#xff0c;边上的权值为w,ww(u)^w(v),问你满足权值范围和边权值的方案数。 思路&#xff1a;我们把树根的值设为0&#xff…

【SPOJ2666】QTree4【链分治】

传送门 题意&#xff1a;给一棵带边权的树&#xff0c;每个点开始时为白色&#xff0c;维护两种操作&#xff1a; 1.改变一个点的颜色&#xff08;白变黑&#xff0c;黑变白&#xff09; 2.询问最远的两个白点之间的距离 树分治国集论文 链分治的本质其实就是树链剖分。它们…

2021牛客第五场 I.Interval Queries-回滚莫队

如图&#xff1a;一道很裸的回滚莫队&#xff0c;注意加入的操作和回滚的操作就好了。 #include <cstdio> #include <cstring> #include <string> #include <cmath> #include <iostream> #include <algorithm> #include <queue> #i…

Linux下搭建asp.net运行环境

最近有个项目&#xff0c;是在Windows平台下开发的&#xff0c;需要把 asp.net web应用移植到 CentOS下&#xff0c;甚是头疼&#xff1b;翻阅资料&#xff0c;发现Jexus是个可行的方案&#xff0c;下面是官方对Jexus的定义&#xff1a;什么是JexusJexus是一款Linux平台上的高性…

Codeforces Round #736 (Div. 2)E. The Three Little Pigs-长除法求多项式系数

https://codeforces.com/contest/1549/problem/E 题目大意自己悟吧&#xff0c;不解释了。 这题在赛时我想的是fft&#xff0c;但是题解说不用fft就能求出多项式的系数&#xff0c;fft在这题时间复杂度有点高。 使用多项式长除法就能求出系数&#xff0c;k^x的系数就是抓x个p…

在 Asp.Net Core WebAPI 中防御跨站请求伪造攻击

什么是跨站请求伪造跨站请求伪造&#xff08;英语&#xff1a;Cross-site request forgery&#xff09;&#xff0c;也被称为 one-click attack 或者 session riding&#xff0c;通常缩写为 CSRF 或者 XSRF&#xff0c; 是一种挟制用户在当前已登录的Web应用程序上执行非本意的…

HDU2021多校第五天 1009 Array-树状数组求二阶前缀和

https://acm.hdu.edu.cn/showproblem.php?pid7020 题目大意&#xff1a;求出区间众数严格大于区间一半大小的子区间的数量。 思路&#xff1a;先对每一个数求出对应的位置排列&#xff0c;S[j]代表1-有多少个i&#xff0c;易得s[r]-s[l]>r-l-(s[r]-s[l])的时候子区间{l1&a…