x=黄金比率；z=1000；p[k_]：=p[k]=素数[k]；
t=表[Max[表[NextPrime[x*p[k]，-1]/p[k]、{k、1、n}]]、{n、1、z}]；
d=删除重复项[t]；tL=选择[d，#>0&]（*下部POBA*）
t=表[Min[Table[NextPrime[x*p[k]]/p[k]，{k，1，n}]]，{n，1，z}]；
d=删除重复项[t]；tU=选择[d，#>0&]（*上部POBA*）
v=排序[Union[tL，tU]，Abs[#1-x]>Abs[#2-x]&]；
b=分母[v]；s=选择[Range[Length[b]]，b[[#]]==最小值[Drop[b，#-1]]&]；
