[Pie]

Buenas. Estoy intentando calcular la constante de Brun (la suma de los inversos de los primos gemelos) en Pascal. No puedo calcular muy allá ya que uso un compilador online (con límite de tiempo), pero aún así el resultado que me da me parece que queda bastante lejos de lo esperado.

Para primos < 10 millones me da:

1.738357043917..

Es normal y converge de forma taaan lenta o puede deberse a algún error de programación?

He sumado el 5 (1/5) dos veces como en la wiki (tenía dudas sobre esto XD), y analizado sólo los números de la forma \( 6n \pm{1} \) para ganar algo de tiempo (o no perderlo más bien) pero de 10 millones no paso (y me hace pensar que incluso con el ordenador me costaría llegar a sacar las primeras cifras XD)

Salu2

[Luis Fuentes]

Hola

 Si, la convergencia es muy lenta. Mira este enlace, en especial su sección 3.2:

http://numbers.computation.free.fr/Constants/Primes/twin.html

Saludos.

[Pie]

Muchas gracias Luís. No sabía que usaban esa extrapolación para calcularla (ya me extrañaba, podría tirarme la vida entera calculándola de forma directa XD)

Saludos.

[feriva]

Buenas. Estoy intentando calcular la constante de Brun (la suma de los inversos de los primos gemelos) en Pascal. No puedo calcular muy allá ya que uso un compilador online (con límite de tiempo), pero aún así el resultado que me da me parece que queda bastante lejos de lo esperado.

Para primos < 10 millones me da:

1.738357043917..

A mí con Python, hasta el millón, me da 1.710776930804... (y está bien hecho, lo he ido comprobando).

Saludos.

[Pie]

Buenas. Estoy intentando calcular la constante de Brun (la suma de los inversos de los primos gemelos) en Pascal. No puedo calcular muy allá ya que uso un compilador online (con límite de tiempo), pero aún así el resultado que me da me parece que queda bastante lejos de lo esperado.

Para primos < 10 millones me da:

1.738357043917..

A mí con Python, hasta el millón, me da 1.710776930804... (y está bien hecho, lo he ido comprobando).

Saludos.

Sí, a mí también me da eso hasta el millón. También da lo mismo en el enlace que puso Luís así que no era fallo de programación.

Saludos.

[Pie]

Por cierto, aunque ya no tenga que ver con el cálculo de esta constante (imposible de esta forma), estaba pensando que también es una pérdida de tiempo analizar los números de la forma \( 6n \pm{1} \) acabados en 4 o 6, ya que siempre habrá algún múltiplo de 5 por ahí.

Supongo que con alguna condición tipo si 6n mod 10 = 4 o 6 entonces nada, se podría solventar. Lo digo por si a alguien le interesa calcular los primos gemelos así sin más y ganar algo más de tiempo. Aunque imagino que esto es bastante obvio y conocido jeje

PD: Si alguien conoce o se le ocurre algún truquito más que no se corte.

Saludos.