前言

\(HE\)沾\(BJ\)的光成功滚回家里了...这堆最大子段和的题抠了半天,然而各位\(dalao\)们都已经去做概率了...先%为敬。

引流之主:老姚的博客

最大M子段和 V1

思路

  • 最简单的ver.数据范围在5000以内,可以考虑暴力一点的做法\(O(n^3)\),定义\(dp\)状态\(dp[i][j]\),递推式子:
\[dp[i][j]=max\{dp[i-1][j],dp[k][j-1]\}+a[i]\ (j-1\le k<i)
\]
  • 其中\(i\)表示序列中前\(i\)个元素且包含第\(i\)个元素,而\(j\)则表示划分了\(j\)个不同的字段。而对于\(a[i]\)来说,要么和最后一个字段合并,要么单独成为一个字段,所以就可以得出动态转移方程。
  • 此时我们可以看出,对于第二维\(j\),只与\(j\)和\(j-1\)有关,所以我们可以互换,使\(dp[i][j]\)表示把序列前\(j\)个元素分成\(i\)个子段且包含\(a[j]\)的最大子段和。此时的状态就只和上一行和左边的状态有关了,所以我们就可以把效率压到\(O(m*n)\)。

    此时的动态转移方程:
\[dp[i][j]=max\{dp[i][j-1],dp[i-1][k]\}+a[j]\ (i-1\le k<j)
\]
  • 还有一个要特判的就是\(i=j\)的时候,见代码即可。

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e4+5;
const ll INF=1e15;
ll a[maxn],dp[2][maxn];
int n,m; int main(){
scanf("%d%d",&n,&m); for(int i=1;i<=n;i++)
scanf("%lld",&a[i]); int k=1; for(int i=1;i<=m;i++,k=k^1){
ll temp=-INF;
for(int j=i;j<=n;j++){
temp=max(temp,dp[k^1][j-1]);
if(i==j)
dp[k][j]=dp[k^1][j-1]+a[j];
else
dp[k][j]=max(temp,dp[k][j-1])+a[j];
}
} ll ans=-INF;
for(int i=m;i<=n;i++)
ans=max(ans,dp[k^1][i]);//注意此处是k^1记录的才是第m行的状态
printf("%lld\n",ans); return 0;
}

最大M子段和 V3

思路

  • 为什么先说V3因为V3是环形的边界处理比V2还要简单。对于V1的DP虽然效率是\(O(m*n)\)但是显然从V2开始就已经有m接近于n的数据了,所以最差的效率为\(O(n^2)\),所以DP显然是不行了。这时候用可撤销贪心。
  • 可撤销贪心简单来说就是贪心时候的悔棋操作,可能最一开始选择的为最大的贡献\(w[i]\),但是第二轮选择可能会发现选择\(w[i-1]\)和\(w[i+1]\)更优,那么此时的答案就应加上\(w[i-1]+w[i+1]-w[i]\),相当于建立了一个新的节点,同时其相邻关系也会发生变化。
  • 对于此题,我们可以把相邻的相同符号的数都加在一起(0加在哪里都一样甚至可以选择直接把0删了),这样就得到了一个正负交替的环形序列,而每个数均代表一个区间,此时设序列中正数个数为\(cnt\),正数的和为\(ans\)。
    • 若\(cnt\le m\),则依照题目要求直接输出\(ans\)。
    • 否则需要删减\(cnt-m\)个区间。
  • 考虑如何删减:
    • 可以选择直接删除一个正数。则每次删除最小的正数。
    • 也可以将一个负数和相邻两个正数合并,即将两个正数区间合并成一个。这时选择最大的负数和两边合并。
  • 但是此时不能随便删除,因为可能通过多次合并减少区间得到的答案会更优,盗老姚的数据:
    • 10 -4 3 -4 8 -100

    • 若先删除最小的正数3,则第二步要选择最大的负数-4,但是此时的3已经被删除了无法合并,因此我们要用到可撤销贪心思想。
  • 实现方法:
    • 将修改后的序列中的正数改为负数,压入大根堆(网上好多题解都说用绝对值,但实际上正变负更方便一点,如果用pair的话都不用重载运算符了)。
    • 可撤销贪心执行\(cnt-m\)次。

