是否有一種舒適的方法可以在PARI/GP中生成n個組合中的k個組合？

Question

假設我在PARI/GP中有n個元素的向量。

我想在n個元素中生成k個組合。

例如，如果載體是[3,7,11,14,18]，且k = 3時，輸出應該是

[3,7,11] 
[3,7,14] 
[3,7,18] 
[3,11,14] 
[3,11,18] 
[3,14,18] 
[7,11,14] 
[7,11,18] 
[7,14,18] 
[11,14,18] 

是否有PARI/GP這樣做的命令，或我必須編寫函數嗎？

Answer 1

不幸的是，PARI沒有一個內置的命令來完成所有你需要的東西。

功能forvec（標誌= 2）適合非常好的組合問題。所以你的功能可以如下：

subsets(A, k) = { 
    my (lst = List()); 
    forvec(v = vector(k, i, [1, #A]), listput(lst, vecextract(A, v)), 2); 

    Vec(lst) 
}; 

subsets([3,7,11,14,18], 3) 
gp> [[3,7,11], [3,7,14], [3,7,18], [3,11,14], [3,11,18], [3,14,18], [7,11,14], [7,11,18], 
    [7,14,18], [11,14,18]]