indivisible-legacy/src/main.c

70 lines
1.4 KiB
C
Raw Normal View History

#include <stdio.h>
#include <stdlib.h>
#include <signal.h>
#include <stdbool.h>
#include <gmp.h>
#include "list.h"
#include "types.h"
2016-12-09 23:50:23 +00:00
#include "optimizers.h"
static bool run;
void leave();
int main(void) {
printf("Indivisible v0.2\n");
// Quit on ^C by setting `run = false'
run = true;
signal(SIGINT, leave);
// Primes we've found
List primes;
initList(&primes);
// The number we're going to be testing for
mpz_t num;
mpz_init(num);
// Add 2, a known prime to this list
mpz_set_ui(num, 2);
addToList(&primes, num);
if(mpz_out_str(stdout, 10, num) == 0) {
fprintf(stderr, "Could not print to `stdout'!\n");
exit(1);
}
printf("\n");
mpz_add_ui(num, num, 1);
do {
// Loop through found primes
for(ulli i = 0; i < primes.size; ++i) {
// If `num' is divisible by a prime then go to the next number
2016-12-10 10:20:01 +00:00
if(mpz_divisible_p(num, primes.list[i]) != 0)
goto nextPrime;
}
// `num' is a prime so we add it to the list and print it
addToList(&primes, num);
2016-12-09 23:50:23 +00:00
if(unlikely(mpz_out_str(stdout, 10, num) == 0)) {
fprintf(stderr, "Could not print to `stdout'!\n");
exit(1);
}
printf("\n");
nextPrime:
// Add 2 (skip even numbers since they're all divisible by 2)
mpz_add_ui(num, num, 2);
2016-12-10 00:15:52 +00:00
} while(likely(run));
2016-12-10 10:20:01 +00:00
mpz_clear(num);
// Deinitialize the list
deInitList(&primes);
return 0;
}
void leave() {
printf("Exiting...\n");
run = false;
}