views:

310

answers:

4

Hi, I'm trying to find the sum of the prime numbers < 2,000,000. This is my solution in java but I can't seem get the correct answer. Please give some input on what could be wrong and general advice on the code is appreciated.

Printing 'sum' gives: 1308111344, which is incorrect.

edit: Thanks for all the help. Changed int to long and < to <= and it worked flawlessly, except for being an inefficient way of finding prime numbers :)

/*
The sum of the primes below 10 is 2 + 3 + 5 + 7 = 17.
Find the sum of all the primes below two million.
*/
class Helper{
 public void run(){
  Integer sum = 0;
  for(int i = 2; i < 2000000; i++){
   if(isPrime(i))
    sum += i;   
  }
  System.out.println(sum);
 }

 private boolean isPrime(int nr){
  if(nr == 2)
   return true;
  else if(nr == 1)
   return false;
  if(nr % 2 == 0)
   return false;

  for(int i = 3; i < Math.sqrt(nr); i += 2){
   if(nr % i == 0)
    return false;
  }  
  return true;  
 }
}   
class Problem{
 public static void main(String[] args){
  Helper p = new Helper();
p.run();
}
}
+7  A: 

The result will be too large to fit into an integer, so you are getting an overflow. Try using a BigInteger or a long instead. In this case a long is enough.

Mark Byers
That's not the problem, or at least not the only problem.
Pointy
A long would be better [faster] than a BigInteger. It's obviously large enough since 1 + 2 + 3 + ... + 2000000 = 2000001000000 which is smaller than the max of a long.
Charles
@Pointy: Of the two mistakes this one is the one that gives the largest magnitude error. But both mistakes cause the answer to be wrong. Unfortunately there is no way to accept two answers so if you want you can copy my answer into yours and accept the glory. ;)
Mark Byers
Yes you're right; I wasn't sure whether it'd squeeze in but another 1/2 cup of coffee and I did the basic math. When I did those Project Euler problems I used Erlang, so I never had to worry about integer sizes (all Erlang integers are like "BigInteger") :-)
Pointy
+6  A: 

You're treating as prime those numbers that are only divisible by their square root (like 25). Instead of i < Math.sqrt(nr) try i <= Math.sqrt(nr).

That's a really inefficient way to find primes, incidentally.

Pointy
Whoa why the drive-by downvote? The algorithm is in fact wrong, for the reason I state here.
Pointy
Wrong. He's looping from 3 to the square root of the number.If the number is, let's say 25 - you don't need to check anything higher than 5. If there is divisor higher than 5 - then the other divisor will be lower than five, and would've been checked in the beginning of the loop.
Max
He's **not** looping to the square root!! The comparison he uses is `<` not `<=`.
Pointy
Good lord, people, **read his code**.
Pointy
Ah, you are right, I slightly misunderstood your response. Sorry, inverted my vote :D
Max
+3  A: 

Your isPrime doesn't work for squares. isPrime(9) returns true.

Pete Kirkham
+1  A: 

As already stated errors were two:

  • you used an int that is not big enough to hold that sum.. you should have used a long
  • you used < instead that <=, and it was a wrong guard for the cycle

Apart from that what you are doing is really inefficient, without going too deep inside this class of algorithms (like Miller-Rabin test) I would suggest you to take a look to the Sieve of Eratosthenes.. a really old approach that teaches how to treat a complex problem in a simple manner to improve elegance and efficiency with a trade-off of memory.

It's really cleaver: it keeps track of a boolean value for every prime up to your 2 millions that asserts if that number is prime or not. Then starting from the first prime it excludes all the successive numbers that are obtained by multiplying the prime it is analyzing for another number. Of couse more it goes and less numbers it will have to check (since it already excluded them)

Code is fair simple (just wrote it on the fly, didn't check it):

    boolean[] numbers = new boolean[2000000];
    long sum = 0;

    for (int i = 0; i < numbers.length; ++i)
        numbers[i] = true;

    for (int i = 2; i < numbers.length; ++i)
        if (!numbers[i])
            continue;
        else {
            int j = i + i;
            while (j < 2000000) {                   
                numbers[j] = false;
                j += i;
            }           
        }

    for (int i = 2; i < 2000000; ++i)
        sum += numbers[i] ? i : 0;

    System.out.println(sum);

Of course this approach is still unsuitable for high numbers (because it has to find all the previous primes anyway and because of memory) but it's a good example for starters to think about problems..

Jack
http://stackoverflow.com/questions/1042902/most-elegant-way-to-generate-prime-numbers/1043247#1043247
starblue
using a bitset is a really nice idea, I would like to have the good old JIT java compiler to do these optimization by itself :)
Jack
Also note that it starts eliminating multiples at `i * i`, all multiples below that must have been eliminated by a smaller factor. This is a big gain.
starblue