bzoj4671: 异或图

Posted ccz181078

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj4671: 异或图相关的知识,希望对你有一定的参考价值。

Description

定义两个结点数相同的图 G1 与图 G2 的异或为一个新的图 G, 其中如果 (u, v) 在 G1 与
G2 中的出现次数之和为 1, 那么边 (u, v) 在 G 中, 否则这条边不在 G 中.
现在给定 s 个结点数相同的图 G1...s, 设 S = {G1, G2, . . . , Gs}, 请问 S 有多少个子集的异
或为一个连通图?

Input

第一行为一个整数s, 表图的个数.
接下来每一个二进制串, 第 i 行的二进制串为 gi, 其中 gi 是原图通过以下伪代码转化得
到的. 图的结点从 1 开始编号, 下面设结点数为 n.
Algorithm 1 Print a graph G = (V, E)
for i = 1 to n do
for j = i + 1 to n do
if G contains edge (i, j) then
print 1
else
print 0
end if
end for
end for
 2 ≤ n ≤ 10,1 ≤ s ≤ 60.

Output

输出一行一个整数, 表示方案数
 
解异或方程组可以求出每种连通情况(用最小表示法表示连通性,要求不同标号的点严格不连通,相同标号的点任意)的方案数
容斥一下可以得到答案
#include<cstdio>
#include<cstring>
typedef unsigned long long u64;
char s[75];
int n=1,m,v[75],xp=0;
u64 ans=0,fac[75],x[75],e[75];
void dfs(int w,int mx=0){
    if(w==n){
        xp=0;
        for(int i=0,p=0;i<n;++i)
        for(int j=i+1;j<n;++j,++p)if(v[i]!=v[j])x[xp++]=e[p];
        int t=0;
        while(xp<m)x[xp++]=0;
        for(int i=0,w=0;i<m;++i){
            if(~x[w]>>i&1){
                for(int j=w+1;j<xp;++j)if(x[j]>>i&1){
                    x[w]^=x[j];
                    break;
                }
            }
            if(x[w]>>i&1){
                for(int j=w+1;j<xp;++j)if(x[j]>>i&1)x[j]^=x[w];
                ++w;
            }else ++t;
        }
        if(mx&1)ans+=fac[mx-1]<<t;
        else ans-=fac[mx-1]<<t;
        return;
    }
    for(int i=0;i<=mx;++i){
        v[w]=i;
        dfs(w+1,mx+(i==mx));
    }
}
int main(){
    scanf("%d",&m);
    for(int t=0,l;t<m;++t){
        scanf("%s",s);
        l=strlen(s);
        while(n*(n-1)/2<l)++n;
        for(int i=0,p=0;i<n;++i)
        for(int j=i+1;j<n;++j,++p)if(s[p]==1)e[p]|=1llu<<t;
    }
    for(int i=fac[0]=1;i<=n;++i)fac[i]=fac[i-1]*i;
    dfs(0);
    printf("%lld",ans);
    return 0;
}

 

以上是关于bzoj4671: 异或图的主要内容,如果未能解决你的问题,请参考以下文章

bzoj 4671: 异或图

bzoj4671: 异或图

BZOJ4671 异或图

bzoj 4671 异或图

bzoj4671: 异或图

bzoj 4671 异或图 —— 容斥+斯特林反演+线性基