洛谷 - P4717 【模板】快速莫比乌斯/沃尔什变换 (FMT/FWT)

题目链接：点击查看

题目分析：看不懂原理我爪巴了

代码：

// Problem: P4717 【模板】快速莫比乌斯/沃尔什变换 (FMT/FWT)
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P4717
// Memory Limit: 250 MB
// Time Limit: 1000 ms
//
// Powered by CP Editor (https://cpeditor.org)// #pragma GCC optimize(2)
// #pragma GCC optimize("Ofast","inline","-ffast-math")
// #pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<cassert>
#include<bitset>
#include<list>
#include<unordered_map>
#define lowbit(x) (x&-x)
using namespace std;
typedef long long LL;
typedef unsigned long long ull;
template<typename T>
inline void read(T &x)
{T f=1;x=0;char ch=getchar();while(0==isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}while(0!=isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();x*=f;
}
template<typename T>
inline void write(T x)
{if(x<0){x=~(x-1);putchar('-');}if(x>9)write(x/10);putchar(x%10+'0');
}
const int inf=0x3f3f3f3f;
const int N=1e6+100;
const int mod=998244353;
int len;
LL a[N],b[N],A[N],B[N];
void FWTor(LL *f,int x)
{for(int mid=1;(mid<<1)<=len;mid<<=1){int R=mid<<1;for(int i=0;i<len;i+=R)for(int j=0;j<mid;j++)f[i+j+mid]=(f[i+j+mid]+f[i+j]*x+mod)%mod;}
}
void FWTand(LL *f,int x)
{for(int mid=1;(mid<<1)<=len;mid<<=1){int R=mid<<1;for(int i=0;i<len;i+=R)for(int j=0;j<mid;j++)f[i+j]=(f[i+j]+f[i+j+mid]*x+mod)%mod; }
}
void FWTxor(LL *f,int x)
{for(int mid=1;(mid<<1)<=len;mid<<=1){int R=mid<<1;for(int i=0;i<len;i+=R)for(int j=0;j<mid;j++){f[i+j]=(f[i+j]+f[i+j+mid])%mod;f[i+j+mid]=(f[i+j]-f[i+j+mid]+mod-f[i+j+mid]+mod)%mod;f[i+j]=f[i+j]*x%mod;f[i+j+mid]=f[i+j+mid]*x%mod;}}
}
void copy() {for(int i=0;i<len;i++) {A[i]=a[i],B[i]=b[i];}
}
void merge() {for(int i=0;i<len;i++) {A[i]=A[i]*B[i]%mod;}
}
void print() {for(int i=0;i<len;i++) {printf("%lld ",A[i]);}puts("");
}
int main()
{#ifndef ONLINE_JUDGE
//  freopen("data.in.txt","r",stdin);
//  freopen("data.out.txt","w",stdout);
#endif
//  ios::sync_with_stdio(false);read(len);len=1<<len;for(int i=0;i<len;i++) {read(a[i]);}for(int i=0;i<len;i++) {read(b[i]);}copy(),FWTor(A,1),FWTor(B,1),merge(),FWTor(A,-1),print();copy(),FWTand(A,1),FWTand(B,1),merge(),FWTand(A,-1),print();copy(),FWTxor(A,1),FWTxor(B,1),merge(),FWTxor(A,(mod+1)>>1),print();return 0;
}

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

专题·树链剖分【including 洛谷·【模板】树链剖分
初见安~~~终于学会了树剖~~~ [兴奋]当初机房的大佬在学树剖的时候我反复强调过:"学树剖没有前途的!!!" 恩.真香. 一.重链与重儿子所谓树剖--树链剖分,就是赋予一个链的 ...
洛谷·【模板】点分树 | 震波【including 点分树
初见安-这里是传送门:洛谷P6329 [模板]点分树 | 震波一.点分树其实你会点分治的话,点分树就是把点分治时的重心提出来重新连城一棵树. 比如当前点是u,求出子树v的重心root后将root与 ...
洛谷.4897.[模板]最小割树(Dinic)
题目链接最小割树模板.具体见:https://www.cnblogs.com/SovietPower/p/9734013.html. ISAP不知为啥T成0分了.. Dinic: //1566ms ...
强连通分量：洛谷P3387 模板：缩点
传送门顾名思义,模板awa #include <cstdio> #include <cstring> #include <cmath> #include < ...
【后缀数组】洛谷P3809模板题
题目背景这是一道模板题. 题目描述读入一个长度为 n n n 的由大小写英文字母或数字组成的字符串,请把这个字符串的所有非空后缀按字典序从小到大排序,然后按顺序输出后缀的第一个字符在原串中的位置. ...
洛谷 P1919 模板】A*B Problem升级版（FFT快速傅里叶）
https://www.luogu.com.cn/problem/P1919 题目背景本题数据已加强,请使用 FFT/NTT,不要再交 Python 代码浪费评测资源. 题目描述给你两个正整数 a ...
洛谷 | P1226 【快速幂】
快速幂题目链接题目描述输入b,p,k的值,求b^p mod k的值.其中b,p,k*k为长整型数. 输入输出格式时空限制时间:1000ms 空间:128MB 代码 #include < ...
洛谷 p3372 模板-线段树 1
题目描述如题,已知一个数列,你需要进行下面两种操作: 1.将某区间每一个数加上x 2.求出某区间每一个数的和输入输出格式输入格式: 第一行包含两个整数N.M,分别表示该数列数字的个数和操作的总个 ...
洛谷.4245.[模板]任意模数NTT(MTT/三模数NTT)
题目链接三模数\(NTT\): 就是多模数\(NTT\)最后\(CRT\)一下...下面两篇讲的都挺明白的. https://blog.csdn.net/kscla/article/details/ ...

洛谷 - P4717 【模板】快速莫比乌斯/沃尔什变换 (FMT/FWT)

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

最新文章

热门文章