强连通分量模板+例题（acwing368/洛谷3275）

模板：

//模板
int scc_cnt,id[N],_size[N];
int dfn[N],low[N];
int stk[N],top;
bool in_stk[N];
int timestamp;
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],low[j]);
        }
    }
    if(dfn[u]==low[u]){
        int y;
        scc_cnt++;
        do{
            y=stk[top--];
            in_stk[y]=false;
            id[y]=scc_cnt;
            _size[scc_cnt]++;
        }while(u!=y);
    }
}

例题：

银河中的恒星浩如烟海，但是我们只关注那些最亮的恒星。

我们用一个正整数来表示恒星的亮度，数值越大则恒星就越亮，恒星的亮度最暗是 1。

现在对于 N 颗我们关注的恒星，有 M 对亮度之间的相对关系已经判明。

你的任务就是求出这 N𝑁 颗恒星的亮度值总和至少有多大。

输入格式

第一行给出两个整数 N 和 M。

之后 M 行，每行三个整数 T,A,B,表示一对恒星 (A,B) 之间的亮度关系。恒星的编号从 1 开始。

如果 T=1，说明 A 和 B 亮度相等。
如果 T=2，说明 A 的亮度小于 B 的亮度。
如果 T=3，说明 A 的亮度不小于 B 的亮度。
如果 T=4，说明 A 的亮度大于 B 的亮度。
如果 T=5，说明 A 的亮度不大于 B 的亮度。

输出格式

输出一个整数表示结果。

若无解，则输出 −1。

数据范围

N≤100000,M≤100000

输入样例：

5 7 
1 1 2 
2 3 2 
4 4 1 
3 4 5 
5 4 5 
2 3 5 
4 5 1 

输出样例：

11

ac代码 ：

这里有两种思路正解是强连通分量线性o（n）过，但acwing数据很弱可以差分约束spfa过，不过洛谷hack数据过不了，这里两种写法都给出来。

//差分约束
#define _CRT_SECURE_NO_WARNINGS
#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#include<stack>
using namespace std;
typedef long long LL;
const int N = 1e5 + 10,M=3*N;
int h[M], w[M], ne[M], e[M];
int d[N];
int cnt[N];
int idx;
int n, m;
bool st[N];
void add(int a, int b, int c)
{
	e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx++;
}

bool spfa()
{
	memset(d, -0x3f, sizeof(d));
	stack<int> q;
	q.push(0);
	d[0] = 0;
	st[0] = true;
	while (q.size())
	{
		int t = q.top();
		q.pop();
		st[t] = false;
		for (int i = h[t]; i != -1; i = ne[i])
		{
			int j = e[i];
			if (d[j] < d[t] + w[i]) {
				d[j] = d[t] + w[i];
				cnt[j] = cnt[t] + 1;
				if (cnt[j] > n) {
					return false;
				}
				if (st[j])
				{
					continue;
				}
				q.push(j);
				st[j] = true;
			}
		}
	}
	return true;
}
int main()
{
	scanf("%d%d", &n, &m);
	memset(h, -1, sizeof(h));
	for (int i = 1; i <= m; i++)
	{
		int a, b, c;
		scanf("%d%d%d", &a, &b, &c);
		if (a == 1) {
			add(c, b, 0);
			add(b, c, 0);
		}
		else if (a == 2)
		{
			add(b, c, 1);
		}
		else if (a == 3)
		{
			add(c, b, 0);
		}
		else if (a == 4)
		{
			add(c, b, 1);
		}
		else {
			add(b, c, 0);
		}
	}
	for (int i = 1; i <= n; i++)
	{
		add(0, i, 1);
	}
	LL sum = 0;
	if (spfa())
	{
		for (int i = 1; i <= n; i++)
		{
			sum += d[i];
			
		}
		printf("%lld", sum);
	}
	else {
		printf("-1");
		
	}
}

//强连通分量
#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10,M=6*N;
int h[N],hs[N],w[M],ne[M],e[M],idx;
int stk[N],scc_cnt,se[N],top,id[N];
int dfn[N],low[N];
int timestamp;
bool in_stk[N];
int d[N];
int n,m;

void add(int h[],int a,int b,int c)
{
    w[idx]=c,e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}

void targan(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]){
            targan(j);
            low[u]=min(low[u],low[j]);
        }else if(in_stk[j])
        {
            low[u]=min(low[u],low[j]);
        }
    }
    if(dfn[u]==low[u])
    {
        scc_cnt++;
        int y;
        do{
            y=stk[top--];
            in_stk[y]=false;
            id[y]=scc_cnt;
            se[scc_cnt]++;
        }while(y!=u);
        
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    memset(h,-1,sizeof(h));
    memset(hs,-1,sizeof(hs));
    for(int i=1;i<=m;i++)
    {
        int t,a,b;
        scanf("%d%d%d",&t,&a,&b);
        if(t==1)
        {
            add(h,a,b,0);
            add(h,b,a,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 if(t==5)
        {
            add(h,a,b,0);
        }
    }
    for(int i=1;i<=n;i++){
        add(h,0,i,1);
    }
    targan(0);
    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];
            int b=id[k];
            if(a==b){
                if(w[j]){
                    printf("-1");
                    return 0;
                }
            }else {
                add(hs,a,b,w[j]);
            }
        }
    }
    for(int i=scc_cnt;i;i--){
        for(int j=hs[i];j!=-1;j=ne[j]){
            int k=e[j];
            d[k]=max(d[k],d[i]+w[j]);
        }
    }
    long long ans=0;
    for(int i=1;i<=scc_cnt;i++)
    {
        ans+=(long long)d[i]*se[i];
    }
    printf("%lld",ans);
    return 0;
}