En primer lugar, quería señalar que la mayor parte de lo que has escrito no es correcto. t(x)∼xlogxt(x)∼xlogx no ∼x . Esto se debe a que si ∑p≤x∼li(x) entonces esperamos que al dividir por logp , voy a conseguir li(x)/log(x) y luego t(x)∼li(x) . Además, los cálculos están muy equivocados. Comprobando en Matlab, lo tengo: r(900000)=9000090.94… y por lo tanto r(900000)−900000=90.94… No has eliminado los primos correctamente. Ahora vamos a probar todo esto en su totalidad.
Algunas pruebas rigurosas: Vamos a repasar todo en detalle. Las dos series que está considerando son
t(x)=logx∑p≤x1logp y r(x)=logx∑1<c≤x1logc, where c>1 is composite. Observe que r(x)+t(x)=∑2≤n≤x1logn por lo que sólo tenemos que evaluar ∑2≤n≤x1logn y ∑p≤x1logp . Primero
Lema 1: ∑2≤n≤x1logn=li(x)−C+O(1logx) donde la constante C viene dada por C=∫∞2{t}t(logt)2dt.
Prueba: Podemos escribirlo como una integral de Riemann-Stieltjes: ∑2≤n≤x1logn=∫x21logtd[t]=∫x21logtdt−∫x21logtd{t} donde [t] y {t} son las partes mínima y fraccionaria de t respectivamente. Utilizando la integración por partes y la definición de li(x) esto es li(x)−{x}logx−∫x2{t}(logt)2dt.
Observación: Aunque no he hecho el cálculo yo mismo, creo que la constante C puede limpiarse en términos de otras constantes conocidas.
Lema 2: Tenemos que ∑p≤x1logp=li(x)−xlogx+O(xe−c√logx).
Prueba: Si θ(x)=∑p≤xlogp podemos escribir ∑p≤x1logp=∫x21(logt)2dθ(t)=∫x21(logt)2dt+∫x21(logt)2d(θ(t)−t). Para el primer término, observe que por integración por partes ∫x21logtdt=xlogx+∫x21(logt)2dt. Para el segundo término podemos utilizar la integración por partes junto con el teorema de los números primos que establece que θ(t)−t=O(xe−c√logx). Entonces ∫x21(logt)2d(θ(t)−t)=(θ(x)−x)(logx)2+2∫x2θ(t)−tt(logt)3dt+O(1) =O(xe−c√logx), y el lema sigue.
Consecuencias: Juntando estos dos lemas, encontramos que: t(x)=li(x)logx−x+O(xe−c√logx), r(x)+t(x)=li(x)logx−Clogx+O(xe−c√logx). Tenemos que poner el término de error extra ya que estamos restando r(x) y esto consumirá el Clogx ya que es mayor. Por lo tanto, r(x)=x+O(xe−c√logx) y r(x)−x=O(xe−c√logx).
Espero que eso ayude,