模板
void tarjan(int u){
dfn[u]=low[u]=++timestamp;
stk[++top]=u,in_stk[u]=true;
for(int i=h[u];i!=-1;i=ne[i]){
int j=e[i];
if(!dfn[j]){
tarjan(j);
low[u]=min(low[u],low[j]);
}
else if(in_stk[j])low[u]=min(low[u],dfn[j]);
}
if(dfn[u]==low[u]){
++ scc_cnt;
int y;
do{
y=stk[top--];
in_stk[y]=false;
id[y]=scc_cnt;
si[scc_cnt]++;
}while(y!=u);
}
}
二.练习题
1.受欢迎的牛(模板题)
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 10010,M=50010;
int h[N],e[M],ne[M],idx=0;
int scc_cnt=0,top=0,timestamps=0;
int stk[N],dfn[N],low[N],dout[N];
bool in_stk[N];
int id[N],nums_size[N];
int n,m;
void add(int a,int b)
{
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void tarjan(int u)
{
dfn[u]=low[u]=++timestamps;
stk[++top]=u;
in_stk[u]=true;
for(int i=h[u];i!=-1;i=ne[i])
{
int j=e[i];
if(!dfn[j])
{
tarjan(j);
low[u]=min(low[u],low[j]);
}
else if(in_stk[j])
{
low[u]=min(low[u],dfn[j]);
}
}
if(dfn[u]==low[u])
{
++scc_cnt;
int y;
do{
y=stk[top--];
in_stk[y]=false;
id[y]=scc_cnt;
nums_size[scc_cnt]++;
}while(y!=u);
}
}
int main()
{
cin>>n>>m;
memset(h,-1,sizeof h);
while (m -- ){
int a,b;
cin>>a>>b;
add(a,b);
}
for(int i=1;i<=n;i++)
{
if(!dfn[i])tarjan(i);
}
for(int i=1;i<=n;i++)
{
for(int j=h[i];j!=-1;j=ne[j])
{
int k=e[j];
int a=id[i],b=id[k];
if(a!=b)
{
dout[a]++;
}
}
}
int zeros=0,res=0;
for(int i=1;i<=scc_cnt;i++)
{
if(!dout[i])
{
zeros++;
res+=nums_size[i];
if(zeros>1)
{
res=0;
break;
}
}
}
cout<<res<<endl;
}
2.学校网络(模板题)
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 110,M=10010;
int h[N],e[M],ne[M],idx=0;
int dfn[N],low[N],stk[N];
int n,m,scc_cnt=0,timestamp=0,top=0;
bool in_stk[N];
int id[N];
int din[N],dout[N];
void add(int a,int b){
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void tarjan(int u)
{
dfn[u]=low[u]=++timestamp;
stk[++top]=u;
in_stk[u]=true;
for(int i=h[u];i!=-1;i=ne[i])
{
int j=e[i];
if(!dfn[j])
{
tarjan(j);
low[u]=min(low[j],low[u]);
}
else if(in_stk[j])
{
low[u]=min(dfn[j],low[u]);
}
}
if(dfn[u]==low[u])
{
++scc_cnt;
int y;
do{
y=stk[top--];
in_stk[y]=false;
id[y]=scc_cnt;
}while(y!=u);
}
}
int main()
{
cin>>n;
memset(h,-1,sizeof h);
int x;
for(int i=1;i<=n;i++)
{
while(cin>>x,x)add(i,x);
}
for(int i=1;i<=n;i++)
{
if(!dfn[i])tarjan(i);
}
for(int i=1;i<=n;i++)
{
for(int j=h[i];j!=-1;j=ne[j])
{
int k=e[j];
int a=id[i],b=id[k];
if(a!=b)
{
din[b]++;
dout[a]++;
}
}
}
int res1=0,res2=0;
for(int i=1;i<=scc_cnt;i++)
{
if(!din[i])res1++;
if(!dout[i])res2++;
}
cout<<res1<<endl;
if(scc_cnt==1)cout<<0<<endl;
else
{
cout<<max(res1,res2)<<endl;
}
}
3.最大半连通子图(tarjan+dp+topsort+判重)
#include <iostream>
#include <cstring>
#include <algorithm>
#include <unordered_set>
using namespace std;
typedef long long LL;
typedef pair<int, int> PII;
const int N = 100010,M=2000010;
int mod;
int h[N],hs[N],e[M],ne[M],idx=0;
int dfn[N],low[N],nums[N],id[N],stk[N];
bool in_stk[N];
int scc_cnt=0,times=0,top=0;
int n,m;
int f[N],g[N];
void add(int h[],int a,int b)
{
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void tarjan(int u)
{
dfn[u]=low[u]=++times;
stk[++top]=u;
in_stk[u]=true;
for(int i=h[u];i!=-1;i=ne[i])
{
int j=e[i];
if(!dfn[j])
{
tarjan(j);
low[u]=min(low[u],low[j]);
}
else if(in_stk[j])low[u]=min(low[u],dfn[j]);
}
if(dfn[u]==low[u])
{
++scc_cnt;
int y;
do{
y=stk[top--];
in_stk[y]=false;
id[y]=scc_cnt;
nums[scc_cnt]++;
}while(u!=y);
}
}
int main()
{
memset(h, -1, sizeof h);
memset(hs,-1,sizeof hs);
scanf("%d%d%d",&n,&m,&mod);
while(m--){
int a,b;
scanf("%d%d",&a,&b);
add(h,a,b);
}
for(int i=1;i<=n;i++)
if(!dfn[i])
tarjan(i);
unordered_set<LL>S;
for(int i=1;i<=n;i++)
for(int j=h[i];j!=-1;j=ne[j]){
int k=e[j];
int a=id[i],b=id[k];
LL hash=a*1000000ll+b;
if(a!=b&&!S.count(hash)){
add(hs,a,b);
S.insert(hash);
}
}
for(int i=scc_cnt;i;i--){
if(!f[i]){
f[i]=nums[i];
g[i]=1;
}
for(int j=hs[i];j!=-1;j=ne[j]){
int k=e[j];
if(f[k]<f[i]+nums[k]){
f[k]=f[i]+nums[k];
g[k]=g[i];
}
else if(f[k]==f[i]+nums[k])
g[k]=(g[k]+g[i])%mod;
}
}
int maxf=0,sum=0;
for(int i=1;i<=scc_cnt;i++)
if(f[i]>maxf){
maxf=f[i];
sum=g[i];
}
else if(f[i]==maxf)sum=(sum+g[i])%mod;
cout<<maxf<<endl;
cout<<sum<<endl;
}
4.银河(可差分约束但spfa可能tle,tarjan+topsort为线性时间)
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
typedef pair<int, int> PII;
const int N = 100010,M = 600010;
int n,m;
int h[N],hs[N],e[M],ne[M],w[M],idx=0;
int dfn[N],low[N],timestamp;
int stk[N],top;
bool in_stk[N];
int id[N],scc_cnt,nums[N];
int dist[N];
void add(int h[],int a,int b,int c)
{
e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}
void tarjan(int u)
{
dfn[u]=low[u]=++timestamp;
stk[++top]=u,in_stk[u]=true;
for(int i=h[u];i!=-1;i=ne[i])
{
int j=e[i];
if(!dfn[j])
{
tarjan(j);
low[u]=min(low[u],low[j]);
}
else if(in_stk[j])
{
low[u]=min(low[u],dfn[j]);
}
}
if(dfn[u]==low[u])
{
++scc_cnt;
int y;
do{
y=stk[top--];
in_stk[y]=false;
id[y]=scc_cnt;
nums[scc_cnt]++;
}while(y!=u);
}
}
int main()
{
cin>>n>>m;
memset(h, -1, sizeof h);
memset(hs,-1,sizeof hs);
//连接超级源点
for(int i=1;i<=n;i++)add(h,0,i,1);
while (m -- )
{
int t,a,b;
cin>>t>>a>>b;
if(t==1)add(h,b,a,0),add(h,a,b,0);
else if(t==2)add(h,a,b,1);
else if(t==3)add(h,b,a,0);
else if(t==4)add(h,b,a,1);
else add(h,a,b,0);
}
tarjan(0);
bool success=true;
for(int i=0;i<=n;i++)
{
for(int j=h[i];j!=-1;j=ne[j])
{
int k=e[j];
int a=id[i],b=id[k];
if(a==b)
{
if(w[j]>0)
{
success=false;
break;
}
}
else add(hs,a,b,w[j]);
}
if(!success)break;
}
if(!success)cout<<-1<<endl;
else
{
for(int i=scc_cnt;i;i--)
{
for(int j=hs[i];j!=-1;j=ne[j])
{
int k=e[j];
dist[k]=max(dist[k],dist[i]+w[j]);
}
}
LL res=0;
for(int i=1;i<=scc_cnt;i++)res+=(LL)dist[i]*nums[i];
cout<<res<<endl;
}
}