2019-07-16 23:09:53 +00:00
|
|
|
"""
|
2020-10-09 03:03:23 +00:00
|
|
|
Problem 14: https://projecteuler.net/problem=14
|
|
|
|
|
2019-07-16 23:09:53 +00:00
|
|
|
Collatz conjecture: start with any positive integer n. Next term obtained from
|
|
|
|
the previous term as follows:
|
|
|
|
|
|
|
|
If the previous term is even, the next term is one half the previous term.
|
|
|
|
If the previous term is odd, the next term is 3 times the previous term plus 1.
|
|
|
|
The conjecture states the sequence will always reach 1 regardless of starting
|
|
|
|
n.
|
|
|
|
|
|
|
|
Problem Statement:
|
|
|
|
The following iterative sequence is defined for the set of positive integers:
|
|
|
|
|
|
|
|
n → n/2 (n is even)
|
|
|
|
n → 3n + 1 (n is odd)
|
|
|
|
|
|
|
|
Using the rule above and starting with 13, we generate the following sequence:
|
|
|
|
|
|
|
|
13 → 40 → 20 → 10 → 5 → 16 → 8 → 4 → 2 → 1
|
|
|
|
|
|
|
|
It can be seen that this sequence (starting at 13 and finishing at 1) contains
|
|
|
|
10 terms. Although it has not been proved yet (Collatz Problem), it is thought
|
|
|
|
that all starting numbers finish at 1.
|
|
|
|
|
|
|
|
Which starting number, under one million, produces the longest chain?
|
|
|
|
"""
|
2021-09-07 11:37:03 +00:00
|
|
|
from __future__ import annotations
|
2019-10-05 05:14:13 +00:00
|
|
|
|
|
|
|
|
2021-09-23 18:55:18 +00:00
|
|
|
def collatz_sequence_length(n: int) -> int:
|
|
|
|
"""Returns the Collatz sequence length for n."""
|
|
|
|
sequence_length = 1
|
2019-07-16 23:09:53 +00:00
|
|
|
while n != 1:
|
|
|
|
if n % 2 == 0:
|
|
|
|
n //= 2
|
|
|
|
else:
|
|
|
|
n = 3 * n + 1
|
2021-09-23 18:55:18 +00:00
|
|
|
sequence_length += 1
|
|
|
|
return sequence_length
|
2019-07-16 23:09:53 +00:00
|
|
|
|
|
|
|
|
2020-10-09 03:03:23 +00:00
|
|
|
def solution(n: int = 1000000) -> int:
|
2019-07-16 23:09:53 +00:00
|
|
|
"""Returns the number under n that generates the longest Collatz sequence.
|
|
|
|
|
2019-07-18 22:34:29 +00:00
|
|
|
# The code below has been commented due to slow execution affecting Travis.
|
|
|
|
# >>> solution(1000000)
|
2020-10-09 03:03:23 +00:00
|
|
|
# 837799
|
2019-07-16 23:09:53 +00:00
|
|
|
>>> solution(200)
|
2020-10-09 03:03:23 +00:00
|
|
|
171
|
2019-07-16 23:09:53 +00:00
|
|
|
>>> solution(5000)
|
2020-10-09 03:03:23 +00:00
|
|
|
3711
|
2019-07-16 23:09:53 +00:00
|
|
|
>>> solution(15000)
|
2020-10-09 03:03:23 +00:00
|
|
|
13255
|
2019-07-16 23:09:53 +00:00
|
|
|
"""
|
|
|
|
|
2021-09-23 18:55:18 +00:00
|
|
|
result = max((collatz_sequence_length(i), i) for i in range(1, n))
|
2020-10-09 03:03:23 +00:00
|
|
|
return result[1]
|
2019-07-16 23:09:53 +00:00
|
|
|
|
|
|
|
|
|
|
|
if __name__ == "__main__":
|
2020-10-09 03:03:23 +00:00
|
|
|
print(solution(int(input().strip())))
|