猫猫cpu的缓存(NW)

在做 dp 的时候经常在某个细节处死扣,总的说还是没把握住题目的特征性和 dp 的整体架构。既然能一眼看出来是 dp 题,当根据题目特征列出个常规的 dp 式子 ,之后才有优化的资本。这是一般 dp 题的流程,重点在列出像样的 dp 式子,写一大坨不可用也不可优化的极端、特殊式子是没有意义的,要在把握 dp 的线性的基础上尝试状态转移,才能将千奇百怪的 dp 题化归成一套有迹可循的系统和模型

读原题面,叽里呱啦写了千把字,结果最后一行给了下面这段话(


简而言之:

给你 m 个 1 到 n 之间的整数,你要找到若干个大小为固定的 k 的闭区间,使得所有这些数都在你找到的某个区间内。你需要最小化这些区间的并集的大小,并输出此大小。本题里区间/区间并集的大小,被定义为,这个区间/区间并集里整数的个数。
1 ≤ k ≤ n ≤ 1 0 9 , 1 ≤ m ≤ 3 ⋅ 1 0 5 , p i ≤ n 1\le k\le n\le 10^9,1\le m\le 3\cdot 10^5,p_i\le n 1≤k≤n≤109,1≤m≤3⋅105,pi≤n 且互不相同。


考场做法

考场上做 dp 题时脑子很乱,到最后实在不行推出来也推出一个很乱的式子,只过了两个样例,显然地暴毙,挂光光了,说明一开始思路的走向就不对。以后开头的推理要更加谨慎,开头思路歪后面很难跳出来。

这是考场的代码, f i , 0 / 1 f_{i,0/1} fi,0/1 表示 i i i 当开头(结尾)时囊括前 i i i 个点的最小值。自己看了都绷不住。

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=3e5+5;
int n,k,m,p[N],f[N][2],l[N];
int main(){
    freopen("cpu.in","r",stdin);
    freopen("cpu.out","w",stdout);
    scanf("%d%d%d",&n,&k,&m);
    for(int i=1;i<=m;i++) scanf("%d",p+i);
    sort(p+1,p+m+1);
    if(n-p[1]+1<=k||m==1) return printf("%d",k),0;
    l[0]=1;
    for(int i=1,L=1;i<=m;i++){
        while (p[i]-p[L]>=k) L++;
        l[i]=L;
    }
    // for(int i=1;i<=m;i++) printf("%d ",l[i]);puts("");
    memset(f,0x3f,sizeof(f));
    f[1][0]=k;if(p[1]>=k) f[1][1]=k;
    for(int i=2;i<=m;i++){
        if(n-p[i]+1>=k) f[i][0]=min(f[i-1][1]+k,f[l[i-1]][0]+max(0,p[i]-p[l[i-1]]));
        if(p[i]>=k) f[i][1]=min(f[l[i]-1][1]+k,f[l[l[i]-1]][0]+max(0,p[i]-p[l[l[i]-1]]-k+1));
        // printf("%d %d\n",f[i][0],f[i][1]);
    }
    printf("%d",f[m][1]);
    return 0;
}

正解

注意题目特征,张超老师说想出第一个部分分就接近正解了。第一个部分分 O ( 2 n ⋅ n ) O(2^n\cdot n) O(2n⋅n),重要在于怎么判断一个序列是否合法:
m m m 个点包含是必然的,注意到题目的特征性,只要保证每个选下的区间长度都大于等于 k k k 即可。

类似的题目还有[USACO18JAN] Stamp Painting G

对 p p p 排序是显然的。

所以我们设 f i f_i fi 表示把前 i i i 个要囊括的点囊括需要的区间长度最小值,由于两个距离小于 k k k 的 p i , p j p_i,p_j pi,pj 转移方法不同,可以先考虑距离大于 k k k 的 p i , p j p_i,p_j pi,pj 如何转移,可以列出:
f i = min ⁡ j = 1 p o s ( f j − 1 + p i − p j + 1 ) f_i=\min_{j=1}^{pos}(f_{j-1}+p_i-p_j+1) fi=minj=1pos(fj−1+pi−pj+1),其中 p o s pos pos 是 i i i 往前第一个满足 p i − p j ≥ k p_i-p_j\ge k pi−pj≥k 的位置。

抽离无关量,得 f i = p i + 1 + min ⁡ j = 1 p o s ( f j − 1 − p j ) f_i=p_i+1+\min_{j=1}^{pos}(f_{j-1}-p_j) fi=pi+1+minj=1pos(fj−1−pj)

显然 p o s pos pos 单调不降,故 i i i 单增过程中 j j j 的范围只增不减,尺取不断取 min ⁡ \min min 即可。

对于距离小于 k k k 的两点,只会在 p i p_i pi 和 p p o s + 1 p_{pos+1} ppos+1 时对答案产生贡献,增量为 k k k。

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
const int N=3e5+5;
int n,k,m,p[N],f[N],w=0,mn=1e9;
int main(){
    freopen("cpu.in","r",stdin);
    freopen("cpu.out","w",stdout);
    scanf("%d%d%d",&n,&k,&m);
    for(int i=1;i<=m;i++) scanf("%d",p+i);
    sort(p+1,p+m+1);
    for(int i=1;i<=m;i++){
        while (p[i]-p[w+1]>=k) mn=min(mn,f[w]-p[++w]);
        f[i]=min(p[i]+1+mn,f[w]+k);
    }
    printf("%d",f[m]);
}

总的来说,要把握住 dp 的常态化思路,抓题目特征,列可行性强的状态,状态转移,最后考虑优化。

相关推荐
Aaron15884 分钟前
全频段 SDR干扰源模块解决方案(星链干扰、LORA无人机干扰)
人工智能·算法·fpga开发·硬件架构·硬件工程·无人机·信息与通信
AI科技星8 分钟前
全域数学·球面拓扑微扰标准系数η=0.01 应用详解(典籍正式版)
人工智能·算法·数学建模·数据挖掘·机器人
逻辑君23 分钟前
物理学研究报告【20260001】
人工智能·算法
AI科技星32 分钟前
算法联盟·全域数学公理体系下黑洞标量毛发与LVK引力波O4全维理论、求导、证明、计算、验证、分析
人工智能·线性代数·算法·架构·学习方法·量子计算
谙弆悕博士32 分钟前
【附C语言源码】C语言 栈结构 实现及其扩展操作
c语言·开发语言·数据结构·算法·链表·指针·
YuanDaima204833 分钟前
图论基础原理与题目说明
数据结构·人工智能·python·算法·图论·手撕代码
AI人工智能+电脑小能手34 分钟前
【大白话说Java面试题 第53题】【JVM篇】第13题:JVM采用什么算法判断一个对象是否需要被回收?
java·jvm·算法·面试
小赵不会秃头38 分钟前
数据结构Day 06:线性结构、库操作及 Makefile 完整学习笔记
java·linux·数据结构·算法·面试
Shan12051 小时前
在C++中尝试封装为函数
开发语言·c++·算法
郝学胜-神的一滴1 小时前
干货版《算法导论》04:渐近复杂度与序列接口实战
java·开发语言·数据结构·c++·python·算法