f[n_]:=加号@@选择[Divisors[n],Sqrt[n/2]<=#<Sqrt[2]&];选择[Range[175],f[#]!=0 &]
(*与西格玛的对称表示有关*)
(*m<=k<=n*的数字k的奇数部分的子序列)
(*使用Wilson的Mathematica程序(见上文),我验证了数字k<=10000时两者的相等性*)
a071562[m_,n_]:=选择[Range[m,n],OddQ[Length[a237270[#]]&]
a071562[1,160](*数据*)
(*使用奇数除数属性实现*)
evenExp[n_]:=第一个[NestWhile[{#[[1]]+1,#[2]]/2}&,{0,n},EvenQ[Last[#]]&]]
oddSRQ[n_]:=模块[{e=2^evenExp[n],楼层[(Sqrt[8n+1]-1)/2]},选择[Divisors[n/e],#<=r&&2 e#>r&]={}]
a071562D[m_,n_]:=选择[Range[m,n],oddSRQ]
|