又是一个强连通分量,首先我们要明确,加了边之后的图是一个什么图,一定存在两个部分,假设为x部分和y部分,那样的话,只要只存在从x出发到y得边就能保证这个图不是一个强连通图,同时边最多。x部分和y部分一定是完全图,因为我们要保证边尽量多嘛,那么这个图总有多少条边,x部分一共有x*(x-1)条边,y部分一共有y*(y-1)条边,然后x和y之间就是x个点,每个点往y中所有的点发射一条边,一共有x*y条边,那么整理以后所有边数=n^2-n-x*y,同时x+y=n,根据基本不等式,x和y相差越大,那么乘积也就越小,前面的所有都是定值,那么明确了一个目的,x部分和y部分必须要保证某一个部分的点越小越好。
首先,x部分和y部分一定是强连通分量,那么好说了,只有入度为0或者出度为0的强连通才能够当x部分或者y部分,我们就统计所有入度或者出度为0的强连通分量,找到点数最小的一个当成x部分,剩下所有的点当成y部分,公式计算出所有的边减去m就是答案。
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <stack>
#include <iostream>
using namespace std;
typedef long long LL;
const int maxn = 100010;
int pre[maxn], sccno[maxn], scc_cnt, dfs_clock, lowlink[maxn], idg[maxn], odg[maxn];
vector<int> G[maxn], scc[maxn];
stack<int> S;
void dfs(int u){
pre[u] = lowlink[u] = ++dfs_clock;
S.push(u);
for (int i = 0; i < G[u].size(); i++){
int v = G[u][i];
if (!pre[v]){
dfs(v);
lowlink[u] = min(lowlink[u], lowlink[v]);
}
else if (!sccno[v])
lowlink[u] = min(lowlink[u], pre[v]);
}
if (lowlink[u] == pre[u]){
scc_cnt += 1;
scc[scc_cnt].clear();
for (;;){
int x = S.top(); S.pop();
if (sccno[x] != scc_cnt) scc[scc_cnt].push_back(x);
sccno[x] = scc_cnt;
if (x == u) break;
}
}
}
void tarjan(int n){
memset(sccno, 0, sizeof(sccno));
memset(pre, 0, sizeof(pre));
dfs_clock = scc_cnt = 0;
for (int i = 0; i < n; i++)
if (!pre[i]) dfs(i);
}
int cs, n, m;
int main(){
//freopen("in.txt", "r", stdin);
scanf("%d", &cs);
for (int k = 1; k <= cs; k++){
for (int i = 0; i <= n; i++)
G[i].clear();
scanf("%d%d", &n, &m);
for (int i = 0; i < m; i++){
int x, y;
scanf("%d%d", &x, &y);
x -= 1; y -= 1;
G[x].push_back(y);
}
tarjan(n);
printf("Case %d: ", k);
if (scc_cnt == 1){
printf("-1\n");
continue;
}
memset(idg, 0, sizeof(idg));
memset(odg, 0, sizeof(odg));
for (int u = 0; u < n; u++)
for (int i = 0; i < G[u].size(); i++){
int v = G[u][i];
if (sccno[u] != sccno[v]){
idg[sccno[v]] += 1;
odg[sccno[u]] += 1;
}
}
int minn = 0x7fffffff;
for (int i = 1; i <= scc_cnt; i++){
int tmp = scc[i].size();
if (idg[i] == 0 || odg[i] == 0) minn = min(minn, tmp);
}
int x = minn, y = n - minn;
LL ans = (LL)n * (LL)n - n - (LL)x * (LL)y - m;
cout << ans << endl;
}
return 0;
}