记得在求逆的时候，对 $x^n$ 取模而不是 $x^{\frac{n}{2}}$。

void PolySqrt(vector<int>& a){
    if(a.size()==1u) return;
    vector<int> Sq((a.size()+1)>>1);
    for(int i=0;i<(int)Sq.size();i++) Sq[i]=a[i];
    PolySqrt(Sq);
    vector<int> Iv=Sq;
    Iv.resize(a.size());/*here*/
    PolyInv(Iv);
    sqrtconvol(a,Sq,Iv);
}