div.rs 15 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462
  1. use ibig::{
  2. ibig,
  3. ops::{DivEuclid, DivRem, DivRemEuclid, RemEuclid},
  4. ubig, IBig,
  5. };
  6. #[test]
  7. fn test_div_rem_ubig() {
  8. let test_cases = [
  9. (ubig!(331), ubig!(10), ubig!(33), ubig!(1)),
  10. (ubig!(17),
  11. ubig!(_0x987987123984798abbcc213789723948792138479837492837498cc),
  12. ubig!(0),
  13. ubig!(17)
  14. ),
  15. (ubig!(_0x987987123984798abbcc213789723948792138479837492837498cc),
  16. ubig!(_0x1234),
  17. ubig!(_0x86054c502f0a4e43e2d0de91f1029d251ce67bbdb88dc3edbb40),
  18. ubig!(_0xfcc)
  19. ),
  20. (ubig!(_0x987987123984798abbcc213789723948792138479837492837498cc),
  21. ubig!(_0x1000),
  22. ubig!(_0x987987123984798abbcc21378972394879213847983749283749),
  23. ubig!(_0x8cc)
  24. ),
  25. (ubig!(_0x987987123984798abbcc213789723948792138479837492837498cc),
  26. ubig!(_0xf234567812345678),
  27. ubig!(_0xa128cfb49d0d746cc0295e163c343aafbffbfa8),
  28. ubig!(_0x9068d997bb10520c)
  29. ),
  30. (ubig!(_0x987987123984798abbcc213789723948792138479837492837498ce),
  31. ubig!(_0x987987123984798abbcc213789723948792138479837492837498cc),
  32. ubig!(1),
  33. ubig!(2),
  34. ),
  35. // Special case for division (64-bit words): top 2 / top 1 overflows.
  36. (
  37. ubig!(_0xffffffffffffffffffffffffffffffff00000000000000000000000000000000),
  38. ubig!(_0xffffffffffffffffffffffffffffffff0000000000000001),
  39. ubig!(_0xffffffffffffffff),
  40. ubig!(_0xfffffffffffffffffffffffffffffffe0000000000000001),
  41. ),
  42. // Random 500-bit by random 250-bit.
  43. (
  44. ubig!(_0x2b8f1bb75f1ca5bf3400549a663d503d298da7f53942cd3c5c6a1bc50598d091e8ca30896413783e9b001572e28808c4dc9598bdd17ef3ce35b40e0368b60),
  45. ubig!(_0x3e880309f5e48d145337aae47694a74f2860db8e49665f03978f1b11665dc80),
  46. ubig!(_0xb254145d6f736c22ed5fca6a41f4c883a59fc32c638710758bb50fa532b31f),
  47. ubig!(_0x195afda8e35e347c65ed01409c73d1c820ed78a87e83cf6cfdad1a25fb357e0),
  48. ),
  49. (
  50. ubig!(_0x3e880309f5e48d145337aae47694a74f2860db8e49665f03978f1b11665dc80),
  51. ubig!(_0x2b8f1bb75f1ca5bf3400549a663d503d298da7f53942cd3c5c6a1bc50598d091e8ca30896413783e9b001572e28808c4dc9598bdd17ef3ce35b40e0368b60),
  52. ubig!(0),
  53. ubig!(_0x3e880309f5e48d145337aae47694a74f2860db8e49665f03978f1b11665dc80),
  54. ),
  55. // 3^300 - 1 by 3^150
  56. (
  57. ubig!(_0xb39cfff485a5dbf4d6aae030b91bfb0ec6bba389cd8d7f85bba3985c19c5e24e40c543a123c6e028a873e9e3874e1b4623a44be39b34e67dc5c2670),
  58. ubig!(_0x359ba2b98ca11d6864a331b45ae7114c01ffbdcf60cc16e692fb63c6e219),
  59. ubig!(_0x359ba2b98ca11d6864a331b45ae7114c01ffbdcf60cc16e692fb63c6e218),
  60. ubig!(_0x359ba2b98ca11d6864a331b45ae7114c01ffbdcf60cc16e692fb63c6e218),
  61. ),
  62. // 7^70-1 by 7^35
  63. (
  64. ubig!(_0x16dc8782276b9f7addf9768f33c8007ce903866a4546c1a190),
  65. ubig!(_0x4c8077a58a0a8cb7c24960e57),
  66. ubig!(_0x4c8077a58a0a8cb7c24960e56),
  67. ubig!(_0x4c8077a58a0a8cb7c24960e56),
  68. ),
  69. // 2^20480-1 by 2^5120-1
  70. (
  71. (ubig!(1) << 20480) - ubig!(1),
  72. (ubig!(1) << 5120) - ubig!(1),
  73. ubig!(1) + (ubig!(1) << 5120) + (ubig!(1) << 10240) + (ubig!(1) << 15360),
  74. ubig!(0)
  75. ),
  76. // 2^20480-1 by 2^15360-1
  77. (
  78. (ubig!(1) << 20480) - ubig!(1),
  79. (ubig!(1) << 15360) - ubig!(1),
  80. ubig!(1) << 5120,
  81. (ubig!(1) << 5120) - ubig!(1)
  82. ),
  83. // 2^19000-1 by 2^5000-1
  84. (
  85. (ubig!(1) << 19000) - ubig!(1),
  86. (ubig!(1) << 5000) - ubig!(1),
  87. (ubig!(1) << 14000) + (ubig!(1) << 9000) + (ubig!(1) << 4000),
  88. (ubig!(1) << 4000) - ubig!(1)
  89. ),
  90. ];
  91. for (a, b, q, r) in &test_cases {
  92. let qr = (q.clone(), r.clone());
  93. assert_eq!(a / b, *q);
  94. assert_eq!(a.clone() / b, *q);
  95. assert_eq!(a / b.clone(), *q);
  96. assert_eq!(a.clone() / b.clone(), *q);
  97. let mut x = a.clone();
  98. x /= b;
  99. assert_eq!(x, *q);
  100. let mut x = a.clone();
  101. x /= b.clone();
  102. assert_eq!(x, *q);
  103. assert_eq!(a % b, *r);
  104. assert_eq!(a.clone() % b, *r);
  105. assert_eq!(a % b.clone(), *r);
  106. assert_eq!(a.clone() % b.clone(), *r);
  107. let mut x = a.clone();
  108. x %= b;
  109. assert_eq!(x, *r);
  110. let mut x = a.clone();
  111. x %= b.clone();
  112. assert_eq!(x, *r);
  113. assert_eq!(a.div_rem(b), qr);
  114. assert_eq!(a.clone().div_rem(b), qr);
  115. assert_eq!(a.div_rem(b.clone()), qr);
  116. assert_eq!(a.clone().div_rem(b.clone()), qr);
  117. assert_eq!(a.div_euclid(b), *q);
  118. assert_eq!(a.clone().div_euclid(b), *q);
  119. assert_eq!(a.div_euclid(b.clone()), *q);
  120. assert_eq!(a.clone().div_euclid(b.clone()), *q);
  121. assert_eq!(a.rem_euclid(b), *r);
  122. assert_eq!(a.clone().rem_euclid(b), *r);
  123. assert_eq!(a.rem_euclid(b.clone()), *r);
  124. assert_eq!(a.clone().rem_euclid(b.clone()), *r);
  125. assert_eq!(a.div_rem_euclid(b), qr);
  126. assert_eq!(a.clone().div_rem_euclid(b), qr);
  127. assert_eq!(a.div_rem_euclid(b.clone()), qr);
  128. assert_eq!(a.clone().div_rem_euclid(b.clone()), qr);
  129. }
  130. }
  131. #[test]
  132. #[should_panic]
  133. fn test_divide_by_0_ubig() {
  134. let _ = ubig!(5) / ubig!(0);
  135. }
  136. #[test]
  137. fn test_div_rem_ibig() {
  138. for a in -20i8..=20i8 {
  139. for b in -20i8..=20i8 {
  140. if b == 0 {
  141. continue;
  142. }
  143. let a_big: IBig = a.into();
  144. let b_big: IBig = b.into();
  145. let q: IBig = (a / b).into();
  146. let r: IBig = (a % b).into();
  147. let qr = (q.clone(), r.clone());
  148. assert_eq!(a_big.clone() / b_big.clone(), q);
  149. assert_eq!(&a_big / b_big.clone(), q);
  150. assert_eq!(a_big.clone() / &b_big, q);
  151. assert_eq!(&a_big / &b_big, q);
  152. let mut x = a_big.clone();
  153. x /= b_big.clone();
  154. assert_eq!(x, q);
  155. let mut x = a_big.clone();
  156. x /= &b_big;
  157. assert_eq!(x, q);
  158. assert_eq!(a_big.clone() % b_big.clone(), r);
  159. assert_eq!(&a_big % b_big.clone(), r);
  160. assert_eq!(a_big.clone() % &b_big, r);
  161. assert_eq!(&a_big % &b_big, r);
  162. let mut x = a_big.clone();
  163. x %= b_big.clone();
  164. assert_eq!(x, r);
  165. let mut x = a_big.clone();
  166. x %= &b_big;
  167. assert_eq!(x, r);
  168. assert_eq!(a_big.clone().div_rem(b_big.clone()), qr);
  169. assert_eq!((&a_big).div_rem(b_big.clone()), qr);
  170. assert_eq!(a_big.clone().div_rem(&b_big), qr);
  171. assert_eq!((&a_big).div_rem(&b_big), qr);
  172. }
  173. }
  174. }
  175. #[test]
  176. fn test_div_rem_euclid_ibig() {
  177. for a in -20i8..=20i8 {
  178. for b in -20i8..=20i8 {
  179. if b == 0 {
  180. continue;
  181. }
  182. let a_big: IBig = a.into();
  183. let b_big: IBig = b.into();
  184. let q: IBig = a.div_euclid(b).into();
  185. let r: IBig = a.rem_euclid(b).into();
  186. let qr = (q.clone(), r.clone());
  187. assert_eq!(a_big.clone().div_euclid(b_big.clone()), q);
  188. assert_eq!((&a_big).div_euclid(b_big.clone()), q);
  189. assert_eq!(a_big.clone().div_euclid(&b_big), q);
  190. assert_eq!((&a_big).div_euclid(&b_big), q);
  191. assert_eq!(a_big.clone().rem_euclid(b_big.clone()), r);
  192. assert_eq!((&a_big).rem_euclid(b_big.clone()), r);
  193. assert_eq!(a_big.clone().rem_euclid(&b_big), r);
  194. assert_eq!((&a_big).rem_euclid(&b_big), r);
  195. assert_eq!(a_big.clone().div_rem_euclid(b_big.clone()), qr);
  196. assert_eq!((&a_big).div_rem_euclid(b_big.clone()), qr);
  197. assert_eq!(a_big.clone().div_rem_euclid(&b_big), qr);
  198. assert_eq!((&a_big).div_rem_euclid(&b_big), qr);
  199. }
  200. }
  201. }
  202. #[test]
  203. #[should_panic]
  204. fn test_divide_by_0_ibig() {
  205. let _ = ibig!(5) / ibig!(0);
  206. }
  207. #[test]
  208. #[allow(clippy::op_ref)]
  209. fn test_div_rem_ubig_unsigned() {
  210. assert_eq!(ubig!(23) / 10u8, ubig!(2));
  211. assert_eq!(ubig!(23) / &10u8, ubig!(2));
  212. assert_eq!(&ubig!(23) / 10u8, ubig!(2));
  213. assert_eq!(&ubig!(23) / &10u8, ubig!(2));
  214. let mut x = ubig!(23);
  215. x /= 10u8;
  216. assert_eq!(x, ubig!(2));
  217. let mut x = ubig!(23);
  218. x /= &10u8;
  219. assert_eq!(x, ubig!(2));
  220. assert_eq!(ubig!(23) % 10u8, 3u8);
  221. assert_eq!(ubig!(23) % &10u8, 3u8);
  222. assert_eq!(&ubig!(23) % 10u8, 3u8);
  223. assert_eq!(&ubig!(23) % &10u8, 3u8);
  224. let mut x = ubig!(23);
  225. x %= 10u8;
  226. assert_eq!(x, ubig!(3));
  227. let mut x = ubig!(23);
  228. x %= &10u8;
  229. assert_eq!(x, ubig!(3));
  230. assert_eq!(ubig!(23).div_rem(10u8), (ubig!(2), 3u8));
  231. assert_eq!(ubig!(23).div_rem(&10u8), (ubig!(2), 3u8));
  232. assert_eq!((&ubig!(23)).div_rem(10u8), (ubig!(2), 3u8));
  233. assert_eq!((&ubig!(23)).div_rem(&10u8), (ubig!(2), 3u8));
  234. }
  235. #[test]
  236. fn test_div_rem_euclid_ubig_unsigned() {
  237. assert_eq!(ubig!(23).div_euclid(10u8), ubig!(2));
  238. assert_eq!(ubig!(23).div_euclid(&10u8), ubig!(2));
  239. assert_eq!((&ubig!(23)).div_euclid(10u8), ubig!(2));
  240. assert_eq!((&ubig!(23)).div_euclid(&10u8), ubig!(2));
  241. assert_eq!(ubig!(23).rem_euclid(10u8), 3u8);
  242. assert_eq!(ubig!(23).rem_euclid(&10u8), 3u8);
  243. assert_eq!((&ubig!(23)).rem_euclid(10u8), 3u8);
  244. assert_eq!((&ubig!(23)).rem_euclid(&10u8), 3u8);
  245. assert_eq!(ubig!(23).div_rem_euclid(10u8), (ubig!(2), 3u8));
  246. assert_eq!(ubig!(23).div_rem_euclid(&10u8), (ubig!(2), 3u8));
  247. assert_eq!((&ubig!(23)).div_rem_euclid(10u8), (ubig!(2), 3u8));
  248. assert_eq!((&ubig!(23)).div_rem_euclid(&10u8), (ubig!(2), 3u8));
  249. }
  250. #[test]
  251. #[allow(clippy::op_ref)]
  252. fn test_div_rem_ubig_signed() {
  253. assert_eq!(ubig!(23) / 10, ubig!(2));
  254. assert_eq!(ubig!(23) / &10, ubig!(2));
  255. assert_eq!(&ubig!(23) / 10, ubig!(2));
  256. assert_eq!(&ubig!(23) / &10, ubig!(2));
  257. let mut x = ubig!(23);
  258. x /= 10;
  259. assert_eq!(x, ubig!(2));
  260. let mut x = ubig!(23);
  261. x /= &10;
  262. assert_eq!(x, ubig!(2));
  263. assert_eq!(ubig!(23) % 10, 3);
  264. assert_eq!(ubig!(23) % &10, 3);
  265. assert_eq!(&ubig!(23) % 10, 3);
  266. assert_eq!(&ubig!(23) % &10, 3);
  267. let mut x = ubig!(23);
  268. x %= 10;
  269. assert_eq!(x, ubig!(3));
  270. let mut x = ubig!(23);
  271. x %= &10;
  272. assert_eq!(x, ubig!(3));
  273. assert_eq!(ubig!(23).div_rem(10), (ubig!(2), 3));
  274. assert_eq!(ubig!(23).div_rem(&10), (ubig!(2), 3));
  275. assert_eq!((&ubig!(23)).div_rem(10), (ubig!(2), 3));
  276. assert_eq!((&ubig!(23)).div_rem(&10), (ubig!(2), 3));
  277. assert_eq!(ubig!(3) / (-100), ubig!(0));
  278. assert_eq!(ubig!(21) % (-10), 1);
  279. assert_eq!(ubig!(1).div_rem(-10), (ubig!(0), 1));
  280. }
  281. #[test]
  282. #[should_panic]
  283. fn test_div_ubig_signed_overflow() {
  284. let _ = ubig!(100) / (-10);
  285. }
  286. #[test]
  287. #[should_panic]
  288. fn test_div_rem_ubig_signed_overflow() {
  289. let _ = ubig!(101).div_rem(-10);
  290. }
  291. #[test]
  292. fn test_div_rem_euclid_ubig_signed() {
  293. assert_eq!(ubig!(23).div_euclid(10), ubig!(2));
  294. assert_eq!(ubig!(23).div_euclid(&10), ubig!(2));
  295. assert_eq!((&ubig!(23)).div_euclid(10), ubig!(2));
  296. assert_eq!((&ubig!(23)).div_euclid(&10), ubig!(2));
  297. assert_eq!(ubig!(23).rem_euclid(10), 3);
  298. assert_eq!(ubig!(23).rem_euclid(&10), 3);
  299. assert_eq!((&ubig!(23)).rem_euclid(10), 3);
  300. assert_eq!((&ubig!(23)).rem_euclid(&10), 3);
  301. assert_eq!(ubig!(23).div_rem_euclid(10), (ubig!(2), 3));
  302. assert_eq!(ubig!(23).div_rem_euclid(&10), (ubig!(2), 3));
  303. assert_eq!((&ubig!(23)).div_rem_euclid(10), (ubig!(2), 3));
  304. assert_eq!((&ubig!(23)).div_rem_euclid(&10), (ubig!(2), 3));
  305. assert_eq!(ubig!(1).div_euclid(-10), ubig!(0));
  306. assert_eq!(ubig!(21).rem_euclid(-10), 1);
  307. assert_eq!(ubig!(1).div_rem_euclid(-10), (ubig!(0), 1));
  308. }
  309. #[test]
  310. #[should_panic]
  311. fn test_div_euclid_ubig_signed_overflow() {
  312. let _ = ubig!(11).div_euclid(-10);
  313. }
  314. #[test]
  315. #[should_panic]
  316. fn test_div_rem_euclid_ubig_signed_overflow() {
  317. let _ = ubig!(11).div_rem_euclid(-10);
  318. }
  319. #[test]
  320. #[allow(clippy::op_ref)]
  321. fn test_div_rem_ibig_unsigned() {
  322. assert_eq!(ibig!(-23) / 10u8, ibig!(-2));
  323. assert_eq!(ibig!(-23) / &10u8, ibig!(-2));
  324. assert_eq!(&ibig!(-23) / 10u8, ibig!(-2));
  325. assert_eq!(&ibig!(-23) / &10u8, ibig!(-2));
  326. let mut x = ibig!(-23);
  327. x /= 10u8;
  328. assert_eq!(x, ibig!(-2));
  329. let mut x = ibig!(-23);
  330. x /= &10u8;
  331. assert_eq!(x, ibig!(-2));
  332. assert_eq!(ibig!(-23) % 10u8, ibig!(-3));
  333. assert_eq!(ibig!(-23) % &10u8, ibig!(-3));
  334. assert_eq!(&ibig!(-23) % 10u8, ibig!(-3));
  335. assert_eq!(&ibig!(-23) % &10u8, ibig!(-3));
  336. let mut x = ibig!(-23);
  337. x %= 10u8;
  338. assert_eq!(x, ibig!(-3));
  339. let mut x = ibig!(-23);
  340. x %= &10u8;
  341. assert_eq!(x, ibig!(-3));
  342. assert_eq!(ibig!(-23).div_rem(10u8), (ibig!(-2), ibig!(-3)));
  343. assert_eq!(ibig!(-23).div_rem(&10u8), (ibig!(-2), ibig!(-3)));
  344. assert_eq!((&ibig!(-23)).div_rem(10u8), (ibig!(-2), ibig!(-3)));
  345. assert_eq!((&ibig!(-23)).div_rem(&10u8), (ibig!(-2), ibig!(-3)));
  346. }
  347. #[test]
  348. fn test_div_rem_euclid_ibig_unsigned() {
  349. assert_eq!(ibig!(-23).div_euclid(10u8), ibig!(-3));
  350. assert_eq!(ibig!(-23).div_euclid(&10u8), ibig!(-3));
  351. assert_eq!((&ibig!(-23)).div_euclid(10u8), ibig!(-3));
  352. assert_eq!((&ibig!(-23)).div_euclid(&10u8), ibig!(-3));
  353. assert_eq!(ibig!(-23).rem_euclid(10u8), 7u8);
  354. assert_eq!(ibig!(-23).rem_euclid(&10u8), 7u8);
  355. assert_eq!((&ibig!(-23)).rem_euclid(10u8), 7u8);
  356. assert_eq!((&ibig!(-23)).rem_euclid(&10u8), 7u8);
  357. assert_eq!(ibig!(-23).div_rem_euclid(10u8), (ibig!(-3), 7u8));
  358. assert_eq!(ibig!(-23).div_rem_euclid(&10u8), (ibig!(-3), 7u8));
  359. assert_eq!((&ibig!(-23)).div_rem_euclid(10u8), (ibig!(-3), 7u8));
  360. assert_eq!((&ibig!(-23)).div_rem_euclid(&10u8), (ibig!(-3), 7u8));
  361. }
  362. #[test]
  363. #[allow(clippy::op_ref)]
  364. fn test_div_rem_ibig_signed() {
  365. assert_eq!(ibig!(-23) / (-10), ibig!(2));
  366. assert_eq!(ibig!(-23) / &(-10), ibig!(2));
  367. assert_eq!(&ibig!(-23) / (-10), ibig!(2));
  368. assert_eq!(&ibig!(-23) / &(-10), ibig!(2));
  369. let mut x = ibig!(-23);
  370. x /= -10;
  371. assert_eq!(x, ibig!(2));
  372. let mut x = ibig!(-23);
  373. x /= &(-10);
  374. assert_eq!(x, ibig!(2));
  375. assert_eq!(ibig!(-23) % (-10), -3);
  376. assert_eq!(ibig!(-23) % &(-10), -3);
  377. assert_eq!(&ibig!(-23) % (-10), -3);
  378. assert_eq!(&ibig!(-23) % &(-10), -3);
  379. let mut x = ibig!(-23);
  380. x %= -10;
  381. assert_eq!(x, ibig!(-3));
  382. let mut x = ibig!(-23);
  383. x %= &(-10);
  384. assert_eq!(x, ibig!(-3));
  385. assert_eq!(ibig!(-23).div_rem(-10), (ibig!(2), -3));
  386. assert_eq!(ibig!(-23).div_rem(&(-10)), (ibig!(2), -3));
  387. assert_eq!((&ibig!(-23)).div_rem(-10), (ibig!(2), -3));
  388. assert_eq!((&ibig!(-23)).div_rem(&(-10)), (ibig!(2), -3));
  389. }
  390. #[test]
  391. fn test_div_rem_euclid_ibig_signed() {
  392. assert_eq!(ibig!(-23).div_euclid(-10), ibig!(3));
  393. assert_eq!(ibig!(-23).div_euclid(&(-10)), ibig!(3));
  394. assert_eq!((&ibig!(-23)).div_euclid(-10), ibig!(3));
  395. assert_eq!((&ibig!(-23)).div_euclid(&(-10)), ibig!(3));
  396. assert_eq!(ibig!(-23).rem_euclid(-10), 7);
  397. assert_eq!(ibig!(-23).rem_euclid(&(-10)), 7);
  398. assert_eq!((&ibig!(-23)).rem_euclid(-10), 7);
  399. assert_eq!((&ibig!(-23)).rem_euclid(&(-10)), 7);
  400. assert_eq!(ibig!(-23).div_rem_euclid(-10), (ibig!(3), 7));
  401. assert_eq!(ibig!(-23).div_rem_euclid(&(-10)), (ibig!(3), 7));
  402. assert_eq!((&ibig!(-23)).div_rem_euclid(-10), (ibig!(3), 7));
  403. assert_eq!((&ibig!(-23)).div_rem_euclid(&(-10)), (ibig!(3), 7));
  404. }