fuxino/project-euler-solutions
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Correct problem 9 solution

Browse Source
This commit is contained in:
daniele 2019-09-29 18:03:11 +02:00
parent 3a09f3d4f5
commit 2cdbca922a
Signed by: fuxino
GPG Key ID: 6FE25B4A3EE16FDA
2 changed files with 64 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

42
C/p009.c
View File

@ -11,10 +11,11 @@
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include "projecteuler.h"
int main(int argc, char **argv)
{
int a, b, c, m, n, found = 0;
int a, b, c, tmpa, tmpb, tmpc, i, m, n, found = 0;
double elapsed;
struct timespec start, end;
@ -24,16 +25,39 @@ int main(int argc, char **argv)
* Euclid's formula, until the one where a+b+c=1000 is found.*/
for(m = 2; found == 0; m++)
{
for(n = 1; n < m; n++)
for(n = 1; n < m && !found; n++)
{
a = m * m - n * n;
b = 2 * m * n;
c = m * m + n * n;
if(a + b + c == 1000)
if(gcd(m, n) == 1 && (m % 2 == 0 && n % 2 != 0 || m % 2 != 0 && n % 2 == 0))
{
found = 1;
break;
a = m * m - n * n;
b = 2 * m * n;
c = m * m + n * n;
if(a + b + c == 1000)
{
found = 1;
break;
}
i = 2;
do
{
tmpa = a * i;
tmpb = b * i;
tmpc = c * i;
if(tmpa + tmpb + tmpc == 1000)
{
a = tmpa;
b = tmpb;
c = tmpc;
found = 1;
break;
}
i++;
}while(tmpa + tmpb + tmpc < 1000);
}
}
}

37
Python/p009.py
View File

@ -10,6 +10,8 @@
#
# Find the product abc.
from math import gcd
from timeit import default_timer
def main():
@ -23,14 +25,37 @@ def main():
# Euclid's formula, until the one where a+b+c=1000 is found.
while not found:
for n in range(1, m):
a = m * m - n * n
b = 2 * m * n
c = m * m + n * n
if a + b + c == 1000:
found = 1
if found == 1:
break
if gcd(m, n) == 1 and (m % 2 == 0 and n % 2 != 0 or m % 2 != 0 and n % 2 == 0):
a = m * m - n * n
b = 2 * m * n
c = m * m + n * n
if a + b + c == 1000:
found = 1
break
i = 2
while True:
tmpa = a * i
tmpb = b * i
tmpc = c * i
if tmpa + tmpb + tmpc == 1000:
a = tmpa
b = tmpb
c = tmpc
found = 1
break
if tmpa + tmpb + tmpc > 1000:
break
i = i + 1
m = m + 1
end = default_timer()