# read("transforms") ; # read("transforms3") ; # # nimprodP2 := proc(nlog2::integer,mlog2::integer) # option remember; # local a2,b2,a,j ; # if ANDnos(nlog2,mlog2) = 0 then # return 2^(nlog2+mlog2) ; # else # a2 := convert(max(nlog2,mlog2),base,2) ; # b2 := convert(min(nlog2,mlog2),base,2) ; # a := 2^XORnos(nlog2,mlog2) ; # for j from 0 to nops(b2)-1 do # if op(j+1,a2) = 1 and op(j+1,b2)= 1 then # a := A051775(a,3*2^(2^j)/2 ) ; # end if; # end do: # a ; # end if; # end proc: # # A051775 := proc(n::integer,m::integer) # option remember; # local a,nbin,mbin,nlog2,mlog2 ; # if n*m = 0 then # return 0 ; # elif n = 1 or m=1 then # return n*m; # end if; # a := 0 ; # nbin := convert(n,base,2) ; # mbin := convert(m,base,2) ; # for nlog2 from 0 to nops(nbin)-1 do # if op(nlog2+1,nbin) > 0 then # for mlog2 from 0 to nops(mbin)-1 do # if op(mlog2+1,mbin) > 0 then # a := nimsum(a,nimprodP2(nlog2,mlog2)) ; # end if; # end do; # end if ; # end do; # return a; # end proc: # A058734 := proc(n) # A051775(n,n+1) ; # end proc: # L := [seq(A058734(n),n=0..1000)] ; # # LISTTOBFILE("b058734.txt",L,0) ; 0 0 1 2 2 1 3 12 4 2 5 8 6 3 7 15 8 5 9 11 10 4 11 14 12 7 13 1 14 6 15 240 16 8 17 42 18 9 19 100 20 10 21 32 22 11 23 231 24 13 25 35 26 12 27 102 28 15 29 41 30 14 31 124 32 20 33 38 34 21 35 168 36 22 37 44 38 23 39 107 40 17 41 47 42 16 43 170 44 19 45 37 46 18 47 184 48 28 49 14 50 29 51 192 52 30 53 4 54 31 55 131 56 25 57 7 58 24 59 194 60 27 61 13 62 26 63 214 64 39 65 165 66 38 67 203 68 37 69 175 70 36 71 120 72 34 73 172 74 35 75 201 76 32 77 166 78 33 79 204 80 47 81 141 82 46 83 163 84 45 85 135 86 44 87 144 88 42 89 132 90 43 91 161 92 40 93 142 94 41 95 205 96 51 97 129 98 50 99 111 100 49 101 139 102 48 103 28 104 54 105 136 106 55 107 109 108 52 109 130 110 53 111 132 112 59 113 169 114 58 115 7 116 57 117 163 118 56 119 244 120 62 121 160 122 63 123 5 124 60 125 170 126 61 127 210 128 94 129 156 130 95 131 34 132 92 133 150 134 93 135 241 136 91 137 149 138 90 139 32 140 89 141 159 142 88 143 131 144 86 145 180 146 87 147 74 148 84 149 190 150 85 151 25 152 83 153 189 154 82 155 72 156 81 157 183 158 80 159 201 160 74 161 184 162 75 163 134 164 72 165 178 166 73 167 149 168 79 169 177 170 78 171 132 172 77 173 187 174 76 175 203 176 66 177 144 178 67 179 238 180 64 181 154 182 65 183 125 184 71 185 153 186 70 187 236 188 69 189 147 190 68 191 218 192 121 193 59 194 120 195 229 196 123 197 49 198 122 199 134 200 124 201 50 202 125 203 231 204 126 205 56 206 127 207 191 208 113 209 19 210 112 211 141 212 115 213 25 214 114 215 110 216 116 217 26 218 117 219 143 220 118 221 16 222 119 223 120 224 109 225 31 226 108 227 65 228 111 229 21 230 110 231 226 232 104 233 22 234 105 235 67 236 106 237 28 238 107 239 247 240 101 241 55 242 100 243 41 244 103 245 61 246 102 247 10 248 96 249 62 250 97 251 43 252 98 253 52 254 99 255 65280 256 128 257 642 258 129 259 1676 260 130 261 648 262 131 263 3727 264 133 265 651 266 132 267 1678 268 135 269 641 270 134 271 7792 272 136 273 682 274 137 275 1764 276 138 277 672 278 139 279 3687 280 141 281 675 282 140 283 1766 284 143 285 681 286 142 287 16124 288 148 289 678 290 149 291 1576 292 150 293 684 294 151 295 3819 296 145 297 687 298 144 299 1578 300 147 301 677 302 146 303 7736 304 156 305 654 306 157 307 1600 308 158 309 644 310 159 311 3587 312 153 313 647 314 152 315 1602 316 155 317 653 318 154 319 32342 320 167 321 549 322 166 323 1611 324 165 325 559 326 164 327 3832 328 162 329 556 330 163 331 1609 332 160 333 550 334 161 335 7756 336 175 337 525 338 174 339 1571 340 173 341 519 342 172 343 3600 344 170 345 516 346 171 347 1569 348 168 349 526 350 169 351 15949 352 179 353 513 354 178 355 1775 356 177 357 523 358 176 359 3740 360 182 361 520 362 183 363 1773 364 180 365 514 366 181 367 7684 368 187 369 553 370 186 371 1671 372 185 373 547 374 184 375 3700 376 190 377 544 378 191 379 1669 380 188 381 554 382 189 383 65106 384 222 385 540 386 223 387 1698 388 220 389 534 390 221 391 3697 392 219 393 533 394 218 395 1696 396 217 397 543 398 216 399 7683 400 214 401 564 402 215 403 1738 404 212 405 574 406 213 407 3737 408 211 409 573 410 210 411 1736 412 209 413 567 414 208 415 15945 416 202 417 568 418 203 419 1542 420 200 421 562 422 201 423 3605 424 207 425 561 426 206 427 1540 428 205 429 571 430 204 431 7755 432 194 433 528 434 195 435 1646 436 192 437 538 438 193 439 3837 440 199 441 537 442 198 443 1644 444 197 445 531 446 196 447 32346 448 249 449 699 450 248 451 1637 452 251 453 689 454 250 455 3590 456 252 457 690 458 253 459 1639 460 254 461 696 462 255 463 7743 464 241 465 659 466 240 467 1549 468 243 469 665 470 242 471 3822 472 244 473 666 474 245 475 1551 476 246 477 656 478 247 479 16120 480 237 481 671 482 236 483 1729 484 239 485 661 486 238 487 3682 488 232 489 662 490 233 491 1731 492 234 493 668 494 235 495 7799 496 229 497 695 498 228 499 1705 500 231 501 701 502 230 503 3722 504 224 505 702 506 225 507 1707 508 226 509 692 510 227 511 22464 512 320 513 578 514 321 515 2636 516 322 517 584 518 323 519 1615 520 325 521 587 522 324 523 2638 524 327 525 577 526 326 527 9904 528 328 529 618 530 329 531 2596 532 330 533 608 534 331 535 1703 536 333 537 611 538 332 539 2598 540 335 541 617 542 334 543 5692 544 340 545 614 546 341 547 2792 548 342 549 620 550 343 551 1579 552 337 553 623 554 336 555 2794 556 339 557 613 558 338 559 9976 560 348 561 590 562 349 563 2688 564 350 565 580 566 351 567 1731 568 345 569 583 570 344 571 2690 572 347 573 589 574 346 575 38550 576 359 577 741 578 358 579 2699 580 357 581 751 582 356 583 1592 584 354 585 748 586 355 587 2697 588 352 589 742 590 353 591 9868 592 367 593 717 594 366 595 2787 596 365 597 711 598 364 599 1744 600 362 601 708 602 363 603 2785 604 360 605 718 606 361 607 5773 608 371 609 705 610 370 611 2607 612 369 613 715 614 368 615 1628 616 374 617 712 618 375 619 2605 620 372 621 706 622 373 623 9924 624 379 625 745 626 378 627 2631 628 377 629 739 630 376 631 1716 632 382 633 736 634 383 635 2629 636 380 637 746 638 381 639 22162 640 286 641 732 642 287 643 2658 644 284 645 726 646 285 647 1713 648 283 649 725 650 282 651 2656 652 281 653 735 654 280 655 9923 656 278 657 756 658 279 659 2570 660 276 661 766 662 277 663 1625 664 275 665 765 666 274 667 2568 668 273 669 759 670 272 671 5769 672 266 673 760 674 267 675 2758 676 264 677 754 678 265 679 1749 680 271 681 753 682 270 683 2756 684 269 685 763 686 268 687 9867 688 258 689 720 690 259 691 2734 692 256 693 730 694 257 695 1597 696 263 697 729 698 262 699 2732 700 261 701 723 702 260 703 38554 704 313 705 635 706 312 707 2725 708 315 709 625 710 314 711 1734 712 316 713 626 714 317 715 2727 716 318 717 632 718 319 719 9983 720 305 721 595 722 304 723 2765 724 307 725 601 726 306 727 1582 728 308 729 602 730 309 731 2767 732 310 733 592 734 311 735 5688 736 301 737 607 738 300 739 2561 740 303 741 597 742 302 743 1698 744 296 745 598 746 297 747 2563 748 298 749 604 750 299 751 9911 752 293 753 631 754 292 755 2665 756 295 757 637 758 294 759 1610 760 288 761 638 762 289 763 2667 764 290 765 628 766 291 767 43904 768 448 769 194 770 449 771 3276 772 450 773 200 774 451 775 2255 776 453 777 203 778 452 779 3278 780 455 781 193 782 454 783 14384 784 456 785 234 786 457 787 3236 788 458 789 224 790 459 791 2087 792 461 793 227 794 460 795 3238 796 463 797 233 798 462 799 10428 800 468 801 230 802 469 803 3176 804 470 805 236 806 471 807 2219 808 465 809 239 810 464 811 3178 812 467 813 229 814 466 815 14456 816 476 817 206 818 477 819 3072 820 478 821 196 822 479 823 2115 824 473 825 199 826 472 827 3074 828 475 829 205 830 474 831 59414 832 487 833 101 834 486 835 3083 836 485 837 111 838 484 839 2232 840 482 841 108 842 483 843 3081 844 480 845 102 846 481 847 14348 848 495 849 77 850 494 851 3171 852 493 853 71 854 492 855 2128 856 490 857 68 858 491 859 3169 860 488 861 78 862 489 863 10253 864 499 865 65 866 498 867 3247 868 497 869 75 870 496 871 2268 872 502 873 72 874 503 875 3245 876 500 877 66 878 501 879 14404 880 507 881 105 882 506 883 3271 884 505 885 99 886 504 887 2100 888 510 889 96 890 511 891 3269 892 508 893 106 894 509 895 43026 896 414 897 92 898 415 899 3298 900 412 901 86 902 413 903 2097 904 411 905 85 906 410 907 3296 908 409 909 95 910 408 911 14403 912 406 913 116 914 407 915 3210 916 404 917 126 918 405 919 2265 920 403 921 125 922 402 923 3208 924 401 925 119 926 400 927 10249 928 394 929 120 930 395 931 3142 932 392 933 114 934 393 935 2133 936 399 937 113 938 398 939 3140 940 397 941 123 942 396 943 14347 944 386 945 80 946 387 947 3118 948 384 949 90 950 385 951 2237 952 391 953 89 954 390 955 3116 956 389 957 83 958 388 959 59418 960 441 961 251 962 440 963 3109 964 443 965 241 966 442 967 2118 968 444 969 242 970 445 971 3111 972 446 973 248 974 447 975 14463 976 433 977 211 978 432 979 3149 980 435 981 217 982 434 983 2222 984 436 985 218 986 437 987 3151 988 438 989 208 990 439 991 10424 992 429 993 223 994 428 995 3201 996 431 997 213 998 430 999 2082 1000 424