# This b-file was computed by Antti Karttunen, 2019-10-27 # in GP/PARI CALCULATOR Version 2.9.4 (released) amd64 running linux (x86-64/GMP-6.1.2 kernel) 64-bit version, compiled: Dec 19 2017, gcc version 7.3.0 (Ubuntu 7.3.0-1ubuntu1), threading engine: pthread # with the following program: # # default(parisizemax,2^31); # up_to = 1081; \\ = binomial(46+1,2) # A257993(n) = { for(i=1,oo,if(n%prime(i),return(i))); } # A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); }; # A328578(n) = A257993(A276086(A276086(n))); # memoA328631sq = Map(); # A328631sq(n, k) = { my(v=0); if(!mapisdefined(memoA328631sq,[n,k-1],&v),if(1==k, v=-1, v = A328631sq(n, k-1))); for(i=1+v,oo,if(A328578(i)==n,mapput(memoA328631sq,[n,k],i); return(i))); }; # A328631list(up_to) = { my(v = vector(up_to), i=0); for(a=1,oo, for(col=1,a, i++; if(i > up_to, return(v)); v[i] = A328631sq(col,(a-(col-1))))); (v); }; # v328631 = A328631list(up_to); # A328631(n) = v328631[n]; # for(n=1,up_to,write("b328631.txt", n, " ", A328631(n))); 1 1 2 3 3 0 4 5 5 12 6 2 7 7 8 24 9 6 10 4 11 9 12 30 13 18 14 8 15 10 16 11 17 42 18 34 19 16 20 14 21 20 22 13 23 54 24 36 25 32 26 38 27 22 28 28 29 15 30 60 31 48 32 64 33 58 34 26 35 50 36 82 37 17 38 72 39 66 40 152 41 62 42 40 43 52 44 88 45 116 46 19 47 84 48 78 49 184 50 112 51 44 52 56 53 106 54 118 55 148 56 21 57 90 58 96 59 210 60 166 61 46 62 74 63 110 64 140 65 178 66 208 67 23 68 102 69 108 70 242 71 176 72 68 73 76 74 128 75 142 76 196 77 412 78 418 79 25 80 114 81 122 82 274 83 214 84 70 85 80 86 130 87 146 88 200 89 416 90 598 91 838 92 27 93 120 94 126 95 362 96 218 97 86 98 94 99 134 100 160 101 202 102 562 103 626 104 1016 105 1048 106 29 107 132 108 138 109 394 110 226 111 92 112 98 113 158 114 164 115 206 116 568 117 628 118 1036 119 1228 120 1466 121 31 122 144 123 154 124 440 125 240 126 220 127 100 128 182 129 170 130 328 131 586 132 778 133 1040 134 1252 135 1468 136 1858 137 33 138 150 139 156 140 448 141 360 142 224 143 104 144 262 145 188 146 352 147 590 148 800 149 1042 150 1256 151 1648 152 1882 153 2276 154 35 155 162 156 168 157 452 158 650 159 238 160 124 161 266 162 190 163 374 164 592 165 802 166 1046 167 1258 168 1670 169 1886 170 2296 171 4612 172 37 173 174 174 186 175 484 176 658 177 248 178 136 179 268 180 194 181 376 182 596 183 826 184 1192 185 1406 186 1676 187 1888 188 2300 189 4616 190 9028 191 39 192 180 193 198 194 572 195 660 196 270 197 172 198 284 199 292 200 380 201 614 202 830 203 1196 204 1408 205 1678 206 2060 207 2302 208 4618 209 9208 210 11546 211 41 212 192 213 212 214 604 215 780 216 272 217 230 218 290 219 296 220 382 221 616 222 832 223 1198 224 1430 225 1798 226 2062 227 2306 228 6688 229 9232 230 13648 231 16162 232 43 233 204 234 222 235 634 236 844 237 286 238 232 239 304 240 298 241 506 242 620 243 956 244 1216 245 1436 246 1826 247 2066 248 4198 249 6712 250 9236 251 13822 252 16168 253 18478 254 45 255 216 256 234 257 638 258 848 259 356 260 236 261 308 262 316 263 508 264 718 265 958 266 1220 267 1438 268 1828 269 2068 270 4376 271 6716 272 9238 273 13826 274 18238 275 20780 276 23098 277 47 278 228 279 244 280 646 281 856 282 370 283 250 284 310 285 320 286 530 287 746 288 976 289 1222 290 1450 291 1846 292 2086 293 4378 294 6868 295 11128 296 13828 297 18262 298 20782 299 25196 300 27688 301 49 302 246 303 252 304 662 305 1080 306 424 307 254 308 314 309 322 310 532 311 748 312 1000 313 1226 314 1456 315 1850 316 2090 317 4400 318 6892 319 11306 320 13846 321 18266 322 20786 323 25198 324 27712 325 30028 326 51 327 258 328 264 329 694 330 1200 331 428 332 256 333 334 334 326 335 550 336 764 337 1004 338 1238 339 1462 340 1852 341 2092 342 4402 343 6896 344 11330 345 15746 346 18268 347 20788 348 25348 349 27716 350 57748 351 90086 352 53 353 276 354 282 355 782 356 1280 357 500 358 260 359 338 360 340 361 554 362 770 363 1028 364 1240 365 1586 366 1868 367 2096 368 4406 369 6898 370 11332 371 15748 372 18418 373 22678 374 25372 375 27718 376 60026 377 90088 378 147838 379 55 380 288 381 294 382 814 383 1288 384 538 385 278 386 350 387 344 388 560 389 772 390 1030 391 1244 392 1588 393 1874 394 2216 395 4408 396 6914 397 11336 398 15922 399 18440 400 22876 401 25376 402 29816 403 60046 404 117802 405 149938 406 180178 407 57 408 306 409 312 410 872 411 1474 412 544 413 280 414 400 415 346 416 578 417 776 418 1034 419 1250 420 1618 421 1876 422 2218 423 4528 424 6922 425 11338 426 15928 427 18442 428 22880 429 25378 430 29968 431 60050 432 117806 433 150118 434 209996 435 240232 436 59 437 318 438 324 439 904 440 1478 441 584 442 300 443 404 444 368 445 580 446 788 447 1132 448 1370 449 1630 450 1880 451 2272 452 4552 453 6926 454 11482 455 15946 456 18446 457 22886 458 25400 459 29990 460 60052 461 117808 462 175558 463 210172 464 240238 465 297988 466 61 467 336 468 332 469 992 470 1486 471 622 472 302 473 406 474 390 475 602 476 790 477 1136 478 1372 479 1640 480 1978 481 3778 482 4556 483 8788 484 11486 485 15950 486 18448 487 23032 488 25402 489 29996 490 60058 491 119902 492 177658 493 210176 494 267952 495 300268 496 330328 497 63 498 348 499 342 500 1024 501 1500 502 630 503 330 504 410 505 392 506 686 507 794 508 1138 509 1378 510 1658 511 2002 512 3958 513 4576 514 8812 515 11506 516 15952 517 18466 518 23036 519 27292 520 29998 521 85468 522 119906 523 177860 524 210178 525 267958 526 300290 527 360146 528 390358 529 65 530 366 531 354 532 1070 533 1620 534 682 535 358 536 442 537 472 538 688 539 796 540 1156 541 1394 542 1664 543 2006 544 3982 545 4580 546 8816 547 11510 548 16102 549 20158 550 23056 551 27296 552 30022 553 87568 554 119908 555 177862 556 210196 557 270056 558 300292 559 360322 560 390382 561 450238 562 67 563 378 564 364 565 1078 566 1910 567 728 568 386 569 446 570 476 571 706 572 812 573 1160 574 1396 575 1738 576 2024 577 3988 578 4582 579 8968 580 11512 581 16106 582 20338 583 23060 584 27448 585 30026 586 87746 587 120058 588 177866 589 210200 590 270058 591 300296 592 360326 593 390386 594 450440 595 480478 596 69 597 396 598 372 599 1082 600 1918 601 766 602 388 603 460 604 478 605 710 606 898 607 1184 608 1400 609 1762 610 2026 611 4138 612 4586 613 8998 614 11516 615 16124 616 20516 617 23062 618 27472 619 55228 620 87748 621 120086 622 177868 623 210206 624 270232 625 325708 626 360346 627 418106 628 450442 629 508198 630 1021016 631 71 632 408 633 384 634 1114 635 1920 636 836 637 398 638 464 639 496 640 712 641 920 642 1186 643 1402 644 1766 645 2030 646 4160 647 4588 648 9014 649 11530 650 16126 651 20518 652 23068 653 27476 654 55408 655 87766 656 120088 657 179962 658 210208 659 270236 660 327986 661 360356 662 420178 663 475858 664 510296 665 1021018 666 2039728 667 73 668 420 669 402 670 1202 671 2040 672 896 673 430 674 466 675 502 676 716 677 922 678 1210 679 1420 680 1768 681 2032 682 4162 683 4600 684 9022 685 11534 686 16130 687 20540 688 23084 689 27478 690 55432 691 87770 692 120104 693 179966 694 235616 695 270238 696 328010 697 360358 698 420202 699 477958 700 510298 701 1501496 702 2042008 703 3033028 704 75 705 432 706 414 707 1234 708 2104 709 910 710 434 711 470 712 518 713 730 714 944 715 1312 716 1424 717 1790 718 2048 719 4166 720 4604 721 9026 722 11536 723 16136 724 20542 725 23086 726 27496 727 55436 728 87772 729 120106 730 179968 731 235618 732 270256 733 328012 734 385768 735 420380 736 478136 737 510472 738 1529188 739 2042032 740 3062848 741 3573566 742 77 743 444 744 422 745 1264 746 2108 747 934 748 436 749 488 750 520 751 734 752 946 753 1316 754 1426 755 1792 756 2050 757 4168 758 6296 759 9146 760 11540 761 16144 762 20566 763 24778 764 27500 765 55438 766 87776 767 120110 768 180116 769 237716 770 295468 771 328016 772 387862 773 420386 774 478138 775 510476 776 1529212 777 2042036 778 3063026 779 4081762 780 4594558 781 79 782 450 783 426 784 1268 785 2116 786 980 787 458 788 490 789 524 790 736 791 950 792 1318 793 1442 794 1796 795 2054 796 4180 797 6298 798 9148 799 13406 800 16148 801 20570 802 24986 803 27502 804 57328 805 87778 806 120112 807 180136 808 237718 809 295648 810 328018 811 387866 812 420388 813 478162 814 510478 815 1529216 816 2520208 817 3063028 818 4081766 819 4594588 820 5105098 821 81 822 462 823 438 824 1276 825 2312 826 1018 827 482 828 494 829 526 830 740 831 952 832 1342 833 1526 834 1810 835 2056 836 4184 837 6472 838 9170 839 13408 840 16150 841 20572 842 24988 843 27506 844 57508 845 89666 846 120118 847 180140 848 237868 849 295672 850 329908 851 387868 852 420404 853 478166 854 510496 855 1529218 856 2522308 857 3063046 858 4081768 859 5075038 860 5585578 861 6606598 862 83 863 474 864 454 865 1292 866 2314 867 1124 868 556 869 512 870 548 871 754 872 964 873 1346 874 1528 875 1816 876 2072 877 4186 878 6496 879 9172 880 13426 881 16154 882 20576 883 25108 884 27508 885 57532 886 89668 887 145498 888 180142 889 237896 890 295676 891 330106 892 388018 893 420410 894 478168 895 510500 896 1531288 897 2522516 898 3063056 899 4083862 900 5075062 901 5615578 902 6634318 903 7147136 904 85 905 480 906 456 907 1324 908 2318 909 1162 910 566 911 514 912 652 913 758 914 968 915 1358 916 1550 917 1838 918 2078 919 4190 920 6500 921 9176 922 13430 923 16160 924 20692 925 25132 926 27656 927 57536 928 89842 929 145526 930 180146 931 237916 932 297772 933 330110 934 388042 935 420416 936 480236 937 510502 938 1531316 939 2522518 940 3513508 941 4083866 942 5075066 943 5615602 944 6636418 945 7147138 946 8168152 947 87 948 492 949 468 950 1412 951 2386 952 1208 953 724 954 536 955 656 956 760 957 970 958 1360 959 1552 960 1840 961 2080 962 4192 963 6502 964 9178 965 13432 966 17842 967 20698 968 25136 969 27658 970 57538 971 89846 972 145528 973 180160 974 237922 975 297776 976 330112 977 388046 978 420418 979 480266 980 510506 981 1531318 982 2547928 983 3543326 984 4083868 985 5100478 986 5615608 987 6636592 988 7627612 989 8648428 990 8678668 991 89 992 504 993 486 994 1444 995 2396 996 1230 997 742 998 608 999 674 1000 810 1001 974 1002 1364 1003 1556 1004 1844 1005 2156 1006 4196 1007 6506 1008 9190 1009 13436 1010 17846 1011 20716 1012 25138 1013 27674 1014 57686 1015 89866 1016 147626 1017 180164 1018 237926 1019 297778 1020 330118 1021 388048 1022 445618 1023 480416 1024 988678 1025 1531492 1026 2550028 1027 3543328 1028 4084042 1029 5102578 1030 6066058 1031 6636598 1032 7627616 1033 8648608 1034 9159148 1035 9699688 1036 91 1037 510 1038 498 1039 1470 1040 2440 1041 1232 1042 818 1043 610 1044 676 1045 862 1046 998 1047 1384 1048 1570 1049 1862 1050 2180 1051 4316 1052 6508 1053 9202 1054 13438 1055 17848 1056 20720 1057 25156 1058 27680 1059 57688 1060 89870 1061 147628 1062 180166 1063 239816 1064 297956 1065 330266 1066 388066 1067 447718 1068 480418 1069 990952 1070 1531496 1071 2550208 1072 3543502 1073 4084046 1074 5102758 1075 6093778 1076 6636616 1077 7627618 1078 8648632 1079 9186868 1080 18888868 1081 29099062