[HDOJ4635]Strongly connected（强连通分量，缩点）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4635

题意：给一张图，问最多往这张图上加多少条边，使这张图仍然无法成为一个强连通图。

起初是先分析样例可以知道，一个强连通分量里应当加边加至一个完全图，这样对整个图的连通贡献是没有的。然后把每个强连通分量里的边数扩展至num-out-1，但是发现连通分量之间的关系不好处理，所以反着做。

考虑一张有向完全图的边数为n*(n-1)，再去掉现在已经有的边数m，那么不考虑无法成为强连通图这一要求，我们还可以加n*(n-1)-m条边。现在问题变成了最少去掉多少条边，这个图不是一个强连通图。

我们不希望浪费任何一个强连通分量，所以我们最终的情况应当是有一个强连通分量是无法到达其他强连通分量，并且其他强连通分量早已连接成一整个强连通分量。要保证这种边尽量多，也就是保证排除的点尽可能少，我们枚举其中的一个强连通分量（因为单个分量的点数一定小于多个连通分量点数相加后的最小值）。接下来在图上找入度或者出度为0的连通分量，从这种分量出或者入是不会影响连通性的，所以得解。

  1 /*
  2 ━━━━━┒ギリギリ♂ eye！
  3 ┓┏┓┏┓┃キリキリ♂ mind！
  4 ┛┗┛┗┛┃＼○／
  5 ┓┏┓┏┓┃ /
  6 ┛┗┛┗┛┃ノ)
  7 ┓┏┓┏┓┃
  8 ┛┗┛┗┛┃
  9 ┓┏┓┏┓┃
 10 ┛┗┛┗┛┃
 11 ┓┏┓┏┓┃
 12 ┛┗┛┗┛┃
 13 ┓┏┓┏┓┃
 14 ┃┃┃┃┃┃
 15 ┻┻┻┻┻┻
 16 */
 17 #include <algorithm>
 18 #include <iostream>
 19 #include <iomanip>
 20 #include <cstring>
 21 #include <climits>
 22 #include <complex>
 23 #include <fstream>
 24 #include <cassert>
 25 #include <cstdio>
 26 #include <bitset>
 27 #include <vector>
 28 #include <deque>
 29 #include <queue>
 30 #include <stack>
 31 #include <ctime>
 32 #include <set>
 33 #include <map>
 34 #include <cmath>
 35 using namespace std;
 36 #define fr first
 37 #define sc second
 38 #define cl clear
 39 #define BUG puts("here!!!")
 40 #define W(a) while(a--)
 41 #define pb(a) push_back(a)
 42 #define Rint(a) scanf("%d", &a)
 43 #define Rll(a) scanf("%lld", &a)
 44 #define Rs(a) scanf("%s", a)
 45 #define Cin(a) cin >> a
 46 #define FRead() freopen("in", "r", stdin)
 47 #define FWrite() freopen("out", "w", stdout)
 48 #define Rep(i, len) for(int i = 0; i < (len); i++)
 49 #define For(i, a, len) for(int i = (a); i < (len); i++)
 50 #define Cls(a) memset((a), 0, sizeof(a))
 51 #define Clr(a, x) memset((a), (x), sizeof(a))
 52 #define Full(a) memset((a), 0x7f7f, sizeof(a))
 53 #define lp p << 1
 54 #define rp p << 1 | 1
 55 #define pi 3.14159265359
 56 #define RT return
 57 typedef long long LL;
 58 typedef long double LD;
 59 typedef unsigned long long ULL;
 60 typedef pair<int, int> pii;
 61 typedef pair<string, int> psi;
 62 typedef map<string, int> msi;
 63 typedef vector<int> vi;
 64 typedef vector<LL> vl;
 65 typedef vector<vl> vvl;
 66 typedef vector<bool> vb;
 67 
 68 typedef struct Edge {
 69     int u;
 70     int v;
 71     int next;
 72     Edge() { next = -1; }
 73 }Edge;
 74 const int maxn = 100100;
 75 const int maxm = 100100;
 76 int n, m;
 77 int head[maxn], ecnt;
 78 Edge edge[maxm];
 79 int bcnt, dindex;
 80 int dfn[maxn], low[maxn];
 81 int stk[maxn], top;
 82 int belong[maxn], num[maxn];
 83 int in[maxn], out[maxn];
 84 bool instk[maxn];
 85 
 86 void init() {
 87     Clr(head, -1);
 88     Cls(edge); Cls(instk);
 89     Cls(dfn); Cls(low); Cls(num);
 90     Cls(belong); Cls(in); Cls(out);
 91     ecnt = top = bcnt = dindex = 0;
 92 }
 93 
 94 void adde(int uu, int vv) {
 95     edge[ecnt].u = uu;
 96     edge[ecnt].v = vv;
 97     edge[ecnt].next = head[uu];
 98     head[uu] = ecnt++;
 99 }
100 
101 void tarjan(int u) {
102     int v = u;
103     dfn[u] = low[u] = ++dindex;
104     stk[++top] = u;
105     instk[u] = 1;
106     for(int i = head[u]; ~i; i=edge[i].next) {
107         v = edge[i].v;
108         if(!dfn[v]) {
109             tarjan(v);
110             low[u] = min(low[u], low[v]);
111         }
112         else if(instk[v]) low[u] = min(low[u], dfn[v]);
113     }
114     if(dfn[u] == low[u]) {
115         bcnt++;
116         do {
117             v = stk[top--];
118             instk[v] = 0;
119             belong[v] = bcnt;
120         } while(v != u);
121     }
122 }
123 
124 int main() {
125     // FRead();
126     int T, _ = 1, u, v;
127     Rint(T);
128     W(T) {
129         printf("Case %d: ", _++);
130         init();
131         Rint(n); Rint(m);
132         Rep(i, m) {
133             Rint(u); Rint(v);
134             adde(u, v);
135         }
136         For(i, 1, n+1) if(!dfn[i]) tarjan(i);
137         if(bcnt == 1) {
138             printf("-1\n");
139             continue;
140         }
141         For(u, 1, n+1) {
142             num[belong[u]]++;
143             for(int i = head[u]; ~i; i=edge[i].next) {
144                 int v = edge[i].v;
145                 if(belong[u] != belong[v]) {
146                     in[belong[v]]++; out[belong[u]]++;
147                 }
148             }
149         }
150         int sum = n*(n-1)-m;
151         int ret = 0;
152         For(i, 1, bcnt+1) {
153             if(in[i] == 0 || out[i] == 0) {
154                 ret = max(ret, sum-num[i]*(n-num[i]));
155             }
156         }
157         printf("%d\n", ret);
158     }
159     RT 0;
160 }

 

转载于:https://www.cnblogs.com/kirai/p/5523407.html