有向图的强连通分量——银河

银河

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

我们用一个正整数来表示恒星的亮度，数值越大则恒星就越亮，恒星的亮度最暗是 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

题解：

其实这道题和差分约束糖果那道题一样，只不过那道题有可能被卡掉，而tarjan做法可以保证是线性的时间复杂度。
距离看下，让我们求最小值，那么就是缩点过后的最长路。什么时候无解呢？我们可以知道我们的环一定是存在于我们的强连通分量中的。一旦我们强连通分量中有一个值的边权大于0，那么就相当于有正环了。因为其他点可以到达他，他又可以继续到达其他点造成循环。所以强连通分量中的所有边权一定为0，也就是不等式中的相等的那种情况。所以无解的情况也就出来了。如果有解的话那么我就在DAG上面跑一边DP就可以了。

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int N = 100100, M = 600100;
int n, m;
int h[N], hs[N], e[M], ne[M], w[M], idx;
int dfn[N], low[N], timestamp;
bool in_stk[N];
int id[N], scc_cnt, s[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 ++ ;
}
stack<int> st;
void tarjan(int u)
{
    dfn[u] = low[u] = ++ timestamp;
    st.push(u), in_stk[u] = true;
    for (int i = h[u]; ~i; 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=st.top(); st.pop();
            in_stk[y] = false;
            id[y] = scc_cnt;
            s[scc_cnt] ++ ;
        } while (y != u);
    }
}
signed main() {
    scanf("%lld%lld", &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;
        scanf("%lld%lld%lld", &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; 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) puts("-1");
    else {
        for (int i = scc_cnt; i; i--) {
            for (int j = hs[i]; ~j; j = ne[j]) {
                int k = e[j];
                dist[k] = max(dist[k], dist[i] + w[j]);
            }
        }
        long long res = 0;
        for (int i = 1; i <= scc_cnt; i++) res += 1LL* dist[i] * s[i];
        printf("%lld\n", res);
    }
    return 0;
}