>>28に追加
P(x)=x-1-(Li(x^2+x-1)-Li(x^2-1))?

ここで、Li(x)は(補正)対数積分で
Li(x)=∫[t=2,x]1/ln(t)dt