FWT在三种位运算下都满足FWT(a×b)=FWT(a)*FWT(b)
其中or卷积和and卷积还可以通过FMT实现(本质上就是个高维前缀和)

#include<bits/stdc++.h>
#define N 1100000
#define eps 1e-7
#define inf 1e9+7
#define db double
#define ll long long
#define ldb long double
using namespace std;
inline int read()
{char ch=0;int x=0,flag=1;while(!isdigit(ch)){ch=getchar();if(ch=='-')flag=-1;}while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}return x*flag;
}
const ll mo=998244353;
int ksm(int x,int k)
{int ans=1;while(k){if(k&1)ans=1ll*ans*x%mo;k>>=1;x=1ll*x*x%mo;}return ans;
}
void fwt_or(int *f,int n,int flag)
{for(int k=2,kk=1;k<=(1<<n);k<<=1,kk<<=1)for(int i=0;i<(1<<n);i+=k)for(int j=0;j<kk;j++)f[i+j+kk]=(f[i+j+kk]+flag*f[i+j])%mo;
}
void fwt_and(int *f,int n,int flag)
{for(int k=2,kk=1;k<=(1<<n);k<<=1,kk<<=1)for(int i=0;i<(1<<n);i+=k)for(int j=0;j<kk;j++)f[i+j]=(f[i+j]+flag*f[i+j+kk])%mo;
}
void fwt_xor(int *f,int n,int flag)
{for(int k=2,kk=1;k<=(1<<n);k<<=1,kk<<=1)for(int i=0;i<(1<<n);i+=k)for(int j=0;j<kk;j++){int t=f[i+j+kk];f[i+j+kk]=(f[i+j]-t+mo)%mo;f[i+j]=(f[i+j]+t)%mo;}if(flag==-1){int inv=ksm(1<<n,mo-2);for(int i=0;i<(1<<n);i++)f[i]=1ll*f[i]*inv%mo;}
}
int a[N],b[N];
void Or(int n)
{fwt_or(a,n,+1);fwt_or(b,n,+1);for(int i=0;i<(1<<n);i++)a[i]=1ll*a[i]*b[i]%mo;fwt_or(a,n,-1);
}
void And(int n)
{fwt_and(a,n,+1);fwt_and(b,n,+1);for(int i=0;i<(1<<n);i++)a[i]=1ll*a[i]*b[i]%mo;fwt_and(a,n,-1);
}
void Xor(int n)
{fwt_xor(a,n,+1);fwt_xor(b,n,+1);for(int i=0;i<(1<<n);i++)a[i]=1ll*a[i]*b[i]%mo;fwt_xor(a,n,-1);
}
int A[N],B[N];
int main()
{int n=read();for(int i=0;i<(1<<n);i++)A[i]=read();for(int i=0;i<(1<<n);i++)B[i]=read();for(int i=0;i<(1<<n);i++)a[i]=A[i],b[i]=B[i];Or(n);for(int i=0;i<(1<<n);i++)printf("%lld ",(a[i]%mo+mo)%mo);printf("\n");for(int i=0;i<(1<<n);i++)a[i]=A[i],b[i]=B[i];And(n);for(int i=0;i<(1<<n);i++)printf("%lld ",(a[i]%mo+mo)%mo);printf("\n");for(int i=0;i<(1<<n);i++)a[i]=A[i],b[i]=B[i];Xor(n);for(int i=0;i<(1<<n);i++)printf("%lld ",(a[i]%mo+mo)%mo);printf("\n");return 0;
}

转载于:https://www.cnblogs.com/Creed-qwq/p/10456765.html

