ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

FWT整理

2021-08-23 22:33:00  阅读:179  来源: 互联网

标签:UFWT merge int Aj 整理 FWT oplus


FWT

by AmanoKumiko

1.Or

\[A_i=∑_{j|i=i}Aj \]

\[FWT[A]=merge(FWT[A_0],FWT[A_0]+FWT[A_1]) \]

\[UFWT[A]=merge(UFWT[A_0],UFWT[A_1]-UFWT[A_0]) \]


2.And

\[A_i=∑_{j\ and\ i=i}Aj \]

\[FWT[A]=merge(FWT[A_0]+FWT[A_1],FWT[A_1]) \]

\[UFWT[A]=merge(UFWT[A_0]-UFWT[A_1],UFWT[A_1]) \]


3.Xor

\[i\oplus j=popcount(i\ and \ j)\ mod\ 2 \]

\[A_i=∑_{i\oplus j=0}A_j-∑_{i\oplus j=1}A_j \]

\[FWT[A]=merge(FWT[A_0]+FWT[A_1],FWT[A_0]-FWT[A_1]) \]

\[UFWT[A]=merge(\frac{UFWT[A_0]+UFWT[A_1]}{2},\frac{UFWT[A_0]-UFWT[A_1]}{2}) \]


4.code

struct poly{
	int len;
	LL val[N];
	void Or(int opt){
		for(int mid=2;mid<len;mid<<=1){
			for(int i=0;i<len;i+=mid){
				F(j,i,i+mid/2)val[j+mid/2]+=val[j]*opt;
			}
		}
	}
	void And(int opt){
		for(int mid=2;mid<len;mid<<=1){
			for(int i=0;i<len;i+=mid){
				F(j,i,i+mid/2)val[j]+=val[j+mid/2]*opt;
			}
		}
	}
	void Xor(int opt){
		for(int mid=2;mid<len;mid<<=1){
			for(int i=0;i<len;i+=mid){
				F(j,i,i+mid/2){
					LL l=val[j],r=val[j+mid/2];
					val[j]=(l+r)/(opt==1?1:2);val[j+mid/2]=(l-r)/(opt==1?1:2);
				}
			}
		}
	}
}

标签:UFWT,merge,int,Aj,整理,FWT,oplus
来源: https://www.cnblogs.com/AmanoKumiko/p/15178009.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有