[BZOJ]1023 仙人掌图

[BZOJ]1023

题意：

　　给出一颗仙人掌，求仙人掌上两点间最短距离的最大值(即仙人掌的直径)。

题解：

　　树的直径可以用Dp来搞一搞，这个仙人掌其实也类似。
　　如果用$f[i]$表示从$i$向下扩展的最大距离，而且这颗仙人掌上没有环，我们就可以用这样的方法去更新答案：

$$ans = max(ans,f[u] + f[v] + 2)$$
　　其中$u$和$v$是$i$儿子中$f$值最大的和次大的。这里有一个小trick来实现，先用$f[u] + 1$来更新$f[i]$，再用$f[i] + f[v] + 1$来更新$ans$，这样出来一定是最大加次大。
　　接着考虑遇到环的情况。首先我们得知道这里有一个环其次才能想如何去处理，如果在$dfs$完$i$的子树后$i$的一个儿子不再是$i$的儿子滑稽，说明这里有一个环，而且环的顶点是$i$，我们对这个环进行处理。
　　在$dfs$的过程中如果遍历到一个点和其是在同一个环中且$i$不是顶点，就不更新$i$的$f$值。这样我们环里的点的$f$值就是连向了环之外的一条链的最长值。
　　像这样
于是在环上可能更新答案的就是$f[u] + f[v] + Dis(u,v)$，其中$Dis(u,v)$就是$u - v$，$u$和$v$是两个点。这个东西整理一下就是这个$f[u] + u + f[v] - v$，考虑固定$u$，用单调队列维护后面那个东西，把环里的点都跑一遍就行了。最后还要用环里的点更新一下顶点$i$的$f$值，因为$i$上面的点有可能要用到这个环。
　　意思基本就是这样，具体细节还是要看代码。

代码

#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

const int Maxn = 500010;

struct node {int to,next;}e[Maxn];

int deep[Maxn],fa[Maxn],a[Maxn],q[Maxn],low[Maxn],dfn[Maxn],h[Maxn],f[Maxn],n,m,ans,tot,tim;

int read(){
    int x=0; char ch=getchar();
    while (ch<'0' || ch>'9') ch=getchar();
    while (ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }
    return x;
}

void Add(int u,int v){e[++tot].to = v;e[tot].next = h[u];h[u] = tot;}

void Solve(int u,int v){
    int cnt = deep[v] - deep[u] + 1,he = 1,t = 1;q[1] = 1;
    for(int i = v;i != u;i = fa[i])a[cnt--] = f[i];
    a[1] = f[u];
    cnt = deep[v] - deep[u] + 1;
    for(int i = 1;i <= cnt;i++)a[i + cnt] = a[i];
    for(int i = 2;i <= cnt + (cnt >> 1);i++){
        if(i - q[he] > (cnt >> 1))he++;
        ans = max(ans,a[i] + a[q[he]] + i - q[he]);
        while(he <= t && a[q[t]] - q[t] <= a[i] - i)t--;
        q[++t] = i;
    }
    for(int i = 2;i <= cnt;i++)f[u] = max(f[u],a[i] + min(i - 1,cnt - i + 1));
}

void Dfs(int u){
    dfn[u] = low[u] = ++tim;
    for(int i = h[u];i;i = e[i].next){int v = e[i].to;
        if(v == fa[u])continue;
        if(!dfn[v]){
            deep[v] = deep[u] + 1;
            fa[v] = u;Dfs(v);
        }low[u] = min(low[u],low[v]);
        if(low[v] > dfn[u]){ans = max(ans,f[u] + f[v] + 1);f[u] = max(f[u],f[v] + 1);}
    }
    for(int i = h[u];i;i = e[i].next){int v = e[i].to;
        if(fa[v] != u && dfn[v] > dfn[u])Solve(u,v);
    }
}

int main(){
    n = read();m = read();
    for(int i = 1;i <= m;i++){
        int k,last = 0;k = read();
        for(int j = 1;j <= k;j++){
            int now;now = read();
            if(last != 0){Add(last,now);Add(now,last);}last = now;
        }
    }Dfs(1);
    printf("%d\n",ans);
    return 0;
}