如何在PARI / gp中反转矢量

时间:2018-05-29 11:20:54

标签: pari-gp

我需要在PARI / gp中反转一个向量。我找不到内置函数,所以我尝试了这个:

vector(10^4,i,vector(10^4,i,i)[10^4-i+1])

这很慢 - 花了将近四秒钟。所以我试过这个:

Vec(polrecip(Pol(vector(10^6,i,i))))

这更快(大约100毫秒),即使它正在反转100倍长的矢量。不过,这是一个可怕的黑客。是否有一种“正确”的方法可以快速反转PARI / gp中的向量?

1 个答案:

答案 0 :(得分:1)

事实证明,有一个内置函数:Vecrev()和Polrev()与Vec和Pol做同样的事情,但相反。因此,您可以使用以下内容反转任意向量:

 Vecrev(vector(10^6,i,i))