Problem #46 says
It was proposed by Christian Goldbach that every odd composite number can be written as the sum of a prime and twice a square.
9 = 7 + 2×1^(2)
15 = 7 + 2×2^(2)
21 = 3 + 2×3^(2)
25 = 7 + 2×3^(2)
27 = 19 + 2×2^(2)
33 = 31 + 2×1^(2)
It turns out that the conjecture was false.
What is the smallest odd composite that cannot be written as the sum of a prime and twice a square?
Solution below uses the method of generating a list of primes and then using a set list of numbers to square and checking the combination’s to see if we can find one that makes the number we’re looking for.
Solution provided in f# and requires .net 4. runs in ~50 seconds. could be made much faster by rewriting the “Check” function to use a set of recursive functions instead of for loops so we could break out of them.
#light
open System
open System.Diagnostics
open System.Collections.Generic
//get some prime factors
//set the limit and bound
let UpperLimit = 10000.
let UpperBound = 10000
//generate numbers to check
let LNums = ref [2..(int UpperLimit)] in
// loop through all numbers
// keep removing numbers that divide evenly
for div = 2 to UpperBound do
LNums := List.filter(fun num -> (num % div <> 0 || div >= num))
!LNums
done;
let Primes:list<int> = !LNums
//numbers to check
let nums = [13 .. 6000]
let NonPrime =
(Set.ofList nums) - (Set.ofList Primes) |> Set.toList |> List.filter (fun x -> x % 2 <> 0)
Console.WriteLine("Starting check")
let watch = new Stopwatch()
watch.Start()
let Check number =
let mutable RetValue = true
for n in Primes do
for c in [1..100] do
if (n + (2 * (c * c))) = number then
RetValue <- false
RetValue
let N =
NonPrime
|> Seq.filter Check
|> Seq.toArray
watch.Stop()
Console.WriteLine(N.[0])
Console.WriteLine(watch.Elapsed)
Console.ReadKey()
Tags: .net4, f#, programming, project euler
Posted in Project Euler | Comments (0)
Problem #145 says:
Some positive integers n have the property that the sum [ n + reverse(n) ] consists entirely of odd (decimal) digits. For instance, 36 + 63 = 99 and 409 + 904 = 1313. We will call such numbers reversible; so 36, 63, 409, and 904 are reversible. Leading zeroes are not allowed in either n or reverse(n).
There are 120 reversible numbers below one-thousand.
How many reversible numbers are there below one-billion (10^(9))?
The solution below is c# with .Net 4.0. It does a very basic divide and conquer on the one billion numbers to check. This will take a while to run depending on your computer, it ran reasonably fast on a quad core machine. Pretty simple stuff could be made much faster but was written out of boredom.
using System;
using System.Collections;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Linq.Parallel;
using System.Threading;
using System.Threading.Tasks;
using System.Diagnostics;
namespace Net4Test
{
class Program
{
static object LockHandle = new object();
static long Answer = 0;
static void Main(string[] args)
{
Stopwatch sw = new Stopwatch();
sw.Start();
Parallel.For(1L, 1000000000L, i =>
{
if (i % 50000 == 0)
Console.WriteLine(i);
if (Reversable(i))
{
lock (LockHandle)
{
Answer += 1;
}
}
});
Console.WriteLine(Answer);
sw.Stop();
Console.WriteLine(sw.Elapsed);
Console.ReadLine();
}
static bool Reversable(long l)
{
bool RetValue = true;
char[] lchar = l.ToString().ToCharArray();
Array.Reverse(lchar);
if (lchar[0] != '0')
{
long Reversel = (long)Convert.ToDouble(new string(lchar));
char[] AChar = (l + Reversel).ToString().ToCharArray();
foreach (char c in AChar)
{
if (!IsOdd(Convert.ToInt32(c.ToString())))
RetValue = false;
}
}
else
RetValue = false;
return RetValue;
}
static bool IsOdd(int i)
{
if (i == 0)
return false;
if (i % 2 == 0)
return false;
else
return true;
}
}
}
Tags: c#, Math, programming, project euler
Posted in Project Euler | Comments (0)
(updates on the site. pagerank just hit 3 on google. hopefully this will increase traffic a bit. and this post is the first post on the site done in windows 7 (took 6 hours for my pc to update
)
Problem #50 says:
The prime 41, can be written as the sum of six consecutive primes:
41 = 2 + 3 + 5 + 7 + 11 + 13
This is the longest sum of consecutive primes that adds to a prime below one-hundred.
The longest sum of consecutive primes below one-thousand that adds to a prime, contains 21 terms, and is equal to 953.
Which prime, below one-million, can be written as the sum of the most consecutive primes?
The solution below is pretty much brute force. We generate the primes and then add them together and subtract from the total items in the list of primes. then we get the total of all of the primes minus the first item in the primes list and do it over again. there’s lots of room for improvements and when I get the time I’ll post them. but for now this is what I used to get the solution. Solution provided in c#/mono.
using System;
using System.Collections;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Diagnostics;
namespace ProjectEuler
{
class Program
{
static void Main(string[] args)
{
Stopwatch sw = new Stopwatch();
sw.Start();
ArrayList primes = GeneratePrimes(1000000);
primes.Remove(1);
long result = 0;
long maxlength = 0;
Console.WriteLine("starting check");
for (int start = 0;start<primes.Count;start++)
{
long temp = 0;
for (int i = start; i < primes.Count; i++)
{
temp = temp + (long)primes[i];
}
for (int k = primes.Count - 1; k > start; k--)
{
if (temp < 1000000 && primes.Contains(temp) && k - start >= maxlength)
{
maxlength = k -start;
result = temp;
Console.WriteLine(result);
break;
}
temp = temp - (long)primes[k];
if (temp <= start)
break;
}
}
Console.WriteLine(result);
Console.WriteLine(maxlength);
Console.WriteLine(sw.Elapsed);
Console.ReadLine();
}
static ArrayList GeneratePrimes(long num)
{
ArrayList RetValue = new ArrayList();
RetValue.Add((long)2);
RetValue.Add((long)3);
long Stepper = 5;
long Check = 1;
while (Stepper <= num)
{
foreach (long i in RetValue)
{
if (Stepper % i == 0)
{
Check = 0;
break;
}
if (Math.Sqrt(Stepper) < i)
{
break;
}
}
if (Check == 1)
{
//Console.WriteLine(((float)Stepper / (float)num) * 100);
RetValue.Add(Stepper);
}
Check = 1;
Stepper++;
Stepper++;
}
return RetValue;
}
}
}
Tags: c#, csharp, mono, programming, project euler
Posted in Project Euler | Comments (0)