HHHOJ #1247. 「NOIP 2023 模拟赛 20230715 A」1 题解--zhengjun

发布时间 2023-07-16 15:32:10作者: A_zjzj

法老找来的题,说是找了三道 其他模拟赛的 T4 拼成 T1~T3,另外搞了道 T4。

思维好题,但是放在 T1 有点搞心态,但是还好大样例够强,400 没挂。

然而 T3 大样例输出错了,浪费了我 0.5h,差评

首先发现向左走之后向右走是一定不优的,所以最短路的情况只能先向右再向左。

考虑枚举起点 \(s\),分两步转移 \(s\to u\) 的最短距离 \(f_u\)

向右走

转移为 \(f_i=\min\limits_{j=b_i}^{i-1}\{f_j\}+1\)

容易想到 \(O(n\log n)\) 的做法,但是这里考虑 \(O(n)\)

\(g_x=\max\limits_{f_i=x}\{i\}\)

那么 \(f_i=\min\limits_{g_x \ge b_i}\{x\}\)

由于 \(b_i<i\),所以 \(f_i\le f_{i-1}+1\)

下面证明单调性:

\[t=f_{i-1},g_t=i-1 \]

  • 此时

\[g_0,\cdots,g_{t-1}<g_{t} \]

  • 得到

\[g_0,\cdots,g_{t-2},g_{t-1}<g_t\Longrightarrow\\ g_0,\cdots,g_{t-2}<b_{g_{t}}\le g_{t-1}\Longrightarrow\\ g_0,\cdots,g_{t-2}<g_{t-1} \]

  • 归纳得到:

\[\forall x\le f_{i-1},g_{x-1}<g_x \]

所以在 \(i\) 的时候,可以从 \(f_{i-1}+1\) 开始向下枚举 \(f_i\),不合法时直接 break,做到均摊 \(O(n)\)

向左走

可以在原来的 \(f\) 数组上继续 dp。

转移为:\(f'_i=\min\{f_i,\min\limits_{a_j\le i<j}\{f'_j\}\}\),下面考虑后一项的。

一样的套路,记 \(g_x=\min\limits_{f_i=x}\{a_i\}\)

转移变为:\(f'_i=\min\limits_{g_x\le i}\{x\}\)

也能够证明单调性:

  • 如果 \(g_{x-1}\le i,g_{x}>i\),那么 \(f'_{i+1}\le x\)
    • \(f'_{i+1}=x\),那么 \(g_x\le a_{i+1}\le i\),不符合;
    • \(f'_{i+1}<x\),那么由于 \(f'_{i}\le f'_{i+1}+1\le x\),所以无需考虑上界 \(x\)
  • 所以 \(g_{x-1}\le i\Longrightarrow g_x\le i\),单调性得到证明。

所以 \(f'_i\) 直接从 \(f'_{i+1}+1\) 开始向下枚举验证即可。

代码

#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=6e3+10,INF=1e9;
int n,a[N],b[N];
int pos[N],f[N];
int ans;
void solve(int s){
	fill(f+1,f+1+n,INF),f[s]=0;
	fill(pos,pos+1+n,0),pos[0]=s;
	for(int i=s+1,x=0;i<=n;i++){
		x++;
		for(;x>0&&pos[x-1]>=b[i];x--);
		pos[f[i]=x+1]=i;
	}
	fill(pos,pos+1+n,n+1);
	for(int i=n,x=INF;i>=1;i--){
		x=min(x+1,f[i]);
		for(;x>0&&pos[x-1]<=i;x--);
		f[i]=min(f[i],x+1);
		pos[f[i]]=min(pos[f[i]],a[i]);
	}
	for(int i=1;i<=n;i++)ans^=f[i]*(i+s);
}
int main(){
	scanf("%d",&n);
	for(int i=2;i<=n;i++)scanf("%d",&a[i]);
	for(int i=2;i<=n;i++)scanf("%d",&b[i]);
	for(int i=1;i<=n;i++)solve(i);
	cout<<ans;
	return 0;
}