mborm.net

Problem 74 auf projecteuler.netNach längerer Zeit endlich mal wieder eine neue Lösung, dieses Mal für das Problem 74
Auf meinem alten Laptop mit einem i3-2350M braucht das Skript nur knapp über eine Sekunde:


gcc -O3 -march=native -mtune=native -std=c11 -ffast-math Problem74.c -lm -o Problem74 && time ./Problem74

Result: 402

real	0m1,196s

user	0m1,196s

sys	0m0,000s

Und der Code dazu:
Problem 74 in C
C
#include <stdio.h>
#include <math.h>

typedef unsigned long ulong;

static const ulong faclist[] = {1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880};

static inline ulong getDigitFacSum(ulong n) {
    ulong sum;

    for(sum = 0; n > 0; sum += faclist[n % 10], n /= 10);

    return sum;
}

static ulong getChainLength(ulong n) {
    ulong checklist[60], count = 1, i;

    if((*checklist = getDigitFacSum(n)) == n)
        return count;

    do {
        checklist[count] = getDigitFacSum(checklist[count-1]);

        for(i = 0; i < count; i++)
            if(checklist[count] == checklist[i])
                return ++count;
    } while(++count < 60);

    return 0;
}

int main(void) {
    ulong c = 0, i;

    for(i = 1; i < 1000000; i++)
        if(getChainLength(i) == 60)
            c++;

    printf("Result: %lu\n", c);
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

						#include <stdio.h>
#include <math.h>
 
typedef unsigned long ulong;
 
static const ulong faclist[] = {1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880};
 
static inline ulong getDigitFacSum(ulong n) {
    ulong sum;
 
    for(sum = 0; n > 0; sum += faclist[n % 10], n /= 10);
 
    return sum;
}
 
static ulong getChainLength(ulong n) {
    ulong checklist[60], count = 1, i;
 
    if((*checklist = getDigitFacSum(n)) == n)
        return count;
 
    do {
        checklist[count] = getDigitFacSum(checklist[count-1]);
 
        for(i = 0; i < count; i++)
            if(checklist[count] == checklist[i])
                return ++count;
    } while(++count < 60);
 
    return 0;
}
 
int main(void) {
    ulong c = 0, i;
 
    for(i = 1; i < 1000000; i++)
        if(getChainLength(i) == 60)
            c++;
 
    printf("Result: %lu\n", c);
}

					

			
 
Veröffentlicht am 19. August 2017 von admin in C, Programmierung, projecteuler
Problem 50 auf projecteuler.netHier meine Brute Force Lösung zu Problem 50.
Befehl zum Compilieren und Ausführen sowie die Ausgabe mit Zeitmessung:


gcc -O3 -march=native -mtune=native -std=c11 -ffast-math -fopenmp Prob50.c -lm -o Prob50 && time ./Prob50

Result: 997651 with 543 terms

real	0m2.809s

user	0m16.577s

sys	0m0.012s

Und der Code:
C
#include <stdio.h>

#define ARRLENGTH(x) (sizeof(x) / sizeof((x)[0]))

static unsigned int prmlst[78498];

inline static unsigned int isPrime(unsigned int n) {
    if(n == 1)
        return 0;

    else if(n == 2)
        return 1;

    else {
        int i, r = (int) __builtin_sqrt(n);

        for(i = 2; i <= r && n % i != 0; ++i);

        return (r == i - 1) ? 1 : 0;
    }
}

inline static unsigned int getConPrmCnt(unsigned int n) {
    unsigned int sum = 2, prm = 1, cnt = 1, sumcnt = 0;

    while(sum != n) {
        sum = prmlst[sumcnt++];
        prm = 1;
        cnt = sumcnt;

        while(sum < n) {
            sum += prmlst[cnt++];
            prm++;
        }
    }

    return (sum == n) ? prm : 0;
}

int main(void) {
    unsigned int idx = 0, res = 0;

    for(unsigned int i = 0; i < 1000000; i++)
        if(isPrime(i))
            prmlst[idx++] = i;

    #pragma omp parallel for
    for(unsigned int i = 0; i < ARRLENGTH(prmlst); i++) {
        unsigned int tmp = getConPrmCnt(prmlst[i]);

        if(tmp > res) {
            res = tmp;
            idx = prmlst[i];
        }
    }

    printf("Result: %u with %u terms\n", idx, res);
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58

						#include <stdio.h>
 
#define ARRLENGTH(x) (sizeof(x) / sizeof((x)[0]))
 
static unsigned int prmlst[78498];
 
inline static unsigned int isPrime(unsigned int n) {
    if(n == 1)
        return 0;
 
    else if(n == 2)
        return 1;
 
    else {
        int i, r = (int) __builtin_sqrt(n);
 
        for(i = 2; i <= r && n % i != 0; ++i);
 
        return (r == i - 1) ? 1 : 0;
    }
}
 
inline static unsigned int getConPrmCnt(unsigned int n) {
    unsigned int sum = 2, prm = 1, cnt = 1, sumcnt = 0;
 
    while(sum != n) {
        sum = prmlst[sumcnt++];
        prm = 1;
        cnt = sumcnt;
 
        while(sum < n) {
            sum += prmlst[cnt++];
            prm++;
        }
    }
 
    return (sum == n) ? prm : 0;
}
 
int main(void) {
    unsigned int idx = 0, res = 0;
 
    for(unsigned int i = 0; i < 1000000; i++)
        if(isPrime(i))
            prmlst[idx++] = i;
 
    #pragma omp parallel for
    for(unsigned int i = 0; i < ARRLENGTH(prmlst); i++) {
        unsigned int tmp = getConPrmCnt(prmlst[i]);
 
        if(tmp > res) {
            res = tmp;
            idx = prmlst[i];
        }
    }
 
    printf("Result: %u with %u terms\n", idx, res);
}

					

			
 
Veröffentlicht am 21. Mai 2016 von admin in C, Programmierung, projecteuler
Problem 35 auf projecteuler.netHier meine Lösung zu Problem 35. Für alle Primzahlen unter 1 Mio werden die jeweiligen Permutationen untersucht. Sind diese auch alles Primzahlen handelt es sich um eine kreisförmige Primzahl (Circular Prime Number).
Output mit Geschwindigkeit:


Result: 55

real	0m0.598s

user	0m0.717s

sys	0m0.020s

und dazu der Java-Code:
Java
public final class Prob35 {

    private final static boolean isPrime(int n) {
        if(n == 1)
            return false;
 
        if(n == 2)
            return true;
 
        int r = (int) Math.floor(Math.sqrt(n));
        int i;

        for(i = 2; i <= r && n % i != 0; ++i) {}
 
        return (r == i - 1);
    }

    private final static boolean isCircular(int n) {
        String s = Integer.toString(n);

        for (int i = 0; i < s.length(); i++)
            if(!isPrime((Integer.parseInt(s.substring(i) + s.substring(0, i)))))
                return false;

        return true;
    }

    public static void main(String[] argv) {
        int cnt = 0; 

        for(int i = 0; i < 1_000_000; i++)
            if(isPrime(i))
                if(isCircular(i))
                    ++cnt;

        System.out.println("Result: "+cnt);
    }
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

						public final class Prob35 {
 
    private final static boolean isPrime(int n) {
        if(n == 1)
            return false;
 
        if(n == 2)
            return true;
 
        int r = (int) Math.floor(Math.sqrt(n));
        int i;
 
        for(i = 2; i <= r && n % i != 0; ++i) {}
 
        return (r == i - 1);
    }
 
    private final static boolean isCircular(int n) {
        String s = Integer.toString(n);
 
        for (int i = 0; i < s.length(); i++)
            if(!isPrime((Integer.parseInt(s.substring(i) + s.substring(0, i)))))
                return false;
 
        return true;
    }
 
    public static void main(String[] argv) {
        int cnt = 0; 
 
        for(int i = 0; i < 1_000_000; i++)
            if(isPrime(i))
                if(isCircular(i))
                    ++cnt;
 
        System.out.println("Result: "+cnt);
    }
}

					

			
 
Veröffentlicht am 19. Mai 2016 von admin in Java, Programmierung, projecteuler
Problem 44 auf projecteuler.netHier meine Lösung zu Problem 44 auf projecteuler.net. Dank dieser Wiki-Seite konnte man schnell eine Funktion schreiben, die prüft ob eine Zahl eine Pentagonalzahl (oder Fünfeckszahl) ist. Die Lösung habe ich wieder in C geschrieben. Mit folgendem gcc-Befehl lässt sich der Code compilieren und ausführen:


gcc -Wall -O3 -march=native -mtune=native -ffast-math -std=c11 Prob44.c -lm -o Prob44 && time ./Prob44

Das Ergebnis mit der Rechendauer sieht dann so aus:


Result found: 5482660

real	0m0.021s

user	0m0.018s

sys	0m0.004s

Und der Code schaut so aus:
C
#include <stdio.h>

typedef enum { false, true } bool;

bool isPentagonal(long n) {
    double x = (__builtin_sqrt(24 * n + 1) + 1) / 6;
    return (x == (long) x) ? true : false;
}

int main(void) {
    bool flag = true;

    for(int i = 1; flag == true; ++i) {
        int n = i * (3 * i - 1) / 2;

        for(int j = i - 1; j > 0; j--) {
            int d = j * (3 * j - 1) / 2;

            if(isPentagonal(n - d) == true && isPentagonal(n + d) == true) {
                printf("Result found: %d\n", (n - d));
                flag = false;
            }
        }
    }
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25

						#include <stdio.h>
 
typedef enum { false, true } bool;
 
bool isPentagonal(long n) {
    double x = (__builtin_sqrt(24 * n + 1) + 1) / 6;
    return (x == (long) x) ? true : false;
}
 
int main(void) {
    bool flag = true;
 
    for(int i = 1; flag == true; ++i) {
        int n = i * (3 * i - 1) / 2;
 
        for(int j = i - 1; j > 0; j--) {
            int d = j * (3 * j - 1) / 2;
 
            if(isPentagonal(n - d) == true && isPentagonal(n + d) == true) {
                printf("Result found: %d\n", (n - d));
                flag = false;
            }
        }
    }
}

					

			
 
Veröffentlicht am 30. Dezember 2015 von admin in C, Programmierung, projecteuler
Problem 5 auf projecteuler.netHier meine C Lösung zu Problem 5 auf projecteuler.net. Weil ich faul war habe ich einfach auf Brute Force gesetzt und alle Zahlen durchprobiert. Trotz Brute Force braucht das Programm nur 1,8 Sekunden um die Lösung zu berechnen:


gcc -Wall -O3 -march=native -mtune=native -std=c11 Prob5.c -lm -o Prob5 && time ./Prob5

Result: 232792560

real	0m1.801s

user	0m1.800s

sys	0m0.000s

Und hier der Code:
C
#include <stdio.h>

int main(void) {
    for(int i = 1;;i++) {
        int tmp = 0;
        for(int j = 2; j <= 20; j++) {
            if(i % j != 0) {
                tmp = 1;
                break;
            }
        }
        if(tmp == 0) {
            printf("Result: %i\n", i);
            break;
        }
    }
    return 0;
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

						#include <stdio.h>
 
int main(void) {
    for(int i = 1;;i++) {
        int tmp = 0;
        for(int j = 2; j <= 20; j++) {
            if(i % j != 0) {
                tmp = 1;
                break;
            }
        }
        if(tmp == 0) {
            printf("Result: %i\n", i);
            break;
        }
    }
    return 0;
}

					

			
 
Veröffentlicht am 19. Dezember 2015 von admin in C, Programmierung, projecteuler
Problem 3 auf projecteuler.netHier mal wieder eine weitere Project-Euler Lösung. Dieses Mal für das Problem 3 und wieder in der Programmiersprache C.


gcc -Wall -march=native -mtune=native -O3 -std=c11 -ffast-math prob3.c -lm -o prob3 && time ./prob3

Result: 6857

real	0m0.007s

user	0m0.007s

sys	0m0.000s

Hier der Code:
C
#include <stdio.h>

#define NUM 600851475143

typedef enum { false, true } bool;

inline static unsigned long long square(unsigned long long n) {
    return n * n;
}

static bool isPrime(unsigned long long n) {
    if(n == 1)
        return false;

    else if(n == 2)
        return true;

    else {
        unsigned long long i, r = (unsigned long long) __builtin_sqrt(n);

        for(i = 3; i <= r && n % i != 0; i+=2);

        return (r == i - 1) ? true : false;
    }
}

int main(void) {
    unsigned long long max_fac = 0;

    for(unsigned long long x = 1; square(x) <= NUM; ++x)
        if((NUM % x) == 0)
            if(isPrime(x) == true)
                if(x > max_fac)
                    max_fac = x;

    printf("Result: %llu\n", max_fac);
    return 0;
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38

						#include <stdio.h>
 
#define NUM 600851475143
 
typedef enum { false, true } bool;
 
inline static unsigned long long square(unsigned long long n) {
    return n * n;
}
 
static bool isPrime(unsigned long long n) {
    if(n == 1)
        return false;
 
    else if(n == 2)
        return true;
 
    else {
        unsigned long long i, r = (unsigned long long) __builtin_sqrt(n);
 
        for(i = 3; i <= r && n % i != 0; i+=2);
 
        return (r == i - 1) ? true : false;
    }
}
 
int main(void) {
    unsigned long long max_fac = 0;
 
    for(unsigned long long x = 1; square(x) <= NUM; ++x)
        if((NUM % x) == 0)
            if(isPrime(x) == true)
                if(x > max_fac)
                    max_fac = x;
 
    printf("Result: %llu\n", max_fac);
    return 0;
}

					

			
 
Veröffentlicht am 18. Dezember 2015 von admin in C, Programmierung, projecteuler
Problem 493 auf projecteuler.netHier die Lösung in der Programmiersprache C zu Problem 493 auf projecteuler.net
Compiliert und ausgeführt wurde mit folgendem Befehl:


gcc -Wall -O3 -march=native -mtune=native -std=c11 Prob493.c -lm -o Prob493 && time ./Prob493

Als Ergbnis kommt dann folgendes heraus:


Result: 6.818741802

real	0m0.001s

user	0m0.000s

sys	0m0.000s

Der Code sieht so aus:
C
#include <stdio.h>

double factorial(double n) {
    return (n == 1) ? 1 : factorial(n - 1) * n;
}

double binomial(double n, double k) {
    return factorial(n) / (factorial(k) * factorial(n - k));
}

int main(void) {
    double result = 7 * (1 - binomial(60, 20) / binomial(70, 20));
    printf("Result: %.9f", result);

    return 0;
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16

						#include <stdio.h>
 
double factorial(double n) {
    return (n == 1) ? 1 : factorial(n - 1) * n;
}
 
double binomial(double n, double k) {
    return factorial(n) / (factorial(k) * factorial(n - k));
}
 
int main(void) {
    double result = 7 * (1 - binomial(60, 20) / binomial(70, 20));
    printf("Result: %.9f", result);
 
    return 0;
}

					

			
 
Veröffentlicht am 25. Oktober 2015 von admin in C, Programmierung, projecteuler
Problem 6 auf projecteuler.netAls ich heute Langeweile hatte, habe ich einmal die Lösung für das Problem 6 in Bash geschrieben, weil ich ohnehin schon den ganzen Tag an verschiedenen Bash-Skripten saß. Das Ergebnis für das Problem 6 berechnet sich sehr schnell:


25164150

real	0m0.002s

user	0m0.000s

sys	0m0.002s

Der Code schaut so aus:
Shell
#!/usr/bin/env bash

PART1=0
PART2=0
RESULT=0

for i in {1..100}
do
    PART1=$((PART1 + i * i))
    PART2=$((PART2 + i))
done

RESULT=$((PART2 * PART2 - PART1))

echo $RESULT

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15

						#!/usr/bin/env bash
 
PART1=0
PART2=0
RESULT=0
 
for i in {1..100}
do
    PART1=$((PART1 + i * i))
    PART2=$((PART2 + i))
done
 
RESULT=$((PART2 * PART2 - PART1))
 
echo $RESULT

					

			
 
Veröffentlicht am 16. Oktober 2015 von admin in Bash, Programmierung, projecteuler
Problem 102 auf projecteuler.netHier ist mal wieder eine neue Lösung auf projecteuler.net, dieses Mal für Problem 102.
Als erstes wurden 2 Methoden erstellt. Eine für die Flächenberechnung des Dreiecks (https://en.wikipedia.org/wiki/Triangle#Using_coordinates) und eine, die prüft, ob der Ursprung im Dreieck enthalten ist.
Danach wurde nur noch die Textdatei mit den zufällig generierten Werten eingelesen und entsprechend aufgesplittet sodass das Programm die Daten nutzen kann.
Die Berechnung ging ebenfalls wieder sehr schnell:



Result: 228

real	0m0.161s

user	0m0.231s

sys	0m0.022s

Hier der Code in Java:
Java
import java.io.IOException;
import java.nio.charset.StandardCharsets;
import java.nio.file.Files;
import java.nio.file.Paths;
import java.util.Arrays;

public final class Prob102 {

    private final static int areaOfTriangle(int ax, int ay, int bx, int by, int cx, int cy) {
        return Math.abs((ax - cx) * (by - ay) - (ax - bx) * (cy - ay));
    }

    private final static boolean triangleContainsPoint(int ax, int ay, int bx, int by, int cx, int cy, int px, int py) {
        return (areaOfTriangle(ax, ay, bx, by, cx, cy) ==
                areaOfTriangle(ax, ay, bx, by, px, py) +
                areaOfTriangle(ax, ay, px, py, cx, cy) +
                areaOfTriangle(px, py, bx, by, cx, cy));
    }

    public static void main(String[] args) {
        String[] l = null;
        int result = 0;

        try {
            final byte[] f = Files.readAllBytes(Paths.get("triangles.txt"));
            l = (new String(f, StandardCharsets.UTF_8)).split("\\r?\\n");
        }
        catch (IOException e) {
            e.printStackTrace();
        }

        for(int i = 0; i < l.length; i++) {
            int[] C = Arrays.stream(l[i].split(",")).mapToInt(Integer::parseInt).toArray();

            if(triangleContainsPoint(C[0], C[1], C[2], C[3], C[4], C[5], 0, 0))
                result++;
        }

        System.out.println("Result: "+result);
    }
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41

						import java.io.IOException;
import java.nio.charset.StandardCharsets;
import java.nio.file.Files;
import java.nio.file.Paths;
import java.util.Arrays;
 
public final class Prob102 {
 
    private final static int areaOfTriangle(int ax, int ay, int bx, int by, int cx, int cy) {
        return Math.abs((ax - cx) * (by - ay) - (ax - bx) * (cy - ay));
    }
 
    private final static boolean triangleContainsPoint(int ax, int ay, int bx, int by, int cx, int cy, int px, int py) {
        return (areaOfTriangle(ax, ay, bx, by, cx, cy) ==
                areaOfTriangle(ax, ay, bx, by, px, py) +
                areaOfTriangle(ax, ay, px, py, cx, cy) +
                areaOfTriangle(px, py, bx, by, cx, cy));
    }
 
    public static void main(String[] args) {
        String[] l = null;
        int result = 0;
 
        try {
            final byte[] f = Files.readAllBytes(Paths.get("triangles.txt"));
            l = (new String(f, StandardCharsets.UTF_8)).split("\\r?\\n");
        }
        catch (IOException e) {
            e.printStackTrace();
        }
 
        for(int i = 0; i < l.length; i++) {
            int[] C = Arrays.stream(l[i].split(",")).mapToInt(Integer::parseInt).toArray();
 
            if(triangleContainsPoint(C[0], C[1], C[2], C[3], C[4], C[5], 0, 0))
                result++;
        }
 
        System.out.println("Result: "+result);
    }
}

					

			
 
Veröffentlicht am 30. August 2015 von admin in Java, Programmierung, projecteuler
Problem 89 auf projecteuler.netHeute hatte ich mal wieder etwas Zeit ein weiteres Problem auf projecteuler.net zu lösen. Gelöst habe ich dieses Mal Problem 89 in der Programmiersprache Java.
Das Problem ist dank dieser Hilfe sehr einfach: https://projecteuler.net/about=roman_numerals
Man muss also nur „DCCCC“, „LXXXX“, „VIIII“, „CCCC“, „XXXX“ und „IIII“ ersetzen. Danach muss man nur noch die Differenz der Anzahl der Zeichen vor und nach dem Replacen herausfinden.
Das Programm von mir berechnet die Lösung sehr schnell:


Result: 743

real	0m0.163s

user	0m0.170s

sys	0m0.034s

Und so schaut der Code aus:
Java
import java.io.IOException;
import java.nio.charset.StandardCharsets;
import java.nio.file.Files;
import java.nio.file.Paths;

public final class Prob89 {
    private static String romanNumbers;
    private static String[] singleRomanNumbers;

    public static void main(String[] args) {
        try {
            final byte[] encoded = Files.readAllBytes(Paths.get("89.txt"));
            romanNumbers = new String(encoded, StandardCharsets.UTF_8);
        } catch (IOException e) {
            e.printStackTrace();
        }

        singleRomanNumbers = romanNumbers.split("\\r?\\n");

        for (int i = 0; i < singleRomanNumbers.length; i++) {
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("DCCCC", "CM");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("LXXXX", "XC");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("VIIII", "IX");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("CCCC", "CD");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("XXXX", "XL");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("IIII", "IV");
        }

        StringBuilder b = new StringBuilder();

        for (String s : singleRomanNumbers)
            b.append(s);

        String newRomanNumbers = b.toString();

        final int result = romanNumbers.replace("\n", "").length() - newRomanNumbers.length();
        System.out.println("Result: " + result);
    }
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39

						import java.io.IOException;
import java.nio.charset.StandardCharsets;
import java.nio.file.Files;
import java.nio.file.Paths;
 
public final class Prob89 {
    private static String romanNumbers;
    private static String[] singleRomanNumbers;
 
    public static void main(String[] args) {
        try {
            final byte[] encoded = Files.readAllBytes(Paths.get("89.txt"));
            romanNumbers = new String(encoded, StandardCharsets.UTF_8);
        } catch (IOException e) {
            e.printStackTrace();
        }
 
        singleRomanNumbers = romanNumbers.split("\\r?\\n");
 
        for (int i = 0; i < singleRomanNumbers.length; i++) {
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("DCCCC", "CM");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("LXXXX", "XC");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("VIIII", "IX");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("CCCC", "CD");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("XXXX", "XL");
            singleRomanNumbers[i] = singleRomanNumbers[i].replace("IIII", "IV");
        }
 
        StringBuilder b = new StringBuilder();
 
        for (String s : singleRomanNumbers)
            b.append(s);
 
        String newRomanNumbers = b.toString();
 
        final int result = romanNumbers.replace("\n", "").length() - newRomanNumbers.length();
        System.out.println("Result: " + result);
    }
}

					

			
 
Veröffentlicht am 18. Juli 2015 von admin in Java, Programmierung, projecteuler
Problem 55 auf projecteuler.netDieses Wochenende habe ich mich mal wieder hingesetzt und eine weiteres Problem auf projecteuler.net gelöst.
Gelöst habe ich Problem 55 in der Programmiersprache C.
Kompiliert habe ich mit folgendem Befehl:



gcc -O3 -Wall Problem55.c -lm -o Problem55

Das Ergebnis berechnet sich sehr schnell:



Result: 249

real 0m0.019s

user 0m0.019s

sys  0m0.000s

So schaut der Quellcode des Programms aus:
C
#include <stdio.h>

__uint128_t reverseInt(__uint128_t n) {
   __uint128_t r = 0;

   while (n != 0) {
      r *= 10;
      r += n % 10;
      n /= 10;
   }
   
   return r;
}

int palindromic(__uint128_t n) {
    return (n == reverseInt(n)) ? 0 : 1;
}

int main(void) {
    unsigned int res = 0, i, t;
    __uint128_t r;
	
    for(i = 0; i < 10000; i++) {
        r = i + reverseInt(i);

        for(t = 0; t < 50; t++) {
            if(palindromic(r) == 0)
                break;
            else
                r += reverseInt(r);
        }
		
        if(palindromic(r) == 1)
            res++;
    }
	
    printf("Result: %d!\n", res);
    return 0;
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39

						#include <stdio.h>
 
__uint128_t reverseInt(__uint128_t n) {
   __uint128_t r = 0;
 
   while (n != 0) {
      r *= 10;
      r += n % 10;
      n /= 10;
   }
   
   return r;
}
 
int palindromic(__uint128_t n) {
    return (n == reverseInt(n)) ? 0 : 1;
}
 
int main(void) {
    unsigned int res = 0, i, t;
    __uint128_t r;
	
    for(i = 0; i < 10000; i++) {
        r = i + reverseInt(i);
 
        for(t = 0; t < 50; t++) {
            if(palindromic(r) == 0)
                break;
            else
                r += reverseInt(r);
        }
		
        if(palindromic(r) == 1)
            res++;
    }
	
    printf("Result: %d!\n", res);
    return 0;
}

					

			
 
Veröffentlicht am 14. Juni 2015 von admin in C, Programmierung, projecteuler
Problem 26 auf projecteuler.netHier meine erste Python-Lösung auf projecteuler.net für das Problem 26.


Result: 983

real	0m0.246s

user	0m0.241s

sys	0m0.004s

Python
result = [0, 0]

for i in range(1, 1001):
    a = [0, 1000]
    num = 10 % i

    while num != 1 and a[1] > 0:
        num *= 10
        num %= i
        a[0] += 1
        a[1] -= 1

    if a[0] > result[0] and a[1] > 1:
        result[0] = a[0]
        result[1] = i

print("Result: ", result[1])

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

						result = [0, 0]
 
for i in range(1, 1001):
    a = [0, 1000]
    num = 10 % i
 
    while num != 1 and a[1] > 0:
        num *= 10
        num %= i
        a[0] += 1
        a[1] -= 1
 
    if a[0] > result[0] and a[1] > 1:
        result[0] = a[0]
        result[1] = i
 
print("Result: ", result[1])

					

			
 
Veröffentlicht am 3. Mai 2015 von admin in Programmierung, projecteuler, Python
Problem 27 auf projecteuler.netHier meine Lösung zu Problem 27 auf projecteuler.net in der Programmiersprache Java.
Bei der Lösung habe ich auf brute force gesetzt. Trotzdem berechnet sich das Ergebnis recht schnell:


Result: -59231

real    0m0.316s

user    0m0.325s

sys     0m0.008s

Java
package net.mborm.euler.problems;

public class Problem27 {
    private static boolean is_prime(int n) {
        if(n == 1)
            return false;

        if(n == 2)
            return true;

        int r = (int) Math.floor(Math.sqrt(n));
        int i;

        for(i = 2; i <= r; ++i)
            if(n % i == 0)
                break;

        return (r == i - 1);
    }

    public static void main(String[] args) {
        int[] sequence = new int[2];
        int result = 0;

        for(int a = -1000; a <= 1000; a++) {
            for(int b = -1000; b <= 1000; b++) {
                int n = 0;

                while(is_prime(Math.abs(n * n + a * n + b)))
                    n++;

                if(n > result) {
                    result = n;
                    sequence[0] = a;
                    sequence[1] = b;
                }
            }
        }

        System.out.println("Result: "+sequence[0] * sequence[1]);

    }
}

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43

						package net.mborm.euler.problems;
 
public class Problem27 {
    private static boolean is_prime(int n) {
        if(n == 1)
            return false;
 
        if(n == 2)
            return true;
 
        int r = (int) Math.floor(Math.sqrt(n));
        int i;
 
        for(i = 2; i <= r; ++i)
            if(n % i == 0)
                break;
 
        return (r == i - 1);
    }
 
    public static void main(String[] args) {
        int[] sequence = new int[2];
        int result = 0;
 
        for(int a = -1000; a <= 1000; a++) {
            for(int b = -1000; b <= 1000; b++) {
                int n = 0;
 
                while(is_prime(Math.abs(n * n + a * n + b)))
                    n++;
 
                if(n > result) {
                    result = n;
                    sequence[0] = a;
                    sequence[1] = b;
                }
            }
        }
 
        System.out.println("Result: "+sequence[0] * sequence[1]);
 
    }
}

					

			
 
Veröffentlicht am 30. April 2015 von admin in Java, Programmierung, projecteuler
Problem 11 auf projecteuler.netDie Lösung zu Problem 11 habe ich ebenfalls in PHP umgesetzt.
Meiner Meinung nach war das Schwierigste an der Aufgabe sich in verschiedene Richtungen im Array zu bewegen. Wenn das klappt, berechnet sich die Lösung allerdings sehr schnell:


Result: 70600674

real	0m0.012s

user	0m0.008s

sys	0m0.004s

PHP
<?php

$p = 0;

$grid = "08 02 22 97 38 15 00 40 00 75 04 05 07 78 52 12 50 77 91 08 
49 49 99 40 17 81 18 57 60 87 17 40 98 43 69 48 04 56 62 00 
81 49 31 73 55 79 14 29 93 71 40 67 53 88 30 03 49 13 36 65 
52 70 95 23 04 60 11 42 69 24 68 56 01 32 56 71 37 02 36 91 
22 31 16 71 51 67 63 89 41 92 36 54 22 40 40 28 66 33 13 80 
24 47 32 60 99 03 45 02 44 75 33 53 78 36 84 20 35 17 12 50 
32 98 81 28 64 23 67 10 26 38 40 67 59 54 70 66 18 38 64 70 
67 26 20 68 02 62 12 20 95 63 94 39 63 08 40 91 66 49 94 21 
24 55 58 05 66 73 99 26 97 17 78 78 96 83 14 88 34 89 63 72 
21 36 23 09 75 00 76 44 20 45 35 14 00 61 33 97 34 31 33 95 
78 17 53 28 22 75 31 67 15 94 03 80 04 62 16 14 09 53 56 92 
16 39 05 42 96 35 31 47 55 58 88 24 00 17 54 24 36 29 85 57 
86 56 00 48 35 71 89 07 05 44 44 37 44 60 21 58 51 54 17 58 
19 80 81 68 05 94 47 69 28 73 92 13 86 52 17 77 04 89 55 40 
04 52 08 83 97 35 99 16 07 97 57 32 16 26 26 79 33 27 98 66 
88 36 68 87 57 62 20 72 03 46 33 67 46 55 12 32 63 93 53 69 
04 42 16 73 38 25 39 11 24 94 72 18 08 46 29 32 40 62 76 36 
20 69 36 41 72 30 23 88 34 62 99 69 82 67 59 85 74 04 36 16 
20 73 35 29 78 31 90 01 74 31 49 71 48 86 81 16 23 57 05 54 
01 70 54 71 83 51 54 69 16 92 33 48 61 43 52 01 89 19 67 48 ";


$grid = preg_replace('/0(\d)/', " $1", $grid);
$grid = explode("\n", $grid);

for ($i = 0; $i < count($grid); $i++)
  $grid[$i] = preg_split('/[ ]+/', trim($grid[$i]));

for($n = 0; $n <= 19; $n++) {
    for($k = 0; $k <= 16; $k++) {
        $p = max($p, $grid[$n][$k+0] * $grid[$n][$k+1] * $grid[$n][$k+2] * $grid[$n][$k+3]);
        $p = max($p, $grid[$k+0][$n] * $grid[$k+1][$n] * $grid[$k+2][$n] * $grid[$k+3][$n]);
    }
}

for($n = 0; $n <= 16; $n++) {
    for($k = 0; $k <= 16; $k++) {
        $p = max($p, $grid[$n+0][$k+0] * $grid[$n+1][$k+1] * $grid[$n+2][$k+2] * $grid[$n+3][$k+3]);
        $p = max($p, $grid[$n+0][$k+3] * $grid[$n+1][$k+2] * $grid[$n+2][$k+1] * $grid[$n+3][$k+0]);
    }
}

echo "Result: ".$p."\n";

?>

				
					
				1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49

						<?php
 
$p = 0;
 
$grid = "08 02 22 97 38 15 00 40 00 75 04 05 07 78 52 12 50 77 91 08 
49 49 99 40 17 81 18 57 60 87 17 40 98 43 69 48 04 56 62 00 
81 49 31 73 55 79 14 29 93 71 40 67 53 88 30 03 49 13 36 65 
52 70 95 23 04 60 11 42 69 24 68 56 01 32 56 71 37 02 36 91 
22 31 16 71 51 67 63 89 41 92 36 54 22 40 40 28 66 33 13 80 
24 47 32 60 99 03 45 02 44 75 33 53 78 36 84 20 35 17 12 50 
32 98 81 28 64 23 67 10 26 38 40 67 59 54 70 66 18 38 64 70 
67 26 20 68 02 62 12 20 95 63 94 39 63 08 40 91 66 49 94 21 
24 55 58 05 66 73 99 26 97 17 78 78 96 83 14 88 34 89 63 72 
21 36 23 09 75 00 76 44 20 45 35 14 00 61 33 97 34 31 33 95 
78 17 53 28 22 75 31 67 15 94 03 80 04 62 16 14 09 53 56 92 
16 39 05 42 96 35 31 47 55 58 88 24 00 17 54 24 36 29 85 57 
86 56 00 48 35 71 89 07 05 44 44 37 44 60 21 58 51 54 17 58 
19 80 81 68 05 94 47 69 28 73 92 13 86 52 17 77 04 89 55 40 
04 52 08 83 97 35 99 16 07 97 57 32 16 26 26 79 33 27 98 66 
88 36 68 87 57 62 20 72 03 46 33 67 46 55 12 32 63 93 53 69 
04 42 16 73 38 25 39 11 24 94 72 18 08 46 29 32 40 62 76 36 
20 69 36 41 72 30 23 88 34 62 99 69 82 67 59 85 74 04 36 16 
20 73 35 29 78 31 90 01 74 31 49 71 48 86 81 16 23 57 05 54 
01 70 54 71 83 51 54 69 16 92 33 48 61 43 52 01 89 19 67 48 ";
 
 
$grid = preg_replace('/0(\d)/', " $1", $grid);
$grid = explode("\n", $grid);
 
for ($i = 0; $i < count($grid); $i++)
  $grid[$i] = preg_split('/[ ]+/', trim($grid[$i]));
 
for($n = 0; $n <= 19; $n++) {
    for($k = 0; $k <= 16; $k++) {
        $p = max($p, $grid[$n][$k+0] * $grid[$n][$k+1] * $grid[$n][$k+2] * $grid[$n][$k+3]);
        $p = max($p, $grid[$k+0][$n] * $grid[$k+1][$n] * $grid[$k+2][$n] * $grid[$k+3][$n]);
    }
}
 
for($n = 0; $n <= 16; $n++) {
    for($k = 0; $k <= 16; $k++) {
        $p = max($p, $grid[$n+0][$k+0] * $grid[$n+1][$k+1] * $grid[$n+2][$k+2] * $grid[$n+3][$k+3]);
        $p = max($p, $grid[$n+0][$k+3] * $grid[$n+1][$k+2] * $grid[$n+2][$k+1] * $grid[$n+3][$k+0]);
    }
}
 
echo "Result: ".$p."\n";
 
?>

					

			
 
Veröffentlicht am 28. April 2015 von admin in PHP, Programmierung, projecteuler
Problem 36 auf projecteuler.netHier meine Lösung zu Problem 36 auf projecteuler.net
Die Lösung für dieses Problem ließ sich extrem einfach in PHP bauen – und das Ganze läuft innerhalb von ~ 1 Sekunde:


872187

real	0m1.049s

user	0m1.048s

sys	0m0.000s

PHP
<?php

$result = 0;

for($i = 0; $i <= 1000000; $i++)
    if(decbin($i) == strrev(decbin($i)) && $i == strrev($i))
        $result += $i;

print $result;

?>

				
					
				1
2
3
4
5
6
7
8
9
10
11

						<?php
 
$result = 0;
 
for($i = 0; $i <= 1000000; $i++)
    if(decbin($i) == strrev(decbin($i)) && $i == strrev($i))
        $result += $i;
 
print $result;
 
?>

					

			
 
Veröffentlicht am 28. April 2015 von admin in PHP, Programmierung, projecteuler
Older posts

		Neueste Beiträge
		
					apache2 TLS-Config für Qualys Lab Bestnote mit Server Hardening
									

					TEDP – The Energy Data Project
									

					GPU Benchmark integriert
									

					Neue Benchmark Version
									

					Eigene API für die luftdaten.info Wetterstation
									


		
KategorienAllgemein
Bash
C
Hardware
IT-Sicherheit
Java
Netzwerk
Perl
PHP
Programmierung
projecteuler
Python
Server
ArchivJanuar 2021
November 2020
Oktober 2019
September 2019
März 2019
Juli 2018
Juni 2018
August 2017
Mai 2016
Dezember 2015
Oktober 2015
August 2015
Juli 2015
Juni 2015
Mai 2015
April 2015
mborm.net

Ein Computerprogramm tut was Du schreibst, nicht was Du willst.

Neueste Beiträge

Kategorien

Archiv
