# 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: # A051776 := proc(n,m) # A051775(n,m) ; # end proc: # L := [] ; # for d from 1 to 60 do # for y from d to 1 by -1 do # L := [op(L), A051776(d-y+1,y)] ; # end do; # end do: # # LISTTOBFILE("b051776.txt",L,1) ; 1 1 2 2 3 2 4 3 5 3 6 3 7 4 8 1 9 1 10 4 11 5 12 8 13 2 14 8 15 5 16 6 17 10 18 12 19 12 20 10 21 6 22 7 23 11 24 15 25 6 26 15 27 11 28 7 29 8 30 9 31 13 32 2 33 2 34 13 35 9 36 8 37 9 38 12 39 14 40 14 41 7 42 14 43 14 44 12 45 9 46 10 47 14 48 4 49 10 50 8 51 8 52 10 53 4 54 14 55 10 56 11 57 15 58 7 59 11 60 13 61 5 62 13 63 11 64 7 65 15 66 11 67 12 68 13 69 5 70 15 71 3 72 3 73 3 74 3 75 15 76 5 77 13 78 12 79 13 80 4 81 6 82 3 83 6 84 7 85 4 86 7 87 6 88 3 89 6 90 4 91 13 92 14 93 6 94 8 95 7 96 9 97 1 98 15 99 15 100 1 101 9 102 7 103 8 104 6 105 14 106 15 107 7 108 11 109 13 110 12 111 12 112 8 113 13 114 8 115 12 116 12 117 13 118 11 119 7 120 15 121 16 122 5 123 9 124 9 125 1 126 10 127 6 128 5 129 5 130 6 131 10 132 1 133 9 134 9 135 5 136 16 137 17 138 32 139 10 140 5 141 4 142 9 143 1 144 1 145 12 146 1 147 1 148 9 149 4 150 5 151 10 152 32 153 17 154 18 155 34 156 48 157 1 158 11 159 15 160 5 161 9 162 11 163 11 164 9 165 5 166 15 167 11 168 1 169 48 170 34 171 18 172 19 173 35 174 51 175 64 176 14 177 2 178 2 179 6 180 2 181 14 182 2 183 6 184 2 185 2 186 14 187 64 188 51 189 35 190 19 191 20 192 33 193 49 194 68 195 80 196 4 197 12 198 14 199 10 200 4 201 4 202 10 203 14 204 12 205 4 206 80 207 68 208 49 209 33 210 20 211 21 212 40 213 50 214 72 215 85 216 96 217 11 218 10 219 3 220 2 221 15 222 2 223 3 224 10 225 11 226 96 227 85 228 72 229 50 230 40 231 21 232 22 233 42 234 60 235 76 236 90 237 102 238 112 239 2 240 4 241 8 242 14 243 14 244 8 245 4 246 2 247 112 248 102 249 90 250 76 251 60 252 42 253 22 254 23 255 43 256 63 257 70 258 95 259 107 260 119 261 128 262 13 263 13 264 5 265 11 266 5 267 13 268 13 269 128 270 119 271 107 272 95 273 70 274 63 275 43 276 23 277 24 278 41 279 61 280 66 281 82 282 109 283 121 284 136 285 144 286 7 287 3 288 7 289 7 290 3 291 7 292 144 293 136 294 121 295 109 296 82 297 66 298 61 299 41 300 24 301 25 302 44 303 62 304 78 305 87 306 110 307 126 308 140 309 153 310 160 311 8 312 15 313 10 314 15 315 8 316 160 317 153 318 140 319 126 320 110 321 87 322 78 323 62 324 44 325 25 326 26 327 46 328 52 329 74 330 88 331 104 332 122 333 132 334 158 335 170 336 176 337 3 338 1 339 1 340 3 341 176 342 170 343 158 344 132 345 122 346 104 347 88 348 74 349 52 350 46 351 26 352 27 353 47 354 55 355 75 356 93 357 101 358 125 359 139 360 151 361 175 362 187 363 192 364 12 365 8 366 12 367 192 368 187 369 175 370 151 371 139 372 125 373 101 374 93 375 75 376 55 377 47 378 27 379 28 380 45 381 53 382 79 383 83 384 99 385 115 386 131 387 159 388 165 389 189 390 204 391 208 392 6 393 6 394 208 395 204 396 189 397 165 398 159 399 131 400 115 401 99 402 83 403 79 404 53 405 45 406 28 407 29 408 36 409 54 410 67 411 86 412 103 413 116 414 135 415 150 416 163 417 182 418 196 419 221 420 224 421 9 422 224 423 221 424 196 425 182 426 163 427 150 428 135 429 116 430 103 431 86 432 67 433 54 434 36 435 29 436 30 437 38 438 56 439 71 440 89 441 97 442 127 443 143 444 145 445 169 446 183 447 200 448 214 449 238 450 240 451 240 452 238 453 214 454 200 455 183 456 169 457 145 458 143 459 127 460 97 461 89 462 71 463 56 464 38 465 30 466 31 467 39 468 59 469 77 470 92 471 108 472 120 473 141 474 152 475 172 476 188 477 205 478 219 479 231 480 255 481 24 482 255 483 231 484 219 485 205 486 188 487 172 488 152 489 141 490 120 491 108 492 92 493 77 494 59 495 39 496 31 497 32 498 37 499 57 500 73 501 81 502 106 503 118 504 133 505 149 506 166 507 186 508 193 509 217 510 233 511 245 512 8 513 8 514 245 515 233 516 217 517 193 518 186 519 166 520 149 521 133 522 118 523 106 524 81 525 73 526 57 527 37 528 32 529 33 530 48 531 58 532 69 533 84 534 105 535 113 536 129 537 156 538 161 539 177 540 201 541 212 542 229 543 250 544 56 545 25 546 56 547 250 548 229 549 212 550 201 551 177 552 161 553 156 554 129 555 113 556 105 557 84 558 69 559 58 560 48 561 33 562 34 563 50 564 16 565 65 566 91 567 111 568 117 569 137 570 155 571 171 572 185 573 197 574 223 575 235 576 241 577 40 578 42 579 42 580 40 581 241 582 235 583 223 584 197 585 185 586 171 587 155 588 137 589 117 590 111 591 91 592 65 593 16 594 50 595 34 596 35 597 51 598 19 599 128 600 94 601 98 602 114 603 134 604 146 605 174 606 178 607 198 608 210 609 226 610 254 611 88 612 59 613 27 614 59 615 88 616 254 617 226 618 210 619 198 620 178 621 174 622 146 623 134 624 114 625 98 626 94 627 128 628 19 629 51 630 35 631 36 632 49 633 17 634 132 635 160 636 100 637 124 638 142 639 154 640 164 641 180 642 202 643 222 644 236 645 244 646 72 647 76 648 9 649 9 650 76 651 72 652 244 653 236 654 222 655 202 656 180 657 164 658 154 659 142 660 124 661 100 662 160 663 132 664 17 665 49 666 36 667 37 668 56 669 18 670 136 671 165 672 176 673 123 674 138 675 147 676 162 677 191 678 194 679 211 680 234 681 251 682 120 683 93 684 112 685 26 686 112 687 93 688 120 689 251 690 234 691 211 692 194 693 191 694 162 695 147 696 138 697 123 698 176 699 165 700 136 701 18 702 56 703 37 704 38 705 58 706 28 707 140 708 170 709 182 710 144 711 130 712 148 713 168 714 190 715 206 716 216 717 228 718 242 719 104 720 110 721 98 722 100 723 100 724 98 725 110 726 104 727 242 728 228 729 216 730 206 731 190 732 168 733 148 734 130 735 144 736 182 737 170 738 140 739 28 740 58 741 38 742 39 743 59 744 31 745 134 746 175 747 187 748 151 749 192 750 157 751 173 752 181 753 203 754 213 755 237 756 253 757 152 758 127 759 83 760 119 761 30 762 119 763 83 764 127 765 152 766 253 767 237 768 213 769 203 770 181 771 173 772 157 773 192 774 151 775 187 776 175 777 134 778 31 779 59 780 39 781 40 782 57 783 29 784 130 785 162 786 189 787 153 788 200 789 224 790 167 791 179 792 199 793 215 794 227 795 247 796 136 797 128 798 65 799 69 800 10 801 10 802 69 803 65 804 128 805 136 806 247 807 227 808 215 809 199 810 179 811 167 812 224 813 200 814 153 815 189 816 162 817 130 818 29 819 57 820 40 821 41 822 60 823 30 824 142 825 167 826 190 827 158 828 204 829 233 830 240 831 184 832 207 833 218 834 239 835 248 836 184 837 145 838 180 839 86 840 54 841 31 842 54 843 86 844 180 845 145 846 184 847 248 848 239 849 218 850 207 851 184 852 240 853 233 854 204 855 158 856 190 857 167 858 142 859 30 860 60 861 41 862 42 863 62 864 20 865 138 866 168 867 184 868 154 869 196 870 238 871 250 872 208 873 195 874 209 875 225 876 243 877 168 878 162 879 166 880 172 881 34 882 32 883 32 884 34 885 172 886 166 887 162 888 168 889 243 890 225 891 209 892 195 893 208 894 250 895 238 896 196 897 154 898 184 899 168 900 138 901 20 902 62 903 42 904 43 905 63 906 23 907 139 908 173 909 181 910 157 911 203 912 231 913 255 914 219 915 64 916 220 917 232 918 252 919 216 920 179 921 151 922 191 923 211 924 53 925 29 926 53 927 211 928 191 929 151 930 179 931 216 932 252 933 232 934 220 935 64 936 219 937 255 938 231 939 203 940 157 941 181 942 173 943 139 944 23 945 63 946 43 947 44 948 61 949 21 950 143 951 163 952 179 953 147 954 195 955 239 956 245 957 221 958 76 959 96 960 230 961 246 962 200 963 196 964 133 965 141 966 199 967 203 968 11 969 11 970 203 971 199 972 141 973 133 974 196 975 200 976 246 977 230 978 96 979 76 980 221 981 245 982 239 983 195 984 147 985 179 986 163 987 143 988 21 989 61 990 44 991 45 992 52 993 22 994 131 995 166 996 183 997 148 998 199 999 230 1000 243 1001 214 1002 68 1003 109 1004 112 1005 249 1006 248 1007 213 1008 252 1009 158 1010 251 1011 222 1012 255 1013 28 1014 255 1015 222 1016 251 1017 158 1018 252 1019 213 1020 248 1021 249 1022 112 1023 109 1024 68 1025 214 1026 243 1027 230 1028 199 1029 148 1030 183 1031 166 1032 131 1033 22 1034 52 1035 45 1036 46 1037 54 1038 24 1039 135 1040 169 1041 177 1042 159 1043 207 1044 225 1045 249 1046 215 1047 72 1048 102 1049 126 1050 80 1051 232 1052 230 1053 238 1054 224 1055 239 1056 225 1057 233 1058 231 1059 231 1060 233 1061 225 1062 239 1063 224 1064 238 1065 230 1066 232 1067 80 1068 126 1069 102 1070 72 1071 215 1072 249 1073 225 1074 207 1075 159 1076 177 1077 169 1078 135 1079 24 1080 54 1081 46 1082 47 1083 55 1084 27 1085 141 1086 172 1087 188 1088 152 1089 205 1090 232 1091 252 1092 220 1093 77 1094 107 1095 119 1096 95 1097 44 1098 247 1099 223 1100 243 1101 149 1102 244 1103 212 1104 240 1105 21 1106 240 1107 212 1108 244 1109 149 1110 243 1111 223 1112 247 1113 44 1114 95 1115 119 1116 107 1117 77 1118 220 1119 252 1120 232 1121 205 1122 152 1123 188 1124 172 1125 141 1126 27 1127 55 1128 47 1129 48 1130 53 1131 25 1132 137 1133 161 1134 186 1135 150 1136 197 1137 229 1138 246 1139 218 1140 65 1141 105 1142 121 1143 85 1144 60 1145 12 1146 205 1147 193 1148 129 1149 137 1150 194 1151 206 1152 13 1153 13 1154 206 1155 194 1156 137 1157 129 1158 193 1159 205 1160 12 1161 60 1162 85 1163 121 1164 105 1165 65 1166 218 1167 246 1168 229 1169 197 1170 150 1171 186 1172 161 1173 137 1174 25 1175 53 1176 48 1177 49 1178 16 1179 26 1180 133 1181 164 1182 185 1183 145 1184 193 1185 236 1186 241 1187 209 1188 73 1189 100 1190 117 1191 90 1192 12 1193 29 1194 28 1195 210 1196 189 1197 156 1198 177 1199 217 1200 57 1201 20 1202 57 1203 217 1204 177 1205 156 1206 189 1207 210 1208 28 1209 29 1210 12 1211 90 1212 117 1213 100 1214 73 1215 209 1216 241 1217 236 1218 193 1219 145 1220 185 1221 164 1222 133 1223 26 1224 16 1225 49 1226 50 1227 18 1228 32 1229 129 1230 171 1231 191 1232 149 1233 201 1234 235 1235 251 1236 217 1237 69 1238 111 1239 123 1240 81 1241 28 1242 46 1243 14 1244 60 1245 169 1246 163 1247 167 1248 173 1249 33 1250 35 1251 35 1252 33 1253 173 1254 167 1255 163 1256 169 1257 60 1258 14 1259 46 1260 28 1261 81 1262 123 1263 111 1264 69 1265 217 1266 251 1267 235 1268 201 1269 149 1270 191 1271 171 1272 129 1273 32 1274 18 1275 50 1276 51 1277 19 1278 35 1279 192 1280 174 1281 178 1282 146 1283 198 1284 226 1285 254 1286 210 1287 70 1288 98 1289 114 1290 94 1291 108 1292 63 1293 63 1294 47 1295 172 1296 182 1297 154 1298 186 1299 94 1300 58 1301 22 1302 58 1303 94 1304 186 1305 154 1306 182 1307 172 1308 47 1309 63 1310 63 1311 108 1312 94 1313 114 1314 98 1315 70 1316 210 1317 254 1318 226 1319 198 1320 146 1321 178 1322 174 1323 192 1324 35 1325 19 1326 51 1327 52 1328 17 1329 33 1330 196 1331 240 1332 180 1333 156 1334 206 1335 234 1336 244 1337 212 1338 74 1339 110 1340 124 1341 84 1342 124 1343 72 1344 45 1345 29 1346 184 1347 140 1348 140 1349 132 1350 70 1351 66 1352 12 1353 12 1354 66 1355 70 1356 132 1357 140 1358 140 1359 184 1360 29 1361 45 1362 72 1363 124 1364 84 1365 124 1366 110 1367 74 1368 212 1369 244 1370 234 1371 206 1372 156 1373 180 1374 240 1375 196 1376 33 1377 17 1378 52 1379 53 1380 24 1381 34 1382 200 1383 245 1384 208 1385 155 1386 202 1387 227 1388 242 1389 223 1390 66 1391 99 1392 122 1393 91 1394 76 1395 89 1396 84 1397 14 1398 132 1399 153 1400 156 1401 147 1402 114 1403 91 1404 122 1405 23 1406 122 1407 91 1408 114 1409 147 1410 156 1411 153 1412 132 1413 14 1414 84 1415 89 1416 76 1417 91 1418 122 1419 99 1420 66 1421 223 1422 242 1423 227 1424 202 1425 155 1426 208 1427 245 1428 200 1429 34 1430 24 1431 53 1432 54 1433 26 1434 44 1435 204 1436 250 1437 214 1438 224 1439 194 1440 228 1441 248 1442 222 1443 78 1444 104 1445 116 1446 82 1447 92 1448 106 1449 70 1450 112 1451 144 1452 166 1453 138 1454 188 1455 106 1456 108 1457 96 1458 102 1459 102 1460 96 1461 108 1462 106 1463 188 1464 138 1465 166 1466 144 1467 112 1468 70 1469 106 1470 92 1471 82 1472 116 1473 104 1474 78 1475 222 1476 248 1477 228 1478 194 1479 224 1480 214 1481 250 1482 204 1483 44 1484 26 1485 54 1486 55 1487 27 1488 47 1489 198 1490 255 1491 219 1492 231 1493 64 1494 237 1495 253 1496 213 1497 75 1498 101 1499 125 1500 93 1501 172 1502 123 1503 119 1504 99 1505 234 1506 179 1507 183 1508 171 1509 236 1510 117 1511 85 1512 125 1513 19 1514 125 1515 85 1516 117 1517 236 1518 171 1519 183 1520 179 1521 234 1522 99 1523 119 1524 123 1525 172 1526 93 1527 125 1528 101 1529 75 1530 213 1531 253 1532 237 1533 64 1534 231 1535 219 1536 255 1537 198 1538 47 1539 27 1540 55 1541 56 1542 25 1543 45 1544 194 1545 242 1546 221 1547 233 1548 72 1549 112 1550 247 1551 211 1552 71 1553 103 1554 115 1555 87 1556 188 1557 132 1558 101 1559 81 1560 254 1561 206 1562 161 1563 149 1564 244 1565 204 1566 79 1567 75 1568 15 1569 15 1570 75 1571 79 1572 204 1573 244 1574 149 1575 161 1576 206 1577 254 1578 81 1579 101 1580 132 1581 188 1582 87 1583 115 1584 103 1585 71 1586 211 1587 247 1588 112 1589 72 1590 233 1591 221 1592 242 1593 194 1594 45 1595 25 1596 56 1597 57 1598 28 1599 46 1600 206 1601 247 1602 222 1603 238 1604 76 1605 121 1606 80 1607 216 1608 79 1609 106 1610 127 1611 88 1612 140 1613 149 1614 144 1615 66 1616 194 1617 219 1618 210 1619 130 1620 192 1621 213 1622 220 1623 80 1624 55 1625 18 1626 55 1627 80 1628 220 1629 213 1630 192 1631 130 1632 210 1633 219 1634 194 1635 66 1636 144 1637 149 1638 140 1639 88 1640 127 1641 106 1642 79 1643 216 1644 80 1645 121 1646 76 1647 238 1648 222 1649 247 1650 206 1651 46 1652 28 1653 57 1654 58 1655 30 1656 36 1657 202 1658 248 1659 216 1660 234 1661 68 1662 126 1663 90 1664 96 1665 67 1666 97 1667 113 1668 83 1669 156 1670 166 1671 130 1672 184 1673 214 1674 228 1675 196 1676 246 1677 216 1678 226 1679 198 1680 252 1681 43 1682 41 1683 41 1684 43 1685 252 1686 198 1687 226 1688 216 1689 246 1690 196 1691 228 1692 214 1693 184 1694 130 1695 166 1696 156 1697 83 1698 113 1699 97 1700 67 1701 96 1702 90 1703 126 1704 68 1705 234 1706 216 1707 248 1708 202 1709 36 1710 30 1711 58 1712 59 1713 31 1714 39 1715 203 1716 253 1717 213 1718 237 1719 75 1720 119 1721 95 1722 107 1723 128 1724 108 1725 120 1726 92 1727 236 1728 183 1729 179 1730 171 1731 39 1732 241 1733 249 1734 225 1735 167 1736 251 1737 243 1738 231 1739 108 1740 52 1741 16 1742 52 1743 108 1744 231 1745 243 1746 251 1747 167 1748 225 1749 249 1750 241 1751 39 1752 171 1753 179 1754 183 1755 236 1756 92 1757 120 1758 108 1759 128 1760 107 1761 95 1762 119 1763 75 1764 237 1765 213 1766 253 1767 203 1768 39 1769 31 1770 59 1771 60 1772 29 1773 37 1774 207 1775 243 1776 211 1777 227 1778 67 1779 127 1780 85 1781 109 1782 140 1783 176 1784 118 1785 86 1786 252 1787 192 1788 161 1789 153 1790 51 1791 15 1792 239 1793 223 1794 191 1795 131 1796 233 1797 209 1798 112 1799 76 1800 14 1801 14 1802 76 1803 112 1804 209 1805 233 1806 131 1807 191 1808 223 1809 239 1810 15 1811 51 1812 153 1813 161 1814 192 1815 252 1816 86 1817 118 1818 176 1819 140 1820 109 1821 85 1822 127 1823 67 1824 227 1825 211 1826 243 1827 207 1828 37 1829 29 1830 60