10scisci~5201411510111Python#Truepri=[Truefor_inrange（lit+1）]p=2while（p*p#pri[p]if（pri[p]==True）:#p10000000lit=1000000010502230pris=np.array（[235711131719232931#gaps=np.diff（pris）#$x$pi（x）（x）[pi（x）approxfrac{x}{lnx}]$x$prirangePythont..（）