【BZOJ 4671】异或图 【斯特林反演】【线性基】【贝尔数复杂度】

传送门

题意:定义两个图的异或的边集为在两张图中恰出现一次的边。给sssnnn个点的图的集合,求异或和为连通图的子集数。

s≤60,n≤10s \leq 60,n \leq 10s60,n10

GiG_iGi表示异或出iii个连通块的子集数,答案就是G1G_1G1

GGG并不好求废话,我们考虑一个看起来很对的求法:

虽然边很多但点很少, 我们可以暴力枚举连通块中的点,复杂度是O(B(n))O(B(n))O(B(n)),BBB为贝尔数,B(10)=21147B(10)=21147B(10)=21147

可以确定的是不同连通块间没有边,但连通块间仍然不好确定,所以不能叫连通块。为了区分,我们把枚举出来的称为「集合」。

FiF_iFi为异或出iii个集合的子集数,即只要求不同集合间没有边,同一集合间没有要求没有要求。

发现GGG可以表示出FFF:因为FFF的每个「集合」一定是若干连通块构成的非空集合,即一个斯特林数

Fm=∑i=mn{im}GiF_m=\sum_{i=m}^n\left\{\begin{matrix}i\\m\end{matrix}\right\}G_iFm=i=mn{im}Gi

斯特林反演一波

Gm=∑i=mn(−1)i−m[im]FiG_m=\sum_{i=m}^n(-1)^{i-m}\left[\begin{matrix}i\\m\end{matrix}\right]F_iGm=i=mn(1)im[im]Fi

答案为

G1=∑i=1n(−1)i−1[i1]Fi=∑i=1n(−1)i−1(i−1)!FiG_1=\sum_{i=1}^n(-1)^{i-1}\left[\begin{matrix}i\\1\end{matrix}\right]F_i \\=\sum_{i=1}^n(-1)^{i-1}(i-1)!F_iG1=i=1n(1)i1[i1]Fi=i=1n(1)i1(i1)!Fi

现在我们只需要求出FFF

枚举出集合后,只考虑所有图中跨集合的边(实现时取个&即可),然后需要求出多少个子集异或和为000

把所有图插进线性基,设最终大小为sizsizsiz,那么自由元个数为s−sizs-sizssiz。首先只含线性基中元素的非空子集不会有贡献,然后对于任意自由元的子集,根据线性基的定义,都可以找到一个线性基的子集把它异或成000