一个不怕死的HACK

  • 老姚的代码
#include <bits/stdc++.h>
typedef long long LL;
const int maxn=2e5+5;
struct Node{
int id;
LL w;
Node(){};
Node(int x,LL y){id=x;w=y;}
bool operator <(const Node &a)const{
return w<a.w;
}
};
int n,m,N=0,L[maxn],R[maxn];
LL ans=0,a[maxn],b[maxn];
std::priority_queue<Node> q;
bool flag[maxn];
void Init(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i){
scanf("%lld",&a[i]);//原始数组
if(!N || (a[i]>=0)!=(b[N]>=0))
b[++N]=a[i];//新数组,上个数和当前数符号不一样节点++
else
b[N]+=a[i];//符号一样,累加
}
if((b[1]>=0)==(b[N]>=0))//因为是环,b[1]和b[N]同号就合并
b[1]+=b[N--];//合并后新数组个数减一,记得N--
}
void Solve(){
Init();
int tot=0;//记录正数个数
for(int i=1;i<=N;++i){//遍历新数组,压入大根堆
L[i]=i-1; R[i]=i+1;//初始化i的左右邻居
if(b[i]>0){//正数累加到答案,然后变负
ans+=b[i]; b[i]=-b[i]; tot++;
}
q.push(Node(i,b[i]));
}
R[N]=1; L[1]=N;//注意是环形
if(m>=tot){//正数个数小于m,则全部选
printf("%lld\n",ans);return;
}
m=tot-m;//正数大于m则合并或删除tot-m个区间
while(m--){
Node t=q.top(); q.pop();
int i=t.id;
if(flag[i]){++m; continue;}
else{//可撤销贪心
ans+=b[i];
flag[L[i]]=1;
flag[R[i]]=1;
b[i]=b[L[i]]+b[R[i]]-b[i];
R[L[L[i]]]=i;
L[R[R[i]]]=i;
q.push(Node(i,b[i]));//压入新点
L[i]=L[L[i]];
R[i]=R[R[i]];
}
}
printf("%lld\n",std::max(0LL,ans));
}
int main(){
Solve();
return 0;
}
  • 在51Nod直接AC了,不得不说51Nod真的都是随机数据。
  • 有个很明显的错误就是如果给出的序列都是一个符号,例如下面这个HACK数据:
    • Input:

      5 4

      1 1 1 1 1

      Output:

      5

    • 然而老姚的代码会输出0,问题在这行:
    • if((b[1]>=0)==(b[N]>=0))
      b[1]+=b[N--];
    • 如果全为一个符号,最后\(N=1\),显然\(b[1]=b[N]\),这样会导致\(N\)变成\(0\)。如果对于全负,答案必输出0没什么影响,但是全正的话就不行了。不过也很好弄,加个特判即可。具体见代码(把老姚的代码改了改增加可读性(雾))。

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=2e5+10; struct Node{
int id;
ll w;
Node(){};
Node(int x,ll y){
id=x;w=y;
}
friend bool operator <(const Node &a,const Node &b){
return a.w<b.w;
}
}; int n,m,N;
int L[maxn],R[maxn];
ll ans;
ll a[maxn],b[maxn];
bool vis[maxn],flag=1;
priority_queue<Node> q; void Init(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%lld",&a[i]); if(a[i]<0)flag=0;
else ans+=a[i]; if(!N || (a[i]>=0)!=(b[N]>=0))
b[++N]=a[i];
else
b[N]+=a[i];
}
if((b[1]>=0)==(b[N]>=0))
b[1]+=b[N--];
} int main(){
Init(); int tot=0;
for(int i=1;i<=N;i++){
L[i]=i-1;R[i]=i+1;
if(b[i]>0){
b[i]=-b[i];
tot++;
}
q.push(Node(i,b[i]));
} R[N]=1;L[1]=N; if(m>=tot||flag){
printf("%lld",ans);
return 0;
} m=tot-m; while(m--){
Node t=q.top();
q.pop(); int i=t.id;
if(vis[i]){
m++;
continue;
}else{
ans+=b[i];
q.push(Node(i,b[i])); vis[L[i]]=1;vis[R[i]]=1;
b[i]=b[L[i]]+b[R[i]]-b[i];
R[L[L[i]]]=i;
L[R[R[i]]]=i;
L[i]=L[L[i]];
R[i]=R[R[i]];
}
} printf("%lld",max(ans,0LL)); return 0;
}

