From 2a53f53d74167353cf434af59392c72a885cb7f6 Mon Sep 17 00:00:00 2001
From: Rob Landley <rob@landley.net>
Date: Fri, 1 Aug 2014 09:08:00 -0500
Subject: 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...
---
 toys/other/factor.c | 77 +++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 77 insertions(+)
 create mode 100644 toys/other/factor.c

(limited to 'toys/other/factor.c')

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);
+  }
+}
-- 
cgit v1.2.3