题目

【BZOJ4671】异或图

很有意思的题

做法

直接处理显然很难,我们考虑范围扩大以求容斥或反演这类的帮助

\(f_i\)表示至少有\(i\)个联通块的方案,形如设立\(i\)个联通块轮廓,联通块内连边随意,联通块与联通块之间无连边

\(g_i\)表示恰好有\(i\)个联通块的方案,形如设立\(i\)个联通块轮廓,在保证内部联通的情况下,外部块与块间无连边

显然:$$f_x=\sum\limits_{i=x}^n\begin{Bmatrix}i\x\end{Bmatrix}g_i$$

根据斯特林反演:$$g_x=\sum\limits_{i=x}^n (-1)^{i-x}\begin{bmatrix}i\x\end{bmatrix}f_i$$

故\(g_1=\sum\limits_{i=1}^n (-1)^{i-1}\begin{bmatrix}i\\1\end{bmatrix}f_i\)

而\(\begin{bmatrix}i\\1\end{bmatrix}\)是阶乘形式:\(\begin{bmatrix}i\\1\end{bmatrix}=(i-1)!\)

化简答案为:\(g_1=\sum\limits_{i=1}^n (-1)^{i-1}(i-1)!f_i\)

考虑\(f_i\)如何求出:状压点所属联通块状态,则我们要选择图集使块与块之间无边,考虑枚举每个图的\(S\)表示点与点之间的连边(不属同一联通块),我们压到线性基里去,\(ele\)表示线性基元素,这些元素是不能选择的(相异),故答案为\(2^{N-ele}\)

Code

#include<bits/stdc++.h>
typedef int LL;
const LL maxn=109;
LL N,n;
LL G[maxn][maxn][maxn],a[maxn];
char s[maxn];
long long ans,p[maxn],S,fac[15];
void Dfs(LL x,LL up){
if(x==n+1){
memset(p,0,sizeof(p)); LL ele(0);
for(LL i=1;i<=N;++i){
S=0; LL tot(0);
for(LL j=1;j<=n;++j)
for(LL k=j+1;k<=n;++k)
if(a[k]!=a[j]){
S|=(1ll<<tot)*G[i][j][k];
++tot;
}
for(LL j=0;j<tot;++j){
if(S&(1ll<<j)){
if(!p[j]){
p[j]=S;
++ele;
break;
}else
S^=p[j];
}
}
}
ans+=1ll*((up&1)?1:-1)*fac[up-1]*(1ll<<N-ele);
return;
}
for(LL i=1;i<=up+1;++i){
a[x]=i;
Dfs(x+1,std::max(up,i));
}
}
int main(){
scanf("%d",&N);
for(LL i=1;i<=N;++i){
scanf(" %s",s+1);
LL len(strlen(s+1));
if(!n){
n=1;
for(;n*(n-1)/2!=len;++n);
}
LL now(0);
for(LL j=1;j<=n;++j) for(LL k=j+1;k<=n;++k) G[i][j][k]=s[++now]-'0';
}
fac[0]=fac[1]=1; for(LL i=2;i<=n;++i) fac[i]=fac[i-1]*i;
Dfs(1,0);
printf("%lld",ans);
return 0;
}

