HDU 6072 Logical Chain（Kosaraju+bitset）

http://acm.hdu.edu.cn/showproblem.php?pid=6072

题意：

给你$n*n$的矩阵，每次修改k条边，让你计算其中能相互到达的点对有多少。

思路：

其实就是求强连通分量，如果一个强连通分量里有n个点，那么这里面的点对就有$n*(n-1)/2$。用Kosaraju计算即可，但是这样是会超时的，还需要用bitset来优化。

__builtin_函数中处理二进制位的函数：

int __builtin_ffs (unsigned x) 返回x中最后一个1是从右往左第几位

int __builtin_popcount (unsigned x) 返回x中1的个数

int __builtin_ctz (unsigned x) 返回x末尾0的个数(x等于0时未定义)

int __builtin_clz (unsigned x) 返回x中前导0的个数(x等于0时未定义)

int __builtin_parity (unsigned x) 返回x中1的奇偶性

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<stack>
#include<queue>
#include<cmath>
#include<map>
#include<set>
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
const int maxn=250+5;

int n, m;
int num;
char s[maxn];

struct Bitset//用unsigned数组实现bitset，为了使用__builtin_ctz()
{
    unsigned v[8];//8个unsigned表示8*32=256位

    void reset()//清零
    {
        for(int i=0;i<8;++i) v[i]=0;
    }
    void set(int x)//把某一位设为1
    {
        v[x>>5]|=1u<<(x&31);
    }
    void flip(int x)//翻转某一位
    {
        v[x>>5]^=1u<<(x&31);
    }
    bool test(int x)//返回某一位是否为1
    {
        return v[x>>5]>>(x&31)&1;
    }
}vis, G[maxn], rG[maxn];

vector<int> S;

void dfs1(int u)
{
    vis.flip(u);
    for(int i=0;i<8;i++)
    {
        while(true)
        {
            unsigned tmp=vis.v[i]&G[u].v[i];  //tmp就计算出了可以访问的点
            if(!tmp)  break;
            dfs1(i<<5|__builtin_ctz(tmp));  //计算出tmp末尾有多少0，有多少0就代表了它是第几位，这儿的话一个点一个点的来
        }
    }
    S.push_back(u);
}

void dfs2(int u)
{
    vis.flip(u);
    ++num;
    for(int i=0;i<8;i++)
    {
        while(true)
        {
            unsigned tmp=vis.v[i]&rG[u].v[i];
            if(!tmp)  break;
            dfs2(i<<5|__builtin_ctz(tmp));
        }
    }
}

void Kosaraju()
{
    S.clear();
    int ans=0;
    for(int i=0;i<n;i++)  vis.set(i);
    for(int i=0;i<n;i++)  if(vis.test(i))  dfs1(i);
    for(int i=0;i<n;i++)  vis.set(i);
    for(int i=n-1;i>=0;i--)
    {
        if(vis.test(S[i]))
        {
            num=0;
            dfs2(S[i]);
            ans+=num*(num-1)/2;
        }
    }
    printf("%d
",ans);
}


int main()
{
    //freopen("in.txt","r",stdin);
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++) G[i].reset(),rG[i].reset();
        for(int i=0;i<n;i++)
        {
            scanf("%s",s);
            for(int j=0;j<n;j++)  if(s[j]=='1')
            {
                G[i].flip(j);
                rG[j].flip(i);
            }
        }
        while(m--)
        {
            int t; scanf("%d",&t);
            while(t--)
            {
                int u,v;
                scanf("%d%d",&u,&v);
                G[u-1].flip(v-1);
                rG[v-1].flip(u-1);
            }
            Kosaraju();
        }
    }
    return 0;
}