[hdu4656] Evaluation - 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&