FTT续写之范德蒙逆矩阵

求范德蒙逆矩阵也是一个大坑,用到的数学知识也很重要,拉格朗日插值法,最重要是代数基本定理写的有点丑陋,代数基本定理可以自己尝试证明下。虽然整个推导流程清晰明了,但是一些细节让人头疼,比如分解时候用到泰勒展开。

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
FTT算法,全称为快速傅里叶变换(Fast Fourier Transform),是一种快速计算傅里叶变换的算法。以下是C++实现的FTT算法程序: ```c++ #include<bits/stdc++.h> using namespace std; const double PI=acos(-1.0); const int MAXN=1<<18; struct Complex{ double real,imag; Complex operator + (const Complex &t) const{ return Complex{real+t.real,imag+t.imag}; } Complex operator - (const Complex &t) const{ return Complex{real-t.real,imag-t.imag}; } Complex operator * (const Complex &t) const{ return Complex{real*t.real-imag*t.imag,real*t.imag+imag*t.real}; } Complex operator / (const double &t) const{ return Complex{real/t,imag/t}; } }a[MAXN],b[MAXN]; int n,m,limit,r[MAXN]; void FFT(Complex *a,int f){ for(int i=0;i<limit;i++) if(i<r[i]) swap(a[i],a[r[i]]); for(int mid=1;mid<limit;mid<<=1){ Complex Wn(cos(PI/mid),sin(PI/mid)*f); for(int j=0;j<limit;j+=mid<<1){ Complex w(1,0); for(int k=0;k<mid;k++,w=w*Wn){ Complex x=a[j+k],y=w*a[j+k+mid]; a[j+k]=x+y,a[j+k+mid]=x-y; } } } if(f==-1) for(int i=0;i<limit;i++) a[i]=a[i]/limit; } int main(){ scanf("%d%d",&n,&m); for(int i=0;i<=n;i++) scanf("%lf",&a[i].real); for(int i=0;i<=m;i++) scanf("%lf",&b[i].real); for(limit=1;limit<n+m+1;limit<<=1); for(int i=0;i<limit;i++) r[i]=(r[i>>1]>>1)|((i&1)?(limit>>1):0); FFT(a,1),FFT(b,1); for(int i=0;i<=limit;i++) a[i]=a[i]*b[i]; FFT(a,-1); for(int i=0;i<=n+m;i++) printf("%d ",(int)(a[i].real+0.5)); return 0; } ``` 其中,结构体 `Complex` 表示复数,包含实部和虚部。`FFT` 函数实现了快速傅里叶变换,其中 `f` 表示进行正变换还是逆变换,正变换时 `f=1`,逆变换时 `f=-1`。主函数中输入了两个多项式 `a` 和 `b`,并对它们进行了乘法运算,最后输出结果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值