From e5888a1ffdae813d7575f5fb02275c6bb07e5199 Mon Sep 17 00:00:00 2001
From: Taru Karttunen <taruti@taruti.net>
Date: Wed, 30 Mar 2011 15:46:40 +0300
Subject: Import sources from 2011-03-30 iso image

---
 sys/src/cmd/factor.c | 96 ++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 96 insertions(+)
 create mode 100755 sys/src/cmd/factor.c

(limited to 'sys/src/cmd/factor.c')

diff --git a/sys/src/cmd/factor.c b/sys/src/cmd/factor.c
new file mode 100755
index 000000000..0232d4e54
--- /dev/null
+++ b/sys/src/cmd/factor.c
@@ -0,0 +1,96 @@
+#include <u.h>
+#include <libc.h>
+#include <bio.h>
+
+#define	whsiz	(sizeof(wheel)/sizeof(wheel[0]))
+
+double	wheel[] =
+{
+	 2,10, 2, 4, 2, 4, 6, 2, 6, 4,
+	 2, 4, 6, 6, 2, 6, 4, 2, 6, 4,
+	 6, 8, 4, 2, 4, 2, 4, 8, 6, 4,
+	 6, 2, 4, 6, 2, 6, 6, 4, 2, 4,
+	 6, 2, 6, 4, 2, 4, 2,10,
+};
+
+Biobuf	bin;
+
+void	factor(double);
+
+void
+main(int argc, char *argv[])
+{
+	double n;
+	int i;
+	char *l;
+
+	if(argc > 1) {
+		for(i=1; i<argc; i++) {
+			n = atof(argv[i]);
+			factor(n);
+		}
+		exits(0);
+	}
+
+	Binit(&bin, 0, OREAD);
+	for(;;) {
+		l = Brdline(&bin, '\n');
+		if(l ==  0)
+			break;
+		n = atof(l);
+		if(n <= 0)
+			break;
+		factor(n);
+	}
+	exits(0);
+}
+
+void
+factor(double n)
+{
+	double quot, d, s;
+	int i;
+
+	print("%.0f\n", n);
+	if(n == 0)
+		return;
+	s = sqrt(n) + 1;
+	while(modf(n/2, &quot) == 0) {
+		print("     2\n");
+		n = quot;
+		s = sqrt(n) + 1;
+	}
+	while(modf(n/3, &quot) == 0) {
+		print("     3\n");
+		n = quot;
+		s = sqrt(n) + 1;
+	}
+	while(modf(n/5, &quot) == 0) {
+		print("     5\n");
+		n = quot;
+		s = sqrt(n) + 1;
+	}
+	while(modf(n/7, &quot) == 0) {
+		print("     7\n");
+		n = quot;
+		s = sqrt(n) + 1;
+	}
+	d = 1;
+	for(i=1;;) {
+		d += wheel[i];
+		while(modf(n/d, &quot) == 0) {
+			print("     %.0f\n", d);
+			n = quot;
+			s = sqrt(n) + 1;
+		}
+		i++;
+		if(i >= whsiz) {
+			i = 0;
+			if(d > s)
+				break;
+		}
+	}
+	if(n > 1)
+		print("     %.0f\n",n);
+	print("\n");
+}
-- 
cgit v1.2.3