题目链接:https://ac.nowcoder.com/acm/contest/881/D

看此博客之前请先参阅吕凯飞的论文《集合幂级数的性质与应用及其快速算法》,论文中很多符号会被本文延用!

题目大意

  给定一个 n * m 的二维矩阵和 k,定义$count(x) = \sum\limits_{i = 1}^{n} \prod\limits_{j = 1}^{m} [v_{i, j} \& x 所表示的二进制位有奇数个一] $,求如下式子:

$$\begin{align*}
\bigoplus\limits_{x = 0}^{2^k - 1} (count(x) * 3^x mod (10^9 + 7))
\end{align*}$$

分析

  首先对于每一个数 x,给它的 k 位二进制位标号,从 1 ~ k,那么每一个数就可以唯一用一个集合 X 来表示,比如 k = 5, x = 10110,那么 X = {2, 3, 5}。
  定义 U 为全集,包含全部 1 ~ k。(为了方便,后面对应字母的大写就代表这个数对应的集合)
  于是我们可以重新定义 count(x) :$count(x) = count(X) = \sum\limits_{i = 1}^{n} \prod\limits_{j = 1}^{m} [V_{i, j} \cap X 有奇数个元素] $。
  进而:$count(X) = \frac{1}{2^m}\sum\limits_{i = 1}^{n} \prod\limits_{j = 1}^{m} (1 - (-1)^{|V_{i, j} \cap X|}) $,其中:$\prod\limits_{j = 1}^{m} (1 - (-1)^{|V_{i, j} \cap X|}) = 1 + \sum\limits_{j = 1}^{m} (-1)^{|X \cap V_{i, j}| + 1} + \sum\limits_{j_1 = 1}^{m} \sum\limits_{j_2 = 1}^{m} [j_1 \neq j_2] (-1)^{|X \cap V_{i, j_1} \cap V_{i, j_2}| + 2} + \dots + (-1)^{|X \cap (\bigcap\limits_{j = 1}^m V_{i, j})| + m}$
   又:$(-1)^{|Y|} * (-1)^{|X \cap T|} = (-1)^{|(X \cap T) \oplus Y|} = (-1)^{|(X \oplus Y) \cap (T \oplus Y)|}$
 
   ????????????????????(求大佬指点QAQ)
 
  所以$count(X) = \frac{1}{2^m}\sum\limits_{T \subseteq 2^U} f_T * (-1)^{T \cap X} = \frac{1}{2^m} * \hat{f_X}$
  于是我们只要先算出$f$,然后通过 FWT 算出所有 count(X) 就好了。
  时间复杂度为$O(n2^m + k2^k)$

代码如下

2019 牛客多校第一场 D Parity of Tuples的更多相关文章

  1. 2019牛客多校第一场 I Points Division(动态规划+线段树)

    2019牛客多校第一场 I Points Division(动态规划+线段树) 传送门:https://ac.nowcoder.com/acm/contest/881/I 题意: 给你n个点,每个点有 ...

  2. 2019牛客多校第一场E ABBA(DP)题解

    链接:https://ac.nowcoder.com/acm/contest/881/E 来源:牛客网 ABBA 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 524288K,其他语 ...

  3. 2019牛客多校第一场A-Equivalent Prefixes

    Equivalent Prefixes 传送门 解题思路 先用单调栈求出两个序列中每一个数左边第一个小于自己的数的下标, 存入a[], b[].然后按照1~n的顺序循环,比较 a[i]和b[i]是否相 ...

  4. 2019牛客多校第一场 A.Equivalent Prefixes

    题目描述 Two arrays u and v each with m distinct elements are called equivalent if and only if RMQ(u,l,r ...

  5. 2019牛客多校第一场 E-ABBA(dp)

    ABBA 题目传送门 解题思路 用dp[i][j]来表示前i+j个字符中,有i个A和j个B的合法情况个数.我们可以让前n个A作为AB的A,因为如果我们用后面的A作为AB的A,我们一定也可以让前面的A对 ...

  6. 【2019牛客多校第一场】XOR

    题意: 给你一个集合A,里边有n个正整数,对于所有A的.满足集合内元素异或和为0的子集S,问你∑|S| n<=1e5,元素<=1e18 首先可以转化问题,不求∑|S|,而是求每个元素属于子 ...

  7. 2019 牛客多校第一场 B Integration

    题目链接:https://ac.nowcoder.com/acm/contest/881/B 题目大意 给定 n 个不同的正整数 ai,求$\frac{1}{\pi}\int_{0}^{\infty} ...

  8. 2019牛客多校第一场E ABBA 贪心 + DP

    题意:问有多少个有(n + m)个A和(n + m)个B的字符串可以凑出n个AB和m个BA. 思路:首先贪心的发现,如果从前往后扫,遇到了一个A,优先把它看成AB的A,B同理.这个贪心策略用邻项交换很 ...

  9. 2019 牛客多校第一场 F Random Point in Triangle

    题目链接:https://ac.nowcoder.com/acm/contest/881/F 题目大意 给定二维平面上 3 个整数表示的点 A,B,C,在三角形 ABC 内随机选一点 P,求期望$E ...

随机推荐

  1. java中获取类加载路径和项目根路径的5种方式分析

    package my; import Java.io.File; import java.io.IOException; import java.net.URL; public class MyUrl ...

  2. DNS域名服务器的搭建

    父域的DNS(svr7): 可以解析父域名下保存的域名地址,即解析.tedu.cn下的域名 一.安装软件包bind.bind-chroot   bind是DNS解析服务需要用到的服务软件包,bind- ...

  3. 关于ios 的耳机线线控

    如题 线控主要就是指的用ios耳机的  中间那个键来自定义一些程序自己的事件(比如玩吃鸡的时候,作为开枪扳机) 而想要控制耳机的音量加减的事件, 有一种方法,就是监听 系统音量的变化,根据音量变化来触 ...

  4. This inspection detects shadowing names defined in outer scopes.

    错误信息:This inspection detects shadowing names defined in outer scopes. 检查到波浪处的单词已在函数外部定义. 解决:使用global ...

  5. [转] JPA 2.0 with EclipseLink - 教程

    原文: http://www.vogella.com/articles/JavaPersistenceAPI/article.html Lars Vogel Version 2.2 Copyright ...

  6. CentOS6 图形界面(gnome)安装

    CentOS6相对于CentOS5的安装有了不少的进步,有不少默认的选项可以选择,如: Desktop :基本的桌面系统,包括常用的桌面软件,如文档查看工具. Minimal Desktop :基本的 ...

  7. maven surefire入门

    一.maven常用命令: mvn compile mvn install mvn test mvn clean mvn dependency:resolve -X   #查看完整的debug信息!!! ...

  8. JAR API

    JAR API包括使用 manifest 文件的类.Manifest类的一个对象表示一个manifest文件. 在代码中创建一个Manifest对象,如下所示: 1 Manifest manifest ...

  9. Linux下NODE配置NODE_PATH变量

    在当我们使用了npm install express -g安装之后的,会安装在全局目录下,但是如果直接用Node来启动程序的话,是无法发现这个包的,原因就在于NODE_PATH这个环境变量设置的不正确 ...

  10. C语言指针函数和函数指针

    #include <stdio.h> char *test(); void test1(); int main() { /********************************* ...