zoukankan      html  css  js  c++  java
  • [HDOJ4635]Strongly connected(强连通分量,缩点)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4635

    题意:给一张图,问最多往这张图上加多少条边,使这张图仍然无法成为一个强连通图。

    起初是先分析样例可以知道,一个强连通分量里应当加边加至一个完全图,这样对整个图的连通贡献是没有的。然后把每个强连通分量里的边数扩展至num-out-1,但是发现连通分量之间的关系不好处理,所以反着做。

    考虑一张有向完全图的边数为n*(n-1),再去掉现在已经有的边数m,那么不考虑无法成为强连通图这一要求,我们还可以加n*(n-1)-m条边。现在问题变成了最少去掉多少条边,这个图不是一个强连通图。

    我们不希望浪费任何一个强连通分量,所以我们最终的情况应当是有一个强连通分量是无法到达其他强连通分量,并且其他强连通分量早已连接成一整个强连通分量。要保证这种边尽量多,也就是保证排除的点尽可能少,我们枚举其中的一个强连通分量(因为单个分量的点数一定小于多个连通分量点数相加后的最小值)。接下来在图上找入度或者出度为0的连通分量,从这种分量出或者入是不会影响连通性的,所以得解。

      1 /*
      2 ━━━━━┒ギリギリ♂ eye!
      3 ┓┏┓┏┓┃キリキリ♂ mind!
      4 ┛┗┛┗┛┃\○/
      5 ┓┏┓┏┓┃ /
      6 ┛┗┛┗┛┃ノ)
      7 ┓┏┓┏┓┃
      8 ┛┗┛┗┛┃
      9 ┓┏┓┏┓┃
     10 ┛┗┛┗┛┃
     11 ┓┏┓┏┓┃
     12 ┛┗┛┗┛┃
     13 ┓┏┓┏┓┃
     14 ┃┃┃┃┃┃
     15 ┻┻┻┻┻┻
     16 */
     17 #include <algorithm>
     18 #include <iostream>
     19 #include <iomanip>
     20 #include <cstring>
     21 #include <climits>
     22 #include <complex>
     23 #include <fstream>
     24 #include <cassert>
     25 #include <cstdio>
     26 #include <bitset>
     27 #include <vector>
     28 #include <deque>
     29 #include <queue>
     30 #include <stack>
     31 #include <ctime>
     32 #include <set>
     33 #include <map>
     34 #include <cmath>
     35 using namespace std;
     36 #define fr first
     37 #define sc second
     38 #define cl clear
     39 #define BUG puts("here!!!")
     40 #define W(a) while(a--)
     41 #define pb(a) push_back(a)
     42 #define Rint(a) scanf("%d", &a)
     43 #define Rll(a) scanf("%lld", &a)
     44 #define Rs(a) scanf("%s", a)
     45 #define Cin(a) cin >> a
     46 #define FRead() freopen("in", "r", stdin)
     47 #define FWrite() freopen("out", "w", stdout)
     48 #define Rep(i, len) for(int i = 0; i < (len); i++)
     49 #define For(i, a, len) for(int i = (a); i < (len); i++)
     50 #define Cls(a) memset((a), 0, sizeof(a))
     51 #define Clr(a, x) memset((a), (x), sizeof(a))
     52 #define Full(a) memset((a), 0x7f7f, sizeof(a))
     53 #define lp p << 1
     54 #define rp p << 1 | 1
     55 #define pi 3.14159265359
     56 #define RT return
     57 typedef long long LL;
     58 typedef long double LD;
     59 typedef unsigned long long ULL;
     60 typedef pair<int, int> pii;
     61 typedef pair<string, int> psi;
     62 typedef map<string, int> msi;
     63 typedef vector<int> vi;
     64 typedef vector<LL> vl;
     65 typedef vector<vl> vvl;
     66 typedef vector<bool> vb;
     67 
     68 typedef struct Edge {
     69     int u;
     70     int v;
     71     int next;
     72     Edge() { next = -1; }
     73 }Edge;
     74 const int maxn = 100100;
     75 const int maxm = 100100;
     76 int n, m;
     77 int head[maxn], ecnt;
     78 Edge edge[maxm];
     79 int bcnt, dindex;
     80 int dfn[maxn], low[maxn];
     81 int stk[maxn], top;
     82 int belong[maxn], num[maxn];
     83 int in[maxn], out[maxn];
     84 bool instk[maxn];
     85 
     86 void init() {
     87     Clr(head, -1);
     88     Cls(edge); Cls(instk);
     89     Cls(dfn); Cls(low); Cls(num);
     90     Cls(belong); Cls(in); Cls(out);
     91     ecnt = top = bcnt = dindex = 0;
     92 }
     93 
     94 void adde(int uu, int vv) {
     95     edge[ecnt].u = uu;
     96     edge[ecnt].v = vv;
     97     edge[ecnt].next = head[uu];
     98     head[uu] = ecnt++;
     99 }
    100 
    101 void tarjan(int u) {
    102     int v = u;
    103     dfn[u] = low[u] = ++dindex;
    104     stk[++top] = u;
    105     instk[u] = 1;
    106     for(int i = head[u]; ~i; i=edge[i].next) {
    107         v = edge[i].v;
    108         if(!dfn[v]) {
    109             tarjan(v);
    110             low[u] = min(low[u], low[v]);
    111         }
    112         else if(instk[v]) low[u] = min(low[u], dfn[v]);
    113     }
    114     if(dfn[u] == low[u]) {
    115         bcnt++;
    116         do {
    117             v = stk[top--];
    118             instk[v] = 0;
    119             belong[v] = bcnt;
    120         } while(v != u);
    121     }
    122 }
    123 
    124 int main() {
    125     // FRead();
    126     int T, _ = 1, u, v;
    127     Rint(T);
    128     W(T) {
    129         printf("Case %d: ", _++);
    130         init();
    131         Rint(n); Rint(m);
    132         Rep(i, m) {
    133             Rint(u); Rint(v);
    134             adde(u, v);
    135         }
    136         For(i, 1, n+1) if(!dfn[i]) tarjan(i);
    137         if(bcnt == 1) {
    138             printf("-1
    ");
    139             continue;
    140         }
    141         For(u, 1, n+1) {
    142             num[belong[u]]++;
    143             for(int i = head[u]; ~i; i=edge[i].next) {
    144                 int v = edge[i].v;
    145                 if(belong[u] != belong[v]) {
    146                     in[belong[v]]++; out[belong[u]]++;
    147                 }
    148             }
    149         }
    150         int sum = n*(n-1)-m;
    151         int ret = 0;
    152         For(i, 1, bcnt+1) {
    153             if(in[i] == 0 || out[i] == 0) {
    154                 ret = max(ret, sum-num[i]*(n-num[i]));
    155             }
    156         }
    157         printf("%d
    ", ret);
    158     }
    159     RT 0;
    160 }
  • 相关阅读:
    解决vmware Invalid memory setting (sched.mem.min)
    PostgreSQL教程
    rpm包安装过程中依赖问题“libc.so.6 is needed by XXX”解决方法
    使用厂商MIB库查找设备OID值 并实施监控的方法
    【交换机】我司交换机上常用的一些MIB以及对应的OID说明
    CentOS6.8-minimal安装gnome桌面 安装NVC远程桌面连接
    LINUX新建和增加SWAP分区
    Caused by: java.lang.ClassNotFoundException: org.aspectj.weaver.reflect.ReflectionWorld$ReflectionWo
    类与对象
    Volley框架源代码分析
  • 原文地址:https://www.cnblogs.com/kirai/p/5523407.html
Copyright © 2011-2022 走看看