hdu 4612 Warm up - 双连通 + 树形DP

/*hdu 4612 Warm up - 双连通 + 树形DP题意很简单通过双连通缩图，构成一棵树，然后求最长链，连接首尾即可。剩下的连通块就是所求https://www.byvoid.com/blog/biconnect该算法是R.Tarjan发明的。对图深度优先搜索，定义DFS(u)为u在搜索树（以下简称为树）中被遍历到的次序号。定义Low(u)为u或u的子树中能通过非父子边追溯到的最早的节点，即DFS序号最小的节点。根据定义，则有：Low(u)=Min { DFS(u) DFS(v) (u,v)为后向边(返祖边) 等价于 DFS(v)<DFS(u)且v不为u的父亲节点 Low(v) (u,v)为树枝边(父子边) }一个顶点u是割点，当且仅当满足(1)或(2) (1) u为树根，且u有多于一个子树。 (2) u不为树根，且满足存在(u,v)为树枝边(或称父子边，即u为v在搜索树中的父亲)，使得DFS(u)<=Low(v)。一条无向边(u,v)是桥，当且仅当(u,v)为树枝边，且满足DFS(u)<Low(v)。*/#pragma comment(linker,"/STACK:102400000,102400000")#include<stdio.h>  #include<string.h>  const int N=200000+10,M=1000000+10;int head[N],low[N],fang[N],du[N];  struct node  {      int v,next;  }e[M*2];//这地儿老郁闷了，没有注意到无向图存边需要的是边数量的2倍，之前老是WA  int yong,n,r,m,dp[N][2],result;  bool visit[M];int min(int a,int b)  {      if(a<b)          return a;      return b;  }  void tarjan(int u)  //之前的通过传递父节点来规避直接走回去的情况存在漏洞，当1-2,2-3,3-4,4-2时，4-2-1就退回去了{      int i,j;      fang[u]=low[u]=yong++;  dp[u][0] = dp[u][1] = 0;    for(i=head[u];i!=-1;i=e[i].next)      {          j=e[i].v;          if(!visit[i >> 1]){  visit[i >> 1] = true;            if(fang[j]==0)//跟强连通一样这往下              {                  tarjan(j);  result += fang[u] < low[j];//统计连通块，比实际数目少一个，就是回溯的时候的最后一个int temp = dp[j][0] + (fang[u] < low[j]);//dp[i][0]第i个节点的子树中最长的链，dp[i][1]第i个节点的子树中第二长的链，                if(temp > dp[u][0]){                    dp[u][1] = dp[u][0];                    dp[u][0] = temp;                }else if(temp > dp[u][1]){                    dp[u][1] = temp;                }                low[u]=min(low[u],low[j]);              }else low[u]=min(low[u],fang[j]);          }      }  }  void adde(int u,int v){e[yong].v=v;e[yong].next=head[u];head[u]=yong++;}inline int max(int a,int b){if(a>b) return a;return b;}int main(){int i,a,b;while(scanf("%d%d",&n,&m),n+m){yong=0;  result=0;        memset(fang,0,sizeof(fang));          memset(low,0,sizeof(low));          memset(head,-1,sizeof(head));  for(i=0;i<m;++i){scanf("%d%d",&a,&b);adde(a,b);adde(b,a);}yong=1;  memset(visit, 0, sizeof(visit));        tarjan(1);int temp = 0;for(i = 1; i <=n; ++ i){            temp = max(temp, dp[i][0] + dp[i][1]);//+的时候没有算当前点所在的块，同时因为result也少算一个，扯平了        }        printf("%d\n", result - temp);}return 0;}