【BZOJ4484】【JSOI2015】最小表示（拓扑排序，bitset）

Description

对于一个N个点（每个点从1到N编号），M条边的有向图，JYY发现，如果从图中删去一些边，那么原图的连通性会发生改变；而也有一些边，删去之后图的连通性并不会发生改变。
JYY想知道，如果想要使得原图任意两点的连通性保持不变，我们最多能删掉多少条边呢？
为了简化一下大家的工作量，这次JYY保证他给定的有向图一定是一个有向无环图。

$N\le30,000,M\le100,000$

---

Solution

先按拓扑序将每个点可以由哪些点到达用bitset存起来，将边反向后再搞一次得到每个点可以到达哪些点。

考虑两个点$u,v$之间的边如果可以被删除，那么一定存在另一条经过其他点的路径。
那么我们找到$u$可以到达的点和可以到达$v$的点的交集判断即可。

---

Code

/**************************************
 * Au: Hany01
 * Prob: [BZOJ4484][JSOI2015] 最小表示
 * Date: Apr 4th, 2018
 * Email: hany01@foxmail.com
**************************************/

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;
typedef pair<int, int> PII;
typedef vector<int> VI;
#define File(a) freopen(a".in", "r", stdin), freopen(a".out", "w", stdout)
#define rep(i, j) for (register int i = 0, i##_end_ = j; i < i##_end_; ++ i)
#define For(i, j ,k) for (register int i = (j), i##_end_ = (k); i <= i##_end_; ++ i)
#define Fordown(i, j, k) for (register int i = (j), i##_end_ = (k); i >= i##_end_; -- i)
#define Set(a, b) memset(a, b, sizeof(a))
#define SZ(a) ((int)(a.size()))
#define ALL(a) a.begin(), a.end()
#define pb(a) push_back(a)
#define mp(a, b) make_pair(a, b)
#define INF (0x3f3f3f3f)
#define INF1 (2139062143)
#define Mod (1000000007)
#define y1 wozenmezhemecaia 
#ifdef hany01
#define debug(...) fprintf(stderr, __VA_ARGS__)
#else
#define debug(...)
#endif

template<typename T> inline bool chkmax(T &a, T b) { return a < b ? a = b, 1 : 0; }
template<typename T> inline bool chkmin(T &a, T b) { return b < a ? a = b, 1 : 0; }

inline int read() {
    register char c_; register int _, __;
    for (_ = 0, __ = 1, c_ = getchar(); !isdigit(c_); c_ = getchar()) if (c_ == '-')  __ = -1;
    for ( ; isdigit(c_); c_ = getchar()) _ = (_ << 1) + (_ << 3) + (c_ ^ 48);
    return _ * __;
}

const int maxn = 30005, maxm = 100005;

int n, deg[2][maxn], v[2][maxm], nex[2][maxm], beg[2][maxn], e[2];
bitset<maxn> bs[2][maxn];

inline void add(int uu, int vv, int id) { v[id][++ e[id]] = vv, nex[id][e[id]] = beg[id][uu], beg[id][uu] = e[id], ++ deg[id][vv]; }

void Work(int id)
{
    static queue<int> q;
    For(i, 1, n) {
        if (!deg[id][i]) q.push(i);
        bs[id][i][i] = 1;
    }
    while (!q.empty()) {
        register int u = q.front(); q.pop();
        for (register int i = beg[id][u]; i; i = nex[id][i]) {
            bs[id][v[id][i]] |= bs[id][u];
            if (!-- deg[id][v[id][i]]) q.push(v[id][i]);
        }
    }
}

int main()
{
#ifdef hany01
    File("bzoj4484");
#endif

    static int cnt = 0, m, uu, vv;

    n = read(), m = read();
    For(i, 1, m) uu = read(), vv = read(), add(uu, vv, 0), add(vv, uu, 1);

    Work(0), Work(1);

    For(u, 1, n) for (register int i = beg[0][u]; i; i = nex[0][i])
        if ((bs[1][u] & bs[0][v[0][i]]).count() > 2) ++ cnt;
    printf("%d\n", cnt);

    return 0;
}