有向图的强连通分量-----tarjan算法

模板

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;
    }
}