题目:点击打开链接
题意:一条路径上必经的边为关键边,现在让你找一条路径,使得其关键边最多,输出最多的数量。
分析:(无向图缩点+图的重构+树的直径)如果一条路径上面有环,那么这个环的任意一条边都不是关键边。所以我们缩点,那么就算在一棵树上找最多的关键边,显然就是求直径。以前只学了有向图缩点的模板,现在发现其实无向图、有向图缩点差不多,无向图要防止回退到已经访问过的点。树的直径的求法大致如下,先任意找一个点,然后找出与它相距最远的一个点,然后再找出与这个找到的点相距最远的点,它们之间的距离就是树的直径,dfs两边就行了,我也不知道怎么证明,书上是怎么说的,感觉也是这样。
题意:一条路径上必经的边为关键边,现在让你找一条路径,使得其关键边最多,输出最多的数量。
分析:(无向图缩点+图的重构+树的直径)如果一条路径上面有环,那么这个环的任意一条边都不是关键边。所以我们缩点,那么就算在一棵树上找最多的关键边,显然就是求直径。以前只学了有向图缩点的模板,现在发现其实无向图、有向图缩点差不多,无向图要防止回退到已经访问过的点。树的直径的求法大致如下,先任意找一个点,然后找出与它相距最远的一个点,然后再找出与这个找到的点相距最远的点,它们之间的距离就是树的直径,dfs两边就行了,我也不知道怎么证明,书上是怎么说的,感觉也是这样。
我的代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<stdlib.h>
#include<vector>
#include<stack>
using namespace std;
const int N = 3e5+5;
stack<int> sta;
vector<int> mp[N],g[N];
int dfn[N],low[N],vis[N],num[N],degree[N];
int n,m,cnt,id,ans;
void init(){
cnt=0;
id=0;
ans=0;
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
memset(vis,0,sizeof(vis));
memset(num,0,sizeof(num));
memset(degree,0,sizeof(degree));
while(!sta.empty())
sta.pop();
for(int i=1;i<=n;i++)
mp[i].clear(),g[i].clear();
}
void tarjan(int x,int fa){
dfn[x]=low[x]=++id;
sta.push(x);
vis[x]=1;
for(int i=0;i<mp[x].size();i++){
int t=mp[x][i];
if(t==fa) continue;
if(!dfn[t]) {
tarjan(t,x);
low[x]=min(low[x],low[t]);
}else if(vis[t]) low[x]=min(low[x],dfn[t]);
}
if(dfn[x]==low[x]){
int tp;
cnt++;
do{
tp=sta.top();
vis[tp]=0;
num[tp]=cnt;
sta.pop();
}while(tp!=x);
}
}
void reb(){
for(int i=1;i<=n;i++)
for(int j=0;j<mp[i].size();j++)
if(num[i]!=num[mp[i][j]]) g[num[i]].push_back(num[mp[i][j]]);
// for(int i=1;i<=cnt;i++,cout<<endl)
// for(int j=0;j<g[i].size();j++)
// cout<<g[i][j]<<" ";
}
int dis[N];
void dfs(int x,int v) {
dis[x]=v,vis[x]=1;
for(int i=0;i<g[x].size();i++)
if(!vis[g[x][i]]) dfs(g[x][i],v+1);
}
void sv() {
for(int i=1;i<=n;i++)
if(!dfn[i])
tarjan(i,0);
reb();
memset(vis,0,sizeof(vis));
dfs(1,0);
int st=1,mx=dis[1];
for(int i=2;i<=cnt;i++)
if(dis[i]>mx) mx=dis[i],st=i;
memset(vis,0,sizeof(vis));
dfs(st,0);
st=1,mx=dis[1];
for(int i=2;i<=cnt;i++)
if(dis[i]>mx) mx=dis[i],st=i;
cout<<mx<<endl;
}
int main(){
while(cin>>n>>m){
init();
for(int i=0;i<m;i++){
int a,b;
cin>>a>>b;
mp[a].push_back(b),mp[b].push_back(a);
}
sv();
}
}
大佬的代码:
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=600010;
int From[maxn],Laxt[maxn],To[maxn<<2],Next[maxn<<2],cnt;
int low[maxn],dfn[maxn],times,q[maxn],head,scc_cnt,scc[maxn];
vector<int>G[maxn];
int dis[maxn],S,T,ans;
void add(int u,int v)
{
Next[++cnt]=Laxt[u]; From[cnt]=u;
Laxt[u]=cnt; To[cnt]=v;
}
void tarjan(int u,int fa)
{
dfn[u]=low[u]=++times;
q[++head]=u;
for(int i=Laxt[u];i;i=Next[i]){
if(To[i]==fa) continue;
if(!dfn[To[i]]) {
tarjan(To[i],u);
low[u]=min(low[u],low[To[i]]);
}
else low[u]=min(low[u],dfn[To[i]]);
}
if(low[u]==dfn[u]){
scc_cnt++;
while(true){
int x=q[head--];
scc[x]=scc_cnt;
if(x==u) break;
}
}
}
void dfs(int u,int f)
{
dis[u]=dis[f]+1;
for(int i=0;i<G[u].size();i++){
if(G[u][i]!=f) dfs(G[u][i],u);
}
}
int main()
{
int N,M,u,v,i,j;
scanf("%d%d",&N,&M);
for(i=1;i<=M;i++){
scanf("%d%d",&u,&v);
add(u,v); add(v,u);
}
tarjan(1,0);
for(i=1;i<=N;i++){
for(j=Laxt[i];j;j=Next[j]){
if(scc[i]!=scc[To[j]])
G[scc[i]].push_back(scc[To[j]]);
}
}
dfs(1,0);
for(i=1;i<=scc_cnt;i++) if(dis[i]>dis[S]) S=i;
dfs(S,0);
for(i=1;i<=scc_cnt;i++) ans=max(ans,dis[i]-1);
printf("%d\n",ans);
return 0;
}