Python/maths/basic_maths.py

81 lines
1.7 KiB
Python
Raw Normal View History

"""Implementation of Basic Math in Python."""
2018-10-19 12:48:28 +00:00
import math
def prime_factors(n: int) -> list:
"""Find Prime Factors.
>>> prime_factors(100)
[2, 2, 5, 5]
"""
2018-10-19 12:48:28 +00:00
pf = []
while n % 2 == 0:
pf.append(2)
n = int(n / 2)
for i in range(3, int(math.sqrt(n)) + 1, 2):
2018-10-19 12:48:28 +00:00
while n % i == 0:
pf.append(i)
n = int(n / i)
if n > 2:
pf.append(n)
return pf
def number_of_divisors(n: int) -> int:
"""Calculate Number of Divisors of an Integer.
>>> number_of_divisors(100)
9
"""
2018-10-19 12:48:28 +00:00
div = 1
temp = 1
while n % 2 == 0:
temp += 1
n = int(n / 2)
div *= temp
for i in range(3, int(math.sqrt(n)) + 1, 2):
2018-10-19 12:48:28 +00:00
temp = 1
while n % i == 0:
temp += 1
n = int(n / i)
div *= temp
2018-10-19 12:48:28 +00:00
return div
def sum_of_divisors(n: int) -> int:
"""Calculate Sum of Divisors.
>>> sum_of_divisors(100)
217
"""
2018-10-19 12:48:28 +00:00
s = 1
temp = 1
while n % 2 == 0:
temp += 1
n = int(n / 2)
if temp > 1:
2019-10-05 05:14:13 +00:00
s *= (2 ** temp - 1) / (2 - 1)
for i in range(3, int(math.sqrt(n)) + 1, 2):
2018-10-19 12:48:28 +00:00
temp = 1
while n % i == 0:
temp += 1
n = int(n / i)
if temp > 1:
2019-10-05 05:14:13 +00:00
s *= (i ** temp - 1) / (i - 1)
return int(s)
2018-10-19 12:48:28 +00:00
def euler_phi(n: int) -> int:
"""Calculate Euler's Phi Function.
>>> euler_phi(100)
40
"""
2018-10-19 12:48:28 +00:00
s = n
for x in set(prime_factors(n)):
s *= (x - 1) / x
return int(s)
2018-10-19 12:48:28 +00:00
if __name__ == "__main__":
print(prime_factors(100))
print(number_of_divisors(100))
print(sum_of_divisors(100))
print(euler_phi(100))