summaryrefslogtreecommitdiff
path: root/euler346.c
blob: dc10d6513b22c0cc0bbb9b83547177eb08bf25bc (plain)
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
#include <stdio.h>
#include <stdlib.h>
#include <assert.h>

#define MAXN 1000000000000
#define NRUMAX 10000000

long long ru[NRUMAX];
int nru = 0;

int comp(const void *a, const void *b)
{
	long long x = *(long long*)a;
	long long y = *(long long*)b;

	if (x < y)
		return -1;
	else if (x > y)
		return 1;
	else 
		return 0;
}

int main()
{
	long long sum_repunit = 1; /* 1 is a strong repunit */
	for (int base = 2; base < 1000000; base++) {
		long long repunit = 1 + base;
		long long nb = base;
		while (1) {
			nb *= base;
			repunit += nb;
			if (repunit >= MAXN)
				break;
			assert(nru < NRUMAX);
			ru[nru] = repunit;
			nru++;
		}
	}

	qsort(ru, nru, sizeof(long long), comp);

	printf("Find %d strong repunits.\n", nru + 1);
	for (int i = 0; i < nru; i++) {
		if (i > 0 && ru[i] == ru[i-1])
			continue;
		sum_repunit += ru[i];
	}
	printf("Their sum is %lld\n", sum_repunit);
}