HDU4635 Strongly connected

点此跳转到题目

又是一个强连通分量，首先我们要明确，加了边之后的图是一个什么图，一定存在两个部分，假设为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;
}