BigNum.swift 8.0 KB

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