zoukankan      html  css  js  c++  java
  • HDU 4635 多校第四场 1004 强联通

    我还有什么好说,还有什么好说。。。。。。

    我是SBSBSBSBSBSBSBSBSBSBSBSBBSBSBSBSBSBSBSBSBS........................

    题意 思路什么的都不写了。。。。。。。。。。。我只是为了记录一下我是SBSBSBSBSBBSBSBSBSBSBSBSBSBSBB

    以后不要再SBSBSBSBSBSBSBSBSBBSBSBSBSBSBSBSBBSBSBSBBSBSBSBBSBSBSBSBSBBSBSBSBSBBSBSBSBSB

    /*******************/


    我还是写一下吧。

    首先缩点,缩点之后就是一些强联通分量,其实我们很容易想到,要使最后的图不是强联通图,那么我们可以枚举每一个强联通的分量,设这个分量是最后不与其他图不强联通的地方,那么除了这个分量,其他分量可以为完全图。 设总点数为n , 该分量点数为n' 。

    那么我们可以计算出除了该点之外的其他分量变成完全图所需要的边数为,( n - n') * (n - n' - 1) - 已有的边数。

    所谓已有的边数就是其他强联通里面已经存在的边,我们要减掉。具体看代码实现吧。

    然后还有一部分就是其实我们枚举出来的这个强联通分量,他也是可以成为一个完全图的,那么增加的边数为(n' * (n' -1)) - 已有的边数。

    最后还有一部分就是完全图到枚举的强联通分量之间的边,这部分的边就是n * n' - (n到n'之间的桥数) 。

    最后更新答案取最大值就可以了。

    我SB的把枚举那个条件判断成out[i] == 0,死都出不来,最后加上in[i] == 0就A掉了。

    我是SB吗?

    #include <iostream>
    #include <cstdio>
    #include <algorithm>
    #include <string>
    #include <cmath>
    #include <cstring>
    #include <queue>
    #include <set>
    #include <vector>
    #include <stack>
    #include <map>
    #include <iomanip>
    #define PI acos(-1.0)
    #define Max 2505
    #define inf 1<<28
    #define LL(x) ( x << 1 )
    #define RR(x) ( x << 1 | 1 )
    #define REP(i,s,t) for( int i = ( s ) ; i <= ( t ) ; ++ i )
    #define ll __int64
    #define mem(a,b) memset(a,b,sizeof(a))
    #define mp(a,b) make_pair(a,b)
    #define PII pair<int,int>
    using namespace std;
    
    
    inline void RD(int &ret) {
        char c;
        do {
            c = getchar();
        } while(c < '0' || c > '9') ;
        ret = c - '0';
        while((c=getchar()) >= '0' && c <= '9')
            ret = ret * 10 + ( c - '0' );
    }
    inline void OT(int a) {
        if(a >= 10)OT(a / 10) ;
        putchar(a % 10 + '0') ;
    }
    #define M 1000005
    #define N 100005
    int n , m ;
    struct kdq {
        int e , next ;
    } ed[M] ,ed1[M] ;
    struct B{
        int s, e ;
    }br[N] ;
    int head[N] , num ;
    int head1[N] ,num1 ;
    void add(int s , int e) {
        ed[num].e = e ;
        ed[num].next = head[s] ;
        head[s] = num ++ ;
    }
    void add1(int s ,int e){
        ed1[num1].e = e ;
        ed1[num1].next = head1[s] ;
        head1[s] = num1 ++ ;
    }
    bool vis[N] ;
    int BB[N] ;
    int dfn[N] ,low[N] , st[N] ,belong[N] ,cnt[N] ,out[N] , in[N] ,edg[N] , edn ;
    int ca , tp , dp ;
    void init() {
        num1 = edn = num = ca = tp = dp = 0 ;
        mem(dfn ,-1) ;
        mem(low ,0) ;
        mem(head, -1) ;
        mem(st, 0) ;
        mem(vis ,0) ;
        mem(belong ,0) ;
        mem(in , 0) ;
        mem(out, 0) ;
        mem(edg, 0) ;
        mem(head1, -1) ;
        mem(cnt ,0) ;
        mem(BB ,0) ;
    }
    void tarjan(int now) {
        dfn[now] = low[now] = dp ++ ;
        st[tp ++ ] = now ;
        vis[now] = 1 ;
        for (int i = head[now] ; i != -1 ; i = ed[i].next ) {
            int v = ed[i].e ;
    
            if(dfn[v] == -1) {
                tarjan(v) ;
                //edn ++ ;
                low[now] = min(low[now],low[v]) ;
            } else if(vis[v]) {
                low[now] = min(low[now] ,dfn[v]) ;
            }
        }
        if(low[now] == dfn[now]) {
            int tt ;
            ca ++ ;
            do {
                tt = st[-- tp] ;
                vis[tt] = 0 ;
                belong[tt] = ca ;
                cnt[ca] ++ ;
            } while(tt != now) ;
        }
    }
    struct PP{
        int num , id ;
    }p[N] ;
    
    bool cmp(const PP& a ,const PP& b){
        return a.num > b.num ;
    }
    int main() {
        int T ;
        cin >> T ;
        int cc = 0 ;
        while( T -- ) {
            cin >> n >> m ;
            init() ;
            for (int i = 0 ; i < m ; i ++ ) {
                int a,  b ;
                RD(a) ;
                RD(b) ;
                add(a , b) ;
            }
            for (int i = 1 ; i <= n ; i ++ ) {
                if(dfn[i] == -1)tarjan(i) ;
            }
            int bn = 0 ;
            for (int i = 1 ; i <= n ; i ++ ) {
                for (int j = head[i] ; ~j ; j = ed[j].next ) {
                    int y = belong[ed[j].e] ;
                    int x = belong[i] ;
                    if(x != y) {
                        br[bn].s = x ;
                        br[bn ++].e = y ;
                        out[x] ++ ;
                        in[y] ++ ;
                    }else {
                        edg[x] ++ ;
                    }
                }
            }
            printf("Case %d: ",++cc) ;
            ll ans = 0 ;
            if(ca == 1){
                puts("-1") ;
            }
            else {
    
                int numB = 0 ;
                for (int i = 0 ; i < bn ; i ++ ){
                    int s = br[i].s ;
                    int e = br[i].e ;
                    BB[e] ++ ;
                }
    
                for (int i = 1 ; i <= ca ; i ++ ){
                    if(out[i] == 0||in[i] == 0){//in[i] == 0 ,没有判断WA到死,我是SBSBSBSBSBSBSBSBSBSBSBSBSBSBSBSBSBSBSBBSBSB
                        int nowP = n - cnt[i] ;//除该分量外其他所有的点数
                        int nowE = m - edg[i] - BB[i];//除该分量以外所有的边数
                        ll sum = (ll)(nowP - 1) * (ll)nowP - nowE ;//连接除i外所有的分量
                        sum += (ll)nowP * (ll)cnt[i] - BB[i] ;//连接完分量之后和i连接
                        sum += (ll)cnt[i] * (ll)(cnt[i] - 1) - edg[i] ;//i分量内所有的点都连起来
                        ans = max(ans ,sum) ;
                    }
                }
                cout << ans << endl;
            }
        }
        return 0 ;
    }
    


  • 相关阅读:
    java中的Runtime类
    java web 三大核心组件Filter
    springBoot 整合mybaits 逆向工程
    java内存区域之程序计数器
    spring mvc 注解整理(一)
    ELK 日志查询分析nginx日志
    npm 淘宝镜像
    两种网页截图功能
    url 编码解码
    对防火墙的理解
  • 原文地址:https://www.cnblogs.com/riskyer/p/3231016.html
Copyright © 2011-2022 走看看