正题

题目链接:https://www.luogu.com.cn/problem/P4717

题目大意

给出两个长度为2n2^n2n的数列A,BA,BA,B求
Cn=∑iorj=nAiBjC_{n}=\sum_{i\ or\ j=n}A_iB_jCn=i or j=n∑AiBj
Cn=∑iandj=nAiBjC_{n}=\sum_{i\ and\ j=n}A_iB_jCn=i and j=n∑AiBj
Cn=∑ixorj=nAiBjC_{n}=\sum_{i\ xor\ j=n}A_iB_jCn=i xor j=n∑AiBj

解题思路

和FFTFFTFFT一样的思路，我们需要将A,BA,BA,B转换为点集表达式相乘再转回去，这里就直接抛结论了。

对于ororor有FWT(A)=(FWT(A),FWT(A)+FWT(B))FWT(A)=(FWT(A),FWT(A)+FWT(B))FWT(A)=(FWT(A),FWT(A)+FWT(B))

IFWT(A)=(IFWT(A),IFWT(A)−IFWT(B))IFWT(A)=(IFWT(A),IFWT(A)-IFWT(B))IFWT(A)=(IFWT(A),IFWT(A)−IFWT(B))
对于andandand有FWT(A)=(FWT(A)+FWT(B),FWT(B))FWT(A)=(FWT(A)+FWT(B),FWT(B))FWT(A)=(FWT(A)+FWT(B),FWT(B))

IFWT(A)=(IFWT(A)−IFWT(B),IFWT(B))IFWT(A)=(IFWT(A)-IFWT(B),IFWT(B))IFWT(A)=(IFWT(A)−IFWT(B),IFWT(B))
对于xorxorxor有FWT(A)=(FWT(A)+FWT(B),FWT(A)−FWT(B))FWT(A)=(FWT(A)+FWT(B),FWT(A)-FWT(B))FWT(A)=(FWT(A)+FWT(B),FWT(A)−FWT(B))

IFWT(A)=(IFWT(A)+IFWT(B)2,IFWT(A)−IFWT(B)2)IFWT(A)=(\frac{IFWT(A)+IFWT(B)}{2},\frac{IFWT(A)-IFWT(B)}{2})IFWT(A)=(2IFWT(A)+IFWT(B),2IFWT(A)−IFWT(B))
因为不用转成虚数所以代码好写很多，时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)

codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll inv2=499122177,N=(1<<18),P=998244353;
ll lg,n,aa[N],bb[N],a[N],b[N],c[3][N];void mul(ll *c,ll *a,ll *b){for(ll i=0;i<n;i++)c[i]=a[i]*b[i]%P;return;
}void FWT_or(ll *f,ll op){for(ll p=2;p<=n;p<<=1){ll len=p>>1;for(ll k=0;k<n;k+=p)for(ll i=k;i<k+len;i++)(f[i+len]+=f[i]*op+P)%=P;}return;
}
void solve_or(ll *c){memcpy(a,aa,sizeof(a));memcpy(b,bb,sizeof(b));FWT_or(a,1);FWT_or(b,1);mul(c,a,b); FWT_or(c,P-1);return;
}void FWT_and(ll *f,ll op){for(ll p=2;p<=n;p<<=1){ll len=p>>1;for(ll k=0;k<n;k+=p)for(ll i=k;i<k+len;i++)(f[i]+=f[i+len]*op+P)%=P;}return;
}
void solve_and(ll *c){memcpy(a,aa,sizeof(a));memcpy(b,bb,sizeof(b));FWT_and(a,1);FWT_and(b,1);mul(c,a,b); FWT_and(c,P-1);return;
}void FWT_xor(ll *f,ll op){for(ll p=2;p<=n;p<<=1){ll len=p>>1;for(ll k=0;k<n;k+=p)for(ll i=k;i<k+len;i++){ll x=f[i],y=f[i+len];f[i]=(x+y)*op%P;f[i+len]=(x-y+P)*op%P;}}return;
}
void solve_xor(ll *c){memcpy(a,aa,sizeof(a));memcpy(b,bb,sizeof(b));FWT_xor(a,1);FWT_xor(b,1);mul(c,a,b); FWT_xor(c,inv2);return;
}signed main()
{scanf("%lld",&lg);n=1<<lg;for(ll i=0;i<n;i++)scanf("%lld",&aa[i]);for(ll i=0;i<n;i++)scanf("%lld",&bb[i]);solve_or(c[0]);solve_and(c[1]);solve_xor(c[2]);for(ll i=0;i<n;i++)printf("%lld ",c[0][i]);putchar('\n');for(ll i=0;i<n;i++)printf("%lld ",c[1][i]);putchar('\n');for(ll i=0;i<n;i++)printf("%lld ",c[2][i]);putchar('\n');return 0;
}

