在FACTOR函数中的PARI / GP中，-1是一个错误或表示其他东西？

Question

首先运行函数b（n）：

? b(n) = lcm(vector(n, i, i))/n

函数c（n）之后：

? c(n)=sum(j=1,n,sum(i=1,n,(-1)^(i+j)/(i+j-1)))

上次运行d（n）：

? d(n)=factor(denominator(c(n))/b(n))~

并测试202

? d(202)

结果是：

%8 =
[3 7 17 19 31 211 223 227 229 233 239 241 251 257 263 269 271 277 281 283 293 
307 311 313 317 331 337 347 349 353 359 367 373 379 383 389 397 401]

[1 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]

什么表示分解结果中的-1？

Answer 1

您尝试考虑有理数。注意，由于type(denominator(c(202))/b(202))，t_FRAC是t_INT而不是denominator(c(202))/b(202) = <some big number>/31。因此-1意味着除数31的功率为-1，这里没有错误。