? f(n) = s=0; for(x=1,n,if(isprime(x^2+1),s++)); s ? f(100) %1 = 19 ? f(1000) %2 = 112 ? f(10000) %3 = 841 ? f(100000) %4 = 6656
4 | 2 | 2 |
6 | 3 | 3 |
8 | 3 | 5 |
10 | 3 | 7 |
12 | 5 | 7 |
? gb(n) = forprime(p=2,n,if(isprime(n-p),return([p,n-p]))); ? gb(4) %7 = [2, 2] ? gb(6) %8 = [3, 3] ? gb(100) %9 = [3, 97] ? gb(1000) %10 = [3, 997] ? gb(570) \\ takes no time at all! %11 = [7, 563]
? t(n) = s=0; forprime(p=2,n,if(isprime(p+2),s++)); s ? t(10^2) %12 = 8 ? t(10^3) %13 = 35 ? t(10^4) %14 = 205 ? t(10^5) %15 = 1224Surely keeps getting bigger!!
As it turns out, these three assertions are all OLD famous extremely difficult unsolved problems! Anyone who proves one of them will be very famous.
Assertion 2 is called ``The Goldbach Conjecture''; Goldbach reformulated it in a letter to Euler in 1742. It's featured in the following recent novel:
The publisher of that novel offers a MILLION dollar prize for the solution to the Goldbach conjecture:
http://www.faber.co.uk/faber/million_dollar.asp?PGE=&ORD=faber&TAG=&CID=The Goldbach conjecture is true for all , see
http://www.informatik.uni-giessen.de/staff/richstein/ca/Goldbach.html
Assertion 3 is the ``Twin Primes Conjecture''. According to
http://perso.wanadoo.fr/yves.gallot/primes/chrrcds.html#twinon May 17, 2001, David Underbakke and Phil Carmody discovered a 32220 digits twin primes record with a set of different programs: . This is the current ``world record''.
With a computer, even if you can't solve one of these ``Grand Challenge'' problems, at least you can perhaps work very hard and prove it for more cases than anybody before you, especially since computers keep getting more powerful. This can be very fun, especially as you search for a more efficient algorithm to extend the computations.