缩完之后合法顺序变少了
把上一节那张有环的图缩点:每个强连通分量捏成一个点,分量之间的边保留(去掉重边和自环)。
n01 里那张 DAG 有 4 种合法顺序。缩点之后的这张图有几种?两个数一起输出:
#include <algorithm>
#include <climits>
#include <iostream>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <string>
#include <utility>
#include <vector>
using namespace std;
using Edges = vector<pair<int, int>>;
const int N = 7;
// 七个构建任务,(a, b) 表示 a 做完才能做 b
const Edges DEP = {{0, 1}, {0, 2}, {1, 3}, {2, 3}, {3, 4}, {4, 5}, {4, 6}};
// 同一张图再加一条 5 → 3,于是有了环
const Edges G = {{0, 1}, {0, 2}, {1, 3}, {2, 3}, {3, 4}, {4, 5}, {5, 3}, {4, 6}};
vector<vector<int>> build(const Edges& edges, int n = N) {
vector<vector<int>> g(n);
for (auto [a, b] : edges) g[a].push_back(b);
for (auto& adj : g) sort(adj.begin(), adj.end());
return g;
}
// 第一遍:用手写栈做 DFS,记下每个点「完成」的先后
vector<int> order_by_finish(const vector<vector<int>>& g) {
vector<bool> seen(N, false);
vector<int> order;
for (int s = 0; s < N; s++) {
if (seen[s]) continue;
vector<pair<int, int>> st = {{s, 0}}; // (点, 下一条要看的边)
seen[s] = true;
while (!st.empty()) {
int u = st.back().first;
int i = st.back().second;
if (i < (int)g[u].size()) {
st.back().second++;
int v = g[u][i];
if (!seen[v]) {
seen[v] = true;
st.push_back({v, 0});
}
} else {
order.push_back(u);
st.pop_back();
}
}
}
return order;
}
// Kosaraju:第二遍在反图上、按完成顺序的逆序收点;返回排好序的分量列表
vector<vector<int>> kosaraju(const Edges& edges, bool second_reversed = true) {
vector<vector<int>> g = build(edges);
Edges rev;
for (auto [a, b] : edges) rev.push_back({b, a});
vector<vector<int>> rg = build(rev);
vector<int> order = order_by_finish(g);
const vector<vector<int>>& second = second_reversed ? rg : g;
vector<bool> seen(N, false);
vector<vector<int>> groups;
for (int k = N - 1; k >= 0; k--) {
int s = order[k];
if (seen[s]) continue;
vector<int> st = {s}, grp;
seen[s] = true;
while (!st.empty()) {
int u = st.back();
st.pop_back();
grp.push_back(u);
for (int v : second[u]) {
if (!seen[v]) {
seen[v] = true;
st.push_back(v);
}
}
}
sort(grp.begin(), grp.end());
groups.push_back(grp);
}
sort(groups.begin(), groups.end());
return groups;
}
vector<int> comp_of(const vector<vector<int>>& gs) {
vector<int> comp(N);
for (int i = 0; i < (int)gs.size(); i++) {
for (int u : gs[i]) comp[u] = i;
}
return comp;
}
// 缩点后的边:两头分量不同才保留,去掉重边,排好序
Edges condense(const Edges& edges, const vector<int>& comp) {
set<pair<int, int>> s;
for (auto [a, b] : edges) {
if (comp[a] != comp[b]) s.insert({comp[a], comp[b]});
}
return Edges(s.begin(), s.end());
}
// 枚举全部排列,数一数有几种不违反任何一条依赖
int topo_count(const Edges& edges, int n = N) {
vector<int> p(n);
iota(p.begin(), p.end(), 0);
int cnt = 0;
do {
vector<int> pos(n);
for (int i = 0; i < n; i++) pos[p[i]] = i;
bool ok = true;
for (auto [a, b] : edges) {
if (pos[a] > pos[b]) ok = false;
}
if (ok) cnt++;
} while (next_permutation(p.begin(), p.end()));
return cnt;
}
int main() {
vector<vector<int>> gs = kosaraju(G);
Edges ce = condense(G, comp_of(gs));
cout << topo_count(DEP) << "/" << topo_count(ce, gs.size()) << endl;
}
(本题用 g++ -std=c++17 -O0 编译。)
全部评论