博客作者分享了在解决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);
扫一扫
540
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
