Aliquot sequence classifications: Difference between revisions

Content added Content deleted

Inline

@@ Line 146: / Line 146: @@
  non-terminating 15355717786080 44534663601120 144940087464480 471714103310688 1130798979186912 2688948041357088 6050151708497568 13613157922639968 35513546724070632 74727605255142168 162658586225561832 353930992506879768 642678347124409032 112510261154846...
 </lang>
+=={{header|jq}}==
+{{works with|jq|1.4}}
+<lang jq># "until" is available in more recent versions of jq
+# than jq 1.4
+def until(cond; next):
+  def _until:
+    if cond then . else (next|_until) end;
+  _until;
+# unordered
+def proper_divisors:
+  . as $n
+  | if $n > 1 then 1,
+      ( range(2; 1 + (sqrt|floor)) as $i
+        | if ($n % $i) == 0 then $i,
+            (($n / $i) | if . == $i then empty else . end)
+          else empty
+          end)
+    else empty
+    end;
+# sum of proper divisors, or 0
+def pdsum:
+  [proper_divisors] | add // 0;
+# input is n
+# maxlen defaults to 16;
+# maxterm defaults to 2^47
+def aliquot(maxlen; maxterm):
+  (maxlen // 15) as $maxlen
+  | (maxterm // 40737488355328) as $maxterm
+  | if . == 0 then "terminating at 0"
+    else
+    # [s, slen, new] = [[n], 1, n]
+    [ [.], 1, .]
+    | until( type == "string" or .[1] > $maxlen or .[2] > $maxterm;
+             .[0] as $s | .[1] as $slen
+             | ($s | .[length-1] | pdsum) as $new
+             | if ($s|index($new)) then
+                 if $s[0] == $new then
+                     if $slen == 1 then "perfect \($s)"
+                     elif $slen == 2 then "amicable: \($s)"
+                     else "sociable of length \($slen): \($s)"
+		     end
+                 elif ($s | .[length-1]) == $new then "aspiring: \($s)"
+                 else "cyclic back to \($new): \($s)"
+		 end
+               elif $new == 0 then "terminating: \($s + [0])"
+               else [ ($s + [$new]), ($slen + 1), $new ]
+               end )
+    | if type == "string" then . else "non-terminating: \(.[0])" end
+    end;
+def task:
+  def pp: "\(.): \(aliquot(null;null))";
+     (range(1; 11) | pp),
+     "",
+     ((11, 12, 28, 496, 220, 1184, 12496, 1264460,
+, 909, 562, 1064, 1488, 15355717786080) | pp);
+task</lang>
+{{out}}
+<lang sh>$ jq -n -r -f aliquot.jq
+: terminating: [1,0]
+: terminating: [2,1,0]
+: terminating: [3,1,0]
+: terminating: [4,3,1,0]
+: terminating: [5,1,0]
+: perfect [6]
+: terminating: [7,1,0]
+: terminating: [8,7,1,0]
+: terminating: [9,4,3,1,0]
+: terminating: [10,8,7,1,0]
+: terminating: [11,1,0]
+: terminating: [12,16,15,9,4,3,1,0]
+: perfect [28]
+: perfect [496]
+: amicable: [220,284]
+: amicable: [1184,1210]
+: sociable of length 5: [12496,14288,15472,14536,14264]
+1264460: sociable of length 4: [1264460,1547860,1727636,1305184]
+: aspiring: [790,650,652,496]
+: aspiring: [909,417,143,25,6]
+: cyclic back to 284: [562,284,220]
+: cyclic back to 1184: [1064,1336,1184,1210]
+: non-terminating: [1488,2480,3472,4464,8432,9424,10416,21328,22320,55056,95728,96720,236592,459792,881392,882384]
+15355717786080: non-terminating: [15355717786080,44534663601120]</lang>
 =={{header|Perl}}==
@@ Line 211: / Line 301: @@
 non-term [1488 2480 3472 4464 8432 9424 10416 21328 22320 55056 95728 96720 236592 459792 881392 882384]
 15355717786080   non-term [15355717786080 44534663601120]</pre>
 =={{header|Perl 6}}==
 <lang perl6>sub propdivsum (\x) {