【模板/经典题型】FWT相关推荐

  1. 【模板/经典题型】并查集维护生成树

    这里的并查集是按秩合并并查集. 这种方法维护生成树的时候可以维护一个点到根的权值. 但是由于合并的时候做了一个类似换根一样的操作,因此这个权值通常只能是异或之类的. //find-union-set ...

  2. 【模板/经典题型】闵可夫斯基和

    闵可夫斯基和,是两个欧几里得空间的点集的和,以德国数学家闵可夫斯基命名. 点集A与B的闵可夫斯基和就是{o|o=a+b},其中a属于A,b属于B. 对于凸包这种特殊的图形,它的闵可夫斯基和有一些较好的 ...

  3. JavaScript中的经典题型(类数组、CSS Sprites、事件委托、经典去重、原型链、闭包、深浅克隆、附带思路流程和源码)

    JavaScript中的经典题型 一.JavaScript中的经典题型 1..如何判断一个数组和类数组? 首先要明白什么是类数组. 类数组:类数组是一个普通对象,他的原型是Object.而真实的数组是 ...

  4. pwn学习总结(五) —— 堆溢出经典题型整理

    pwn学习总结(五) -- 堆溢出经典题型整理 fastbin + 栈溢出 fastbin + 函数构造 fastbin + 堆执行 fastbin + malloc_hook fastbin + 栈 ...

  5. pwn学习总结(三) —— 栈溢出经典题型整理

    pwn学习总结(三) -- 栈溢出经典题型整理 ret2text ret2shellcode rop ret2libc 使用DynELF实现远程libc泄露 ret2syscall ret2libc ...

  6. java 品尝饮料,java细节经典题型

    28. 选项中哪一行代码可以替换题目中//add code here 而不产生编译错误?() [java] view plaincopy 1. public abstract class MyClas ...

  7. mysql中select 的题型_MYSQL经典题型详情解析

    学完了mysql后发现有很多地方不是很明白,于是总结了mysql的经典题型,不论是工作还是面试,我相信还是有一定帮助的. 例题一 ​ 在我的数据库中数据如下(排序有些差别,但是不影响结果) ​ 分析: ...

  8. 高中数学排列组合二项定理经典题型汇编(名师总结)

    今天老师给同学们整理出高中数学二项式定理相关内容,该题在高考中常以选择题或填空题形式出现,难度为容易或中等.下面就一起来看看二项式定理主干知识,以及常考题型~ 一.知识点梳理 二.经典题型汇总 今天的 ...

  9. 高中数学函数奇偶性知识点经典题型解题技巧【提分秘籍】

    作者:vxbomath 今天分享高中数学函数奇偶性的知识点技巧一例题讲解:通过老师的的分享希望能给同学们带来帮助! 知识点: 函数的奇偶性一.函数奇偶性的定义: (1)如果对于函数f(x)定义域内的任 ...

最新文章

  1. Target runtime Apache Tomcat 6.0 is not defined
  2. python内积 卷积
  3. 1400协议是什么和28181区别_舞阳钢铁生产30CrMnSiA国军标、国标、舞钢技术协议有什么区别...
  4. C指针原理(14)-C指针基础
  5. 关于C语言中递归的一点点小问题
  6. java 什么是耦合_什么是耦合、解耦
  7. led显示屏控制卡接线图解_Led显示屏出现花屏是什么原因
  8. 如何在Oracle官网上下载历史版本的JDK
  9. [下载地址] Subclipse 1.10.9(SVN) _附说明
  10. 选择排序 - python实现详解
  11. SQL调优(SQL TUNING)之远程支持完成性能大幅优化
  12. IT Library第4期《备份域升为主域控制器》
  13. 全球十大数据治理解决方案提供商
  14. 通过企业微信Api向企业微信用户发送消息的实现
  15. 远程桌面计算机正在使用,远程桌面使用问题
  16. Win10 更改不了注册表,重启之后恢复原样
  17. mac怎么验机,都应该查什么
  18. 逻辑设计基础_芯片设计--TCAM概述
  19. Structure of a Google Docs document 谷歌文档的结构
  20. 日语基础复习 Day 16

热门文章

  1. Hibernate标识符属性(主键)生成策略全析
  2. 监听ListView滚动到最底部
  3. Lisp语言: 在Windows下搭建CLisp环境
  4. Java动态调用方法
  5. 后台开发经典书籍--深入理解Nginx模块开发与架构
  6. linux命令lscpu
  7. MySQL第3天:MySQL的架构介绍之linux版安装
  8. C语言中:在头文件中使用static定义变量意味着什么?
  9. UART0串口编程(五):串口编程(UART0)之UC/OS(一)UC/OS下的串口发送任务编程
  10. 常见的局域网拓扑结构