传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1023
思路:类似树形DP记录一个f[i]表示最远(因为有环,所以这个定义是有一些限制条件的)
先用点双缩点,每个环的信息可以挂到最高点上
树上的差不多
对于环上的,从一边扫过去,因为dis(i,j)有单调性,用单调队列搞一搞即可
更详细的题解:http://ydcydcy1.blog.163.com/blog/static/21608904020131493113160/
#include
#include
#include
#include
#include
#define abs(a) (a<0?(-(a)):a)
const int maxn=100010,maxm=500010;
typedef unsigned int uint;
using namespace std;
int n,m,low[maxn],dfn[maxn],tim,sta[maxn],top,bnm[maxn],bcnt,ans,q[maxn],head,tail,a[maxn],f[maxn];bool vis[maxm];
vector bel[maxn],scc[maxn];
struct Tgraph{
int pre[maxm],now[maxn],son[maxm],tot;
void add(int a,int b){pre[++tot]=now[a],now[a]=tot,son[tot]=b;}
void ins(int a,int b){add(a,b),add(b,a);}
void dp1(int x,int fa){
int id=x-n,cnt=scc[id].size(),top=0,siz=0;
for (int i=0;icnt/2) head++;
if (head<=tail) ans=max(ans,f[a[i]]+f[a[q[head]]]+i-q[head]);
while (head<=tail&&f[a[q[tail]]]-q[tail]cnt/2) d=cnt-d;
f[x]=max(f[x],f[scc[id][i]]+d);
}
//printf("ans%d\n",ans);
}
void dp2(int x,int fa){
int max1=0,max2=0;
for (int y=now[x];y;y=pre[y]) if (son[y]!=fa){
max2=max(max2,f[son[y]]+(son[y]<=n));
if (max2>max1) swap(max1,max2);
}
f[x]=max1;ans=max(ans,max1+max2);
}
void tree_dp(int x,int fa){
for (int y=now[x];y;y=pre[y]) if (fa!=son[y]) tree_dp(son[y],x);
if (x>n) dp1(x,fa);else dp2(x,fa);
}
}g1,g2;
void tarjan(int x,int fa){
dfn[x]=low[x]=++tim,sta[++top]=x;
for (int y=g1.now[x];y;y=g1.pre[y]) if (!vis[y]){
int v=g1.son[y];vis[y]=vis[y^1]=1;
if (!dfn[v]) tarjan(v,x),low[x]=min(low[x],low[v]);
else low[x]=min(low[x],dfn[v]);
}
if (low[x]==dfn[x]&&fa) top--,g2.ins(x,fa),bnm[x]++,bnm[fa]++;
if (low[x]==dfn[fa]){
int xx;bcnt++;
bel[fa].push_back(bcnt),scc[bcnt].push_back(fa);
do{xx=sta[top--],bel[xx].push_back(bcnt),scc[bcnt].push_back(xx);}while (x!=xx);
}
}
void rebuild(){
/*puts("");for (int i=1;i<=bcnt;i++){printf("i%d ",i);
for (uint j=0;j=2)
for (uint j=0;j