[hdu4656] Evaluation - FFT,CRT

最新推荐文章于 2024-09-04 16:36:43 发布

原创

最新推荐文章于 2024-09-04 16:36:43 发布 · 1.2k 阅读

1 ·

CC 4.0 BY-SA版权

博客作者分享了在解决HDU4656问题时遇到的困难，提到式子推导的复杂性，并提及通过参考叉姐的课件来理解和应用快速傅里叶变换(FFT)及中国剩余定理(CRT)解决问题的过程。

式子真难推…还是看了叉姐的课件

#include"bits/stdc++.h"

using namespace std;
typedef long long LL;

const int Mod=1000003,m1=998244353,m2=1004535809,pwMod=Mod-1;
const LL P=1002772198720536577LL;

inline int CRT(int r1,int r2){
	return (((r1-r2)*334845110LL%m2*m1+r1)%P+P)%P%Mod;
}

int power(int a,int t,int P){
	int r=1;
	while(t){
		if(t&1)r=(LL)r*a%P;
		a=(LL)a*a%P;t>>=1;
	}
	return r;
}

int _wn[25];

void NTT(int*A,int len,int dft,int P){
	int i,j=len>>1,k,l,w,wn,c=0,u,v;
	for(i=1;i<len-1;i++){
		if(i<j)swap(A[i],A[j]);
		for(k=len>>1;(j^=k)<k;k>>=1);
	}
	for(i=0;i<=21;i++)_wn[i]=power(3,P-1>>i,P);