pari-gp - 如何在 PARI/gp 中反转向量

Question

我需要反转 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 中的向量？

Answer 1

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

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