BigNum.swift 8.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292
  1. //
  2. // BigNum.swift
  3. // Swifter
  4. //
  5. // Copyright © 2016 Damian Kołakowski. All rights reserved.
  6. //
  7. import Foundation
  8. public struct BigNum: Equatable, Comparable, CustomStringConvertible {
  9. //
  10. // Big Integer
  11. //
  12. // TODO/Improvments
  13. // - add init(...) for string literal
  14. // - add powmod operation
  15. // - fix '+' and '-' operators for negative numbers
  16. internal var digits = [UInt8]()
  17. internal var signum = 0
  18. internal init(_ digits: [UInt8], _ signum: Int) {
  19. self.digits.append(contentsOf: digits)
  20. self.signum = signum
  21. }
  22. internal init(_ digits: ArraySlice<UInt8>, _ signum: Int) {
  23. self.digits.append(contentsOf: digits)
  24. self.signum = signum
  25. }
  26. public init(_ text: String) {
  27. if text.unicodeScalars.isEmpty {
  28. self.signum = 0
  29. self.digits = [0]
  30. return
  31. }
  32. if text.unicodeScalars.count == 1 && text.unicodeScalars[text.unicodeScalars.startIndex].value == 48 {
  33. self.signum = 0
  34. self.digits = [0]
  35. return
  36. }
  37. self.signum = 1
  38. for c in text.unicodeScalars.reversed() {
  39. switch c.value {
  40. case 48...57:
  41. self.digits.append(UInt8(c.value-48))
  42. case 45:
  43. self.signum = -1
  44. default: break
  45. }
  46. }
  47. }
  48. public var description: String {
  49. return (self.signum < 0 ? "-" : "") + self.digits.reversed().reduce("") { $0 + String($1) }
  50. }
  51. }
  52. infix operator == { precedence 130 }
  53. public func == (_ left: BigNum, right: BigNum) -> Bool {
  54. return (left.digits == right.digits) && (left.signum == right.signum)
  55. }
  56. infix operator != { precedence 130 }
  57. public func != (_ left: BigNum, right: BigNum) -> Bool {
  58. return !(left.digits == right.digits)
  59. }
  60. infix operator + { associativity left precedence 140 }
  61. public func + (_ left: BigNum, _ right: BigNum) -> BigNum {
  62. if right.signum == 0 {
  63. return left
  64. }
  65. if left.signum == 0 {
  66. return right
  67. }
  68. if left.signum < 0 && right.signum < 0 {
  69. return BigNum((BigNum(left.digits, 1) + BigNum(right.digits, 1)).digits, -1)
  70. }
  71. if left.signum > 0 && right.signum < 0 {
  72. return (left - BigNum(right.digits, 1))
  73. }
  74. if left.signum < 0 && right.signum > 0 {
  75. return BigNum((right - BigNum(left.digits, 1)).digits, -1)
  76. }
  77. var result = [UInt8]()
  78. var carry: UInt8 = 0
  79. for i in 0..<max(left.digits.count, right.digits.count) {
  80. let sum = (((i < left.digits.count) ? left.digits[i] : 0) +
  81. ((i < right.digits.count) ? right.digits[i] : 0)) + carry
  82. if sum >= 10 {
  83. carry = sum / 10
  84. result.append(sum - 10)
  85. } else {
  86. carry = 0
  87. result.append(sum)
  88. }
  89. }
  90. if carry != 0 {
  91. result.append(carry)
  92. }
  93. while let last = result.last, last == 0 { result.removeLast() }
  94. return result.isEmpty ? BigNum(result, 0) : BigNum(result, 1)
  95. }
  96. infix operator - { associativity left precedence 140 }
  97. public func - (_ left: BigNum, _ right: BigNum) -> BigNum {
  98. if right.signum == 0 {
  99. return left
  100. }
  101. if left.signum == 0 {
  102. return BigNum(right.digits, right.signum * (-1))
  103. }
  104. if left == right {
  105. return BigNum([0], 0)
  106. }
  107. if right > left {
  108. return BigNum((right - left).digits, -1)
  109. }
  110. if left.signum < 0 && right.signum > 0 {
  111. return BigNum((BigNum(left.digits, 1) + right).digits, -1)
  112. }
  113. var result = [UInt8]()
  114. var carry: Int = 0
  115. for i in 0..<max(right.digits.count, left.digits.count) {
  116. let sub = (Int((i < left.digits.count) ? left.digits[i] : 0) -
  117. Int((i < right.digits.count) ? right.digits[i] : 0)) + carry
  118. if sub < 0 {
  119. carry = -1
  120. result.append(UInt8(10+sub))
  121. } else {
  122. carry = 0
  123. result.append(UInt8(sub))
  124. }
  125. }
  126. while let last = result.last, last == 0 { result.removeLast() }
  127. return result.isEmpty ? BigNum(result, 0) : BigNum(result, 1)
  128. }
  129. infix operator * { associativity left precedence 150 }
  130. public func * (_ left: BigNum, _ right: BigNum) -> BigNum {
  131. if (left.signum == 0) || (right.signum == 0) {
  132. return BigNum([], 0)
  133. }
  134. var mulResults = Array<[UInt8]>()
  135. for i in 0..<left.digits.count {
  136. var row = [UInt8](repeating: 0, count: i)
  137. var carry: UInt8 = 0
  138. for j in 0..<right.digits.count {
  139. let mul = (left.digits[i] * right.digits[j]) + carry
  140. if mul > 9 {
  141. carry = mul / 10
  142. row.append(mul % 10)
  143. } else {
  144. carry = 0
  145. row.append(mul)
  146. }
  147. }
  148. if carry != 0 {
  149. row.append(carry % 10)
  150. if carry > 9 { row.append(carry / 10) }
  151. }
  152. mulResults.append(row)
  153. }
  154. var sum = mulResults.reduce(BigNum([0], 1)) { $0.partialResult + BigNum($0.1, 1) }
  155. sum.signum = left.signum * right.signum
  156. return sum
  157. }
  158. infix operator / { associativity left precedence 150 }
  159. public func / (_ left: BigNum, _ right: BigNum) -> (quotient: BigNum, reminder: BigNum) {
  160. if left < right {
  161. return (BigNum([0], 0), left)
  162. }
  163. if left.digits == right.digits {
  164. return (BigNum([1], left.signum * right.signum), BigNum([0], 0))
  165. }
  166. var shiftIndex = left.digits.count - right.digits.count
  167. var tmp = BigNum(left.digits[shiftIndex..<left.digits.count], 1)
  168. while tmp < right {
  169. shiftIndex = shiftIndex - 1
  170. tmp = BigNum(left.digits[shiftIndex..<left.digits.count], 1)
  171. }
  172. var quotient = [UInt8]()
  173. var rest = BigNum([0], 0)
  174. while shiftIndex >= 0 {
  175. var i: UInt8 = 1
  176. while (right * BigNum([i], 1)) < tmp {
  177. i = i + 1
  178. }
  179. quotient.append(i-1)
  180. rest = tmp - (right * BigNum([i-1], 1))
  181. shiftIndex = shiftIndex - 1
  182. if shiftIndex < 0 {
  183. break
  184. }
  185. tmp = BigNum([left.digits[shiftIndex]] + rest.digits, 1)
  186. }
  187. return (BigNum(quotient.reversed(), left.signum * right.signum), rest)
  188. }
  189. infix operator % { associativity left precedence 150 }
  190. public func % (_ left: BigNum, _ right: BigNum) -> BigNum {
  191. return (left / right).reminder
  192. }
  193. infix operator > { precedence 130 }
  194. public func > (_ left: BigNum, _ right: BigNum) -> Bool {
  195. if left.signum != right.signum {
  196. return left.signum > right.signum
  197. }
  198. if left.digits.count != right.digits.count {
  199. if left.signum < 0 {
  200. return left.digits.count < right.digits.count
  201. } else {
  202. return left.digits.count > right.digits.count
  203. }
  204. }
  205. for i in (0..<left.digits.count).reversed() {
  206. if left.digits[i] != right.digits[i] {
  207. if left.signum < 0 {
  208. return left.digits[i] < right.digits[i]
  209. } else {
  210. return left.digits[i] > right.digits[i]
  211. }
  212. }
  213. }
  214. return false
  215. }
  216. infix operator < { precedence 130 }
  217. public func < (_ left: BigNum, _ right: BigNum) -> Bool {
  218. if left.signum != right.signum {
  219. return left.signum < right.signum
  220. }
  221. if right.digits.count != left.digits.count {
  222. if left.signum < 0 {
  223. return left.digits.count > right.digits.count
  224. } else {
  225. return left.digits.count < right.digits.count
  226. }
  227. }
  228. for i in (0..<left.digits.count).reversed() {
  229. if left.digits[i] != right.digits[i] {
  230. if left.signum < 0 {
  231. return right.digits[i] < left.digits[i]
  232. } else {
  233. return right.digits[i] > left.digits[i]
  234. }
  235. }
  236. }
  237. return false
  238. }
  239. infix operator >= { precedence 130 }
  240. public func >= (_ left: BigNum, _ right: BigNum) -> Bool {
  241. return (left == right) || (left > right)
  242. }
  243. infix operator <= { precedence 130 }
  244. public func <= (_ left: BigNum, _ right: BigNum) -> Bool {
  245. return (left == right) || (left < right)
  246. }