Abundant odd numbers: Difference between revisions
Add MAD
Thundergnat (talk | contribs) m (Use wiki markup for OEIS link) |
Not a robot (talk | contribs) (Add MAD) |
||
Line 2,826:
1000: the proper divisors of 492975 sum to 519361
Above 1e6: the proper divisors of 1000125 sum to 1076547</pre>
=={{header|MAD}}==
<lang MAD> NORMAL MODE IS INTEGER
INTERNAL FUNCTION(ND)
ENTRY TO ODDSUM.
SUM = 1
SQN = SQRT.(ND)
THROUGH CHECK, FOR CN=3, 2, CN.G.SQN
TM = ND/CN
WHENEVER TM*CN.E.ND
SUM = SUM + CN
WHENEVER TM.NE.CN, SUM = SUM + TM
CHECK END OF CONDITIONAL
FUNCTION RETURN SUM
END OF FUNCTION
SEEN = 0
NUM = 1
THROUGH SHOW, FOR NUM=1, 2, SEEN.G.1000
WHENEVER NUM.L.ODDSUM.(NUM)
SEEN = SEEN + 1
WHENEVER SEEN.LE.25 .OR. SEEN.E.1000,
0 PRINT FORMAT OUTFMT,SEEN,NUM,ODDSUM.(NUM)
SHOW END OF CONDITIONAL
BILION THROUGH BILION, FOR NUM=NUM, 2,
0 NUM.G.1000000000 .AND. NUM.L.ODDSUM.(NUM)
PRINT FORMAT HUGENO,NUM,ODDSUM.(NUM)
VECTOR VALUES OUTFMT =
0 $4HNO. ,I4,S1,3HIS ,I6,S1,7HDIVSUM ,I6*$
VECTOR VALUES HUGENO =
0 $25HFIRST ABOVE 1 BILLION IS ,I10,S1,7HDIVSUM ,I10*$
END OF PROGRAM</lang>
{{out}}
<pre>NO. 1 IS 945 DIVSUM 975
NO. 2 IS 1575 DIVSUM 1649
NO. 3 IS 2205 DIVSUM 2241
NO. 4 IS 2835 DIVSUM 2973
NO. 5 IS 3465 DIVSUM 4023
NO. 6 IS 4095 DIVSUM 4641
NO. 7 IS 4725 DIVSUM 5195
NO. 8 IS 5355 DIVSUM 5877
NO. 9 IS 5775 DIVSUM 6129
NO. 10 IS 5985 DIVSUM 6495
NO. 11 IS 6435 DIVSUM 6669
NO. 12 IS 6615 DIVSUM 7065
NO. 13 IS 6825 DIVSUM 7063
NO. 14 IS 7245 DIVSUM 7731
NO. 15 IS 7425 DIVSUM 7455
NO. 16 IS 7875 DIVSUM 8349
NO. 17 IS 8085 DIVSUM 8331
NO. 18 IS 8415 DIVSUM 8433
NO. 19 IS 8505 DIVSUM 8967
NO. 20 IS 8925 DIVSUM 8931
NO. 21 IS 9135 DIVSUM 9585
NO. 22 IS 9555 DIVSUM 9597
NO. 23 IS 9765 DIVSUM 10203
NO. 24 IS 10395 DIVSUM 12645
NO. 25 IS 11025 DIVSUM 11946
NO. 1000 IS 492975 DIVSUM 519361
FIRST ABOVE 1 BILLION IS 1000000575 DIVSUM 1083561009</pre>
=={{header|Maple}}==
| |||