（Python）
从数学导入isqrt
来自sympy import mobius
定义A062503型（n） ：
定义f（x）：对于范围（1，isqrt（x）+1）中的k，返回n-1+x-sum（mobius（k）*（x//k**2））
kmin，kmax=1,2
而f（kmax）>=kmax：
kmax≤1
为True时：
kmid=kmax+kmin>>1
如果f（kmid）<kmid：
kmax=kmid
其他：
kmin=kmid
如果kmax-kmin<=1：
打破
返回kmax**2#柴华武2024年8月19日
