式子真难推…还是看了叉姐的课件
#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&