[uoj84]水题走四方
瞬移后无法区分两者,分开时不妨称下一次瞬移的为分身,即仅允许分身瞬移回本身
考虑本身,即从根节点出发向下移动的一条路径,并称路径上分身曾瞬移到的点为关键节点(包括根节点)
对关键节点dp,定义$f_{k}$表示当前两者均在$k$且$k$子树外所有点均被经过的最短时间
枚举上一个关键节点,分析两者的移动,不难得到转移即
$$
f_{k}=\min_{x为k的祖先}\big(f_{x}+(g_{x}-g_{k})-dep_{x}(h_{x}-h_{k})+\max(dep_{k}-\max(d_{\max},dep_{x}),0)\big)
$$
(其中$dep_{k}$表示$k$的深度$,h_{k}$表示$k$子树中叶子个数$,g_{k}$表示$k$子树中叶子深度和$,d_{\max}$表示$x$子树内除$k$子树外最深的叶子)
另外,初始状态即$f_{1}=1$,最终答案即$\min f_{k}+g_{k}-dep_{k}h_{k}$
若$d_{\max}\ne 0$且$dep_{k}>d_{\max}$,可以理解为其到达深度为$d_{\max}$的点后分身直接瞬移,显然这样不劣
同时,若$d_{\max}=0$,即$x$到$k$的这条链上(除$k$以外)没有其他儿子,直接从$x=fa_{k}$依次转移即可
进一步的,若$x$不是第一个满足$dep_{k}\le d_{\max}$的祖先,那么不妨先转移到该祖先再转移到$k$,代入可得不劣
综上,仅需考虑$x=fa_{k}$(在$fa_{k}$仅有$k$一个儿子时)和第一个满足$dep_{k}\le d_{\max}$的祖先
关于后者,对树长链剖分,注意到其所在长链顶端的父亲一定满足
换言之,每一条长链内是独立的,分别维护一个单调栈执行上述过程即可
需要注意空间限制,由于保证编号,可以使用遍历的方式代替dfs
时间复杂度为$o(n)$,可以通过
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include<bits/stdc++.h> 2 using namespace std; 3 #define N 5000005 4 #define ll long long 5 vector<int>V[N]; 6 int n,n0,st[N],vis[N],fa[N],dep[N],h[N],mx[N],cmx[N],son[N],top[N]; 7 ll ans,g[N],f[N]; 8 char s[N<<1]; 9 void get_val(int k,int x,int mx){ 10 f[k]=min(f[k],f[x]+(g[x]-g[k])-(ll)dep[x]*(h[x]-h[k])+max(dep[k]-max(mx,dep[x]),0)); 11 } 12 int main(){ 13 scanf("%d%s",&n,s+1); 14 for(int i=1;i<=(n<<1);i++) 15 if (s[i]==')')st[0]--; 16 else{ 17 n0++; 18 if (st[0]){ 19 fa[n0]=st[st[0]]; 20 vis[st[st[0]]]=1; 21 } 22 st[++st[0]]=n0; 23 } 24 for(int i=2;i<=n;i++)dep[i]=dep[fa[i]]+1; 25 for(int i=n;i>1;i--){ 26 if (!vis[i])h[i]=1,mx[i]=g[i]=dep[i]; 27 h[fa[i]]+=h[i],g[fa[i]]+=g[i]; 28 int d=mx[i]; 29 if (mx[fa[i]]<d)son[fa[i]]=i,swap(mx[fa[i]],d); 30 cmx[fa[i]]=max(cmx[fa[i]],d); 31 } 32 for(int i=1;i<=n;i++){ 33 if ((i==1)||(son[fa[i]]!=i))top[i]=i; 34 else top[i]=top[fa[i]]; 35 V[top[i]].push_back(i); 36 } 37 memset(f,0x3f,sizeof(f)); 38 f[1]=0; 39 for(int i=1;i<=n;i++) 40 if (!V[i].empty()){ 41 if (i!=1){ 42 for(int j=0;j<V[i].size();j++)get_val(V[i][j],fa[i],mx[fa[i]]); 43 } 44 st[0]=0; 45 for(int j=0;j<V[i].size();j++){ 46 if (j)get_val(V[i][j],V[i][j-1],cmx[V[i][j-1]]); 47 while ((st[0])&&(cmx[st[st[0]]]<dep[V[i][j]]))st[0]--; 48 if (st[0])get_val(V[i][j],st[st[0]],cmx[st[st[0]]]); 49 while ((st[0])&&(cmx[st[st[0]]]<=cmx[V[i][j]]))st[0]--; 50 st[++st[0]]=V[i][j]; 51 } 52 } 53 ans=1e18; 54 for(int i=1;i<=n;i++)ans=min(ans,f[i]+g[i]-(ll)dep[i]*h[i]); 55 printf("%lld\n",ans); 56 return 0; 57 }