Add factor.

I was reading http://www.muppetlabs.com/~breadbox/txt/rsa.html and it mentioned
"factor" and I noticed it was in coreutils. I'm not sure why it's in coreutils,
but it's pretty trivial, so...

2 files changed, 95 insertions, 0 deletions

diff --git a/scripts/test/factor.test b/scripts/test/factor.test
new file mode 100755
index 00000000..a3e4cbf8
--- /dev/null
+++ b/scripts/test/factor.test
@@ -0,0 +1,18 @@
+#!/bin/bash
+
+[ -f testing.sh ] && . testing.sh
+
+#testing "name" "command" "result" "infile" "stdin"
+
+testing "factor -32" "factor -32" "-32: -1 2 2 2 2 2\n" "" ""
+testing "factor 0" "factor 0" "0: 0\n" "" ""
+testing "factor 1" "factor 1" "1: 1\n" "" ""
+testing "factor 2" "factor 2" "2: 2\n" "" ""
+testing "factor 3" "factor 3" "3: 3\n" "" ""
+testing "factor 4" "factor 4" "4: 2 2\n" "" ""
+testing "factor 10000000017" "factor 10000000017" \
+        "10000000017: 3 3 3 7 7 7 1079797\n" "" ""
+testing "factor 10000000018" "factor 10000000018" \
+        "10000000018: 2 131 521 73259\n" "" ""
+testing "factor 10000000019" "factor 10000000019" \
+        "10000000019: 10000000019\n" "" ""
diff --git a/toys/other/factor.c b/toys/other/factor.c
new file mode 100644
index 00000000..e3992a8b
--- /dev/null
+++ b/toys/other/factor.c
@@ -0,0 +1,77 @@
+/* factor.c - Factor integers
+ *
+ * Copyright 2014 Rob Landley <rob@landley.net>
+ *
+ * No standard, but it's in coreutils
+
+USE_FACTOR(NEWTOY(factor, 0, TOYFLAG_USR|TOYFLAG_BIN))
+
+config FACTOR
+  bool "factor"
+  default y
+  help
+    usage: factor NUMBER...
+
+    Factor integers.
+*/
+
+#include "toys.h"
+
+static void factor(char *s)
+{
+  long l, ll;
+
+  l = strtol(s, &s, 0);
+  if (*s) {
+    error_msg("%s: not integer");
+    return;
+  }
+
+  printf("%ld:", l);
+
+  // Negative numbers have -1 as a factor
+  if (l < 0) {
+    printf(" -1");
+    l *= -1;
+  }
+
+  // Deal with 0 and 1 (and 2 since we're here)
+  if (l < 3) {
+    printf(" %ld\n", l);
+    return;
+  }
+
+  // Special case factors of 2
+  while (l && !(l&1)) {
+    printf(" 2");
+    l >>= 1;
+  }
+
+  // test odd numbers.
+  for (ll=3; ;ll += 2) {
+    if (ll*ll>l) {
+      if (l>1) printf(" %ld", l);
+      break;
+    }
+    while (!(l%ll)) {
+      printf(" %ld", ll);
+      l /= ll;
+    }
+  }
+  xputc('\n');
+}
+
+void factor_main(void)
+{
+  if (toys.optc) {
+    char **ss;
+
+    for (ss = toys.optargs; *ss; ss++) factor(*ss);
+  } else for (;;) {
+    char *s = 0;
+    size_t len = 0;
+
+    if (-1 == getline(&s, &len, stdin)) break;
+    factor(s);
+  }
+}