Codeforces908H. New Year and Boolean Bridges【并查集+强联通+FWT】

题目大意：

有一个n个点的有向图。
定义i能到达j时f(i,j)=1，否则f(i,j)=0。
对于每个点对(i,j)，给定以下三个条件中的某一个为真：
(1) f(i,j) and f(j,i)=1；
(2) f(i,j) or f(j,i)=1；
(3) f(i,j) xor f(j,i)=1。
求满足条件时的最小边数。
1<=n<=47。

解题思路：

等价于满足以下两个条件：
(1)满足 $and$ 的点对一定在同一个强连通分量里。
(2)一个强连通分量里不能有 $xor$，否则无解。

将满足 $and$ 的点对用并查集连起来当做一个块。
一些块能合并成一个块当且仅当它们之间没有 $xor$。
大小为 $1$ 的块要跟外面连通显然贡献是 $1$，对于 $size\ge 2$ 的块，最优肯定是用 $size$ 条边连一个环，而所有块还要连成一棵树，所以最后就是最小化$\ge 2$的块数 $cnt$，答案就是 $n+cnt-1$。

注意$\ge 2$的下等于$n/2$，可以用状压了，一下$m$ 表示初始时 $\ge 2$ 的块数。

先dp预处理出每个集合能否合并，将不合法点对的值设为 $1$，那么求子集和，为0即可合并，否则不行。暴力求子集和是$O(3^m)$的，但其实子集并卷积FWT时就是求子集和，所以是$O(m2^m)$的。

暴力就是枚举子集。
令 $f(i)$ 表示 $i$ 代表的块合并后的最小块数。
$f(i)=min\{f(j)+f(i \oplus j)\} (j\subset i)$
时间复杂度$O(3^m)$

令$f(i)$有值当且仅当$i$代表的块能合并成一块。
好比邻接矩阵的 $k$ 次方对应两点间是否有长度为 $k$ 的路径，我们需要找到最小的 $k$ 使得 $k$个 $f$ 进行集合并卷积后 $f(2^m-1)$ 有值，那么 $cnt=k$。
(开始认为应该是子集卷积，后来发现集合并卷积也没错)

一个直观想法是每次做FWT-UFWT后检查，但最坏会做 $m-1$次，不太能过。
其实可以模拟UFWT的过程处理出每个位置的点值对$f(2^m-1)$的系数，先将 $f$ FWT，每次点乘后不用再UFWT回去，而是直接$O(2^m)$算出$f(2^m-1)$即可。

复杂度为$O(n2^{n/2})$

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=50;
const int M=8388610;
int n,m,fmp[N],s[N],idx[N],f[M],g[M],a[M];
char mp[N][N];
int F(int x){return fmp[x]==x?x:fmp[x]=F(fmp[x]);}
void FWT(int a[],int len)
{
    for(int i=1;i<len;i<<=1)
        for(int j=0;j<len;j+=(i<<1))
            for(int k=j;k<j+i;k++)
                f[k+i]+=f[k];
}
void init(int len)
{
    for(int i=0;i<len;i++)a[i]=1;
    for(int i=1;i<len;i<<=1)
        for(int j=0;j<len;j+=(i<<1))
            for(int k=j;k<j+i;k++)
                a[k]*=-1;       
}
int main()
{
    //freopen("lx.in","r",stdin);
    scanf("%d",&n);
    for(int i=0;i<n;i++)scanf("%s",mp[i]),fmp[i]=i;
    for(int i=0;i<n;i++)
        for(int j=0;j<i;j++)
            if(mp[i][j]=='A')fmp[F(i)]=F(j);
    for(int i=0;i<n;i++)++s[F(i)];
    for(int i=0;i<n;i++)if(F(i)==i&&s[i]>1)idx[i]=1<<m++;
    if(!m){cout<<n-1<<'\n';return 0;}
    for(int i=0;i<n;i++)
        for(int j=0;j<i;j++)
            if(mp[i][j]=='X')
            {
                if(F(i)==F(j)){puts("-1");return 0;}
                else if(idx[F(i)]&&idx[F(j)])f[idx[F(i)]|idx[F(j)]]=1;
            }
    int len=1<<m;FWT(f,len);
    for(int i=0;i<len;i++)f[i]=!f[i],g[i]=1;
    FWT(f,len),init(len);
    for(int i=0;i<m;i++)
    {
        int cur=0;
        for(int j=0;j<len;j++)g[j]*=f[j],cur+=a[j]*g[j];
        if(cur){cout<<n+i<<'\n';return 0;}
    }
    return 0;
}