forked from nayuki/Project-Euler-solutions
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathp071.java
More file actions
62 lines (54 loc) · 2.32 KB
/
p071.java
File metadata and controls
62 lines (54 loc) · 2.32 KB
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
59
60
61
62
/*
* Solution to Project Euler problem 71
* Copyright (c) Project Nayuki. All rights reserved.
*
* https://www.nayuki.io/page/project-euler-solutions
* https://github.com/nayuki/Project-Euler-solutions
*/
public final class p071 implements EulerSolution {
public static void main(String[] args) {
System.out.println(new p071().run());
}
/*
* We consider each (integer) denominator d from 1 to 1000000 by brute force.
* For a given d, what is the largest integer n such that n/d < 3/7?
*
* - If d is a multiple of 7, then the integer n' = (d / 7) * 3 satisfies n'/d = 3/7.
* Hence we choose n = n' - 1 = (d / 7) * 3 - 1, so that n/d < 3/7.
* Since (d / 7) * 3 is already an integer, it is equal to floor(d * 3 / 7),
* which will unifie with the next case. Thus n = floor(d * 3 / 7) - 1.
* - Otherwise d is not a multiple of 7, so choosing n = floor(d * 3 / 7)
* will automatically satisfy n/d < 3/7, and be the largest possible n
* due to the definition of the floor function.
*
* When we choose n in this manner, it might not be coprime with d. In other words,
* the simplified form of the fraction n/d might have a denominator smaller than d.
*
* Let's process denominators in ascending order. Each denominator generates a pair
* of integers (n, d) that conceptually represents a fraction, without simplification.
* Whenever the current value of n/d is strictly larger than the previously saved value,
* we save this current value of (n, d).
*
* If we handle denominators in this way - starting from 1, counting up consecutively -
* then it is guaranteed that our final saved pair (n, d) is in lowest terms. This is
* because if (n, d) is not in lowest terms, then its reduced form (n', d') would have
* been saved when the smaller denominator d' was processed, and because n/d is
* not larger than n'/d' (they are equal), the saved value would not be overwritten.
* Hence in this entire computation we can avoid explicitly simplifying any fraction at all.
*/
private static final int LIMIT = 1000000;
public String run() {
int maxN = 0;
int maxD = 1;
for (int d = 1; d <= LIMIT; d++) {
int n = d * 3 / 7;
if (d % 7 == 0)
n--;
if ((long)n * maxD > (long)maxN * d) { // n/d > maxN/maxD
maxN = n;
maxD = d;
}
}
return Integer.toString(maxN);
}
}