test_pow.py 6.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171
  1. import math
  2. import unittest
  3. class PowTest(unittest.TestCase):
  4. def powtest(self, type):
  5. if type != float:
  6. for i in range(-1000, 1000):
  7. self.assertEqual(pow(type(i), 0), 1)
  8. self.assertEqual(pow(type(i), 1), type(i))
  9. self.assertEqual(pow(type(0), 1), type(0))
  10. self.assertEqual(pow(type(1), 1), type(1))
  11. for i in range(-100, 100):
  12. self.assertEqual(pow(type(i), 3), i*i*i)
  13. pow2 = 1
  14. for i in range(0, 31):
  15. self.assertEqual(pow(2, i), pow2)
  16. if i != 30 : pow2 = pow2*2
  17. for othertype in (int,):
  18. for i in list(range(-10, 0)) + list(range(1, 10)):
  19. ii = type(i)
  20. for j in range(1, 11):
  21. jj = -othertype(j)
  22. pow(ii, jj)
  23. for othertype in int, float:
  24. for i in range(1, 100):
  25. zero = type(0)
  26. exp = -othertype(i/10.0)
  27. if exp == 0:
  28. continue
  29. self.assertRaises(ZeroDivisionError, pow, zero, exp)
  30. il, ih = -20, 20
  31. jl, jh = -5, 5
  32. kl, kh = -10, 10
  33. asseq = self.assertEqual
  34. if type == float:
  35. il = 1
  36. asseq = self.assertAlmostEqual
  37. elif type == int:
  38. jl = 0
  39. elif type == int:
  40. jl, jh = 0, 15
  41. for i in range(il, ih+1):
  42. for j in range(jl, jh+1):
  43. for k in range(kl, kh+1):
  44. if k != 0:
  45. if type == float or j < 0:
  46. self.assertRaises(TypeError, pow, type(i), j, k)
  47. continue
  48. asseq(
  49. pow(type(i),j,k),
  50. pow(type(i),j)% type(k)
  51. )
  52. def test_powint(self):
  53. self.powtest(int)
  54. def test_powfloat(self):
  55. self.powtest(float)
  56. def test_other(self):
  57. # Other tests-- not very systematic
  58. self.assertEqual(pow(3,3) % 8, pow(3,3,8))
  59. self.assertEqual(pow(3,3) % -8, pow(3,3,-8))
  60. self.assertEqual(pow(3,2) % -2, pow(3,2,-2))
  61. self.assertEqual(pow(-3,3) % 8, pow(-3,3,8))
  62. self.assertEqual(pow(-3,3) % -8, pow(-3,3,-8))
  63. self.assertEqual(pow(5,2) % -8, pow(5,2,-8))
  64. self.assertEqual(pow(3,3) % 8, pow(3,3,8))
  65. self.assertEqual(pow(3,3) % -8, pow(3,3,-8))
  66. self.assertEqual(pow(3,2) % -2, pow(3,2,-2))
  67. self.assertEqual(pow(-3,3) % 8, pow(-3,3,8))
  68. self.assertEqual(pow(-3,3) % -8, pow(-3,3,-8))
  69. self.assertEqual(pow(5,2) % -8, pow(5,2,-8))
  70. for i in range(-10, 11):
  71. for j in range(0, 6):
  72. for k in range(-7, 11):
  73. if j >= 0 and k != 0:
  74. self.assertEqual(
  75. pow(i,j) % k,
  76. pow(i,j,k)
  77. )
  78. if j >= 0 and k != 0:
  79. self.assertEqual(
  80. pow(int(i),j) % k,
  81. pow(int(i),j,k)
  82. )
  83. def test_big_exp(self):
  84. import random
  85. self.assertEqual(pow(2, 50000), 1 << 50000)
  86. # Randomized modular tests, checking the identities
  87. # a**(b1 + b2) == a**b1 * a**b2
  88. # a**(b1 * b2) == (a**b1)**b2
  89. prime = 1000000000039 # for speed, relatively small prime modulus
  90. for i in range(10):
  91. a = random.randrange(1000, 1000000)
  92. bpower = random.randrange(1000, 50000)
  93. b = random.randrange(1 << (bpower - 1), 1 << bpower)
  94. b1 = random.randrange(1, b)
  95. b2 = b - b1
  96. got1 = pow(a, b, prime)
  97. got2 = pow(a, b1, prime) * pow(a, b2, prime) % prime
  98. if got1 != got2:
  99. self.fail(f"{a=:x} {b1=:x} {b2=:x} {got1=:x} {got2=:x}")
  100. got3 = pow(a, b1 * b2, prime)
  101. got4 = pow(pow(a, b1, prime), b2, prime)
  102. if got3 != got4:
  103. self.fail(f"{a=:x} {b1=:x} {b2=:x} {got3=:x} {got4=:x}")
  104. def test_bug643260(self):
  105. class TestRpow:
  106. def __rpow__(self, other):
  107. return None
  108. None ** TestRpow() # Won't fail when __rpow__ invoked. SF bug #643260.
  109. def test_bug705231(self):
  110. # -1.0 raised to an integer should never blow up. It did if the
  111. # platform pow() was buggy, and Python didn't worm around it.
  112. eq = self.assertEqual
  113. a = -1.0
  114. # The next two tests can still fail if the platform floor()
  115. # function doesn't treat all large inputs as integers
  116. # test_math should also fail if that is happening
  117. eq(pow(a, 1.23e167), 1.0)
  118. eq(pow(a, -1.23e167), 1.0)
  119. for b in range(-10, 11):
  120. eq(pow(a, float(b)), b & 1 and -1.0 or 1.0)
  121. for n in range(0, 100):
  122. fiveto = float(5 ** n)
  123. # For small n, fiveto will be odd. Eventually we run out of
  124. # mantissa bits, though, and thereafer fiveto will be even.
  125. expected = fiveto % 2.0 and -1.0 or 1.0
  126. eq(pow(a, fiveto), expected)
  127. eq(pow(a, -fiveto), expected)
  128. eq(expected, 1.0) # else we didn't push fiveto to evenness
  129. def test_negative_exponent(self):
  130. for a in range(-50, 50):
  131. for m in range(-50, 50):
  132. with self.subTest(a=a, m=m):
  133. if m != 0 and math.gcd(a, m) == 1:
  134. # Exponent -1 should give an inverse, with the
  135. # same sign as m.
  136. inv = pow(a, -1, m)
  137. self.assertEqual(inv, inv % m)
  138. self.assertEqual((inv * a - 1) % m, 0)
  139. # Larger exponents
  140. self.assertEqual(pow(a, -2, m), pow(inv, 2, m))
  141. self.assertEqual(pow(a, -3, m), pow(inv, 3, m))
  142. self.assertEqual(pow(a, -1001, m), pow(inv, 1001, m))
  143. else:
  144. with self.assertRaises(ValueError):
  145. pow(a, -1, m)
  146. with self.assertRaises(ValueError):
  147. pow(a, -2, m)
  148. with self.assertRaises(ValueError):
  149. pow(a, -1001, m)
  150. if __name__ == "__main__":
  151. unittest.main()