我需要在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中的向量?
答案 0 :(得分:1)
事实证明,有一个内置函数:Vecrev()和Polrev()与Vec和Pol做同样的事情,但相反。因此,您可以使用以下内容反转任意向量:
Vecrev(vector(10^6,i,i))