真的分治fft

以前学的分治fft f[i]=sigma(f[i-x]*g[x]),其中g[x]已知

那么我们可以用cdq分治来做（l,mid 对mid+1,t的影响）

而现在的$f[i]=sum(f(i-x)*f(x))$

我们如果沿用刚才的方法 会发现有$f(t-h)$这一项

而$t-h>mid$是有可能的

所以我们要在后续处理这件事情

先将$f[l,mid]*f[l,mid]$乘起来

如果$t-h<h$ 还要算$f[1,t-h]*f[h,mid]$ 注意还要乘2

注意多次用fft 每次还原a,b数组 因为那个n是要比m大的，所以只清空到m的0是有问题的

1e5的数据开O2要2s 预料之中了

#include <bits/stdc++.h>
using namespace std;
#define rint register int
#define IL inline
#define rep(i,h,t) for (rint i=h;i<=t;i++)
#define dep(i,t,h) for (rint i=t;i>=h;i--)
struct cp{
  double a,b;
};
const int N=3e5+;
int n,m,l,r[N],f[N];
int a[N],b[N],w[N],G=;
const int mo=;
const int p=mo;
IL int fsp(rint x,rint y)
{
  rint ans=;
  while (y)
  {
    if (y&) ans=(1ll*ans*x)%mo;
    x=(1ll*x*x)%mo;
    y>>=;
  }
  return ans;
}
void fft(int *a,int o)
{
    for (int i=;i<n;i++)
      if (i>r[i]) swap(a[i],a[r[i]]);
    for (rint i=;i<n;i*=)
    {
        rint wn=fsp(G,(p-)/(i*)); w[]=;
        rep(j,,i-) w[j]=(1ll*w[j-]*wn)%p;
        for (rint j=;j<n;j+=(i*))
        {
            rint *x=a+j,*y=a+i+j;
            for (rint k=;k<i;k++)
            {
                const int t=(1ll*w[k]*y[k])%p;
                y[k]=x[k]-t;
                if (y[k]<) y[k]+=p;
                if (y[k]>p) y[k]-=p;
                x[k]=x[k]+t;
                if (x[k]>p) x[k]-=p;
            }
        }
    }
    if (o==-)
    {
      reverse(&a[],&a[n]);
      for (int i=,inv=fsp(n,p-);i<n;i++)
        a[i]=1ll*a[i]*inv%p;
    }
}
void query()
{
    l=;
    for (n = ; n <= m; n <<= ) l++;
    for (int i=;i<n;i++) r[i]=(r[i/]/)|((i&)<<(l-));
    fft(a,),
    fft(b,);
    for (int i=;i<n;i++) a[i]=1ll*a[i]*b[i]%p;
    fft(a,-);
}
#define mid ((h+t)/2)
IL void js(rint &x,rint y)
{
  x+=y;
  x%=mo;
}
void cdq_fz(int h,int t)
{
  if (h==t) return;
  cdq_fz(h,mid);
  if (t-h<h)
  {
    n=m=t-h+(mid-h+);
    rep(i,,t-h) a[i-]=f[i];
    rep(i,h,mid) b[i-h]=f[i];
    query();
    rep(i,mid+,t)
      if (i-h->=) js(f[i],(*a[i-h-])%mo);
    rep(i,,n) a[i]=b[i]=;
  }
  n=m=*(mid-h+);
  rep(i,h,mid) a[i-h]=b[i-h]=f[i];
  query();
  rep(i,mid+,t)
    if (i-(*h+)>=) js(f[i],a[i-(*h+)]);
  rep(i,,n) a[i]=b[i]=;
  cdq_fz(mid+,t);
}
int main()
{
  freopen("1.in","r",stdin);
  freopen("1.out","w",stdout);
  ios::sync_with_stdio(false);
  int k;
  cin>>k;
  f[]=; f[]=;
  cdq_fz(,k);
  cout<<f[k]<<endl;
  return ;
}

对拍

#include <bits/stdc++.h>
using namespace std;
#define rep(i,h,t) for (int i=h;i<=t;i++)
const int mo=;
int f[];
int main()
{
  freopen("1.in","r",stdin);
  freopen("2.out","w",stdout);
  ios::sync_with_stdio(false);
  int n;
  cin>>n;
  f[]=; f[]=;
  rep(i,,n)
    rep(k,,i-)
      f[i]+=(1ll*f[k]*f[i-k-])%mo,f[i]%=mo;
  cout<<f[n]<<endl;
  return ;
}

还有说多项式右移

f=f*f>>1+g