Anonymous user
Proper divisors: Difference between revisions
m
→version 3: re-arranged some statements for better reading.
m (→version 2: moved some statements.) |
m (→version 3: re-arranged some statements for better reading.) |
||
Line 3,724:
say right(n, max(20, w) ) 'has' center(#, 4) "proper divisors: " q
end /*n*/
m=0 /*M ≡ maximum number of Pdivs (so far).*/
do r=1 for range; q=
if #<m then iterate /*Less then max? Then ignore this #. */
@.#= @.# @. r;
end /*r*/ /* [↑] process 2nd range of integers.*/
say
say m ' is the highest number of proper divisors in range 1──►'range,
", and it's for: " subword(@.m, 3)
Line 3,745 ⟶ 3,744:
Pdivs: procedure; parse arg x 1 z,b; x= abs(x); if x==1 then return '' /*unity?*/
odd= x // 2; if x==0 then return '∞' /*zero ?*/
r= 0; q=
a=1 /* [↓] use all, or only odd #s. ___ */
do j=2 +odd by 1 +odd to r -(r*r==x) /*divide by some integers up to √ X */
| |||