Strongly connected HDU - 4635

题意：

给定一个有向图，求最大可以增加多少条边使得这个仍然不是强连通。

题解：

考虑逆向思维
我们知道一个有向完全图的边数为 n * (n - 1)
比如现在有k(k >= 2)个连通块，其中一个连通块的点数为 x
那么我们将这 x 个点和其他 n - x 个点的边删掉，删掉后的图一定不是强连通的

因此我们只需要考虑让需要删掉的边 x * (n - x) 最少，那最后留下的边就越多
所以我们只需要使 n * (n - 1) - x * (n - x) - m 最大即可

因为n, m都是定值，这其实也就是一个定义域在[1, n]，关于 n/2 对称的二次函数
由于是函数是下凹的，所以我们的 x 取靠近 1 或者 靠近 n 的

其实我们直接取靠近 1 的即可，因为其他连通块的点数一定 <= n - x，和 n 的绝对值一定大于等于 x 和 1 的绝对值

因此直接找到入度或出度为 0 的连通块中点的数量最小的即可

代码：

/*
 * @Author : Nightmare
 */
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define ls 2 * rt
#define rs 2 * rt + 1
#define gcd(a,b) __gcd(a,b)
#define eps 1e-8
#define lowbit(x) (x&(-x))
#define N 200005
#define M 2000005
#define mod 1000000007
#define inf 0x3f3f3f3f
int n, m, cas, cnt, low[N], dfn[N], instack[N], st[N], top;
int a[N], b[N], sc, scc[N], in[N], out[N], cont[N];
vector<int> edge[N];
void add(int a, int b){ edge[a].push_back(b); }
void init(){
    for(int i = 1 ; i <= n ; i ++) edge[i].clear();
    for(int i = 1 ; i <= n ; i ++) low[i] = dfn[i] = instack[i] = in[i] = out[i] = cont[i] = 0;
    sc = cnt = top = 0;
}
void tarjan(int u){
    low[u] = dfn[u] = ++cnt; st[++top] = u; instack[u] = 1;
    for(auto &v : edge[u]){
        if(!dfn[v]) tarjan(v), low[u] = min(low[u], low[v]);
        else if(instack[v]) low[u] = min(low[u], dfn[v]);
    }
    if(low[u] == dfn[u]){
        ++sc; int v;
        do{
            v = st[top --]; instack[v] = 0; scc[v] = sc; cont[sc] ++;
        }while(u != v);
    }
}
void solve(){
    printf("Case %d: ", ++cas);
    scanf("%d %d", &n, &m);
    init();
    for(int i = 1 ; i <= m ; i ++) scanf("%d %d", &a[i], &b[i]), add(a[i], b[i]);
    for(int i = 1 ; i <= n ; i ++) if(!dfn[i]) tarjan(i);
    if(sc == 1){ puts("-1"); return ; }
    for(int i = 1 ; i <= m ; i ++)
        if(scc[a[i]] != scc[b[i]])
            out[scc[a[i]]] ++, in[scc[b[i]]] ++;
    int mi = n;
    // 入度为 0 表示 删去的是 其他所有点到该连通块中的所有点的有向边
    // 出度为 0 表示 删去的是 该连通块中的所有点到其他所有点的有向边
    for(int i = 1 ; i <= sc ; i ++)
        if(!in[i] || !out[i])
            mi = min(mi, cont[i]);
    ll ans = 1ll * n * (n - 1) - m - 1ll * mi * (n - mi);
    printf("%lld\n", ans);
}
signed main(){
#ifndef ONLINE_JUDGE
    freopen("F:\\in.txt", "r", stdin);
#endif
    int T; scanf("%d", &T); while(T--) solve();
#ifndef ONLINE_JUDGE
    cerr << "Time elapsed: " << 1.0 * clock() / CLOCKS_PER_SEC << " s.\n";
#endif
    return 0;
}