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 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 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 92 auf projecteuler.netHier ist meine Lösung zu Problem 92 auf projecteuler.net. Die Lösung ist mittels der Programmiersprache C umgesetzt und kann mittels des gcc mit diesem Befehl kompiliert werden:


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

Auf meinem Rechner braucht das Skript nur eine sehr kurze Zeit zum Berechnen:



real	0m0.639s

user	0m0.638s

sys	0m0.000s

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

typedef unsigned int uint;

__inline static uint square(uint x) {
    return x * x;
}

__inline static uint solve(uint num) {
    register uint res = 0, val = num;

    while(val > 0) {
        res += square(val % 10);
        val /= 10;
    }

    return res;
}

int main() {
    typedef enum { false, true } bool;
    uint result = 0, i = 0, tmp = 0;

    for(i = 10000000; i != 0; i--) {
        tmp = solve(i);
        bool run = true;

        while(run) {
            switch(tmp) {
                case 1: run = false; break;
                case 89: result++; run = false; break;
                default: tmp = solve(tmp);
            }
        }
    }
    printf("Result: %i", result);
    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>
 
typedef unsigned int uint;
 
__inline static uint square(uint x) {
    return x * x;
}
 
__inline static uint solve(uint num) {
    register uint res = 0, val = num;
 
    while(val > 0) {
        res += square(val % 10);
        val /= 10;
    }
 
    return res;
}
 
int main() {
    typedef enum { false, true } bool;
    uint result = 0, i = 0, tmp = 0;
 
    for(i = 10000000; i != 0; i--) {
        tmp = solve(i);
        bool run = true;
 
        while(run) {
            switch(tmp) {
                case 1: run = false; break;
                case 89: result++; run = false; break;
                default: tmp = solve(tmp);
            }
        }
    }
    printf("Result: %i", result);
    return 0;
}

					

			
 
Veröffentlicht am 14. April 2015 von admin in C, Programmierung, projecteuler

		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
