mirror of
https://github.com/TheAlgorithms/Python.git
synced 2025-02-28 20:08:40 +00:00
115 lines
2.9 KiB
Python
115 lines
2.9 KiB
Python
"""
|
|
Calculate the nth Proth number
|
|
Source:
|
|
https://handwiki.org/wiki/Proth_number
|
|
"""
|
|
|
|
import math
|
|
|
|
|
|
def proth(number: int) -> int:
|
|
"""
|
|
:param number: nth number to calculate in the sequence
|
|
:return: the nth number in Proth number
|
|
Note: indexing starts at 1 i.e. proth(1) gives the first Proth number of 3
|
|
>>> proth(6)
|
|
25
|
|
>>> proth(0)
|
|
Traceback (most recent call last):
|
|
...
|
|
ValueError: Input value of [number=0] must be > 0
|
|
>>> proth(-1)
|
|
Traceback (most recent call last):
|
|
...
|
|
ValueError: Input value of [number=-1] must be > 0
|
|
>>> proth(6.0)
|
|
Traceback (most recent call last):
|
|
...
|
|
TypeError: Input value of [number=6.0] must be an integer
|
|
"""
|
|
|
|
if not isinstance(number, int):
|
|
msg = f"Input value of [number={number}] must be an integer"
|
|
raise TypeError(msg)
|
|
|
|
if number < 1:
|
|
msg = f"Input value of [number={number}] must be > 0"
|
|
raise ValueError(msg)
|
|
elif number == 1:
|
|
return 3
|
|
elif number == 2:
|
|
return 5
|
|
else:
|
|
"""
|
|
+1 for binary starting at 0 i.e. 2^0, 2^1, etc.
|
|
+1 to start the sequence at the 3rd Proth number
|
|
Hence, we have a +2 in the below statement
|
|
"""
|
|
block_index = int(math.log(number // 3, 2)) + 2
|
|
|
|
proth_list = [3, 5]
|
|
proth_index = 2
|
|
increment = 3
|
|
for block in range(1, block_index):
|
|
for _ in range(increment):
|
|
proth_list.append(2 ** (block + 1) + proth_list[proth_index - 1])
|
|
proth_index += 1
|
|
increment *= 2
|
|
|
|
return proth_list[number - 1]
|
|
|
|
|
|
def isprothnumber(number: int) -> bool :
|
|
"""
|
|
:param number: nth number to calculate in the sequence
|
|
:return: true if number is a Proth number, false etherwise
|
|
>>> proth(5)
|
|
true
|
|
>>> proth(34)
|
|
false
|
|
>>> proth(-1)
|
|
Traceback (most recent call last):
|
|
...
|
|
ValueError: Input value of [number=-1] must be > 0
|
|
>>> proth(6.0)
|
|
Traceback (most recent call last):
|
|
...
|
|
TypeError: Input value of [number=6.0] must be an integer
|
|
"""
|
|
if number < 1:
|
|
msg = f"Input value of [number={number}] must be > 0"
|
|
raise ValueError(msg)
|
|
|
|
num = number
|
|
num -= 1
|
|
n = 0
|
|
while num%2 == 0 :
|
|
N = N/2
|
|
n += 1
|
|
return num < (2**n)
|
|
|
|
|
|
|
|
|
|
if __name__ == "__main__":
|
|
import doctest
|
|
|
|
doctest.testmod()
|
|
|
|
for number in range(11):
|
|
value = 0
|
|
try:
|
|
value = proth(number)
|
|
except ValueError:
|
|
print(f"ValueError: there is no {number}th Proth number")
|
|
continue
|
|
|
|
print(f"The {number}th Proth number: {value}")
|
|
|
|
listexe = [3, 5, 9, 13, 49, 57, 193, 241, 163, 201]
|
|
|
|
for number in listexe :
|
|
if isprothnumber(number):
|
|
print(f"{number} is a Proth number")
|
|
else :
|
|
print(f"{number} is not a Proth number") |