2019-07-16 23:09:53 +00:00
|
|
|
"""
|
2020-10-25 03:23:16 +00:00
|
|
|
Project Euler Problem 7: https://projecteuler.net/problem=7
|
2020-10-08 05:57:47 +00:00
|
|
|
|
2020-10-25 03:23:16 +00:00
|
|
|
10001st prime
|
2019-07-16 23:09:53 +00:00
|
|
|
|
2020-10-25 03:23:16 +00:00
|
|
|
By listing the first six prime numbers: 2, 3, 5, 7, 11, and 13, we
|
|
|
|
can see that the 6th prime is 13.
|
2019-07-16 23:09:53 +00:00
|
|
|
|
2020-10-25 03:23:16 +00:00
|
|
|
What is the 10001st prime number?
|
|
|
|
|
|
|
|
References:
|
|
|
|
- https://en.wikipedia.org/wiki/Prime_number
|
2019-07-16 23:09:53 +00:00
|
|
|
"""
|
2022-09-14 08:40:04 +00:00
|
|
|
import math
|
2019-10-05 05:14:13 +00:00
|
|
|
|
|
|
|
|
2022-04-08 17:40:45 +00:00
|
|
|
def is_prime(number: int) -> bool:
|
2022-09-14 08:40:04 +00:00
|
|
|
"""Checks to see if a number is a prime in O(sqrt(n)).
|
|
|
|
A number is prime if it has exactly two factors: 1 and itself.
|
|
|
|
Returns boolean representing primality of given number (i.e., if the
|
|
|
|
result is true, then the number is indeed prime else it is not).
|
2020-10-25 03:23:16 +00:00
|
|
|
|
2022-04-08 17:40:45 +00:00
|
|
|
>>> is_prime(2)
|
2020-10-25 03:23:16 +00:00
|
|
|
True
|
2022-09-14 08:40:04 +00:00
|
|
|
>>> is_prime(3)
|
|
|
|
True
|
|
|
|
>>> is_prime(27)
|
2020-10-25 03:23:16 +00:00
|
|
|
False
|
2022-09-14 08:40:04 +00:00
|
|
|
>>> is_prime(2999)
|
2020-10-25 03:23:16 +00:00
|
|
|
True
|
2022-09-14 08:40:04 +00:00
|
|
|
>>> is_prime(0)
|
|
|
|
False
|
|
|
|
>>> is_prime(1)
|
|
|
|
False
|
2020-10-25 03:23:16 +00:00
|
|
|
"""
|
|
|
|
|
2022-09-14 08:40:04 +00:00
|
|
|
if 1 < number < 4:
|
|
|
|
# 2 and 3 are primes
|
|
|
|
return True
|
|
|
|
elif number < 2 or number % 2 == 0 or number % 3 == 0:
|
|
|
|
# Negatives, 0, 1, all even numbers, all multiples of 3 are not primes
|
|
|
|
return False
|
|
|
|
|
|
|
|
# All primes number are in format of 6k +/- 1
|
|
|
|
for i in range(5, int(math.sqrt(number) + 1), 6):
|
|
|
|
if number % i == 0 or number % (i + 2) == 0:
|
2019-07-16 23:09:53 +00:00
|
|
|
return False
|
|
|
|
return True
|
|
|
|
|
|
|
|
|
2020-10-08 05:57:47 +00:00
|
|
|
def solution(nth: int = 10001) -> int:
|
2020-10-25 03:23:16 +00:00
|
|
|
"""
|
|
|
|
Returns the n-th prime number.
|
2019-08-19 13:37:49 +00:00
|
|
|
|
2019-07-16 23:09:53 +00:00
|
|
|
>>> solution(6)
|
|
|
|
13
|
|
|
|
>>> solution(1)
|
|
|
|
2
|
|
|
|
>>> solution(3)
|
|
|
|
5
|
|
|
|
>>> solution(20)
|
|
|
|
71
|
|
|
|
>>> solution(50)
|
|
|
|
229
|
|
|
|
>>> solution(100)
|
|
|
|
541
|
2019-07-18 17:05:14 +00:00
|
|
|
>>> solution(3.4)
|
|
|
|
5
|
|
|
|
>>> solution(0)
|
|
|
|
Traceback (most recent call last):
|
|
|
|
...
|
2020-10-25 03:23:16 +00:00
|
|
|
ValueError: Parameter nth must be greater than or equal to one.
|
2019-07-18 17:05:14 +00:00
|
|
|
>>> solution(-17)
|
|
|
|
Traceback (most recent call last):
|
|
|
|
...
|
2020-10-25 03:23:16 +00:00
|
|
|
ValueError: Parameter nth must be greater than or equal to one.
|
2019-07-18 17:05:14 +00:00
|
|
|
>>> solution([])
|
|
|
|
Traceback (most recent call last):
|
|
|
|
...
|
2020-10-25 03:23:16 +00:00
|
|
|
TypeError: Parameter nth must be int or castable to int.
|
2019-07-18 17:05:14 +00:00
|
|
|
>>> solution("asd")
|
|
|
|
Traceback (most recent call last):
|
|
|
|
...
|
2020-10-25 03:23:16 +00:00
|
|
|
TypeError: Parameter nth must be int or castable to int.
|
2019-07-16 23:09:53 +00:00
|
|
|
"""
|
2020-10-25 03:23:16 +00:00
|
|
|
|
2019-07-18 17:05:14 +00:00
|
|
|
try:
|
2020-10-08 05:57:47 +00:00
|
|
|
nth = int(nth)
|
2020-05-22 06:10:11 +00:00
|
|
|
except (TypeError, ValueError):
|
2020-10-25 03:23:16 +00:00
|
|
|
raise TypeError("Parameter nth must be int or castable to int.") from None
|
2020-10-08 05:57:47 +00:00
|
|
|
if nth <= 0:
|
2020-10-25 03:23:16 +00:00
|
|
|
raise ValueError("Parameter nth must be greater than or equal to one.")
|
2021-10-11 16:33:44 +00:00
|
|
|
primes: list[int] = []
|
2019-07-16 23:09:53 +00:00
|
|
|
num = 2
|
2020-10-08 05:57:47 +00:00
|
|
|
while len(primes) < nth:
|
2022-04-08 17:40:45 +00:00
|
|
|
if is_prime(num):
|
2019-07-16 23:09:53 +00:00
|
|
|
primes.append(num)
|
|
|
|
num += 1
|
|
|
|
else:
|
|
|
|
num += 1
|
|
|
|
return primes[len(primes) - 1]
|
|
|
|
|
|
|
|
|
|
|
|
if __name__ == "__main__":
|
2020-10-25 03:23:16 +00:00
|
|
|
print(f"{solution() = }")
|