diff --git a/sub_crates/sobol/build.rs b/sub_crates/sobol/build.rs index b7af9c7..8eabc12 100644 --- a/sub_crates/sobol/build.rs +++ b/sub_crates/sobol/build.rs @@ -7,7 +7,7 @@ use std::{env, fs::File, io::Write, path::Path}; const NUM_DIMENSIONS: usize = 1024; /// What file to generate the numbers from. -const DIRECTION_NUMBERS_TEXT: &str = include_str!("direction_numbers/joe-kuo-other-2.1024.txt"); +const DIRECTION_NUMBERS_TEXT: &str = include_str!("direction_numbers/new-joe-kuo-5.1024.txt"); fn main() { let out_dir = env::var("OUT_DIR").unwrap(); diff --git a/sub_crates/sobol/direction_numbers/joe-kuo-cessen-3.1024.txt b/sub_crates/sobol/direction_numbers/joe-kuo-cessen-3.1024.txt deleted file mode 100644 index 6ac613c..0000000 --- a/sub_crates/sobol/direction_numbers/joe-kuo-cessen-3.1024.txt +++ /dev/null @@ -1,1029 +0,0 @@ -# d s a m_i - -# Reordered to improve the 2d projections of -# adjecent dimensions. -2 1 0 1 -3 2 1 1 1 -5 3 2 1 3 1 -4 3 1 1 1 1 -10 5 7 1 1 7 7 9 -9 5 13 1 3 1 1 13 -7 4 4 1 1 3 7 -6 4 1 1 3 7 1 -14 6 1 1 3 7 9 1 9 -11 5 14 1 1 5 15 21 -16 6 13 1 1 5 5 31 9 -18 6 19 1 1 7 15 1 13 -8 5 2 1 3 1 7 23 -12 5 11 1 3 5 3 23 -13 5 4 1 1 1 15 27 -17 6 22 1 3 1 13 3 49 -25 7 56 1 3 7 9 5 35 105 -20 7 1 1 1 1 5 13 51 45 -26 7 14 1 1 1 1 19 21 69 -24 7 7 1 3 1 11 27 3 23 -19 6 25 1 3 3 13 29 37 -15 6 16 1 1 7 15 31 37 -21 7 32 1 3 1 7 19 23 51 -22 7 4 1 3 7 1 27 45 61 -23 7 8 1 3 7 3 9 51 1 -32 7 31 1 1 5 5 13 39 109 -28 7 19 1 1 3 13 15 17 25 -29 7 50 1 1 1 13 19 47 19 -37 7 59 1 1 5 5 13 5 19 -30 7 21 1 1 7 1 31 23 127 -33 7 62 1 3 7 15 23 19 127 -36 7 55 1 3 3 3 3 35 83 -31 7 42 1 1 7 1 31 59 123 -34 7 37 1 1 5 3 1 47 15 -35 7 41 1 1 3 13 23 1 61 -46 8 52 1 1 1 3 9 39 109 17 -40 8 21 1 3 1 9 5 63 33 27 -27 7 28 1 3 3 5 11 61 41 -38 8 14 1 1 5 9 13 51 25 117 -39 8 56 1 1 5 11 19 33 111 71 -54 9 8 1 3 1 9 27 41 83 63 89 -50 8 97 1 3 3 5 21 7 73 91 -47 8 67 1 1 1 5 17 17 47 153 -48 8 70 1 1 1 3 31 57 43 17 -44 8 49 1 1 1 1 13 1 109 29 -43 8 47 1 3 1 9 9 23 45 43 -41 8 22 1 3 1 5 19 57 15 181 -42 8 38 1 3 7 5 5 41 115 15 -58 9 25 1 3 5 15 13 7 89 151 391 -55 9 13 1 3 7 11 9 51 59 89 11 -57 9 22 1 3 7 5 29 33 83 229 287 -56 9 16 1 3 5 13 29 21 55 65 45 -53 8 122 1 3 7 5 7 15 39 247 -45 8 50 1 1 5 13 31 29 95 111 -61 9 52 1 1 1 15 1 23 111 175 125 -59 9 44 1 3 3 9 29 3 111 177 475 -52 8 115 1 1 7 3 17 55 37 189 -60 9 47 1 3 1 1 5 11 23 33 337 -49 8 84 1 3 7 3 29 17 91 153 -51 8 103 1 1 7 11 13 9 83 213 - -# Original order. -62 9 55 1 3 7 15 27 47 31 103 245 -63 9 59 1 3 5 7 7 3 127 187 245 -64 9 62 1 1 1 13 13 33 33 249 381 -65 9 67 1 1 3 11 27 45 95 31 419 -66 9 74 1 1 1 13 17 53 35 197 455 -67 9 81 1 1 5 7 21 43 93 237 467 -68 9 82 1 3 5 15 1 19 43 245 193 -69 9 87 1 1 7 11 21 21 23 147 429 -70 9 91 1 3 3 7 27 33 99 127 97 -71 9 94 1 1 5 11 13 21 101 73 17 -72 9 103 1 3 1 13 3 41 7 237 311 -73 9 104 1 3 7 1 1 43 115 215 499 -74 9 109 1 3 1 13 1 9 113 121 473 -75 9 122 1 1 5 5 21 7 5 171 41 -76 9 124 1 3 7 15 31 11 87 89 495 -77 9 137 1 1 5 1 9 9 19 161 325 -78 9 138 1 3 3 11 15 35 75 99 331 -79 9 143 1 3 1 1 9 27 103 93 187 -80 9 145 1 1 3 11 19 43 123 169 325 -81 9 152 1 1 7 9 17 59 65 139 89 -82 9 157 1 1 5 15 19 57 29 151 131 -83 9 167 1 1 1 1 23 49 75 85 365 -84 9 173 1 1 5 13 11 45 31 245 311 -85 9 176 1 3 7 5 5 13 23 1 461 -86 9 181 1 3 1 11 17 23 19 43 393 -87 9 182 1 1 3 5 27 9 37 105 491 -88 9 185 1 3 3 1 15 51 107 231 39 -89 9 191 1 1 1 5 21 57 59 143 47 -90 9 194 1 1 7 7 23 45 29 33 403 -91 9 199 1 3 7 5 31 7 91 5 355 -92 9 218 1 3 5 1 25 7 15 191 45 -93 9 220 1 1 7 5 9 33 63 75 45 -94 9 227 1 3 5 7 27 39 83 121 495 -95 9 229 1 1 7 5 1 25 71 31 33 -96 9 230 1 1 7 13 13 63 85 203 491 -97 9 234 1 3 3 11 21 39 25 213 423 -98 9 236 1 3 5 11 3 7 115 3 191 -99 9 241 1 1 1 9 19 1 51 247 301 -100 9 244 1 1 1 7 13 15 45 243 453 -101 9 253 1 1 1 7 25 1 113 153 105 -102 10 4 1 3 5 11 19 15 3 245 443 651 -103 10 13 1 3 7 15 13 35 5 167 233 9 -104 10 19 1 3 5 11 21 63 95 45 411 189 -105 10 22 1 1 7 7 31 11 81 87 15 677 -106 10 50 1 3 3 9 13 21 27 215 459 923 -107 10 55 1 1 1 9 13 49 23 7 441 31 -108 10 64 1 3 1 7 27 19 5 71 391 159 -109 10 69 1 3 7 9 9 55 77 95 353 519 -110 10 98 1 3 1 1 23 51 51 113 33 725 -111 10 107 1 3 3 1 21 3 107 223 261 913 -112 10 115 1 1 1 7 9 31 83 219 285 695 -113 10 121 1 1 3 15 1 1 23 49 39 485 -114 10 127 1 1 1 9 5 61 71 121 427 83 -115 10 134 1 1 3 3 5 45 111 227 261 957 -116 10 140 1 3 7 7 17 51 61 85 183 831 -117 10 145 1 3 3 5 13 35 27 135 189 743 -118 10 152 1 3 1 5 7 5 7 149 333 287 -119 10 158 1 3 1 13 19 57 127 201 365 359 -120 10 161 1 1 5 9 27 23 79 173 227 299 -121 10 171 1 1 7 7 9 35 123 1 289 43 -122 10 181 1 1 5 15 15 17 57 175 97 585 -123 10 194 1 1 3 13 9 23 123 45 463 137 -124 10 199 1 3 1 9 19 3 81 35 221 75 -125 10 203 1 3 7 13 31 41 53 207 253 957 -126 10 208 1 1 1 5 1 35 51 67 289 417 -127 10 227 1 1 1 11 27 21 121 235 49 3 -128 10 242 1 3 7 5 15 59 25 43 141 965 -129 10 251 1 1 5 3 7 23 121 161 223 245 -130 10 253 1 3 7 15 5 5 49 235 329 777 -131 10 265 1 3 1 11 7 49 41 187 221 639 -132 10 266 1 1 5 9 1 11 41 39 319 203 -133 10 274 1 3 3 3 21 61 23 233 397 251 -134 10 283 1 1 1 1 27 49 87 251 371 973 -135 10 289 1 1 5 11 11 43 57 59 153 695 -136 10 295 1 1 3 7 29 57 87 17 129 511 -137 10 301 1 3 7 7 5 39 71 27 59 49 -138 10 316 1 3 3 5 1 49 65 87 453 515 -139 10 319 1 3 3 9 19 3 3 167 85 653 -140 10 324 1 1 7 15 7 51 67 37 245 943 -141 10 346 1 1 5 7 1 7 21 5 135 999 -142 10 352 1 3 3 5 19 1 7 179 405 589 -143 10 361 1 1 3 15 3 3 43 79 501 287 -144 10 367 1 1 3 9 21 39 59 181 99 737 -145 10 382 1 3 7 1 3 11 91 33 273 651 -146 10 395 1 1 5 11 23 47 97 253 465 741 -147 10 398 1 1 5 1 7 45 33 203 297 193 -148 10 400 1 1 5 3 31 55 63 41 243 953 -149 10 412 1 3 5 9 9 53 51 71 281 447 -150 10 419 1 3 3 9 5 63 77 89 213 527 -151 10 422 1 1 3 11 31 15 3 91 321 997 -152 10 426 1 3 5 7 31 27 117 179 389 645 -153 10 428 1 3 3 9 21 37 41 49 143 401 -154 10 433 1 3 1 3 31 57 77 129 49 245 -155 10 446 1 3 1 15 27 7 23 69 177 1009 -156 10 454 1 3 3 3 11 9 101 141 313 739 -157 10 457 1 1 3 13 27 31 27 129 139 645 -158 10 472 1 3 7 5 7 11 99 145 313 641 -159 10 493 1 3 1 13 15 3 55 183 153 743 -160 10 505 1 3 5 1 23 63 53 237 267 333 -161 10 508 1 3 1 9 5 27 45 49 9 911 -162 11 2 1 1 3 9 11 51 97 231 295 877 353 -163 11 11 1 3 7 3 15 55 125 241 149 417 449 -164 11 21 1 1 7 1 7 31 89 91 259 819 849 -165 11 22 1 3 1 11 3 9 79 77 195 39 937 -166 11 35 1 1 3 15 27 37 123 233 221 661 323 -167 11 49 1 3 7 7 31 59 121 71 19 723 1713 -168 11 50 1 1 5 5 23 43 75 201 7 565 293 -169 11 56 1 1 1 3 13 39 115 209 239 439 1623 -170 11 61 1 3 7 5 29 39 69 241 71 465 77 -171 11 70 1 1 1 13 17 1 91 67 491 969 1709 -172 11 74 1 3 1 3 7 37 17 185 365 571 757 -173 11 79 1 1 5 1 29 47 87 65 171 893 1659 -174 11 84 1 3 5 5 13 5 21 151 91 801 1415 -175 11 88 1 3 3 13 27 35 81 49 417 315 1039 -176 11 103 1 3 3 7 13 11 63 19 35 659 83 -177 11 104 1 1 1 15 11 55 43 181 45 657 681 -178 11 112 1 3 7 1 23 9 73 55 143 731 401 -179 11 115 1 1 3 1 5 49 125 1 113 41 1733 -180 11 117 1 1 3 5 21 63 45 119 67 499 141 -181 11 122 1 3 5 11 29 39 45 93 415 515 413 -182 11 134 1 1 7 15 5 3 79 201 5 473 55 -183 11 137 1 1 3 7 25 31 55 61 175 919 2047 -184 11 146 1 1 7 3 31 61 15 147 409 249 1677 -185 11 148 1 3 3 9 11 37 9 65 383 689 111 -186 11 157 1 1 5 1 31 23 93 15 327 973 765 -187 11 158 1 3 7 15 25 17 115 245 205 169 13 -188 11 162 1 1 5 15 15 27 109 39 289 135 693 -189 11 164 1 1 7 3 23 1 119 199 1 661 1741 -190 11 168 1 3 5 13 7 29 3 55 369 903 715 -191 11 173 1 1 5 1 29 19 31 7 165 1019 467 -192 11 185 1 1 7 11 27 33 105 27 271 55 475 -193 11 186 1 1 5 1 27 9 63 157 507 491 1525 -194 11 191 1 3 3 9 1 41 71 239 471 147 505 -195 11 193 1 3 3 1 25 11 1 35 127 583 2043 -196 11 199 1 3 7 15 1 13 119 183 89 555 1739 -197 11 213 1 1 1 3 3 45 61 163 357 679 641 -198 11 214 1 3 5 13 23 27 57 125 393 917 939 -199 11 220 1 1 3 9 9 21 85 105 111 537 733 -200 11 227 1 1 3 15 9 61 43 103 65 227 315 -201 11 236 1 1 7 9 9 7 31 219 293 441 559 -202 11 242 1 1 3 1 23 21 25 229 303 537 1655 -203 11 251 1 3 7 7 31 5 21 177 435 899 1939 -204 11 256 1 3 3 3 11 43 77 131 353 11 157 -205 11 259 1 3 3 3 5 45 13 39 459 261 237 -206 11 265 1 3 5 7 23 1 55 63 9 755 589 -207 11 266 1 1 7 5 3 11 59 167 145 459 617 -208 11 276 1 3 1 13 31 49 53 235 1 37 1373 -209 11 292 1 1 7 11 27 63 21 59 183 297 15 -210 11 304 1 3 7 15 3 61 57 205 359 799 1591 -211 11 310 1 1 1 3 29 29 121 193 405 161 909 -212 11 316 1 1 5 5 15 43 19 87 337 595 551 -213 11 319 1 3 3 7 27 23 109 127 107 793 787 -214 11 322 1 1 3 5 9 11 63 21 363 403 981 -215 11 328 1 1 3 5 1 39 43 71 5 151 1535 -216 11 334 1 1 7 3 1 33 65 249 81 889 1591 -217 11 339 1 3 5 11 15 53 117 47 13 215 1953 -218 11 341 1 1 3 13 15 13 49 209 137 403 1373 -219 11 345 1 1 3 15 25 13 123 193 275 235 1975 -220 11 346 1 3 1 3 21 55 109 115 411 413 2023 -221 11 362 1 3 3 1 5 7 53 161 221 529 489 -222 11 367 1 1 3 5 25 27 103 31 7 103 343 -223 11 372 1 1 7 11 19 47 45 201 499 255 1501 -224 11 375 1 3 1 1 5 7 79 113 445 515 287 -225 11 376 1 1 7 11 7 35 87 145 115 773 1415 -226 11 381 1 1 1 13 15 25 83 143 347 473 415 -227 11 385 1 3 1 11 23 45 79 245 81 835 589 -228 11 388 1 1 1 13 25 49 17 15 255 421 1195 -229 11 392 1 1 3 1 13 31 103 189 43 389 1497 -230 11 409 1 3 5 7 3 63 21 151 339 525 455 -231 11 415 1 3 3 13 23 29 101 25 461 631 2017 -232 11 416 1 1 5 13 7 57 11 181 215 911 463 -233 11 421 1 3 7 1 1 45 41 101 263 835 1985 -234 11 428 1 3 3 3 25 13 127 245 137 709 993 -235 11 431 1 3 7 3 5 19 23 55 51 367 1513 -236 11 434 1 1 5 3 13 39 101 119 119 835 875 -237 11 439 1 1 7 13 17 27 103 31 223 461 1267 -238 11 446 1 3 5 13 29 29 29 189 13 647 1931 -239 11 451 1 3 5 15 9 27 127 43 341 939 1981 -240 11 453 1 1 7 11 27 27 25 163 83 709 141 -241 11 457 1 3 7 15 31 45 119 141 181 231 411 -242 11 458 1 1 5 9 31 25 113 243 435 55 765 -243 11 471 1 1 7 9 27 5 63 23 133 767 1963 -244 11 475 1 3 7 9 9 55 97 159 371 1013 2013 -245 11 478 1 1 3 7 3 17 13 215 15 281 291 -246 11 484 1 3 5 9 29 7 71 9 181 763 1111 -247 11 493 1 3 7 13 23 55 79 99 165 447 1947 -248 11 494 1 3 7 7 29 59 25 41 313 665 1245 -249 11 499 1 1 3 3 5 29 71 71 475 781 567 -250 11 502 1 1 1 3 19 53 61 215 475 885 459 -251 11 517 1 3 3 13 5 31 95 137 119 249 975 -252 11 518 1 1 1 15 21 51 71 21 193 709 1167 -253 11 524 1 3 7 9 29 61 45 153 351 867 995 -254 11 527 1 1 1 15 7 37 65 51 45 497 1825 -255 11 555 1 3 3 1 23 59 11 127 501 351 1525 -256 11 560 1 1 3 15 31 29 47 105 179 159 991 -257 11 565 1 3 3 15 27 23 93 163 169 1003 1369 -258 11 569 1 3 1 5 31 51 85 221 87 941 317 -259 11 578 1 3 7 15 21 19 71 205 385 371 1401 -260 11 580 1 1 5 11 1 51 31 185 407 233 1171 -261 11 587 1 3 3 9 19 11 93 253 373 679 69 -262 11 589 1 1 7 15 23 45 87 61 113 541 1701 -263 11 590 1 1 5 7 7 61 55 167 345 633 1699 -264 11 601 1 3 7 7 1 23 121 33 507 361 1945 -265 11 607 1 3 7 1 19 47 43 101 67 57 1685 -266 11 611 1 3 3 3 25 63 53 149 483 551 1519 -267 11 614 1 3 1 3 27 1 41 213 57 897 121 -268 11 617 1 3 1 5 19 13 49 81 239 455 957 -269 11 618 1 1 3 7 15 23 89 171 319 117 2041 -270 11 625 1 1 3 13 3 13 41 15 145 659 973 -271 11 628 1 3 7 1 15 17 69 53 245 179 159 -272 11 635 1 3 3 5 5 3 113 183 219 469 1169 -273 11 641 1 1 5 15 13 45 89 87 371 267 1643 -274 11 647 1 3 5 5 1 25 3 247 309 401 2027 -275 11 654 1 3 1 11 13 7 79 17 297 679 1699 -276 11 659 1 3 1 5 13 19 123 255 343 393 1635 -277 11 662 1 3 3 15 13 13 37 229 483 283 1423 -278 11 672 1 3 1 1 7 37 65 67 77 247 1539 -279 11 675 1 1 3 5 25 39 37 19 231 835 797 -280 11 682 1 1 1 11 27 53 1 217 211 7 801 -281 11 684 1 3 7 1 15 55 69 77 491 61 55 -282 11 689 1 1 1 1 29 21 47 167 281 105 1363 -283 11 695 1 3 1 15 15 55 57 207 479 709 867 -284 11 696 1 1 1 9 7 13 115 237 47 127 1849 -285 11 713 1 1 7 9 7 35 83 63 155 441 1105 -286 11 719 1 3 3 9 17 39 19 231 491 593 1997 -287 11 724 1 3 3 13 17 19 65 21 389 391 483 -288 11 733 1 1 5 7 17 7 97 173 413 917 1691 -289 11 734 1 1 3 3 29 49 3 127 337 713 1057 -290 11 740 1 3 7 5 31 23 31 119 161 51 537 -291 11 747 1 1 5 1 29 25 15 187 143 927 715 -292 11 749 1 1 7 11 15 47 55 103 109 465 1439 -293 11 752 1 1 3 13 31 47 115 177 355 163 833 -294 11 755 1 3 5 5 31 13 15 163 57 937 1553 -295 11 762 1 3 5 15 7 29 47 139 231 507 595 -296 11 770 1 3 1 15 31 13 49 253 29 109 727 -297 11 782 1 1 5 15 21 1 107 89 367 825 1297 -298 11 784 1 3 3 9 31 21 31 247 375 767 1761 -299 11 787 1 1 7 3 31 33 119 119 409 955 1411 -300 11 789 1 1 7 1 17 57 101 37 147 367 265 -301 11 793 1 3 3 3 25 53 91 43 351 1021 1243 -302 11 796 1 3 3 5 11 47 107 209 251 765 1089 -303 11 803 1 1 7 5 9 13 23 153 367 585 315 -304 11 805 1 3 3 15 23 57 83 173 413 555 159 -305 11 810 1 3 7 15 7 59 31 233 67 767 75 -306 11 815 1 3 5 3 17 1 35 37 445 847 1853 -307 11 824 1 3 5 13 3 1 75 113 83 467 1057 -308 11 829 1 1 7 1 7 15 125 31 341 133 285 -309 11 830 1 1 5 5 15 11 97 3 59 777 1125 -310 11 832 1 3 7 9 17 55 117 255 429 239 1071 -311 11 841 1 1 7 3 27 37 31 221 385 877 1767 -312 11 847 1 1 5 1 5 31 43 23 77 497 1605 -313 11 849 1 3 5 15 15 37 103 227 33 185 481 -314 11 861 1 3 3 11 29 61 67 53 257 609 1581 -315 11 871 1 1 3 15 31 25 123 59 395 107 587 -316 11 878 1 3 1 7 29 19 51 135 503 71 33 -317 11 889 1 3 1 15 21 31 23 85 17 227 657 -318 11 892 1 1 1 3 9 63 61 205 61 699 599 -319 11 901 1 3 1 15 1 39 45 133 1 391 255 -320 11 908 1 3 5 13 17 13 107 245 55 619 1085 -321 11 920 1 3 3 9 5 25 5 49 509 661 1225 -322 11 923 1 3 7 9 29 61 13 247 341 257 1551 -323 11 942 1 3 7 9 5 39 123 221 391 501 527 -324 11 949 1 3 1 7 21 23 65 233 181 405 1297 -325 11 950 1 1 1 3 23 23 15 253 37 679 1455 -326 11 954 1 3 3 3 3 31 89 105 11 1015 1837 -327 11 961 1 3 1 7 15 59 105 179 477 923 1189 -328 11 968 1 3 5 9 1 9 21 71 145 409 1235 -329 11 971 1 1 7 13 9 43 85 241 497 683 1169 -330 11 973 1 3 3 13 29 41 23 51 1 563 905 -331 11 979 1 3 5 9 25 51 1 71 25 195 363 -332 11 982 1 3 7 7 31 49 13 255 471 653 1127 -333 11 986 1 3 3 7 7 61 63 21 413 817 1369 -334 11 998 1 3 5 1 7 53 43 167 173 261 633 -335 11 1001 1 1 1 13 21 59 49 135 351 187 447 -336 11 1010 1 1 3 7 5 35 103 5 247 555 1439 -337 11 1012 1 1 3 1 3 11 5 153 507 21 1331 -338 12 41 1 3 7 15 25 11 29 51 353 447 337 3621 -339 12 52 1 3 5 13 31 51 103 69 311 645 1373 1155 -340 12 61 1 1 7 7 1 55 101 203 503 711 131 937 -341 12 62 1 3 3 11 7 5 75 63 217 643 871 3511 -342 12 76 1 1 3 15 21 23 75 3 113 203 13 2727 -343 12 104 1 3 5 11 7 1 27 39 411 195 607 1705 -344 12 117 1 1 7 11 23 21 11 113 437 947 1897 1249 -345 12 131 1 1 7 5 23 45 41 169 33 981 1331 2013 -346 12 143 1 3 1 7 11 5 11 121 239 269 1759 829 -347 12 145 1 1 1 9 31 51 77 61 19 121 969 3627 -348 12 157 1 3 1 13 23 45 39 183 21 787 967 1911 -349 12 167 1 1 1 15 21 19 91 231 325 531 473 2543 -350 12 171 1 3 5 3 7 25 91 193 387 181 493 1225 -351 12 176 1 1 7 7 23 21 55 221 283 995 793 485 -352 12 181 1 1 1 15 29 61 1 33 59 725 1053 55 -353 12 194 1 3 3 5 15 43 85 87 425 39 1559 1403 -354 12 217 1 1 5 9 3 7 65 3 177 491 121 2755 -355 12 236 1 1 3 3 15 37 67 145 295 421 1849 2151 -356 12 239 1 1 3 15 17 41 109 125 509 407 1927 2401 -357 12 262 1 1 7 5 9 39 69 1 259 821 1837 947 -358 12 283 1 1 7 9 25 55 127 185 143 211 1615 1949 -359 12 286 1 3 3 13 15 25 25 231 197 329 1319 853 -360 12 307 1 3 7 11 29 1 23 131 473 497 1147 3779 -361 12 313 1 3 1 15 21 23 67 243 191 353 1091 2121 -362 12 319 1 1 7 11 27 15 41 183 235 219 303 795 -363 12 348 1 3 7 15 3 21 19 187 365 453 1069 1983 -364 12 352 1 1 7 1 27 11 95 117 433 453 1839 3863 -365 12 357 1 3 3 5 29 63 17 243 247 609 1999 785 -366 12 391 1 1 1 7 21 17 57 99 227 217 585 397 -367 12 398 1 1 1 13 13 41 75 153 355 175 313 2237 -368 12 400 1 3 7 9 27 63 1 33 383 507 1315 1691 -369 12 412 1 1 3 3 17 33 65 235 395 317 1769 2937 -370 12 415 1 3 1 5 31 19 25 15 373 823 2003 1309 -371 12 422 1 3 3 5 3 33 43 225 445 181 325 797 -372 12 440 1 1 3 1 11 19 49 129 73 983 621 2381 -373 12 460 1 3 3 15 1 33 109 33 363 925 589 3955 -374 12 465 1 1 3 3 11 21 23 99 159 949 977 3103 -375 12 468 1 1 5 3 3 15 41 45 325 703 907 2467 -376 12 515 1 1 7 15 31 11 71 225 199 521 933 2777 -377 12 536 1 3 5 9 13 55 23 199 293 297 2031 3815 -378 12 539 1 1 3 13 25 5 25 229 301 301 1165 2229 -379 12 551 1 1 1 11 1 39 25 211 233 233 355 3131 -380 12 558 1 3 3 9 31 37 113 17 163 547 877 333 -381 12 563 1 1 7 1 31 7 117 175 191 191 1185 2637 -382 12 570 1 1 3 11 27 59 33 79 293 657 1899 1547 -383 12 595 1 3 1 7 29 43 89 227 315 971 1175 693 -384 12 598 1 1 1 5 9 43 17 37 321 907 441 21 -385 12 617 1 3 1 15 25 47 21 11 385 353 1103 3445 -386 12 647 1 1 7 9 21 63 9 147 79 627 471 1683 -387 12 654 1 3 7 7 1 13 29 127 3 629 1213 2981 -388 12 678 1 1 5 9 7 41 41 235 131 379 23 1251 -389 12 713 1 3 5 3 17 31 101 193 271 189 1923 1849 -390 12 738 1 1 1 5 21 43 35 209 359 185 1921 2773 -391 12 747 1 1 7 5 31 5 41 127 415 151 1027 2171 -392 12 750 1 1 1 3 31 53 127 37 273 415 1657 467 -393 12 757 1 3 7 3 29 43 77 35 347 37 1265 1211 -394 12 772 1 1 1 1 3 31 89 57 405 491 1685 2727 -395 12 803 1 3 5 15 29 9 23 79 253 161 273 4001 -396 12 810 1 3 3 13 29 49 73 145 203 215 423 3293 -397 12 812 1 1 5 15 29 43 49 99 335 1017 315 2761 -398 12 850 1 3 1 9 15 39 87 243 331 1015 1771 1303 -399 12 862 1 1 1 5 5 1 15 19 277 699 1957 373 -400 12 906 1 1 7 1 7 13 41 211 167 705 1011 2097 -401 12 908 1 1 7 9 3 25 127 75 271 817 629 2027 -402 12 929 1 3 1 5 31 49 17 113 71 233 443 3005 -403 12 930 1 1 1 1 1 23 87 29 121 485 227 2981 -404 12 954 1 1 1 1 19 3 127 31 57 567 375 1121 -405 12 964 1 1 1 9 15 49 67 49 473 957 1955 29 -406 12 982 1 1 1 1 15 49 71 241 43 611 953 4013 -407 12 985 1 3 3 3 13 57 45 65 5 101 511 3317 -408 12 991 1 3 3 7 25 9 97 35 495 791 223 1555 -409 12 992 1 1 3 1 25 37 73 61 31 525 855 695 -410 12 1067 1 3 5 11 15 45 11 53 393 557 349 2507 -411 12 1070 1 1 5 11 1 43 77 215 149 81 1663 987 -412 12 1096 1 3 5 5 13 55 45 79 305 65 1721 605 -413 12 1099 1 1 1 3 15 19 103 249 253 935 1933 3809 -414 12 1116 1 1 1 11 1 43 119 233 79 539 1397 1669 -415 12 1143 1 1 7 15 9 27 29 125 335 573 1721 551 -416 12 1165 1 1 5 11 29 33 121 47 351 857 1155 2939 -417 12 1178 1 1 1 13 21 13 107 179 135 597 39 2129 -418 12 1184 1 3 5 15 11 5 101 227 457 773 1369 2229 -419 12 1202 1 1 7 3 9 47 53 47 315 221 369 2991 -420 12 1213 1 3 7 9 7 1 117 229 193 105 1591 875 -421 12 1221 1 3 1 13 5 49 47 95 431 37 2037 3793 -422 12 1240 1 1 3 11 17 39 91 225 109 375 1431 3689 -423 12 1246 1 3 5 7 27 39 23 215 55 573 965 321 -424 12 1252 1 1 3 13 1 51 23 181 399 471 429 2409 -425 12 1255 1 1 7 9 17 21 43 177 327 789 1581 3505 -426 12 1267 1 1 7 7 9 39 59 187 203 655 545 2107 -427 12 1293 1 1 7 3 15 13 83 189 213 575 255 1475 -428 12 1301 1 1 1 3 23 27 101 67 299 1003 1351 1799 -429 12 1305 1 3 3 9 13 29 83 251 151 607 2047 447 -430 12 1332 1 3 1 5 5 7 23 145 241 437 1215 489 -431 12 1349 1 3 1 11 3 63 55 53 187 179 1137 3923 -432 12 1384 1 3 7 11 31 29 127 37 509 835 1693 965 -433 12 1392 1 1 1 9 29 1 11 61 459 415 317 2981 -434 12 1402 1 3 7 3 1 61 31 77 55 301 1333 2803 -435 12 1413 1 1 7 13 3 1 117 185 385 771 409 3667 -436 12 1417 1 1 5 13 19 3 43 53 41 111 1617 95 -437 12 1423 1 3 7 13 25 25 125 91 85 99 461 3999 -438 12 1451 1 3 5 7 7 29 15 155 237 443 1187 623 -439 12 1480 1 3 3 1 11 3 35 169 449 917 141 1525 -440 12 1491 1 1 5 7 17 11 85 31 165 35 1069 3395 -441 12 1503 1 3 7 15 11 21 89 61 459 211 959 2411 -442 12 1504 1 1 5 11 25 35 1 137 491 955 1287 3955 -443 12 1513 1 3 3 11 13 47 105 69 119 455 769 1665 -444 12 1538 1 3 3 13 17 9 27 25 445 611 855 1273 -445 12 1544 1 1 3 9 3 5 107 89 443 991 727 871 -446 12 1547 1 1 5 3 19 43 93 163 193 395 99 1717 -447 12 1555 1 1 1 7 23 43 121 205 249 255 1763 1867 -448 12 1574 1 1 3 3 1 45 33 191 475 503 761 1845 -449 12 1603 1 1 3 3 31 39 11 203 9 183 407 3787 -450 12 1615 1 1 5 11 11 3 65 175 131 757 1717 1003 -451 12 1618 1 1 1 15 13 61 29 123 269 531 257 1399 -452 12 1629 1 3 5 5 13 23 23 1 293 921 79 479 -453 12 1634 1 1 3 13 9 55 21 15 47 685 829 2679 -454 12 1636 1 3 5 11 1 39 97 59 37 787 991 263 -455 12 1639 1 3 5 15 31 3 45 159 219 161 1025 665 -456 12 1657 1 3 7 15 11 51 51 19 181 125 1677 2125 -457 12 1667 1 1 3 5 27 23 79 161 97 581 405 2239 -458 12 1681 1 3 1 11 17 1 97 33 31 771 707 1053 -459 12 1697 1 1 1 11 11 51 113 47 51 827 385 3211 -460 12 1704 1 3 3 7 17 3 99 57 59 199 1359 2459 -461 12 1709 1 1 5 11 1 7 19 37 53 435 1095 495 -462 12 1722 1 1 5 13 19 5 95 109 245 737 1163 33 -463 12 1730 1 3 5 3 1 17 79 195 153 857 603 2905 -464 12 1732 1 3 5 5 31 51 71 103 121 5 931 217 -465 12 1802 1 1 5 7 31 39 11 9 223 365 1317 449 -466 12 1804 1 3 5 3 17 19 85 121 161 75 527 3941 -467 12 1815 1 3 5 9 17 53 41 255 207 251 925 2395 -468 12 1826 1 3 3 3 13 5 109 125 293 165 1715 1017 -469 12 1832 1 3 7 1 23 27 45 141 267 187 611 59 -470 12 1843 1 3 1 9 11 9 81 239 469 905 525 1549 -471 12 1849 1 1 1 15 13 21 71 87 199 389 139 3499 -472 12 1863 1 3 7 5 21 41 107 195 369 391 271 3351 -473 12 1905 1 3 3 1 17 19 75 219 415 241 1441 967 -474 12 1928 1 1 1 15 7 7 79 213 391 831 1649 3741 -475 12 1933 1 3 5 3 29 57 99 207 243 915 1969 2499 -476 12 1939 1 1 3 13 1 3 27 211 271 657 1251 3333 -477 12 1976 1 3 3 1 3 31 49 53 273 837 1771 1405 -478 12 1996 1 3 7 7 1 33 125 221 179 601 565 2709 -479 12 2013 1 1 1 7 31 19 81 173 413 849 579 2487 -480 12 2014 1 1 7 11 17 23 97 145 383 1015 547 483 -481 12 2020 1 3 1 3 11 63 41 175 43 283 1441 2695 -482 13 13 1 1 7 9 5 23 61 161 149 49 1851 1281 4615 -483 13 19 1 1 5 7 23 47 111 213 317 107 849 1563 6863 -484 13 26 1 1 7 9 31 53 115 233 497 23 1395 327 551 -485 13 41 1 1 7 13 21 51 17 109 29 1003 1221 2175 4645 -486 13 50 1 3 3 15 7 57 103 131 63 217 607 2311 1237 -487 13 55 1 3 5 3 19 27 107 223 15 149 1777 2211 2241 -488 13 69 1 1 1 1 5 37 87 59 435 37 1771 1675 6241 -489 13 70 1 3 1 7 29 51 83 41 11 545 995 1625 4883 -490 13 79 1 1 1 5 15 43 57 111 55 921 567 1767 1849 -491 13 82 1 1 5 15 5 49 3 23 307 691 1727 2297 6705 -492 13 87 1 3 7 5 31 1 111 53 499 1017 1889 1955 1327 -493 13 93 1 3 7 11 17 35 71 223 149 743 319 3157 5961 -494 13 94 1 1 1 3 21 31 27 79 353 753 341 2121 5455 -495 13 97 1 1 5 3 21 29 27 91 111 359 27 2383 7149 -496 13 100 1 3 7 9 5 51 115 65 275 211 1519 2247 2065 -497 13 112 1 1 1 15 27 7 21 173 79 357 693 873 879 -498 13 121 1 3 1 7 25 47 53 123 491 707 1723 1139 3757 -499 13 134 1 1 7 5 15 55 103 17 17 345 1477 3423 7527 -500 13 138 1 3 1 13 1 3 55 85 251 389 631 117 4879 -501 13 148 1 1 7 1 29 33 17 3 303 299 465 1709 1111 -502 13 151 1 1 7 5 3 29 111 249 383 611 859 1555 7387 -503 13 157 1 3 5 5 29 41 33 63 395 53 351 247 5317 -504 13 161 1 1 5 7 11 63 121 131 249 451 619 2395 4575 -505 13 179 1 3 7 3 3 43 37 211 385 289 791 735 4683 -506 13 181 1 1 3 11 7 13 53 7 497 771 903 119 6601 -507 13 188 1 3 7 1 3 61 117 237 427 335 1749 3105 6493 -508 13 196 1 3 5 7 11 55 121 93 3 953 433 2663 1633 -509 13 203 1 1 7 1 7 51 43 31 31 819 791 1071 569 -510 13 206 1 3 5 11 17 37 71 97 69 681 1133 2425 3217 -511 13 223 1 1 1 15 9 15 31 131 113 567 593 1781 6715 -512 13 224 1 1 7 9 29 25 39 159 495 507 1253 409 6723 -513 13 227 1 1 1 5 13 29 25 197 489 853 87 3173 6389 -514 13 230 1 1 7 11 11 7 27 227 347 777 893 1007 1953 -515 13 239 1 1 3 7 31 9 61 39 251 1 1297 2861 2593 -516 13 241 1 1 5 1 15 45 3 147 143 279 2045 649 5579 -517 13 248 1 1 5 5 9 51 3 249 253 533 1305 1749 1287 -518 13 253 1 3 7 9 13 15 85 165 183 943 41 4061 5845 -519 13 268 1 3 3 1 23 31 49 247 287 83 1833 1249 2789 -520 13 274 1 3 3 7 7 53 125 211 145 837 515 1971 6201 -521 13 283 1 3 3 11 13 11 87 11 33 355 1049 3053 4801 -522 13 286 1 3 5 13 23 25 63 221 385 217 687 455 4083 -523 13 289 1 1 5 1 3 21 119 79 207 693 925 1519 3981 -524 13 301 1 1 3 11 21 35 21 21 297 461 793 1679 1285 -525 13 302 1 3 5 13 31 23 11 129 189 933 283 3435 4593 -526 13 316 1 3 1 1 13 11 113 35 199 149 181 3183 6205 -527 13 319 1 3 1 9 7 59 105 161 63 821 13 2561 755 -528 13 324 1 1 1 7 27 3 87 243 279 13 1341 2383 6005 -529 13 331 1 3 3 13 27 25 91 205 275 663 1875 2399 4179 -530 13 333 1 1 7 15 27 25 51 235 391 951 759 285 3443 -531 13 345 1 3 3 5 1 33 57 107 137 851 1937 1003 2427 -532 13 351 1 1 5 11 27 3 27 81 303 151 433 1893 587 -533 13 358 1 1 1 9 25 45 113 173 471 191 1031 2119 4853 -534 13 375 1 1 7 11 27 23 107 229 485 197 1925 1809 7775 -535 13 379 1 1 1 3 31 61 91 47 3 643 1487 2563 3893 -536 13 381 1 1 7 15 19 27 97 109 283 67 335 3707 1909 -537 13 386 1 1 1 11 23 59 101 249 121 919 1385 2119 7389 -538 13 403 1 1 7 11 7 29 33 111 243 617 125 3887 727 -539 13 405 1 1 3 1 7 41 91 177 315 555 1653 1725 8069 -540 13 419 1 1 5 15 3 63 53 31 177 451 141 3727 1695 -541 13 426 1 1 3 15 25 37 17 229 75 705 327 1367 807 -542 13 428 1 3 5 5 27 1 35 97 171 167 1077 25 5497 -543 13 439 1 1 7 7 11 17 105 41 53 389 941 3273 967 -544 13 440 1 3 3 15 11 61 125 155 361 229 1981 1647 5077 -545 13 446 1 3 1 15 5 11 35 203 425 509 993 1861 1977 -546 13 451 1 3 5 15 25 7 55 89 101 939 7 3535 7543 -547 13 454 1 1 3 1 19 45 117 249 261 529 1299 597 7207 -548 13 458 1 3 5 13 17 47 73 231 353 941 325 115 4005 -549 13 465 1 3 1 9 5 29 11 161 461 793 813 3755 7017 -550 13 468 1 3 3 15 23 21 83 47 143 279 1009 3315 3801 -551 13 472 1 1 7 9 21 55 31 129 325 59 701 2227 1491 -552 13 475 1 1 7 1 21 41 87 151 269 885 795 4087 2919 -553 13 477 1 1 3 11 7 19 103 45 461 109 1921 2497 4969 -554 13 496 1 1 3 7 7 49 9 185 95 277 1325 2393 7585 -555 13 502 1 3 3 13 13 1 71 77 245 865 1693 2187 2589 -556 13 508 1 3 1 15 19 59 17 189 101 511 23 677 4091 -557 13 517 1 3 5 7 9 45 45 233 47 973 795 2339 595 -558 13 521 1 3 1 7 15 35 79 53 429 299 389 237 2261 -559 13 527 1 3 1 1 9 35 27 169 85 755 145 91 257 -560 13 530 1 3 7 11 25 33 41 151 11 401 1335 2111 3309 -561 13 532 1 1 7 3 17 1 61 253 289 917 1157 1513 4783 -562 13 542 1 1 5 7 23 51 59 171 161 815 969 1199 4137 -563 13 552 1 3 3 7 7 27 81 175 89 145 19 1469 3775 -564 13 555 1 3 1 7 7 33 95 179 177 375 1077 3711 2903 -565 13 560 1 3 1 15 15 27 127 245 307 949 1037 3753 4993 -566 13 566 1 3 3 7 29 15 59 39 9 331 1877 3565 6483 -567 13 575 1 1 5 3 1 41 29 185 351 837 747 3263 749 -568 13 577 1 1 7 1 19 47 1 239 489 79 767 1123 8145 -569 13 589 1 3 7 15 7 39 23 239 465 981 185 837 7837 -570 13 590 1 3 5 5 23 15 21 221 259 663 185 1203 2857 -571 13 602 1 1 3 5 31 43 107 199 117 161 403 2227 7969 -572 13 607 1 1 3 5 5 5 105 153 349 601 815 939 3547 -573 13 608 1 1 7 15 23 33 85 51 309 63 1799 1599 1711 -574 13 611 1 1 7 3 15 41 71 169 87 381 1695 3189 3343 -575 13 613 1 1 7 11 15 57 81 107 57 773 111 2741 6399 -576 13 625 1 1 5 15 5 55 65 35 387 797 1093 4039 3041 -577 13 644 1 1 5 3 29 5 87 235 503 539 483 2785 4211 -578 13 651 1 1 3 5 11 53 101 83 61 97 411 2063 2951 -579 13 654 1 1 1 5 25 33 45 253 183 251 1323 1891 5143 -580 13 656 1 1 5 11 7 61 67 201 465 851 1 1533 2839 -581 13 662 1 3 3 7 15 3 51 9 425 671 223 2203 481 -582 13 668 1 1 1 7 25 5 45 153 479 437 221 1085 455 -583 13 681 1 3 7 1 7 53 65 35 153 423 989 2779 4787 -584 13 682 1 3 1 3 29 43 7 219 47 431 627 3735 383 -585 13 689 1 1 5 13 31 49 121 29 235 325 1247 1765 7173 -586 13 696 1 1 1 9 25 15 73 167 191 215 259 111 2289 -587 13 699 1 3 7 9 23 1 119 15 155 153 1419 667 6309 -588 13 707 1 1 5 13 7 53 91 161 445 283 571 3623 5737 -589 13 709 1 3 3 5 3 39 29 47 393 871 345 805 8167 -590 13 714 1 3 1 5 21 9 1 169 471 693 53 967 4145 -591 13 716 1 3 7 11 17 19 123 133 511 15 675 3935 5465 -592 13 719 1 3 5 13 25 51 63 87 503 61 907 105 7599 -593 13 727 1 3 5 13 27 3 67 189 193 131 1227 227 1975 -594 13 734 1 3 5 5 31 17 89 187 391 575 1715 2551 5371 -595 13 738 1 3 1 9 11 63 9 251 399 703 1557 107 675 -596 13 743 1 3 5 15 29 55 15 65 51 811 1387 3815 7973 -597 13 747 1 3 1 13 19 33 87 7 183 97 2011 3849 3015 -598 13 757 1 1 3 11 29 57 27 185 457 181 1341 2033 2181 -599 13 769 1 3 7 5 25 19 19 147 101 73 499 897 3053 -600 13 770 1 3 5 5 29 29 109 5 197 581 1829 3679 283 -601 13 776 1 3 5 1 27 63 71 131 309 193 1663 3009 1041 -602 13 790 1 3 3 5 19 29 95 173 195 205 207 3597 2223 -603 13 799 1 1 5 1 3 45 127 251 175 663 415 1789 6577 -604 13 805 1 3 1 9 1 51 73 165 377 525 381 2829 4619 -605 13 809 1 3 1 15 23 19 79 115 297 75 1505 1407 5865 -606 13 812 1 1 1 5 27 33 93 123 449 747 2027 51 4507 -607 13 820 1 3 1 7 13 39 49 111 101 247 1375 3577 6769 -608 13 827 1 1 7 9 19 49 35 139 57 527 1943 1781 1839 -609 13 829 1 3 7 7 3 7 67 163 289 139 1415 2149 7501 -610 13 835 1 3 5 5 13 47 21 47 3 105 473 2013 2443 -611 13 841 1 1 3 13 31 25 21 75 505 795 1311 2467 5541 -612 13 844 1 1 5 1 25 15 21 113 427 249 293 1089 1589 -613 13 856 1 1 7 5 17 11 3 211 219 37 781 1299 3669 -614 13 859 1 3 1 15 23 3 55 59 273 947 1647 3111 2615 -615 13 862 1 1 3 7 21 33 81 253 267 499 299 4015 1669 -616 13 865 1 3 7 5 5 17 75 205 53 237 1791 3671 4205 -617 13 885 1 3 1 5 3 59 77 39 295 649 1219 3971 7675 -618 13 890 1 3 5 15 1 1 13 207 445 91 571 869 359 -619 13 905 1 3 5 5 15 35 79 175 313 745 1689 1499 5899 -620 13 916 1 3 7 7 19 51 57 131 491 301 1405 3789 6081 -621 13 925 1 1 3 9 27 53 55 193 413 535 715 2823 2063 -622 13 935 1 1 7 13 31 47 47 99 441 125 999 3575 8007 -623 13 939 1 1 1 11 17 45 67 27 329 267 773 565 5559 -624 13 942 1 1 7 11 31 23 125 117 373 523 571 1385 3841 -625 13 949 1 3 7 15 3 15 63 91 23 97 311 3045 2025 -626 13 953 1 3 7 5 29 37 105 63 511 35 537 4047 7449 -627 13 956 1 3 1 13 5 9 29 51 269 395 1737 2097 199 -628 13 961 1 3 7 15 23 33 71 111 497 369 1265 2937 1883 -629 13 968 1 1 3 3 3 49 117 243 57 83 3 3319 5711 -630 13 976 1 1 5 11 25 41 29 123 405 623 1537 3467 699 -631 13 988 1 3 7 13 17 5 87 103 501 581 713 1535 4637 -632 13 995 1 3 3 3 25 59 91 241 163 537 1249 3427 6647 -633 13 997 1 3 7 1 27 7 49 219 51 989 599 2593 6051 -634 13 1007 1 1 7 9 17 19 31 23 199 729 819 287 4689 -635 13 1015 1 1 1 11 13 49 97 229 411 29 1815 3761 3849 -636 13 1016 1 3 1 1 23 35 15 33 3 865 995 1263 3779 -637 13 1027 1 1 7 7 17 1 103 117 455 93 599 1075 4389 -638 13 1036 1 3 7 5 27 3 67 117 101 443 1773 3731 629 -639 13 1039 1 3 1 11 9 51 63 217 273 425 1889 1113 3799 -640 13 1041 1 1 3 3 25 13 37 89 429 493 637 3287 6733 -641 13 1048 1 1 1 9 17 29 45 89 141 633 243 4035 1835 -642 13 1053 1 1 1 9 15 41 61 239 349 3 1669 1767 6257 -643 13 1054 1 3 5 11 21 41 17 65 495 879 1891 667 3807 -644 13 1058 1 1 5 11 7 25 101 231 159 287 1125 3859 4331 -645 13 1075 1 1 1 9 13 59 23 85 393 611 1213 7 5707 -646 13 1082 1 3 5 15 13 33 61 255 163 465 641 703 6391 -647 13 1090 1 1 3 3 29 41 25 19 247 973 1815 235 7215 -648 13 1109 1 1 3 11 17 43 83 105 493 663 1395 1079 745 -649 13 1110 1 1 3 13 17 53 5 55 31 459 673 1253 5601 -650 13 1119 1 3 7 7 1 49 85 73 265 555 43 2235 6787 -651 13 1126 1 1 7 5 15 3 43 201 289 755 1605 1927 6179 -652 13 1130 1 1 7 9 15 5 101 19 17 879 1155 3845 3225 -653 13 1135 1 3 5 15 15 63 103 133 209 121 1865 4039 3335 -654 13 1137 1 3 7 11 15 51 87 241 455 431 1111 1319 1959 -655 13 1140 1 3 7 5 21 47 67 251 91 575 481 3269 1473 -656 13 1149 1 3 1 1 23 5 37 169 105 803 1765 1823 5081 -657 13 1156 1 1 7 9 1 37 55 175 431 77 171 3759 1227 -658 13 1159 1 3 7 9 3 61 111 139 87 649 243 1109 1303 -659 13 1160 1 3 1 11 17 25 9 117 227 821 1219 2567 5851 -660 13 1165 1 1 7 1 19 53 89 57 29 885 275 2289 4341 -661 13 1173 1 3 1 5 5 1 43 237 235 667 1433 2497 5535 -662 13 1178 1 1 1 11 21 9 29 165 153 157 1209 817 7819 -663 13 1183 1 3 1 5 1 37 33 9 493 497 1469 357 5929 -664 13 1184 1 3 3 9 25 55 5 157 397 461 537 1085 6737 -665 13 1189 1 3 3 13 3 61 83 5 5 243 1745 845 1965 -666 13 1194 1 3 1 3 31 63 3 75 355 799 417 791 2763 -667 13 1211 1 3 3 11 5 31 99 163 69 699 809 1651 7249 -668 13 1214 1 1 7 3 29 53 51 185 379 743 611 963 5631 -669 13 1216 1 1 5 11 21 15 59 75 121 247 307 2483 2861 -670 13 1225 1 1 1 5 21 5 33 43 225 675 213 3081 5047 -671 13 1231 1 3 3 15 29 45 19 19 15 605 1531 1305 7725 -672 13 1239 1 3 1 5 19 21 73 131 37 667 1953 1211 7499 -673 13 1243 1 3 1 5 13 3 105 151 131 965 1953 335 4991 -674 13 1246 1 3 1 13 23 37 81 147 353 125 141 4017 7543 -675 13 1249 1 1 1 7 15 13 25 107 431 327 1069 155 1729 -676 13 1259 1 1 5 1 19 3 51 113 349 845 1157 1771 4923 -677 13 1273 1 1 3 7 17 19 47 249 215 959 659 3989 3719 -678 13 1274 1 3 5 5 27 57 41 29 153 1007 1209 1373 2847 -679 13 1281 1 3 5 7 5 21 23 67 127 755 211 2517 8161 -680 13 1287 1 1 7 15 1 5 65 197 435 275 473 1333 3669 -681 13 1294 1 3 7 1 17 19 61 11 425 717 351 2323 6575 -682 13 1296 1 1 7 3 29 13 115 19 235 371 81 3107 3103 -683 13 1305 1 1 3 1 23 19 57 179 7 405 217 3159 6675 -684 13 1306 1 3 3 3 19 39 121 195 205 489 1959 3617 799 -685 13 1318 1 3 5 9 19 35 115 3 17 201 671 1847 597 -686 13 1332 1 3 1 3 7 1 121 77 163 769 977 1223 4215 -687 13 1335 1 1 1 15 9 13 33 57 271 733 1699 3083 813 -688 13 1336 1 1 7 7 15 49 101 183 339 129 1865 2697 3131 -689 13 1341 1 3 1 3 9 39 33 141 459 971 1621 3893 4393 -690 13 1342 1 1 3 3 27 61 53 79 375 613 663 2325 4707 -691 13 1362 1 1 7 1 27 59 59 179 173 333 649 1053 555 -692 13 1364 1 3 3 13 15 55 25 91 45 739 1769 1341 2007 -693 13 1368 1 1 3 15 9 59 27 45 133 941 1037 1559 7145 -694 13 1378 1 3 5 1 13 5 3 167 235 609 1023 1967 771 -695 13 1387 1 1 1 15 19 15 45 143 309 697 337 3399 3599 -696 13 1389 1 1 5 3 9 1 19 41 39 1003 997 1357 5849 -697 13 1397 1 3 3 7 29 37 9 165 9 855 337 2363 2633 -698 13 1401 1 3 5 1 23 19 97 55 411 623 287 3365 7675 -699 13 1408 1 3 3 3 11 1 47 81 377 85 859 3717 1979 -700 13 1418 1 1 3 1 11 3 125 43 471 127 341 2943 231 -701 13 1425 1 1 5 15 29 21 101 193 389 727 1397 2799 4387 -702 13 1426 1 3 5 15 27 19 13 195 105 5 985 767 4573 -703 13 1431 1 1 5 3 13 29 33 123 107 353 723 3821 5425 -704 13 1435 1 3 5 11 23 11 69 219 15 735 1015 431 6511 -705 13 1441 1 3 1 1 27 49 79 91 201 365 255 1309 7165 -706 13 1444 1 3 1 13 21 61 39 111 71 137 551 3237 6681 -707 13 1462 1 3 5 13 19 21 121 227 381 677 421 449 5375 -708 13 1471 1 1 7 5 1 23 95 199 509 197 829 941 2269 -709 13 1474 1 3 7 5 21 61 119 59 97 841 453 2439 7345 -710 13 1483 1 1 1 3 23 21 19 245 281 841 273 1967 2129 -711 13 1485 1 3 3 5 7 43 19 253 239 421 461 4023 6981 -712 13 1494 1 3 1 11 25 3 109 107 55 781 1109 101 889 -713 13 1497 1 3 7 13 1 21 27 3 303 557 1381 407 1849 -714 13 1516 1 1 3 13 3 61 121 119 453 137 2039 3333 333 -715 13 1522 1 3 5 1 21 9 113 9 413 301 1795 2827 1837 -716 13 1534 1 3 5 15 13 31 119 107 445 277 1555 2121 2173 -717 13 1543 1 1 1 5 9 5 5 97 191 949 1603 2181 6879 -718 13 1552 1 3 3 9 5 47 77 43 483 137 61 1811 5309 -719 13 1557 1 3 7 11 25 63 107 207 23 231 1685 4009 391 -720 13 1558 1 3 7 9 3 35 119 83 193 269 289 3051 8137 -721 13 1567 1 3 3 11 25 25 17 243 195 503 1797 1913 749 -722 13 1568 1 3 1 13 5 35 101 57 185 881 2005 3493 2715 -723 13 1574 1 3 1 11 25 51 113 177 175 855 1727 1901 8181 -724 13 1592 1 1 3 9 9 63 109 233 205 357 499 3367 1189 -725 13 1605 1 3 1 15 17 11 35 135 501 319 1465 2095 5269 -726 13 1606 1 1 3 9 1 27 99 201 67 269 135 3891 4067 -727 13 1610 1 1 1 9 15 37 83 135 3 645 499 2597 2705 -728 13 1617 1 3 1 7 17 33 89 51 295 809 743 27 2599 -729 13 1623 1 1 7 1 5 55 11 97 177 515 1651 1183 5803 -730 13 1630 1 1 3 11 21 3 39 249 55 891 1681 709 7833 -731 13 1634 1 3 5 13 11 21 81 21 99 223 355 2015 7195 -732 13 1640 1 3 7 11 25 23 9 131 417 583 1927 277 2923 -733 13 1643 1 1 5 5 1 43 35 87 471 499 1411 2867 2477 -734 13 1648 1 1 5 9 25 43 83 103 183 527 483 935 7097 -735 13 1651 1 1 5 7 31 13 49 169 135 885 1367 2627 5917 -736 13 1653 1 3 7 5 1 39 97 167 95 745 1067 805 2885 -737 13 1670 1 3 7 15 31 41 1 147 15 795 323 309 5075 -738 13 1676 1 3 7 5 23 3 65 143 183 671 993 547 5971 -739 13 1684 1 3 7 1 21 59 29 45 367 363 1707 541 7045 -740 13 1687 1 1 7 15 11 25 77 33 409 495 1713 1479 2857 -741 13 1691 1 1 3 1 11 59 113 217 313 79 379 627 3917 -742 13 1693 1 1 1 15 31 61 73 211 419 965 1817 3965 1425 -743 13 1698 1 1 7 13 9 9 61 119 223 387 595 191 1723 -744 13 1709 1 3 5 3 25 49 83 175 329 1003 729 1399 6553 -745 13 1715 1 1 5 9 5 37 31 43 365 403 1485 1621 3753 -746 13 1722 1 3 5 7 17 5 125 219 311 1 639 1939 4367 -747 13 1732 1 3 5 3 21 63 119 31 151 325 1063 2477 1223 -748 13 1735 1 1 5 13 15 23 1 141 105 21 747 161 3555 -749 13 1747 1 3 7 1 31 21 11 27 267 205 1569 765 1357 -750 13 1749 1 1 3 9 23 59 93 51 81 367 1545 1161 3481 -751 13 1754 1 3 7 11 15 51 43 205 237 617 1075 1051 5285 -752 13 1777 1 3 7 3 11 5 95 203 495 661 1119 1095 3391 -753 13 1784 1 3 3 1 25 39 59 159 273 101 1887 417 3423 -754 13 1790 1 3 5 13 23 45 13 87 139 313 1527 3739 7009 -755 13 1795 1 1 5 1 23 43 127 149 63 823 1075 3123 3079 -756 13 1801 1 1 7 3 11 1 29 127 313 971 1363 1765 6439 -757 13 1802 1 1 3 13 1 19 61 125 345 51 273 2783 5255 -758 13 1812 1 1 1 7 29 33 21 239 351 481 795 2865 2235 -759 13 1828 1 3 3 3 19 31 41 125 327 169 837 753 191 -760 13 1831 1 3 7 15 5 15 5 109 333 127 2021 441 6307 -761 13 1837 1 3 7 9 23 39 79 235 27 573 1989 1401 7597 -762 13 1838 1 1 7 3 23 29 119 53 105 609 1139 321 1111 -763 13 1840 1 1 7 13 13 43 3 211 389 563 1047 1373 4669 -764 13 1845 1 3 3 1 9 7 1 135 455 539 1681 1851 4585 -765 13 1863 1 3 1 13 19 45 65 109 391 563 753 2405 2623 -766 13 1864 1 3 1 13 31 41 29 153 507 559 877 125 4827 -767 13 1867 1 1 7 3 15 53 39 1 143 351 505 1213 4829 -768 13 1870 1 1 5 5 27 49 39 101 501 239 1233 1695 6929 -769 13 1877 1 1 5 9 13 29 103 185 343 943 1709 801 1227 -770 13 1881 1 1 1 3 3 1 63 23 199 53 1629 325 7609 -771 13 1884 1 3 7 1 11 47 61 171 71 493 2035 869 5925 -772 13 1903 1 1 7 7 5 3 95 163 89 395 1033 2051 5589 -773 13 1917 1 3 5 9 3 21 29 143 333 805 267 4077 383 -774 13 1918 1 1 3 13 9 17 87 163 107 1011 1099 3395 987 -775 13 1922 1 3 1 1 5 51 71 249 441 513 791 1831 2757 -776 13 1924 1 1 3 13 23 7 101 211 289 1001 1623 3239 1685 -777 13 1928 1 1 5 15 31 53 91 245 143 237 557 3099 3737 -778 13 1931 1 3 1 1 5 27 51 173 321 85 1385 9 1887 -779 13 1951 1 3 1 15 11 17 49 97 121 189 2031 565 7513 -780 13 1952 1 1 5 5 1 3 55 157 27 783 1543 2645 1079 -781 13 1957 1 1 5 11 27 9 83 155 17 147 1561 3213 4483 -782 13 1958 1 3 5 15 23 63 127 97 459 869 671 3287 995 -783 13 1964 1 3 3 1 5 1 89 229 269 493 709 2719 4803 -784 13 1967 1 1 3 3 27 5 85 55 491 781 711 3193 2051 -785 13 1970 1 1 5 13 21 9 31 107 381 613 991 2289 325 -786 13 1972 1 3 7 11 11 47 115 103 101 945 1993 1345 2031 -787 13 1994 1 3 7 11 27 59 123 73 85 975 397 2619 2095 -788 13 2002 1 1 1 13 19 31 87 45 39 43 1041 2015 5859 -789 13 2007 1 1 3 15 5 19 47 193 293 961 113 3013 889 -790 13 2008 1 1 1 3 31 61 61 209 493 785 1825 3779 8095 -791 13 2023 1 1 3 9 17 51 103 255 133 165 281 3945 2779 -792 13 2030 1 1 7 3 9 33 105 29 449 231 1393 1071 5043 -793 13 2035 1 1 7 5 29 35 115 1 467 969 343 2291 5099 -794 13 2038 1 1 3 5 13 21 99 151 13 861 189 1709 4943 -795 13 2042 1 1 7 15 5 41 53 89 235 515 1819 3907 4861 -796 13 2047 1 3 5 1 7 15 75 65 247 581 853 445 7623 -797 13 2051 1 1 3 5 25 29 27 73 359 729 1991 1075 8107 -798 13 2058 1 1 3 1 5 35 83 119 385 657 1085 2707 4403 -799 13 2060 1 3 5 5 25 55 51 205 257 887 827 2795 925 -800 13 2071 1 1 3 1 5 29 113 195 439 241 127 2783 2577 -801 13 2084 1 3 1 11 15 47 7 149 81 897 235 3935 6981 -802 13 2087 1 3 3 1 27 5 69 125 377 829 1549 3025 2319 -803 13 2099 1 1 7 15 15 19 3 205 291 209 1817 2439 1001 -804 13 2108 1 1 7 9 21 7 101 17 13 615 137 1821 1387 -805 13 2111 1 3 3 11 21 27 61 133 183 699 245 527 471 -806 13 2120 1 3 3 13 5 25 111 69 417 271 1545 2729 475 -807 13 2128 1 3 5 15 11 47 93 93 311 775 789 1163 1563 -808 13 2138 1 3 5 1 13 63 27 113 449 971 1021 73 1229 -809 13 2143 1 3 3 15 17 37 103 13 91 823 1767 1295 6941 -810 13 2144 1 1 3 9 15 59 125 91 195 409 65 2873 1921 -811 13 2153 1 3 1 7 17 31 21 251 145 385 1385 1505 5735 -812 13 2156 1 3 7 11 17 19 85 85 81 583 889 669 5187 -813 13 2162 1 1 5 9 23 53 97 217 129 517 271 2967 3937 -814 13 2167 1 1 1 13 3 45 49 193 35 671 1039 2841 4547 -815 13 2178 1 3 1 15 21 29 5 91 23 247 1701 1585 1493 -816 13 2183 1 3 3 11 15 31 21 165 153 607 221 2781 6159 -817 13 2202 1 3 7 15 15 27 107 25 457 995 849 877 407 -818 13 2211 1 3 1 15 7 17 37 31 203 585 369 609 2509 -819 13 2214 1 3 1 11 17 5 17 71 453 171 1317 1985 5817 -820 13 2223 1 1 5 5 23 39 123 179 49 19 1337 1051 5289 -821 13 2225 1 3 1 13 17 43 47 139 21 851 1253 1939 299 -822 13 2232 1 1 7 1 3 35 57 115 481 653 1975 1319 5225 -823 13 2237 1 3 7 15 5 13 3 249 21 791 1349 1961 6975 -824 13 2257 1 1 5 11 17 3 97 103 201 313 939 693 5397 -825 13 2260 1 1 7 9 29 59 35 91 135 19 1505 1715 3335 -826 13 2267 1 3 5 1 17 53 59 217 1 305 1257 3773 3509 -827 13 2274 1 3 7 7 5 29 47 245 487 525 1409 245 6427 -828 13 2276 1 3 5 15 25 39 49 153 461 1003 1527 2883 3499 -829 13 2285 1 1 7 13 23 11 15 163 9 77 1853 1679 1023 -830 13 2288 1 3 7 13 3 55 101 229 327 537 79 3947 6259 -831 13 2293 1 3 7 5 9 39 3 151 439 181 1473 3613 7425 -832 13 2294 1 3 7 11 21 5 53 11 39 455 1131 393 1525 -833 13 2297 1 1 1 7 5 45 49 229 173 567 1285 2317 2207 -834 13 2303 1 1 1 7 5 33 53 241 477 161 1607 1253 4497 -835 13 2308 1 1 7 15 11 3 107 187 249 529 905 2883 7777 -836 13 2311 1 1 5 7 17 13 123 7 449 627 1679 569 4145 -837 13 2318 1 1 1 1 3 29 57 219 209 51 1475 2989 7373 -838 13 2323 1 3 3 5 19 33 95 241 257 465 1305 189 269 -839 13 2332 1 1 3 1 13 11 67 53 157 693 1583 739 3827 -840 13 2341 1 1 3 9 17 9 5 7 421 425 1851 3847 5147 -841 13 2345 1 1 3 3 23 47 123 215 11 643 1543 1859 4073 -842 13 2348 1 1 1 11 5 35 79 55 225 773 1717 2891 563 -843 13 2354 1 1 1 15 5 17 115 3 107 33 1897 3677 4313 -844 13 2368 1 1 5 7 5 57 75 65 177 269 733 3837 1211 -845 13 2377 1 1 5 15 25 45 125 159 501 17 693 1083 1801 -846 13 2380 1 3 5 3 15 17 17 159 421 449 83 1367 7719 -847 13 2383 1 3 7 3 7 31 121 43 491 501 335 2731 1959 -848 13 2388 1 1 1 11 7 1 43 187 437 113 1205 1255 1319 -849 13 2395 1 1 5 15 15 49 103 139 187 697 907 1229 3205 -850 13 2397 1 1 5 13 9 19 25 89 345 229 1325 859 527 -851 13 2401 1 1 1 3 3 63 105 11 407 131 679 3069 4153 -852 13 2411 1 1 3 11 21 27 113 59 409 845 309 2639 129 -853 13 2413 1 1 5 15 21 47 55 237 201 103 331 3531 4965 -854 13 2419 1 1 5 11 9 11 85 49 209 333 1993 511 7347 -855 13 2435 1 1 5 13 25 49 109 241 297 677 313 657 4503 -856 13 2442 1 1 1 15 23 57 61 117 461 559 837 2781 709 -857 13 2455 1 1 3 3 31 13 109 167 389 1021 1437 31 811 -858 13 2472 1 1 1 9 21 17 19 61 241 445 1935 3161 4061 -859 13 2478 1 3 5 5 17 5 59 117 241 913 1315 1597 713 -860 13 2490 1 3 1 15 5 41 91 191 43 49 1491 4051 2317 -861 13 2507 1 1 1 3 13 33 53 39 341 341 617 2627 421 -862 13 2509 1 3 3 15 23 25 47 251 405 659 721 229 2269 -863 13 2517 1 1 3 3 25 19 123 53 59 805 2033 1693 6223 -864 13 2524 1 1 1 3 13 45 67 67 199 925 145 775 7923 -865 13 2528 1 1 5 7 17 23 49 37 185 253 1743 2771 6899 -866 13 2531 1 3 5 5 15 53 75 185 17 379 1091 4029 5399 -867 13 2538 1 1 1 9 13 37 101 33 499 723 57 467 7227 -868 13 2545 1 1 7 1 21 49 17 49 417 895 1099 3199 7059 -869 13 2546 1 1 5 13 27 35 49 35 213 799 197 3599 4135 -870 13 2555 1 3 5 3 31 33 61 127 87 343 845 1699 3897 -871 13 2557 1 3 1 11 7 59 19 161 119 377 1637 2025 849 -872 13 2564 1 1 3 3 3 51 3 67 77 357 295 2383 7423 -873 13 2573 1 3 1 3 23 25 67 29 433 49 901 1641 3021 -874 13 2579 1 3 7 3 25 59 15 181 193 809 1089 423 1629 -875 13 2592 1 3 5 1 3 19 53 45 277 411 455 9 3651 -876 13 2598 1 3 7 15 13 33 37 15 489 275 1113 3371 3915 -877 13 2607 1 1 5 11 3 3 97 47 471 713 1975 1679 7901 -878 13 2612 1 1 3 15 31 47 13 33 111 607 85 4049 2937 -879 13 2619 1 1 5 1 7 61 17 179 301 245 67 2823 4121 -880 13 2621 1 3 7 15 7 25 97 73 19 899 1357 2041 6493 -881 13 2627 1 1 1 3 3 63 117 231 125 527 1795 1807 7331 -882 13 2633 1 1 5 11 17 35 61 161 415 687 1015 1239 6925 -883 13 2636 1 3 3 5 15 19 1 5 209 545 1341 4019 6177 -884 13 2642 1 3 5 15 13 3 7 39 63 727 1033 2551 449 -885 13 2654 1 1 1 5 11 35 39 143 361 659 1407 2399 6421 -886 13 2660 1 1 7 13 3 23 21 161 153 595 1847 1339 1887 -887 13 2669 1 3 3 9 3 43 75 45 433 811 1045 2773 2415 -888 13 2675 1 1 3 5 23 29 73 167 211 403 539 13 4433 -889 13 2684 1 3 3 13 19 25 113 181 343 349 931 115 5687 -890 13 2694 1 3 1 11 25 39 23 91 359 221 37 575 3035 -891 13 2703 1 1 7 1 19 61 99 81 167 53 1493 2623 5797 -892 13 2706 1 1 3 5 13 13 105 53 449 1003 1405 1403 2557 -893 13 2712 1 1 3 15 1 37 51 69 13 5 1865 3585 3409 -894 13 2715 1 1 7 13 15 31 61 53 439 665 57 445 435 -895 13 2722 1 3 1 3 1 27 1 155 65 191 449 2185 1787 -896 13 2727 1 1 3 11 13 59 15 43 435 477 823 717 299 -897 13 2734 1 1 3 1 23 9 85 155 127 543 1081 221 7357 -898 13 2742 1 3 3 1 7 35 77 9 79 283 449 2729 4193 -899 13 2745 1 1 7 1 9 15 97 99 363 513 1063 3051 1529 -900 13 2751 1 1 5 3 1 17 101 19 425 377 1369 3071 6877 -901 13 2766 1 1 5 1 23 63 11 7 29 717 963 2205 7487 -902 13 2768 1 3 1 13 1 15 43 151 59 315 691 2765 7723 -903 13 2780 1 3 1 5 5 31 17 81 235 419 123 2983 1797 -904 13 2790 1 1 1 11 29 17 123 115 91 531 1393 3537 6387 -905 13 2794 1 3 5 11 31 5 107 61 463 507 1937 4075 4727 -906 13 2796 1 3 7 1 9 27 23 97 125 511 239 1031 153 -907 13 2801 1 3 5 9 13 55 53 195 265 217 357 567 1613 -908 13 2804 1 1 5 15 5 51 51 235 63 393 1915 2495 7067 -909 13 2807 1 1 7 1 15 35 97 65 115 819 1419 2313 6029 -910 13 2816 1 1 5 13 11 35 45 149 45 617 2001 1629 5611 -911 13 2821 1 3 5 11 15 17 39 155 143 473 1269 2537 2571 -912 13 2831 1 1 3 3 1 7 57 227 151 731 1871 245 7605 -913 13 2834 1 1 1 13 21 1 19 189 299 483 513 2235 3079 -914 13 2839 1 1 5 15 13 53 91 55 339 161 273 2759 7251 -915 13 2845 1 3 1 1 17 47 29 107 371 437 715 739 2195 -916 13 2852 1 3 5 11 27 1 35 105 443 623 1401 2683 4327 -917 13 2856 1 1 5 7 13 23 101 235 369 951 1983 1021 7835 -918 13 2861 1 1 1 5 21 57 45 67 89 209 1105 2161 2245 -919 13 2873 1 1 7 15 7 33 69 129 181 231 961 1621 7763 -920 13 2874 1 3 1 5 1 31 105 123 413 989 1311 199 5175 -921 13 2888 1 3 5 13 15 9 123 215 463 197 537 531 2331 -922 13 2893 1 1 3 13 13 45 61 189 91 445 1719 3335 5269 -923 13 2894 1 3 1 11 19 41 65 143 199 177 1695 3647 7169 -924 13 2902 1 1 3 13 13 25 37 89 175 551 1535 1401 1823 -925 13 2917 1 3 3 1 21 15 31 101 329 907 731 1493 451 -926 13 2921 1 1 3 9 5 61 85 247 289 547 1031 3567 1769 -927 13 2922 1 3 1 7 25 35 99 43 25 855 1451 617 2615 -928 13 2929 1 1 1 7 25 9 77 39 127 327 1501 57 1485 -929 13 2935 1 3 5 5 9 13 13 239 125 473 785 2853 5165 -930 13 2946 1 1 1 3 27 43 117 63 175 487 1505 157 6043 -931 13 2951 1 1 5 5 25 43 97 231 389 191 1795 2659 7409 -932 13 2957 1 1 7 7 11 25 47 131 363 311 523 2773 8015 -933 13 2960 1 1 5 11 25 61 123 83 265 995 613 1893 6341 -934 13 2966 1 1 3 3 27 45 41 143 173 1009 1963 4009 3675 -935 13 2972 1 3 3 9 5 11 95 253 415 307 501 3117 4621 -936 13 2976 1 3 5 1 31 41 95 39 473 565 1389 4049 5673 -937 13 2979 1 1 3 11 3 13 87 101 413 769 1223 549 4967 -938 13 2985 1 3 1 3 13 41 29 187 491 197 1685 3859 3799 -939 13 3000 1 3 1 3 17 55 91 131 429 899 29 1193 7323 -940 13 3003 1 1 3 1 11 17 107 113 71 185 1057 897 249 -941 13 3013 1 1 3 5 29 27 87 33 273 405 547 763 2195 -942 13 3018 1 3 1 11 7 43 39 103 475 515 1105 767 3881 -943 13 3020 1 3 3 1 13 21 93 141 407 1001 1939 2747 6479 -944 13 3025 1 3 1 15 1 31 15 69 313 547 987 3569 2795 -945 13 3042 1 3 5 7 1 17 105 113 253 453 1219 2769 197 -946 13 3047 1 3 7 7 29 55 77 169 227 471 1385 2689 7995 -947 13 3048 1 3 7 15 7 15 113 133 205 271 263 3141 5007 -948 13 3051 1 1 5 7 7 57 121 209 319 733 1689 2817 6621 -949 13 3054 1 3 7 9 23 37 101 229 269 429 157 515 6975 -950 13 3056 1 1 7 7 23 31 83 237 299 333 959 1467 419 -951 13 3065 1 3 7 3 23 45 17 33 293 679 625 3651 4205 -952 13 3073 1 3 1 13 17 19 87 237 219 119 841 1923 8141 -953 13 3074 1 1 1 15 15 49 21 189 287 759 695 913 5167 -954 13 3083 1 1 3 15 1 9 125 5 67 603 1691 2863 1299 -955 13 3086 1 3 7 3 5 35 105 27 51 889 1469 3203 6997 -956 13 3091 1 3 5 7 1 53 77 1 313 475 1515 2233 399 -957 13 3097 1 3 5 7 3 31 43 17 425 673 1817 1939 8027 -958 13 3109 1 1 5 5 31 31 119 179 157 635 1 315 1179 -959 13 3116 1 1 5 5 1 7 25 15 311 489 829 1929 623 -960 13 3124 1 1 5 11 1 17 43 205 347 713 2029 977 4345 -961 13 3128 1 1 3 9 17 3 109 253 427 333 203 1609 1851 -962 13 3153 1 3 1 9 21 55 107 203 483 869 1417 587 3965 -963 13 3160 1 3 3 13 27 21 103 109 313 797 723 683 349 -964 13 3165 1 3 5 3 1 35 43 19 367 487 703 2211 1199 -965 13 3172 1 3 1 3 19 33 13 135 329 621 843 407 6523 -966 13 3175 1 1 7 5 13 1 45 133 27 267 1975 2761 5061 -967 13 3184 1 1 3 3 19 59 13 35 489 391 1595 2983 3203 -968 13 3193 1 3 3 7 19 27 125 55 121 733 1437 51 4029 -969 13 3196 1 3 1 15 11 7 73 187 209 127 1259 1733 1897 -970 13 3200 1 3 7 1 27 19 87 189 91 113 237 2073 419 -971 13 3203 1 3 5 3 15 27 107 77 467 953 1317 877 6911 -972 13 3205 1 3 7 13 7 17 107 157 439 871 651 797 2095 -973 13 3209 1 1 3 9 13 19 19 75 13 267 435 1067 6603 -974 13 3224 1 1 7 13 25 7 81 207 149 587 1561 2025 2355 -975 13 3239 1 3 3 5 25 25 79 249 371 305 305 53 1051 -976 13 3251 1 1 3 9 19 63 57 237 39 349 1163 2443 1739 -977 13 3254 1 1 5 5 23 57 61 21 115 107 303 1451 2955 -978 13 3265 1 3 1 3 17 19 121 123 25 367 963 545 2885 -979 13 3266 1 3 3 5 15 57 123 1 193 861 1001 1457 4977 -980 13 3275 1 3 5 9 21 31 37 83 361 491 1933 407 1923 -981 13 3280 1 3 7 5 9 47 33 43 239 837 799 3731 5883 -982 13 3283 1 3 5 5 7 39 49 201 339 783 1937 1473 391 -983 13 3286 1 3 5 13 3 63 39 91 63 645 929 1543 577 -984 13 3301 1 3 3 3 9 59 107 129 455 259 1523 107 285 -985 13 3302 1 1 3 9 23 21 89 13 413 87 1703 1611 5667 -986 13 3305 1 1 3 9 21 9 83 187 159 757 1111 3619 3907 -987 13 3319 1 3 3 9 21 35 27 57 357 659 1243 2725 5769 -988 13 3323 1 1 7 1 3 7 75 145 119 1019 1169 953 2343 -989 13 3326 1 3 3 11 19 15 121 135 187 569 671 131 479 -990 13 3331 1 3 5 1 21 27 63 187 367 181 1885 2513 4499 -991 13 3348 1 3 3 5 27 13 107 183 345 855 123 1491 867 -992 13 3351 1 3 7 9 5 41 43 47 277 657 1227 2265 1777 -993 13 3358 1 1 5 5 31 37 119 5 395 399 787 2869 829 -994 13 3368 1 1 1 5 3 55 43 117 191 577 1185 3407 4673 -995 13 3374 1 3 5 11 15 63 103 127 505 53 297 3411 3919 -996 13 3376 1 3 7 13 15 57 45 17 253 261 767 1951 6599 -997 13 3379 1 3 1 13 3 55 11 113 223 403 249 633 3265 -998 13 3385 1 1 1 5 7 49 47 47 467 633 1991 361 8099 -999 13 3386 1 1 3 5 31 45 87 117 483 331 687 3393 7377 -1000 13 3396 1 1 7 9 29 3 11 147 147 631 1131 2651 3141 -1001 13 3420 1 3 7 9 27 49 41 207 305 403 889 2947 5775 -1002 13 3423 1 3 1 13 1 63 25 127 401 375 727 2375 2945 -1003 13 3430 1 3 3 15 7 19 65 29 383 809 279 2081 1993 -1004 13 3433 1 1 7 7 29 9 39 159 457 435 947 2309 2219 -1005 13 3434 1 1 1 11 15 31 51 147 167 113 145 125 6997 -1006 13 3439 1 1 5 1 27 21 123 249 21 487 331 35 3101 -1007 13 3442 1 3 7 3 5 7 127 173 115 607 1061 2923 6869 -1008 13 3444 1 3 1 5 17 55 13 15 95 365 1419 1111 7561 -1009 13 3453 1 3 3 15 23 45 115 123 477 335 85 431 6893 -1010 13 3464 1 1 5 1 5 49 95 19 429 549 713 661 167 -1011 13 3477 1 1 3 15 3 13 127 189 385 649 67 499 5355 -1012 13 3478 1 1 1 5 9 31 17 21 419 655 287 1173 1023 -1013 13 3482 1 1 5 1 17 11 19 17 195 1015 1621 1197 3497 -1014 13 3487 1 3 1 5 29 55 85 43 381 59 149 391 2109 -1015 13 3497 1 1 7 9 11 43 115 57 211 463 787 3665 5619 -1016 13 3500 1 3 1 7 23 5 55 93 117 925 941 91 6221 -1017 13 3505 1 1 1 13 29 35 89 7 75 759 735 3011 1279 -1018 13 3506 1 1 1 15 13 17 13 175 167 965 1069 2215 6065 -1019 13 3511 1 1 3 15 1 11 109 205 501 419 1071 543 5377 -1020 13 3512 1 1 3 5 11 1 123 147 87 289 325 3693 4939 -1021 13 3515 1 1 3 1 9 21 29 215 365 995 1727 2515 5741 -1022 13 3525 1 1 3 5 3 7 5 147 413 505 661 3071 2521 -1023 13 3532 1 3 7 5 5 41 117 145 455 379 601 455 5377 -1024 13 3538 1 1 7 15 11 45 27 37 457 935 595 3883 1491 diff --git a/sub_crates/sobol/src/lib.rs b/sub_crates/sobol/src/lib.rs index 0b07290..c48d018 100644 --- a/sub_crates/sobol/src/lib.rs +++ b/sub_crates/sobol/src/lib.rs @@ -54,24 +54,33 @@ pub fn sample_owen(dimension: u32, index: u32, scramble: u32) -> f32 { u32_to_0_1_f32(owen_scramble_u32(sobol_u32(dimension, index), scramble)) } -/// Same as `sample()` except applies both Owen scrambling and -/// Cranley-Patterson rotation using the given scramble parameter. +/// Same as `sample_owen()` except uses a slower more full version of +/// Owen scrambling. /// -/// For the technically curious: this first does Owen scrambling, and then -/// Cranley-Patterson. If it were done the other way around it would -/// invalidate the Owen scrambling. -/// -/// To get proper scrambling and rotation, you need to use a different scramble -/// value for each dimension, and those values should be generated more-or-less -/// randomly. For example, using a 32-bit hash of the dimension parameter -/// works well. +/// This is mainly intended to help validate the faster Owen scrambling, +/// and likely shouldn't be used for real things. It is significantly +/// slower. #[inline] -pub fn sample_owen_cranley(dimension: u32, index: u32, scramble: u32) -> f32 { - // Reusing the same scramble parameter for both the Owen scrambling and - // the Cranely-Patterson rotation actually works fine, because the Owen - // scrambling is implemented as a sort of hash, so they don't end up being - // correlated. - u32_to_0_1_f32(owen_scramble_u32(sobol_u32(dimension, index), scramble).wrapping_add(scramble)) +pub fn sample_owen_slow(dimension: u32, index: u32, scramble: u32) -> f32 { + let mut n = sobol_u32(dimension, index); + n ^= scramble; + for i in 0..16 { + let mask = (1 << (31 - i)) - 1; + let hash = { + let mut hash = n & (!mask); + let seed = scramble + i; + let perms = [0x29aaaaa7, 0x736caf6f, 0x54aad35b, 0x2ab35aaa]; + for p in perms.iter().cycle().take(6) { + hash = hash.wrapping_mul(*p); + hash ^= hash.wrapping_shr(16); + hash ^= seed; + } + hash + }; + n ^= hash & mask; + } + + u32_to_0_1_f32(n) } //----------------------------------------------------------------------