GGG的过程中可能会爆ll,但反正最终答案在范围内,溢出就溢出吧。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
using namespace std;
typedef long long ll;
int s,n;
char t[100];
ll h[100],bas[100],f[100],ans;
int col[100];
inline void insert(ll x)
{for (int i=60;i>=0;i--)if (x&(1ll<<i))if (bas[i]) x^=bas[i];else return (void)(bas[i]=x);
}
void dfs(int k,int c)
{if (k>n){ll d=0;for (int i=1,cnt=0;i<n;i++)for (int j=i+1;j<=n;j++,cnt++)if (col[i]!=col[j])d|=1ll<<cnt;memset(bas,0,sizeof(bas));for (int i=1;i<=s;i++) insert(h[i]&d);int cnt=0;for (int i=0;i<=60;i++) cnt+=(bas[i]>0);f[c]+=1ll<<s-cnt;return;}for (col[k]=1;col[k]<=c;col[k]++,dfs(k+1,c));col[k]++,dfs(k+1,c+1);
}
int main()
{scanf("%d",&s);for (int i=1;i<=s;i++){scanf("%s",t);for (int j=strlen(t)-1;j>=0;j--) h[i]=(h[i]<<1)^(t[j]^48);}int len=strlen(t);for (n=1;n*(n-1)/2<len;n++);dfs(1,0);ll fac=1;for (int i=1;i<=n;fac*=i,i++) ans+=((i&1)? 1:-1)*fac*f[i];printf("%lld\n",ans);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/315086.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

.NET Core IdentityServer4实战 第Ⅳ章-集成密码登陆模式

回顾下ClientCredentials模式&#xff0c;在ReSourceApi中定义了我们公开服务&#xff0c;第三方网站想要去访问ReSourceApi则需要在身份验证服务中获取toekn&#xff0c;根据token的内容&#xff0c;硬编码去访问公开服务&#xff08;ResApi&#xff09;,这个还是非常简单的&a…

Codeforces Round #656 (Div. 3) F. Removing Leaves 贪心 + 模拟

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先有一个贪心策略就是每次都找一个叶子节点最多的点&#xff0c;让后删掉他的kkk个叶子节点&#xff0c;现在我们就来考虑如何模拟这个过程。 我们整一个vector<set<int>>ve…

.NetCore中三种注入生命周期的思考

.NetCore彻底诠释了“万物皆可注入”这句话的含义&#xff0c;在.NetCore中到处可见注入的使用。因此core中也提供了三种注入方式的生命周期使用&#xff0c;分别是&#xff1a;AddTransient&#xff1a;每次请求&#xff0c;都获取一个新的实例。即使同一个请求获取多次也会是…

兰州大学第一届 飞马杯 体育课排队 二分 + 最大流 + 输出路径

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 读懂题就会辣&#xff0c;经典模型了&#xff0c;二分时间&#xff0c;让后将其转换成二分图&#xff0c;左边是人&#xff0c;右边是位置&#xff0c;能在规定时间到的连边&#xff0c;跑…

从零开始实现ASP.NET Core MVC的插件式开发(一) - 使用Application Part动态加载控制器和视图...

如果你使用过一些开源CMS的话&#xff0c;肯定会用过其中的的插件化功能&#xff0c;用户可以通过启用或者上传插件包的方式动态添加一些功能&#xff0c;那么在ASP.NET Core MVC中如何实现插件化开发呢&#xff0c;下面我们来探究一下。本系列只是笔者的一些尝试&#xff0c;并…

感性理解Berlekamp-Massey算法

引入 BM算法主要解决的是根据数列求最短线性齐次递推式的问题在OI中主要辅助打表使用 即&#xff1a;已知FFF&#xff0c;求序列AAA使得 Fn∑i1mAiFn−i(n>m)F_n\sum_{i1}^mA_iF_{n-i} \quad(n>m)Fn​i1∑m​Ai​Fn−i​(n>m) 其中mmm尽量小 算法流程 文中的所有…

兰州大学第一届 飞马杯 ★★飞马祝福语★★ 线段树维护dp(动态dp)

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个串&#xff0c;每次将区间都修改为某一个字母&#xff0c;问最终包含多少个FeiMaFeiMaFeiMa子序列。 思路&#xff1a; 首先暴力修改肯定是不行的&#xff0c;复杂度nqnqnq。 如果没有修改操作&am…

兰州大学第一届 飞马杯 ★★快乐苹果树★★ 树链剖分 + 懒标记 + 树状数组

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 第一次听说树链剖分能在fa[top[i]]fa[top[i]]fa[top[i]]的地方加懒标记&#xff0c;学到了学到了。 首先不能被题目吓住&#xff0c;这个题目仔细剖析一下不难发现一些性质&#xff1a; 以…

【BZOJ3328】PYXFIB【矩阵快速幂】【单位根反演】【二项式定理】

传送门 题意&#xff1a; ∑i0⌊nk⌋(nik)Fik\sum_{i0}^{\lfloor\frac nk\rfloor}\binom n{ik}F_{ik}i0∑⌊kn​⌋​(ikn​)Fik​ FFF为斐波拉契数列 n≤1e18,k≤2e4,p≤1e9n\leq 1e18,k\leq2e4,p\leq1e9n≤1e18,k≤2e4,p≤1e9且为质数且模kkk余111 显然就是求 ∑i0n[k∣i](ni…

.NET分布式框架 | Orleans 知多少

引言公司物联网项目集成Orleans以支持高并发的分布式业务&#xff0c;对于Orleans也是第一次接触&#xff0c;本文就分享下个人对Orleans的理解。这里先抛出自己的观点&#xff1a;Orleans 是一个支持有状态云生应用/服务水平伸缩的基于Virtual Actor 模型的.NET分布式框架。下…

nowcoder 牛牛的最大兴趣组 质因子 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先nnn很小的话可以暴力连边&#xff0c;让后染个色求一个颜色最多的即可。但是这个题显然不行&#xff0c;由于是三次方&#xff0c;所以考虑质因子入手。 首先很容易就能想到将所有的数…

ASP.NET Core on K8S学习初探(3)部署API到K8S

“ 终于可以部署ASP.NET Core到K8S中了...”在上一篇《基本概念快速一览》中&#xff0c;我们把基本的一些概念快速地简单地不求甚解地过了一下&#xff0c;本篇开始我们会将ASP.NET Core WebAPI部署到K8S&#xff0c;从而结束初探的旅程。01—准备一个WebAPI这里准备一个空的A…

K-D Tree学习笔记

引入 K-D Tree 是一种处理高维空间的数据结构。 支持O(nk−1k)O(n^{\frac {k-1}k})O(nkk−1​)查询给定超矩形内的点的信息&#xff0c; kkk 为维数。 可以用替罪羊树的思想实现动态插入。 但其实用的最多的是在错误的复杂度内查询奇奇怪怪的点信息。 构建 K-D Tree 是平…

Codeforces Round #633 B. Edge Weight Assignment 结论题 + dp

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 思路&#xff1a; 首先考虑最小值&#xff0c;如果从一个叶子结点出发到任意叶子的距离都为偶数&#xff0c;那么只需要一个值就可以满足条件。如果有奇数的&#xff0c;考虑111 ^ 222 ^ 303030&#xff0c;…

ASP.NET Core Web API中使用Swagger

本节导航Swagger介绍在ASP.NET CORE 中的使用swagger在软件开发中,管理和测试API是一件重要而富有挑战性的工作。在我之前的文章《研发团队,请管好你的API文档》 也专门阐述了通过文档管理工具,来保证API文档和代码的一致性,这样更加有助于团队的协作。以往我们总是通过第三方平…

【ZJOI2015】幻想乡战略游戏【点分树】【带权重心】

题意&#xff1a;nnn个点带边权的树&#xff0c;动态修改点权viv_ivi​&#xff0c;最小化 钦定一个点xxx 后 ∑idist(x,i)∗vi\sum\limits_{i} dist(x,i)*v_ii∑​dist(x,i)∗vi​的值。 n,q≤105n,q \leq10^5n,q≤105&#xff0c;度数不超过202020 限制度数的树上的一些诡异…

P4559 [JSOI2018]列队 主席树

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你nnn个学生以及其位置&#xff0c;mmm个询问&#xff0c;每次询问[l,r][l,r][l,r]的学生跑到[k,kr−l][k,kr-l][k,kr−l]的位置的最小总的移动距离。当然不同的学生不能站在同一位置。 n,m≤5e5,1≤ai,k≤…

C# 8.x 先睹为快

尽管 C# 8.0 还需要几个月的时间才能发布&#xff0c;但是 C# 8.x 的特性计划已经开始了。该计划中的特性有一些是新提出的&#xff0c;而另外一些则是之前在 C# 8 中考虑过的。和往常一样&#xff0c;这个特性列表仍可能会发生变化。如之前的一个报告所讨论的那样&#xff0c;…

【WC2014】紫荆花之恋【替罪羊思想】【动态点分树】【替罪羊树】

传送门 题意&#xff1a;有一棵开始时没有结点的树&#xff0c;nnn次询问&#xff0c;每次新加一点并给定父结点、到父亲的距离、参数rir_iri​,并询问满足dist(u,v)≤rurvdist(u,v)\leq r_ur_vdist(u,v)≤ru​rv​的点对(u,v)(u,v)(u,v)的对数。 n≤105n\leq 10^5n≤105&…

NuGet是什么?理解与使用

如果你了解python&#xff0c;那么它类似pip。如果你了解nodejs&#xff0c;那么它类似npm。如果你了解ruby&#xff0c;那么它类似gem。对&#xff0c;它就是一个包&#xff08;package&#xff09;管理平台&#xff0c;确切的说是 .net平台的包管理工具&#xff0c;它提供了一…