//生成OEIS A367566 b文件的Magma程序//Jon E.Schoenfield，2023年11月24日[1..10000]do中的n如果IsEven（n），则触头：=1；其他的触头：=2；结束条件：；如果（n gt 1）和（n mod 3 ne 0），则触头*：=3；结束条件：；P： =[]；对于[-1,1]do中的jt： =n+j；如果t大于0，则F： =因子分解（t）；对于f do中的fp： =f[1]；如果p gt 3，则B： =p*（p-1）div 2；如果[p-1，B-（p-1）]中的n mod B，则触头*：=p；结束条件：；结束条件：；结束；结束条件：；结束；n、 触头；结束；