P4717-[模板]快速莫比乌斯/沃尔什变换(FMT/FWT)相关推荐

洛谷 - P4717 【模板】快速莫比乌斯/沃尔什变换 (FMT/FWT)
题目链接:点击查看题目分析:看不懂原理我爪巴了代码: // Problem: P4717 [模板]快速莫比乌斯/沃尔什变换 (FMT/FWT) // Contest: Luogu // URL: ...
FFT/NTT/FMT/FWT题目
快速傅里叶变换FFT hdu4609 3-idiots P4173 残缺的字符串 Hash Function [AH2017/HNOI2017]礼物 CF993E Nikita and Order S ...
模板：快速莫比乌斯变换（FMT）+快速沃尔什变换（FWT）（多项式）
文章目录前言解析 OR 定义变换: 逆变换代码 AND 代码 XOR 定义变换逆变换代码所谓快速沃尔什变换,就是快速的沃尔玛什锦专柜变换 (逃) 前言正常卷积的定义:ck=∑i+j= ...
[模板] 快速沃尔什变换
快速沃尔什变换快速沃尔什变换是求这样的式子: 对序列 \(A\), \(B\), 求序列 \(C\), 使得 \[ C_{i}=\sum_{j \oplus k} A_{j} B_{k} \] 其中 ...
快速沃尔什变换（FWT）介绍
快速沃尔什变换(FWT)介绍众所周知,FFT可以有效地解决如下的问题 A=(a1,-,an)B=(b1,...,bn)C=(c1,...,cn)ck=∑i+j=kai∗bjA=(a_1,\dots, ...
T4模板：MVC中用T4模板快速生成代码
T4模板快速生成代码: 以快速生Dal文件为例,下面为T4模板文件的内容 <#@ template debug="false" hostspecific="true ...
php快速搭建后台,基于thinkphp的后台管理系统模板快速搭建，thinkphp后台模板_PHP教程...
基于thinkphp的后台管理系统模板快速搭建,thinkphp后台模板当我们在搭建网站的时候,后端开发人员在编写后台的管理系统的时候,往往会因为缺少一个合适的后台管理系统的模板,而必须去重新编写一 ...
使用模板快速制作H5网页（仅供参考）
#[转载] 使用模板快速制作H5网页(仅供参考) 教程来源: 零基础30分钟学会制作网站--如何使用HTML网页模板(更新P2字幕版) 作者:SolidifieRay 1.HTML整体框架: < ...
SpringBoot集成beetl模板快速入门
SpringBoot集成beetl模板快速入门本文引导 SpringBoot集成beetl模板快速入门创建SpringBoot工程(idea) 新建工程修改maven本地仓库位置配置Sprin ...
快速沃尔什变换（FWT）讲解解决集合卷积的方法
由于某次CSDN更新,文章样式已经丢失,大量公式无法正常显示,本博客也不再维护. 请到快速沃尔什变换FWT查看,有更加完整的公式讲解. 能看到这篇博客的人,一定知道FWT是干什么的.(什么?你不知道? ...

P4717-[模板]快速莫比乌斯/沃尔什变换(FMT/FWT)

正题

题目大意

解题思路

codecodecode

P4717-[模板]快速莫比乌斯/沃尔什变换(FMT/FWT)相关推荐

最新文章

热门文章