最大M子段和 V2

  • 最后就是V2了,一开始尝试特判\(0\)和\(N\)处的边界但发现巨难写,后来发现其实只需要加上一个权值为负无穷的\(0\)号节点就行了...

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=5e4+10; struct Node{
int id;
ll w;
Node(){};
Node(int x,ll y){
id=x;w=y;
}
friend bool operator <(const Node &a,const Node &b){
return a.w<b.w;
}
}; int n,m,N;
int L[maxn],R[maxn];
ll ans;
ll a[maxn],b[maxn];
bool vis[maxn],flag=1;
priority_queue<Node> q; void Init(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%lld",&a[i]);
if(a[i]<0)flag=0;
else ans+=a[i];
if(!N && a[i]<0)continue; if(!N || (a[i]>=0)!=(b[N]>=0))
b[++N]=a[i];
else
b[N]+=a[i];
} while(b[N]<0)b[N--]=0;
} int main(){
Init(); int tot=0;
for(int i=1;i<=N;i++){
L[i]=i-1;R[i]=i+1;
if(b[i]>0){
b[i]=-b[i];
tot++;
}
q.push(Node(i,b[i]));
} R[N]=0;L[0]=N;R[0]=1;
q.push(Node(0,b[0]=-1e15)); if(m>=tot||flag){
printf("%lld",ans);
return 0;
} m=tot-m; while(m--){
Node t=q.top();
q.pop(); int i=t.id;
if(vis[i]){
m++;
continue;
}else{
ans+=b[i];
q.push(Node(i,b[i])); vis[L[i]]=1;vis[R[i]]=1;
b[i]=b[L[i]]+b[R[i]]-b[i];
R[L[L[i]]]=i;
L[R[R[i]]]=i;
L[i]=L[L[i]];
R[i]=R[R[i]];
}
} printf("%lld",max(ans,0LL)); return 0;
}

