bzoj1023 [SHOI2008]cactus仙人掌图 & poj3567 Cactus Reloaded——求仙人掌直径

Posted 2020-12-03 zinn

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1023

　   　http://poj.org/problem?id=3567

仙人掌！直接模仿 lyd 的代码；

大概就是 tarjan 找环 + 单调队列优化 dp，然后缩环成链继续递归；

直接模仿着写的，感觉好妙啊；

不太明白边为什么要开成点数的4倍。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int const maxn=5e4+5;
int n,m,hd[maxn],ct,fa[maxn],f[maxn],a[maxn<<1],ans,q[maxn<<1],low[maxn],dfn[maxn],tim;
struct N{
    int to,nxt,w;
    N(int t=0,int n=0,int w=0):to(t),nxt(n),w(w) {}
}ed[maxn<<2];//<<2
void add(int x,int y){ed[++ct]=N(y,hd[x]); hd[x]=ct;}
void dp(int x,int y)
{
    int tp=0,i,l,r;
    for(i=y;i!=x;i=fa[i])a[++tp]=f[i];
    for(i=1,a[++tp]=f[x];i<m;i++)a[i+tp]=a[i];
    l=r=1;q[1]=1; int p=tp/2;
    for(int i=2;i<=tp+p;i++)//+p!
    {
        while(l<=r&&i-q[l]>p)l++;
        ans=max(ans,a[q[l]]+a[i]+i-q[l]);
        while(l<=r&&a[q[r]]+i-q[r]<=a[i])r--;//<=
        q[++r]=i;//
    }
    for(int i=1;i<tp;i++)//<
        f[x]=max(f[x],a[i]+min(i,tp-i));//min!!
}
void tarjan(int x)
{
    dfn[x]=low[x]=++tim;
    for(int i=hd[x],u;i;i=ed[i].nxt)
    {
        if((u=ed[i].to)==fa[x])continue;
        if(!dfn[u])
        {
            fa[u]=x; tarjan(u);
            low[x]=min(low[x],low[u]);
            if(dfn[x]<low[u])//不在一个环上 
            {
                ans=max(ans,f[x]+f[u]+1);
                f[x]=max(f[x],f[u]+1);
            }
        }
        else low[x]=min(low[x],dfn[u]);
    }
    for(int i=hd[x],u;i;i=ed[i].nxt)
        if(fa[u=ed[i].to]!=x&&dfn[x]<dfn[u])dp(x,u);//&&dfn[x]<dfn[u]  则x是入环点，y是环终点 
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1,x,y,k;i<=m;i++)
    {
        scanf("%d%d",&k,&x);
        for(int j=1;j<k;j++)
        {
            scanf("%d",&y);
            add(x,y); add(y,x);
            x=y;
        }
    }
    tarjan(1);
    printf("%d",ans);
    return 0;
}