Codeforces Round #688 (Div. 2) F. Even Harder(dp+后缀最小值优化)

题目

T(T<=500)组数据，每次给出一个长为n(2<=n<=3e3)的数组，

第i(0<=i<n)个位置是ai，表示i可以跳到[i+1,i+ai]的任意整点位置

特别地，如果ai的值是0，就代表在此点不可跳

现在要修改一些点的a值，使得从0跳到n-1只有唯一一条路径，求最小修改数

保证输入时一定有至少一条路径可达，且sumn<=3e3

题解

改值显然改成0是最优的，相当于切断了这个点的所有后继，

dp[i][j]表示这条路径倒数第一个在i，倒数第二个在j的最小代价 

首先，直接删掉[j+1,i-1]里所有能通往>=i的地方的路径，也就是判[j+1,i-1]的点能不能直接到i

dp[i][j]能从dp[j][k]转移而来，枚举k是O(n^3)的，显然不大行

注意到想要的是k+a[k]<i(因为倒数第三个点不能通往倒数第一个点)的dp[j][k]的最小值

这个东西再开个数组now[][]存下来

更新now数组则用当前状态dp[i][j]更新now[i][j+a[j]]，相当于若干个单点更新，

最后i的状态都扫完的时候，做一个后缀最小值即可，

那使用的时候使用now[j][i-1]即可，即代表了dp[j][k+a[k]<=i-1]的dp最小值

心得

自己推出的一道div2的F的dp，纪念一下

dp题/cf题也往往会有一些性质，瞪是瞪不出来做法的

还是要自己多推一推，不要那么着急去看答案

尝试构造一些自己见过的模型，或是尽量尝试一些能转移的状态

不用刻意求一步写到位，可以后续再慢慢尝试优化

后续

思路来源：https://blog.csdn.net/qq_30361651/article/details/110677008

①发现del[i][j]可以现推出来，不用预处理来求，

②发现dp[i][j]根本就没有用，可以直接被优化干掉，

用now[i][j]的值滚来滚去即可，最后输出now[n-1][n-1]的值

强行赋一个意义的话，就是now[i][j](j>=i)表示路径的倒数第一个点在i，倒数第二个点在k，满足k+a[k]<=j时的最小删点数

只是感觉显得失去了一点计数意义，让这个定义变得非常抽象

代码

#include<bits/stdc++.h>
using namespace std;
const int N=3e3+10,INF=0x3f3f3f3f;
int t,n,dp[N][N],a[N];
int del[N][N],now[N][N];
int main(){
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		dp[0][0]=0;
		for(int i=0;i<n;++i){
			scanf("%d",&a[i]);
			for(int j=0;j<i;++j){
				del[j][i]=0;
			}
		}
		for(int i=0;i<n;++i){
			del[i][i]=0;
			if(i)dp[i][i]=INF;
			for(int j=i-1;j>=0;--j){
				del[j][i]=del[j+1][i]+(j+a[j]>=i);
				dp[i][j]=INF;
			}
			for(int j=i;j<n;++j){
				now[i][j]=INF;
			}
		}
		for(int i=1;i<n;++i){
			if(a[0]<i)break;
			dp[i][0]=min(dp[i][0],del[1][i]);
			now[i][a[0]]=min(now[i][a[0]],dp[i][0]);
			for(int j=a[0]+1;j<n;++j){
				now[i][j]=min(now[i][j-1],now[i][j]);
			}
		}
		for(int i=2;i<n;++i){
			for(int j=1;j<i;++j){
				if(j+a[j]<i)continue;
				//dp[i][j]=min(dp[i][j],min(dp[j][k<i && k+a[k]<i])+del[j+1][i]);
				dp[i][j]=min(dp[i][j],now[j][i-1]+del[j+1][i]);
				now[i][j+a[j]]=min(now[i][j+a[j]],dp[i][j]);
			}
			for(int j=i+1;j<n;++j){
				now[i][j]=min(now[i][j-1],now[i][j]);
			}
		}
		int ans=INF;
		for(int i=n-2;i>=0;--i){
			ans=min(ans,dp[n-1][i]);
		}
		printf("%d\n",ans); 
	}
	return 0;
}