def sumfac(n): return sum([i for i in range(2,n+1) if ((1.*n/i)==int(n/i))]) #return sum([i for i in range(2,n+1) if (n/i)%1==0]) def abundant(n): return [i for i in range(1,n) if sumfac(i)>i] #f= open('abund.dat', 'w') for i in range(1000): #print>>f,i+1,sumfac(i) print i+1,sumfac(i)
Standard input is empty
1 0 2 0 3 2 4 3 5 6 6 5 7 11 8 7 9 14 10 12 11 17 12 11 13 27 14 13 15 23 16 23 17 30 18 17 19 38 20 19 21 41 22 31 23 35 24 23 25 59 26 30 27 41 28 39 29 55 30 29 31 71 32 31 33 62 34 47 35 53 36 47 37 90 38 37 39 59 40 55 41 89 42 41 43 95 44 43 45 83 46 77 47 71 48 47 49 123 50 56 51 92 52 71 53 97 54 53 55 119 56 71 57 119 58 79 59 89 60 59 61 167 62 61 63 95 64 103 65 126 66 83 67 143 68 67 69 125 70 95 71 143 72 71 73 194 74 73 75 113 76 123 77 139 78 95 79 167 80 79 81 185 82 120 83 125 84 83 85 223 86 107 87 131 88 119 89 179 90 89 91 233 92 111 93 167 94 127 95 143 96 119 97 251 98 97 99 170 100 155 101 216 102 101 103 215 104 103 105 209 106 191 107 161 108 107 109 279 110 109 111 215 112 151 113 247 114 113 115 239 116 143 117 209 118 181 119 179 120 143 121 359 122 132 123 185 124 167 125 223 126 155 127 311 128 127 129 254 130 175 131 251 132 131 133 335 134 159 135 203 136 239 137 269 138 137 139 287 140 139 141 335 142 191 143 215 144 167 145 402 146 179 147 221 148 227 149 265 150 149 151 371 152 151 153 299 154 233 155 287 156 191 157 391 158 157 159 239 160 215 161 377 162 191 163 362 164 163 165 293 166 287 167 251 168 167 169 479 170 182 171 323 172 259 173 307 174 173 175 359 176 247 177 371 178 239 179 269 180 179 181 545 182 181 183 335 184 247 185 359 186 227 187 383 188 215 189 335 190 319 191 359 192 191 193 507 194 193 195 293 196 335 197 398 198 197 199 467 200 199 201 464 202 271 203 305 204 239 205 503 206 251 207 311 208 311 209 433 210 239 211 575 212 211 213 377 214 287 215 323 216 263 217 599 218 255 219 329 220 295 221 503 222 251 223 455 224 223 225 503 226 402 227 341 228 227 229 559 230 229 231 431 232 383 233 449 234 233 235 545 236 287 237 419 238 319 239 431 240 239 241 743 242 241 243 398 244 363 245 433 246 341 247 503 248 279 249 479 250 335 251 467 252 251 253 727 254 287 255 383 256 431 257 510 258 257 259 527 260 303 261 587 262 389 263 395 264 263 265 719 266 323 267 479 268 359 269 475 270 269 271 719 272 271 273 557 274 447 275 413 276 371 277 671 278 277 279 419 280 415 281 719 282 281 283 575 284 283 285 503 286 479 287 503 288 335 289 818 290 306 291 539 292 391 293 517 294 293 295 683 296 359 297 569 298 479 299 449 300 335 301 867 302 351 303 455 304 407 305 619 306 371 307 701 308 307 309 671 310 415 311 575 312 311 313 839 314 313 315 473 316 623 317 559 318 317 319 647 320 359 321 761 322 431 323 575 324 359 325 846 326 433 327 491 328 439 329 629 330 383 331 863 332 331 333 587 334 493 335 503 336 407 337 991 338 337 339 548 340 455 341 755 342 383 343 779 344 399 345 659 346 575 347 521 348 347 349 839 350 349 351 743 352 559 353 755 354 353 355 719 356 431 357 629 358 575 359 539 360 359 361 1169 362 380 363 545 364 531 365 783 366 443 367 743 368 367 369 743 370 545 371 683 372 431 373 895 374 373 375 647 376 623 377 719 378 419 379 959 380 379 381 839 382 511 383 575 384 383 385 1019 386 575 387 581 388 571 389 685 390 389 391 1007 392 431 393 854 394 527 395 593 396 479 397 1091 398 397 399 599 400 639 401 960 402 401 403 815 404 447 405 713 406 725 407 719 408 455 409 1079 410 409 411 755 412 551 413 727 414 479 415 935 416 503 417 881 418 559 419 719 420 419 421 1343 422 421 423 635 424 623 425 809 426 557 427 863 428 495 429 755 430 671 431 791 432 431 433 1239 434 433 435 767 436 719 437 769 438 479 439 887 440 439 441 1079 442 740 443 755 444 443 445 1063 446 539 447 671 448 599 449 1015 450 449 451 1208 452 503 453 797 454 607 455 683 456 671 457 1199 458 457 459 689 460 719 461 1007 462 461 463 1151 464 463 465 929 466 767 467 701 468 467 469 1273 470 543 471 863 472 631 473 899 474 527 475 959 476 619 477 1007 478 701 479 719 480 479 481 1511 482 531 483 725 484 767 485 930 486 587 487 1091 488 487 489 929 490 655 491 1025 492 491 493 1175 494 539 495 839 496 935 497 991 498 575 499 1007 500 499 501 1091 502 671 503 755 504 503 505 1559 506 611 507 863 508 731 509 895 510 509 511 1295 512 591 513 1022 514 799 515 773 516 623 517 1231 518 575 519 911 520 695 521 1259 522 521 523 1169 524 523 525 923 526 991 527 791 528 575 529 1487 530 552 531 971 532 779 533 1119 534 587 535 1079 536 647 537 1019 538 719 539 809 540 683 541 1679 542 541 543 815 544 727 545 1133 546 659 547 1343 548 547 549 965 550 805 551 1115 552 599 553 1439 554 639 555 833 556 911 557 979 558 557 559 1247 560 615 561 1487 562 863 563 845 564 563 565 1343 566 683 567 851 568 967 569 1079 570 569 571 1439 572 571 573 1175 574 767 575 1007 576 743 577 1650 578 577 579 920 580 775 581 1259 582 671 583 1175 584 647 585 1109 586 1091 587 881 588 587 589 1595 590 639 591 1079 592 791 593 1177 594 593 595 1439 596 863 597 1049 598 799 599 1007 600 599 601 1859 602 601 603 1055 604 883 605 1063 606 797 607 1223 608 607 609 1259 610 959 611 1115 612 671 613 1637 614 613 615 923 616 1007 617 1439 618 617 619 1247 620 619 621 1343 622 959 623 935 624 719 625 1735 626 780 627 941 628 959 629 1105 630 683 631 1871 632 631 633 1199 634 847 635 953 636 767 637 1511 638 797 639 1079 640 935 641 1529 642 641 643 1295 644 643 645 1343 646 1055 647 1079 648 647 649 1814 650 719 651 1301 652 1023 653 1147 654 653 655 1319 656 791 657 1301 658 961 659 1151 660 659 661 2015 662 661 663 995 664 1007 665 1259 666 959 667 1481 668 719 669 1175 670 895 671 1223 672 743 673 2015 674 673 675 1013 676 1239 677 1280 678 677 679 1367 680 783 681 1619 682 911 683 1151 684 683 685 1819 686 827 687 1199 688 919 689 1363 690 755 691 1727 692 691 693 1217 694 1247 695 1043 696 839 697 1799 698 755 699 1049 700 935 701 1735 702 701 703 1679 704 759 705 1523 706 1151 707 1061 708 815 709 1679 710 709 711 1295 712 1039 713 1349 714 767 715 1727 716 1007 717 1259 718 959 719 1079 720 719 721 2417 722 831 723 1142 724 967 725 1273 726 929 727 1595 728 727 729 1679 730 1092 731 1331 732 791 733 1735 734 733 735 1103 736 1367 737 1511 738 815 739 1637 740 739 741 1595 742 1119 743 1295 744 743 745 1919 746 899 747 1121 748 1091 749 1511 750 863 751 1871 752 751 753 1487 754 1007 755 1259 756 911 757 2239 758 757 759 1139 760 1151 761 1799 762 761 763 1535 764 879 765 1343 766 1403 767 1151 768 839 769 2043 770 769 771 1727 772 1031 773 1357 774 773 775 1715 776 991 777 1469 778 1215 779 1169 780 839 781 2351 782 863 783 1295 784 1199 785 1766 786 947 787 1583 788 787 789 1385 790 1055 791 1439 792 911 793 2339 794 867 795 1193 796 1295 797 1399 798 797 799 1919 800 863 801 1952 802 1169 803 1205 804 887 805 1903 806 1151 807 1343 808 1079 809 1529 810 809 811 2177 812 811 813 1679 814 1087 815 1367 816 983 817 2231 818 879 819 1229 820 1455 821 1763 822 821 823 1655 824 823 825 1559 826 1487 827 1439 828 827 829 2183 830 829 831 1511 832 1111 833 1777 834 1025 835 1679 836 1007 837 1679 838 1279 839 1259 840 839 841 2879 842 870 843 1265 844 1127 845 1483 846 1097 847 1871 848 1063 849 1673 850 1135 851 1673 852 911 853 2015 854 853 855 1487 856 1559 857 1619 858 857 859 2015 860 859 861 1847 862 1343 863 1295 864 863 865 2519 866 1043 867 1301 868 1227 869 1791 870 959 871 2159 872 951 873 1649 874 1273 875 1439 876 1247 877 2071 878 877 879 1319 880 1175 881 2231 882 881 883 2222 884 883 885 1763 886 1439 887 1331 888 887 889 2279 890 1023 891 1619 892 1451 893 1567 894 959 895 1799 896 1079 897 2039 898 1343 899 1349 900 959 901 2820 902 971 903 1511 904 1407 905 1709 906 1091 907 1823 908 907 909 1595 910 1325 911 2015 912 911 913 2479 914 1007 915 1373 916 1487 917 1609 918 1055 919 2159 920 919 921 2159 922 1231 923 1385 924 1007 925 2687 926 1177 927 1391 928 1351 929 1889 930 929 931 2303 932 1139 933 1637 934 1247 935 1403 936 1295 937 2729 938 937 939 1631 940 1255 941 2015 942 941 943 1895 944 1007 945 1859 946 1919 947 1583 948 947 949 2239 950 1035 951 1859 952 1271 953 2159 954 953 955 2105 956 1151 957 1679 958 1439 959 1439 960 1103 961 3047 962 992 963 1595 964 1403 965 1693 966 1163 967 2303 968 967 969 1994 970 1439 971 1763 972 971 973 2547 974 1119 975 1463 976 1735 977 1921 978 977 979 1967 980 1079 981 2393 982 1429 983 1475 984 983 985 2519 986 1187 987 1619 988 1535 989 1959 990 1055 991 2807 992 991 993 2015 994 1327 995 1727 996 1199 997 2351 998 997 999 1499 1000 1519