prime-factoring - 在函数 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 次幂，这里没有错误。