Prince and Princess
神题。。。。。。
本蒟蒻又刷到一个难题。不会做
查了题解后发现这是poj1904的变种。我只想说这题太牛逼了。
给出poj1904的题解:
https://www.cnblogs.com/zxndgv/archive/2011/08/06/2129333.html
其实还是增广路的思想。
由此引发到了强连通分量。
这让我想到了以前做过的一道题,记得是判断最小割中割边的唯一性。好像也是用到了强连通分量。
真的是非常巧妙,要对增广路有十分清晰的认识。同时,还应该充分利用完美匹配的条件。
真的是太强了。
至于这道题,我们因为刚开始不是完美匹配所以需要构造完美匹配:
见题解:https://blog.csdn.net/birdmanqin/article/details/99675178
大佬在这里说:新建虚拟节点使得其变为完美匹配,然后就可以利用完美匹配的性质了。
而这个新建节点,虚拟王子则是喜欢所有真实公主,虚拟公主则是要被所有真实王子喜欢
为什么要这样建立节点呢??
在poj1904中,大佬的分析中说了至关重要的一句话:
我们从xj开始找增广路,如果能找到一条增广路,那么根据完美匹配的性质这条增广路的末端一定是yi
因此我们只需要连接yi->xi就可以利用强连通分量找他了。
那么,来看看我们这里的情况:
假设是这样的一个图:
红色为王子,黑色为公主。
红2为虚构的王子。
我们发现在这种情况下求解出来的并不是正确答案。
因为在虚构节点这种连接方式之下,我们无法找到从2开始的一条增广路最终指向公主1
所以我们应该让虚拟王子只想所有真实公主,才能保证在通过虚拟接待你找增广路时一定能成功。
因为寻你节点毕竟是虚拟的啊,我们不用为其安排配对。
我说的很不好,实在是抱歉。希望能该你一点启发。
然后最后还有一个坑点:就是我们要判断是否有这条边使得王子瑜公主连接,具体见代码
#include<iostream> #include<algorithm> #include<cstdio> #include<queue> #include<stack> #include<vector> using namespace std; const int inf = 2e9; const int max_n = 2100; const int max_m = max_n * max_n; struct edge{ int to, next; }E[max_m << 1]; int head[max_n]; int cnt = 1; void add(int from, int to) { E[cnt].to = to;E[cnt].next = head[from]; head[from] = cnt++; } int leftTo[max_n], rightTo[max_n]; int distleft[max_n], distright[max_n]; int dist; bool seacherpath(int left_tot, int right_tot) { fill(distleft, distleft + max_n, -1); fill(distright, distright + max_n, -1); queue<int> que;dist = inf; for (int i = 1;i <= left_tot;++i) if (leftTo[i] == -1) que.push(i); while (!que.empty()) { int u = que.front();que.pop(); if (distleft[u] > dist)break; for (int i = head[u];i;i = E[i].next) { int v = E[i].to; if (distright[v] != -1)continue; distright[v] = distleft[u] + 1; if (rightTo[v] == -1)dist = distright[v]; else { distleft[rightTo[v]] = distright[v] + 1; que.push(rightTo[v]); } } }return dist != inf; } bool matchpath(int u) { for (int i = head[u];i;i = E[i].next) { int v = E[i].to; if (distright[v] != distleft[u] + 1)continue; if (distright[v] == dist && rightTo[v] != -1)continue; distright[v] = -1; if (rightTo[v] == -1 || matchpath(rightTo[v])) { leftTo[u] = v; rightTo[v] = u; return true; } }return false; }int HK(int left_tot, int right_tot) { fill(leftTo, leftTo + max_n, -1); fill(rightTo, rightTo + max_n, -1); int ans = 0; while (seacherpath(left_tot, right_tot)) { for (int i = 1;i <= left_tot;++i) if (leftTo[i] == -1 && matchpath(i)) ++ans; }return ans; } int dfn[max_n], low[max_n]; int tot0 = 0; stack<int> st; int color[max_n]; int colour = 0; void tarjan(int u) { dfn[u] = low[u] = ++tot0; st.push(u); for (int i = head[u];i;i = E[i].next) { int v = E[i].to; if (!dfn[v]) { tarjan(v); low[u] = min(low[u], low[v]); } else if (color[v] == 0)low[u] = min(low[u], dfn[v]); }if (low[u] != dfn[u])return; colour++; while (st.top() != u) { color[st.top()] = colour; st.pop(); }color[st.top()] = colour;st.pop(); } vector<int> nums[max_n]; bool mp[550][550]; int main() { int T;scanf("%d", &T);int n, m; for (int tcase = 1;tcase <= T;++tcase) { fill(head, head + max_n, 0);cnt = 1; fill(dfn, dfn + max_n, 0); fill(color, color + max_n, 0); for (int i = 0;i < 550;++i)for (int j = 0;j < 550;++j)mp[i][j] = false; colour = tot0 = 0; for (int i = 0;i < max_n;++i)nums[i].clear(); scanf("%d %d", &n, &m); for (int u = 1, len;u <= n;++u) { scanf("%d", &len); for (int i = 1, v;i <= len;++i) { scanf("%d", &v); add(u, v + n); mp[u][v] = true; } }HK(n, n + m); int tot = n + m; for (int u = 1;u <= n;++u) { if (leftTo[u] == -1) { leftTo[u] = ++tot; rightTo[tot] = u; add(leftTo[u], u); for (int i = 1;i <= n;++i) add(i, tot); }else add(leftTo[u], u); }for (int v = n + 1;v <= n + m;++v) { if (rightTo[v] == -1) { rightTo[v] = ++tot; leftTo[tot] = v; add(v, tot); for (int i = n + 1;i <= n + m;++i) add(tot, i); } }for (int i = 1;i <= tot;++i) if (!dfn[i]) tarjan(i); printf("Case #%d:\n", tcase); for (int i = 1;i <= n;++i) { for (int j = n + 1;j <= n + m;++j) if (mp[i][j - n] && color[i] == color[j]) nums[i].push_back(j - n); printf("%d", nums[i].size()); for (int j = 0;j < nums[i].size();++j) printf(" %d", nums[i][j]); printf("\n"); } } }
题单:https://vjudge.net/article/371