题目链接: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 "); 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 ", ret); 158 } 159 RT 0; 160 }