51Nod 最大M子段和系列 V1 V2 V3的更多相关文章

  1. android google map v1 v2 v3 参考

    V1,V2已经不被推荐使用,谷歌强烈推荐使用V3. 本人在选择时着实纠结了良久,现在总结如下: 对于V1,现在已经申请不到API KEY了,所以不要使用这个版本.这个是网址:https://devel ...

  2. 51nod 最大M子段和系列

    1052 最大M子段和 N个整数组成的序列a[1],a[2],a[3],…,a[n],将这N个数划分为互不相交的M个子段,并且这M个子段的和是最大的.如果M >= N个数中正数的个数,那么输出所 ...

  3. 51nod Bash游戏(V1,V2,V3,V4(斐波那契博弈))

    Bash游戏V1 有一堆石子共同拥有N个. A B两个人轮流拿.A先拿.每次最少拿1颗.最多拿K颗.拿到最后1颗石子的人获胜.如果A B都很聪明,拿石子的过程中不会出现失误.给出N和K,问最后谁能赢得 ...

  4. 51nod 最大M子段和系列(1052、1053、1115)

    51nod1052 数据量小,可使用O(N*M)的DPAC,递推公式: dp[i][j]=max(dp[i-1][j-1], dp[i][j-1])+a[j]; dp[i][j]表示前j个数取 i 段 ...

  5. 51Nod 最大公约数之和V1,V2,V3;最小公倍数之和V1,V2,V3

    1040 最大公约数之和 给出一个n,求1-n这n个数,同n的最大公约数的和.比如:n = 6 1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6,加在一起 = 15 输入 1个数N ...

  6. GoogLeNet 之 Inception v1 v2 v3 v4

    论文地址 Inception V1 :Going Deeper with Convolutions Inception-v2 :Batch Normalization: Accelerating De ...

  7. 从Inception v1,v2,v3,v4,RexNeXt到Xception再到MobileNets,ShuffleNet,MobileNetV2

    from:https://blog.csdn.net/qq_14845119/article/details/73648100 Inception v1的网络,主要提出了Inceptionmodule ...

  8. 51nod 1254 最大子段和 V2 ——单调栈

    N个整数组成的序列a[1],a[2],a[3],…,a[n],你可以对数组中的一对元素进行交换,并且交换后求a[1]至a[n]的最大子段和,所能得到的结果是所有交换中最大的.当所给的整数均为负数时和为 ...

  9. 51nod 1254 最大子段和 V2

    N个整数组成的序列a[1],a[2],a[3],…,a[n],你可以对数组中的一对元素进行交换,并且交换后求a[1]至a[n]的最大子段和,所能得到的结果是所有交换中最大的.当所给的整数均为负数时和为 ...

随机推荐

  1. latex 封面

    latex 封面 代码: \begin{titlepage} \heiti \vspace*{64pt} \begin{center} \fontsize{48pt}{0} XX大学\\ \vspac ...

  2. JS基础回顾_Dom

    方法 reverse 将子元素位置颠倒 <!-- log --> <ul id="reverse"> <li>1</li> < ...

  3. oracle之二实例管理及数据库的启动/关闭

    实例管理及数据库的启动/关闭   2.1 实例和参数文件 1.instance 功能:用于管理和访问database.instance在启动阶段读取初始化参数文件(init parameter fil ...

  4. django之安装和项目创建

    dos界面下安装django 自动下载和安装:cmd:pip3 install  django 手动安装: 1.登录django官网下载django 2.下载地址:https://www.django ...

  5. 单链表的前K个的逆序输出

    单链表逆序输出也是常被面试官问到题算法题,所以自己就总结了一下,在此贴出算法,与小伙伴们相互交流. 首先要有三个指针,前两个分别指向首节点,首节点的下一个节点,第三个是临时指针,是为了储存首节点的下一 ...

  6. Vue中你忽略的点

    自定义组件使用 v-model <my-component v-model="data"></my-component> 在组件my-component中, ...

  7. Hadoop学习笔记(一):ubuntu虚拟机下的hadoop伪分布式集群搭建

    hadoop百度百科:https://baike.baidu.com/item/Hadoop/3526507?fr=aladdin hadoop官网:http://hadoop.apache.org/ ...

  8. Python爬取天气预报,Ta的城市开始降温了

    前言 本文的文字及图片来源于网络,仅供学习.交流使用,不具有任何商业用途,如有问题请及时联系我们以作处理. PS:如有需要Python学习资料的小伙伴可以加点击下方链接自行获取 python免费学习资 ...

  9. 【思维】The Four Dimensions of Thinking :长线思维的力量

    "经历过这些苦难之后,我拥抱了一种新的人生哲学,就是更多地关注在那些长期可以获得复利的小收获上,而不是那种频繁的短跑冲刺和精力消耗". 斯坦福教授,著名的心理学家Philip Zi ...

  10. 种子爆破&[GWCTF 2019]枯燥的抽奖

    伪随机数的爆破,种子爆破 做到了一道题,就是有个伪随机数爆破的漏洞,当时尽管是看到了这两个敏感的函数,但是先去看其他的了,没有看到什么漏洞,所以我当时是准备直接强行爆破,之后看到使用伪随机数爆破的方式 ...