1. 程式人生 > >bzoj 1023: [SHOI2008]cactus仙人掌圖【tarjan+dp+單調隊列】

bzoj 1023: [SHOI2008]cactus仙人掌圖【tarjan+dp+單調隊列】

最短距離 esp main 仙人掌 i++ clu pac tarjan 處理

本來想先求出點雙再一個一個處理結果寫了很長發現太麻煩
設f[u]為u點向下的最長鏈
就是再tarjan的過程中,先照常處理,用最長兒子鏈和次長兒子鏈更新按ans,然後處理以這個點為根的環,也就是這個點是dfs第一次到這個環訪問到的點
環用來更新ans的是兒子鏈+到根的一段,這個直接for一邊就行,還有就是一個兒子鏈+環上的一段+另一個兒子鏈,這個把環復制一遍然後單調隊列掃即可
註意距離的定義是最短距離!

#include<iostream>
#include<cstdio>
using namespace std;
const int N=100005;
int n,m,h[N],cnt,fa[N],de[N],f[N],s[N],dfn[N],low[N],dft,ans,a[N],q[N];
struct qwe
{
    int ne,to;
}e[N*10];
int read()
{
    int r=0,f=1;
    char p=getchar();
    while(p>'9'||p<'0')
    {
        if(p=='-')
            f=-1;
        p=getchar();
    }
    while(p>='0'&&p<='9')
    {
        r=r*10+p-48;
        p=getchar();
    }
    return r*f;
}
void add(int u,int v)
{
    cnt++;
    e[cnt].ne=h[u];
    e[cnt].to=v;
    h[u]=cnt;
}
void dfs(int u,int fat)
{
    dfn[u]=low[u]=++dft;
    fa[u]=fat;
    de[u]=de[fat]+1;
    for(int i=h[u];i;i=e[i].ne)
        if(e[i].to!=fat)
        {
            if(!dfn[e[i].to])
            {
                dfs(e[i].to,u);
                low[u]=min(low[u],low[e[i].to]);
            }
            else 
                low[u]=min(low[u],dfn[e[i].to]);
            if(low[e[i].to]>dfn[u])
            {
                ans=max(ans,f[u]+f[e[i].to]+1);
                f[u]=max(f[u],f[e[i].to]+1);
            }
        }
    for(int j=h[u];j;j=e[j].ne)
        if(fa[e[j].to]!=u&&dfn[e[j].to]>dfn[u])
        {
            int tot=de[e[j].to]-de[u]+1,l=1,r=1;
            q[1]=1;
            for(int x=e[j].to;x!=u;x=fa[x])
                a[tot--]=f[x];
            a[tot]=f[u];
            tot=de[e[j].to]-de[u]+1;
            for(int i=1;i<=tot;i++)
                a[i+tot]=a[i];
            for(int i=2;i<=2*tot;i++)
            {
                while(l<=r&&i-q[l]>tot/2)
                    l++;
                if(l<=r)
                    ans=max(ans,a[i]+i+a[q[l]]-q[l]);
                while(l<=r&&a[q[r]]-q[r]<=a[i]-i)
                    r--;
                q[++r]=i;
            }
            for(int i=2;i<=tot;i++)
                f[u]=max(f[u],a[i]+min(i-1,tot-i+1));
        }
}
int main()
{
    n=read(),m=read();
    for(int i=1;i<=m;i++)
    {
        int k=read(),la=read(),x;
        for(int j=2;j<=k;j++)
        {
            x=read();
            add(la,x),add(x,la);
            la=x;
        }
    }
    dfs(1,0);
    printf("%d\n",ans);
    return 0;
}

bzoj 1023: [SHOI2008]cactus仙人掌圖【tarjan+dp+單調隊列】