实现了一下Berlekamp-Massey

//from https://www.cnblogs.com/TSHugh/p/9265155.html 
//在FP中求固定项数数列的线性递推式
//此递推式严格符合数学定义,故可能在末尾出现一些看起来奇怪但是在数学上必要的0,虽然这些0可能在实际操作时没有用处
#pragma GCC optimize("O3")
#include <cstdio>
#include <cstring>
#include <algorithm>
typedef long long LL;
const int N=;
const int P=;
inline int Pow(int x,int y){
  int ret=;
  while(y){
    if(y&)ret=(LL)ret*x%P;
    x=(LL)x*x%P;y>>=;
  }
  return ret;
}
int n,a[N],pos,len,wa[N],cnt,cur[N],tmp[N];
int main(){
  int i,j,x,s;
  scanf("%d",&n);
  for(i=;i<=n;++i)
    scanf("%d",a+i);
  for(i=;i<=n;++i)
    if(a[i])break;
  if(i==n+)return ;
  pos=i;
  len=;
  wa[]=Pow(a[i],P-);
  cnt=i;
  for(++i;i<=n;++i){
    s=;
    for(j=;j<=cnt;++j)
      s=(s+(LL)a[i-j]*cur[j])%P;
    s=(s+P-a[i])%P;
    if(s==)continue;
    memcpy(tmp,cur,cnt+<<);
    for(j=;j<=len;++j)
      cur[i-pos+j-]=(cur[i-pos+j-]+(LL)wa[j]*(P-s))%P;
    wa[]=P-;
    for(j=;j<=cnt;++j)
      wa[j+]=tmp[j];
    x=len;
    len=cnt+;
    s=Pow(s,P-);
    for(j=;j<=len;++j)
      wa[j]=(LL)wa[j]*s%P;
    cnt=std::max(cnt,i-pos+x-);
    pos=i;
  }
  for(i=;i<=cnt;++i)
    printf("%d ",cur[i]);
  puts("");
  return ;
}