在GAP中构造满足特定条件的三元组集合PS的方法及语法问题咨询

首先帮你明确PS的数学定义：
根据你给出的几条条件，PS是下闭的非负整数三元组集合（简单说就是：只要某个三元组在PS里，任何坐标不超过它的非负三元组也都在PS里），同时要满足：

• 边界点$(A-1,0,0)$、$(0,V-1,0)$、$(0,0,R-1)$都属于PS；
• 边界点$(A,0,0)$、$(0,V,0)$、$(0,0,R)$都不属于PS。

其实这个集合可以直接简化表述：所有满足 $p < A$、$q < V$、$r < R$ 的非负整数三元组 $(p,q,r)$。这个集合天然符合下闭性，也完美匹配你给出的所有边界条件。

接下来解决你在GAP代码里遇到的具体问题：

1. 函数拼写错误：SubstractSet → SubtractSet #

你写的SubstractSet是拼写失误，GAP里的正确函数名是SubtractSet（注意是Subtract而非Substract）。这个函数会直接修改第一个集合，把指定元素移除：

PS := [];
for p in [0..10] do
  for q in [0..10] do
    for r in [0..10] do
      AddSet(PS, [p,q,r]);
    od;
  od;
od;
// 移除指定的三个三元组
SubtractSet(PS, [[10,0,0], [0,10,0], [0,0,10]]);
// 查看最终集合
PS;

如果你不想修改原集合，而是生成一个新的集合，可以用Difference函数：

newPS := Difference(PS, [[10,0,0], [0,10,0], [0,0,10]]);

2. if语句的语法问题 #

你碰到的Syntax error: then expected，是因为GAP要求if和then处于同一个逻辑行（或者用分号分隔开）。不过其实完全没必要用if语句来判断，直接构造符合条件的集合会更简洁高效。比如假设这里A=10、V=10、R=10，直接生成所有满足$p < 10$、$q < 10$、$r < 10$的三元组：

PS := [];
A := 10; V := 10; R := 10;
for p in [0..A-1] do
  for q in [0..V-1] do
    for r in [0..R-1] do
      AddSet(PS, [p,q,r]);
    od;
  od;
od;
// 验证条件：检查(A,0,0)不在PS中
[10,0,0] in PS; // 会返回false
// 检查(A-1,0,0)在PS中
[9,0,0] in PS; // 会返回true

3. 更简洁的构造方式 #

你还可以用GAP的Cartesian函数直接生成笛卡尔积，一步到位得到目标集合：

A := 10; V := 10; R := 10;
// 生成所有0到A-1、0到V-1、0到R-1的三元组笛卡尔积
PS := Cartesian([0..A-1], [0..V-1], [0..R-1]);
// 转换为集合类型（如果需要的话）
PS := AsSet(PS);

这种方式不需要后续移除元素，直接生成满足所有条件的集合，效率更高。

备注：内容来源于stack exchange，提问作者Hana Neziri