强连通之HDU4635 Strongly connected

最新推荐文章于 2020-10-21 20:06:56 发布

htw0056

最新推荐文章于 2020-10-21 20:06:56 发布

阅读量243

点赞数

分类专栏：图论

本文链接：https://blog.csdn.net/u013656184/article/details/48878051

版权

图论专栏收录该内容

6 篇文章 0 订阅

订阅专栏

还是参考了别人的思路，有待提高

1.假设该图为完全图，那么一共有n*(n-1)条边

2.已经存在m条边，最多能加n*(n-1)-m条边

3.缩点以后，如果某点的入度或者出度为0，此时，应该保留该点所表示的点数最小的数量k，语文不是很好~其余的点到这个缩点后的点要么全都不指向这个点，要么全都不从这个点指出，即该点永远无法成为强连通图之内的点。最终的结果就是n*(n-1)-m-k*(n-k)

#include <stdio.h>
#include <ctype.h>
#include <string.h>
#include <stdlib.h>
#include <limits.h>
#include <math.h>
#include <algorithm>
#include <stack>
#include <queue>
#include <vector>
#include <map>
#include <set>
#include <string>
#include <sstream>
using namespace std;

#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1

typedef long long LL;
const double pi=4.0*atan(1.0);
const int INF=1<<30;
const int MAXN=100005;
vector<int> g[MAXN];
int dfn[MAXN],lowlink[MAXN],sccno[MAXN],dfs_clock,scc_cnt;
int sum[MAXN],in[MAXN],out[MAXN];
stack<int> s;

void dfs(int u)
{
    dfn[u]=lowlink[u]=++dfs_clock;
    s.push(u);
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        if(!dfn[v])
        {
            dfs(v);
            lowlink[u]=min(lowlink[u],lowlink[v]);
        }
        else if(!sccno[v])
        {
            lowlink[u]=min(lowlink[u],dfn[v]);
        }
    }

    if(lowlink[u]==dfn[u])
    {
        scc_cnt++;
        while(1)
        {
            int x=s.top();
            s.pop();
            sum[scc_cnt]++;
            sccno[x]=scc_cnt;
            if(x==u)
                break;
        }
    }
}

void Tarjan(int n)
{
    dfs_clock=scc_cnt=0;
    memset(sccno,0,sizeof(sccno));
    memset(dfn,0,sizeof(dfn));
    memset(sum,0,sizeof(sum));
    while(!s.empty())
        s.pop();

    for(int i=1;i<=n;i++)//下标1到n
        if(!dfn[i])
            dfs(i);
}

int main()
{
    int i,j,k;
    int T;
    int n,m;
    while(scanf("%d",&T)!=EOF)
    {
        for(int cas=1;cas<=T;cas++)
        {
            scanf("%d%d",&n,&m);
            for(i=1;i<=n;i++)
                g[i].clear();
            for(i=0;i<m;i++)
            {
                scanf("%d%d",&j,&k);
                g[j].push_back(k);
            }
            Tarjan(n);
            if(scc_cnt==1)
            {
                printf("Case %d: -1\n",cas);
                continue;
            }
            memset(in,0,sizeof(in));
            memset(out,0,sizeof(out));
            for(i=1;i<=n;i++)
                for(j=0;j<g[i].size();j++)
                {
                    if(sccno[i]!=sccno[g[i][j]])
                    {
                        out[sccno[i]]++;
                        in[sccno[g[i][j]]]++;
                    }
                }
            int Min=INF;
            for(i=1;i<=scc_cnt;i++)
                if(!in[i]||!out[i])
                    Min=min(Min,sum[i]);
            printf("Case %d: %I64d\n",cas,LL(n)*(n-1)-m-(LL)Min*(n-Min));
        }
        
    }
    return 0;
}