【图论·强连通分量scc】hdu6072Logical Chain

非常有灵性的一道强连通分量。乍一看是道模板题，然而需要用bitset优化时间复杂度至1/64；bitset什么鬼啊喂，很佩服赛场上唯一做出来的大佬orz；

常规部分，学习了scc模板求强连通分量；
第一遍dfs给图上的倒序编号（最后访问到的标最小号）；
第二遍在所有边反向后的图上dfs，从一个强连通分图上的点不可能访问到其他强分图上的点（为什么还没有理解）
cmp[I] = k 表示编号为i的点属于编号为k的强分图；
scc结束；

本题中节点数为t的强分图，对答案贡献t*（t-1）/2；

难点就是bitset 的理解和使用了（紧张）；
首先试了下头文件#include
然而长度为8的bitset<8> b; 只能放8个0/1；
并不是很神奇；遂放弃，看大神的手写版；

这就非常有灵性了。
reset(int x) 好理解，初始化。x=0 b为0，x = 0xff，无限大；
flip(int x) 不好理解。8位，256个数，把x这个数0变1，1变0；
b[x>>5] 每一位表示32个数，x>>5确定x属于哪一位管；
将这一位加上2的x%32次方，表示这个数的存在

#define _CRT_SECURE_NO_WARNINGS#include<iostream>#include<cstdio>#include<vector>#include<cstring>#include<stack>#include<queue>#include<bitset>using namespace std;using std::bitset;typedef long long ll;const int maxv = 255;    struct BIT{        unsigned int b[8];        void reset(int x){ memset(b, x, sizeof(b)); }        void flip(int x){ b[x >> 5] ^= 1U << (x & 31); }        int get(int x){ return (b[x >> 5] >> (x & 31)) & 1; }    }G[maxv], rG[maxv], vis;int V;vector<int> vs;int cmp[maxv];int cnt = 0;inline int ctz(unsigned int s){    int x = 0;    while (!(s & 1))    {        ++x; s >>= 1;    }    return x;}void dfs(int v){    vis.flip(v);    unsigned int s;    for (int i = 0; i < 8; ++i){        while (1){             s = vis.b[i] & G[v].b[i];            if (!s) break;            //这里是什么原理依然没有搞清楚orz；            dfs((i << 5)| ctz(s));          }    }    vs.push_back(v);}void rdfs(int v){    ++cnt;    vis.flip(v);    unsigned int s;    for (int i = 0; i < 8; i++){        while (1){            s = vis.b[i] & rG[v].b[i];            if (!s) break;            rdfs((i << 5 )| ctz(s));        }    }}void scc(){    vis.reset(0xff);vs.clear();    for (int i = 0; i < V; i++){        if (vis.get(i)) dfs(i);    }    vis.reset(0xff);    ll ans = 0;    for (int i = vs.size() - 1; i >= 0; i--){        if (vis.get(vs[i])) {            cnt = 0;            rdfs(vs[i]);            ans += cnt*(cnt - 1) / 2;        }    }    printf("%lld\n", ans);}char c;int t, m;int main(){    scanf("%d", &t);    while (t--){        scanf("%d%d", &V, &m);        for (int i = 0; i <= V; i++){            G[i].reset(0);            rG[i].reset(0);        }        for (int i = 0; i < V; i++){            char g[maxv];            scanf("%s", g);            for (int j = 0; j < V; j++){                if (g[j] == '1') {                    G[i].flip(j);                    rG[j].flip(i);                }            }        }        while (m--){            int k; scanf("%d", &k);            for (int i = 1; i <= k; ++i){                int u, v;        scanf("%d%d", &u, &v);                G[u - 1].flip(v - 1);                rG[v - 1].flip(u - 1);            }            scc();        }    }}