【BZOJ4671】(斯特林反演)的更多相关文章

  1. bzoj4671: 异或图——斯特林反演

    [BZOJ4671]异或图 - xjr01 - 博客园 考虑先算一些限制少的情况 gi表示把n个点的图,划分成i个连通块的方案数 连通块之间不连通很好处理(怎么处理看下边),但是内部必须连通,就很难办 ...

  2. 【BZOJ4671】异或图(斯特林反演)

    [BZOJ4671]异或图(斯特林反演) 题面 BZOJ Description 定义两个结点数相同的图 G1 与图 G2 的异或为一个新的图 G, 其中如果 (u, v) 在 G1 与 G2 中的出 ...

  3. bzoj4671 异或图(斯特林反演,线性基)

    bzoj4671 异或图(斯特林反演,线性基) 祭奠天国的bzoj. 题解时间 首先考虑类似于容斥的东西. 设 $ f_{ i } $ 为至少有 $ i $ 个连通块的方案数, $ g_{ i } $ ...

  4. 【bzoj4671】异或图(容斥+斯特林反演+线性基)

    传送门 题意: 给出\(s,s\leq 60\)张图,每张图都有\(n,n\leq 10\)个点. 现在问有多少个图的子集,满足这些图的边"异或"起来后,这张图为连通图. 思路: ...

  5. Square(斯特林反演)

    题意 给出一个 \(n × m\) 大小的矩形,每个位置可以填上 \([1, c]\) 中的任意一个数,要求填好后任意两行互不等价且任意两列互不等价,两行或两列等价当且仅当对应位置完全相同,求方案数 ...

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

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4671 考虑计算不是连通图的方案,乘上容斥系数来进行容斥. 可以枚举子集划分(复杂度是O(Be ...

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

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4671 首先,考虑容斥,就是设 \( t[i] \) 表示至少有 \( i \) 个连通块的方 ...

  8. TCO14 Wildcard CountTables——斯特林反演

    不知道咕了多长时间的题... 讲了3遍,还是自己搞懂了.. 暂时没有找到题目链接 题意: n×m的网格,每个格子填[1,x]的数,使得不存在两行两列同构. 先保证一个,行相同. 再容斥掉列. 枚举至多 ...

  9. 斯特林数&斯特林反演

    第一类斯特林数 定义 第一类Stirling数\(s(n,m)\),也可记为\(\begin{bmatrix}n\\m\end{bmatrix}\). 第一类Stirling分为无符号第一类Stirl ...

随机推荐

  1. LeetCode 第 231 题 (Power of Two)

    LeetCode 第 231 题 (Power of Two) Given an integer, write a function to determine if it is a power of ...

  2. iOS conformsToProtocol

    - (BOOL)conformsToProtocol:(Protocol *)aProtocol; 是用来检查对象(包括其祖先)是否实现了指定协议类的方法. 今天遇到一个问题,一个类没有实现proto ...

  3. C语言程序设计-同一天生日[综合应用]

    [问题描述] 在一个有200人的大班级中,存在两个人生日相同的概率非常大,现给出每个学生的学号,出生月日,试找出所有生日相同的学生. [输入形式] 第一行为整数n,表示有n个学生,n<=200. ...

  4. [分享]JavaScript Quick Reference Card

    pdf文件:https://files.cnblogs.com/files/MakeView660/JavaScript_Quick_Reference_Card.pdf

  5. Python_selenium之窗口切换

    Python_selenium之窗口切换 1. 运用switch_to.window()方法来进行窗口切换 2. 思路拆分: 浏览器获取百度贴吧网址 点击定位到一个元素,获取当前的句柄 获得所有的句柄 ...

  6. 基于注解的Spring AOP拦截含有泛型的DAO

    出错场景 1.抽象类BaseDao public abstract class BaseDao<T> { public BaseDao() { entityClass = (Class&l ...

  7. angualar入门学习-- 自定义指令 指令编译执行过程

    3个阶段: 一.加载阶段 加载angular.js的源码,找到ng-app确定应用边界范围. 二.编译阶段 compile 查找所有指令,保存在一个列表中 对所有指令按优先级(property属性值) ...

  8. 160809、tomcat中配置多个域名及将tomcat配置成系统服务

    本地测试用的(注意红色部分) 第一步.自己的windows电脑,在c盘中有个hosts文件(搜索一下),做以下修改(其中127.0.0.1是本机地址,192.1638.10.139是我虚拟机中linu ...

  9. JDK源码分析之concurrent包(二) -- 线程池ThreadPoolExecutor

    上一篇我们简单描述了Executor框架的结构,本篇正式开始并发包中部分源码的解读. 我们知道,目前主流的商用虚拟机在线程的实现上可能会有所差别.但不管如何实现,在开启和关闭线程时一定会耗费很多CPU ...

  10. delphi--自带控件(Tcheckbox、TComBobox、TpageControl、TDBGrid,TScrollBox,TEdit)

    一.TCheckBox 使CheckBox选中CheckBox1.Checked := True;选择发生变化 ,OnClick事件if CheckBox1.Checked thenbegin //. ...