1023: [SHOI2008]cactus仙人掌图(DP+单调队列优化)

这道题吗= =首先解决了我多年以来对仙人掌图的疑问，原来这种高大上的东西原来是这个啊= =

然后，看到这种题，首先必须的就是缩点= =

缩点完之后呢，变成在树上找最长路了= =直接树形dp了

那么那些环呢，就是一个环形dp了，可以先把它拆成一条链，然后注意到最长路径=max(f[i]+f[j]-dist(i,j))  拆成链的话dist（i,j）=i-j 然后就发现dist（i,j）有单调性，就可以用单调队列优化了= =

这样写就可以a了= =

ps1：今天发现有人给我留言了真开心QAQ 感觉自己写了这么久还是有人看到的QAQ 继续加油吧！！！

ps2：bzoj的808端口坏了现在上都得改网址真麻烦QAQ

ps3：刷了好久感觉没啥精神了这个月感觉刷不了多少了QAQ 所以众多STOI补番队的成员啊，这个月的占领头版计划就交给你们了QAQ

CODE：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
#define maxn 50100
#define maxm 20000100
struct edges{
 int to,next;
}edge[maxm];
int l,next[maxn];
int addedge(int x,int y){
 edge[++l]=(edges){y,next[x]};next[x]=l;
 edge[++l]=(edges){x,next[y]};next[y]=l;
 return 0;
}
int f[maxn],q[maxn*2],ans,que[maxn],pre[maxn],dep[maxn];
int dp(int y,int x ){
 int l=dep[y]-dep[x]+1;
 int t=y;
 while (l) {
  q[l--]=f[t];
  t=pre[t];
 }
 l=dep[y]-dep[x]+1;
 for (int i=1;i<=l;i++) q[i+l]=q[i];
 int h=0;
 t=1;
 for (int i=2;i<=(l>>1)+1;i++) {
  while (h>=t&&que[h]+q[que[h]]<=i+q[i]) h--;
  que[++h]=i;
 }
 int j=(l>>1)+1;
 for (int i=1;i<=(l<<1)-(l>>1);i++){
  while (h>=t&&que[t]<=i) t++;
  ans=max(ans,que[t]+q[que[t]]+q[i]-i);
  j++;
  while (h>=t&&que[h]+q[que[h]]<=j+q[j]) h--;
  que[++h]=j;
 }
 for (int i=2;i<=l;i++) f[x]=max(f[x],min(i-1,l-i+1)+q[i]);
 return 0;
}
int low[maxn],dfn[maxn],clo;
int dfs(int u,int fa){
 pre[u]=fa;
 dep[u]=dep[fa]+1;
 low[u]=dfn[u]=++clo;
 for (int i=next[u];i;i=edge[i].next)
  if (edge[i].to!=fa) {
   if (!low[edge[i].to]){
    dfs(edge[i].to,u);
    low[u]=min(low[u],low[edge[i].to]);
   }else low[u]=min(low[u],dfn[edge[i].to]);
   if (dfn[u]<low[edge[i].to]) {
    ans=max(ans,f[u]+f[edge[i].to]+1);
    f[u]=max(f[u],f[edge[i].to]+1);
   }
  }
 for (int i=next[u];i;i=edge[i].next)
  if (pre[edge[i].to]!=u&&dfn[u]<dfn[edge[i].to]) dp(edge[i].to,u);
 return 0;
}
int main(){
 int n,m;
 scanf("%d%d",&n,&m);
 for (int i=1;i<=m;i++) {
  int k,x,y;
  scanf("%d",&k);
  scanf("%d",&x);
  for (int j=2;j<=k;j++) {
   scanf("%d",&y);
   addedge(x,y);
   swap(x,y);
  }
 }
 dfs(1,0);
 printf("%d\n",ans); 
 return 0;
}