mborm.net

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

		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.

Problem 44 auf projecteuler.net

Problem 5 auf projecteuler.net

Problem 3 auf projecteuler.net

Neueste Beiträge

Kategorien

Archiv
