题目：

　　　

开方：

不过还是不会二次剩余。

也不知道为什么取了 G(x)-B(x)=0 而不是 G(x)+B(x)=0。

式子是  B(x) = ( A(x) + G²(x) ) / 2*G(x) ，但写的时候这样比较方便：

令 D(x) = b^-1(x) , 则 B(x) =( a(x) * D(x) + b(x) ) / 2；

可以开一个 C(x) 表示 a(x) ，这样就不用动 a 数组了；因为加法和乘法都可以在系数上做，所以也可以不用把 b 数组 ntt 了。

 getinv 的时候可以用 A 和 B 表示 a 和 b ，这样就不用把 a 数组和 b 数组 ntt 了。

　　F(x) = D(x) * F²(x) + 1 ；其中+1是因为D没有常数项，所以不+1的话F也没有常数项，但 f [0]=0。

　　F(x) = ( 1+ sqrt(1-4*D(x)) ) / 2*D(x) ，这里要取减号，因为D没有常数项，如果分子剩下了常数项，就会除出余数，不太对。

　　写成  F(x) = 2 / 2*D(x)*sqrt( 1+4*D(x) )  比较方便。

#include
     
      #include
      
       #include
       
        #include
        
         #define ll long longusing namespace std;const int N=1e5+5,M=N<<2,mod=998244353;int a[M],b[M],A[M],B[M],C[M],D[M],len,r[M],inv2;int rdn(){  int ret=0;bool fx=1;char ch=getchar();  while(ch>'9'||ch<'0'){
   if(ch=='-')fx=0;ch=getchar();}  while(ch>='0'&&ch<='9') ret=ret*10+ch-'0',ch=getchar();  return fx?ret:-ret;}void upd(int &x){x>=mod?x-=mod:0;}int pw(int x,int k){
    int ret=1;while(k){
   if(k&1)ret=(ll)ret*x%mod;x=(ll)x*x%mod;k>>=1;}return ret;}void ntt(int *a,bool fx){  for(int i=0;i
         
          >1;i
          
           >1,a,b); len=n<<1; for(int i=0;i
           
            >1]>>1)+((i&1)?len>>1:0); for(int i=0;i
            
             >1,a,b); for(int i=0;i
             
              >1]>>1)+((i&1)?len>>1:0); ntt(C,0); ntt(D,0); for(int i=0;i