HDU-4635 Strongly connected

在图是简单有向图和图不是强连通图的条件下，问最多可以添加多少边

简单有向图就是没有重边，并且没有自环的有向图，完全有向图就是最大的简单有向图

先按照简单有向图的条件来，已经有了m条边，我们最多能添加ans=n*(n-1)-m条边

然后再按不是强连通图的条件，图中至少应该有一个强连通分量的出度或入度为0

先对给出的图缩点，然后统计每个强连通分量的出度入度，找出出度或入度为0且包含点最少的那个

设分量点数量为minnum，最后可以添加的边数为ans-minnum*(n-minnum)

就是说我们保证选中的这个分量的出度或入度为0，不给它加边，其他的都加上就可以了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
const int MAXN = 100010;//点数
const int MAXM = 100100;//边数
struct Edge
{
	int to,next;
}edge[MAXM];
int head[MAXN],tot;
int Low[MAXN],DFN[MAXN],Stack[MAXN],Belong[MAXN],in[MAXN],out[MAXN],num[MAXN];//Belong数组的值是1~scc
int Index,top;
int scc;//强连通分量的个数
bool Instack[MAXN];
void addedge(int u,int v)
{
	edge[tot].to = v;edge[tot].next = head[u];head[u] = tot++;
}
void Tarjan(int u,int f)
{
	int v;
	Low[u] = DFN[u] = ++Index;
	Stack[top++] = u;
	Instack[u] = true;
	for(int i = head[u];i != -1;i = edge[i].next)
	{
		v = edge[i].to;
		if( !DFN[v] )
		{
			Tarjan(v,u);
			if( Low[u] > Low[v] )Low[u] = Low[v];
		}
		else if(Instack[v] && Low[u] > DFN[v])
		Low[u] = DFN[v];
	}
	if(Low[u] == DFN[u])
	{
		scc++;
		do
		{
			v = Stack[--top];
			Instack[v] = false;
			Belong[v] = scc;
			num[scc]++;
		}while( v != u);
	}
}
void solve(int N)
{
	memset(DFN,0,sizeof(DFN));
	memset(Instack,false,sizeof(Instack));
	memset(num,0,sizeof(num));
	Index = scc = top = 0;
	for(int i = 1;i <= N;i++)
	if(!DFN[i]) Tarjan(i,i);
}
void init()
{
	tot = 0;
	memset(head,-1,sizeof(head));
}
int main()
{
    //freopen("/home/zlwang/test.txt","r",stdin);
    int T,n,m,a,b,kase=0;
    scanf("%d",&T);
    while(T--)
    {
        init();
        scanf("%d%d",&n,&m);
        for(int i=0;i<m;i++)
        {
            scanf("%d%d",&a,&b);
            addedge(a,b);
        }
        solve(n);
        printf("Case %d: ",++kase);
        if(scc==1)
        {
            printf("-1\n");
            continue;
        }
        memset(in,0,sizeof(in));
        memset(out,0,sizeof(out));
        int ans=(n-1)*n-m;
        int minnum=n;
        for(int u=1;u<=n;u++)
        for(int j=head[u];j!=-1;j=edge[j].next)
        {
            int v=edge[j].to;
            if(Belong[u]!=Belong[v])
            {
                in[Belong[v]]++;
                out[Belong[u]]++;
            }
        }
        for(int i=1;i<=scc;i++)
        if(in[i]==0||out[i]==0)
        minnum=min(minnum,num[i]);
        printf("%d\n",ans-minnum*(n-minnum));
    }
    return 0;
}