Meissel–Mertens constant: Difference between revisions

m
Minor improvement to code.
(New post.)
m (Minor improvement to code.)
Line 540:
import java.util.List;
 
/**
* Calculates the Meissel-Mertens constant correct to 9 s.f. in approximately 15 seconds.
*/
public final class MeisselMertensConstant {
Line 550 ⟶ 547:
double sum = 0.0;
for ( double reciprocal : primeReciprocals ) {
sum += reciprocal + Math.log(1.0 - reciprocal);
}
final double constantmeisselMertens = euler + sum;
System.out.println(String.format("%s%.9f", "The Meissel-Mertens constant = ", + constantmeisselMertens));
}
Line 561 ⟶ 558:
sieve.set(2, aLimit + 1);
finalfor ( int squareRooti = (int)2; i <= Math.sqrt(aLimit); i = sieve.nextSetBit(i + 1) ) {
for ( int i = 2; i <= squareRoot; i = sieve.nextSetBit(i + 1) ) {
for ( int j = i * i; j <= aLimit; j += i ) {
sieve.clear(j);
Line 577 ⟶ 573:
 
}
</syntaxhighlight>
{{ out }}
<pre>
The Meissel-Mertens constant = 0.26149721287104916261497213
</pre>
 
908

edits