x=平方[2];z=800;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]]&];