1 bc7cb1a1 2003-11-23 devnull #include <u.h>
2 bc7cb1a1 2003-11-23 devnull #include <libc.h>
3 bc7cb1a1 2003-11-23 devnull #include <bio.h>
5 bc7cb1a1 2003-11-23 devnull #define whsiz (sizeof(wheel)/sizeof(wheel[0]))
7 bc7cb1a1 2003-11-23 devnull double wheel[] =
9 bc7cb1a1 2003-11-23 devnull 2,10, 2, 4, 2, 4, 6, 2, 6, 4,
10 bc7cb1a1 2003-11-23 devnull 2, 4, 6, 6, 2, 6, 4, 2, 6, 4,
11 bc7cb1a1 2003-11-23 devnull 6, 8, 4, 2, 4, 2, 4, 8, 6, 4,
12 bc7cb1a1 2003-11-23 devnull 6, 2, 4, 6, 2, 6, 6, 4, 2, 4,
13 bc7cb1a1 2003-11-23 devnull 6, 2, 6, 4, 2, 4, 2,10,
16 bc7cb1a1 2003-11-23 devnull Biobuf bin;
18 bc7cb1a1 2003-11-23 devnull void factor(double);
21 bc7cb1a1 2003-11-23 devnull main(int argc, char *argv[])
23 bc7cb1a1 2003-11-23 devnull double n;
27 bc7cb1a1 2003-11-23 devnull if(argc > 1) {
28 bc7cb1a1 2003-11-23 devnull for(i=1; i<argc; i++) {
29 bc7cb1a1 2003-11-23 devnull n = atof(argv[i]);
30 bc7cb1a1 2003-11-23 devnull factor(n);
32 bc7cb1a1 2003-11-23 devnull exits(0);
35 bc7cb1a1 2003-11-23 devnull Binit(&bin, 0, OREAD);
36 bc7cb1a1 2003-11-23 devnull for(;;) {
37 bc7cb1a1 2003-11-23 devnull l = Brdline(&bin, '\n');
38 bc7cb1a1 2003-11-23 devnull if(l == 0)
40 bc7cb1a1 2003-11-23 devnull n = atof(l);
41 bc7cb1a1 2003-11-23 devnull if(n <= 0)
43 bc7cb1a1 2003-11-23 devnull factor(n);
45 bc7cb1a1 2003-11-23 devnull exits(0);
49 bc7cb1a1 2003-11-23 devnull factor(double n)
51 bc7cb1a1 2003-11-23 devnull double quot, d, s;
54 bc7cb1a1 2003-11-23 devnull print("%.0f\n", n);
55 bc7cb1a1 2003-11-23 devnull if(n == 0)
57 bc7cb1a1 2003-11-23 devnull s = sqrt(n) + 1;
58 bc7cb1a1 2003-11-23 devnull while(modf(n/2, ") == 0) {
59 bc7cb1a1 2003-11-23 devnull print(" 2\n");
60 bc7cb1a1 2003-11-23 devnull n = quot;
61 bc7cb1a1 2003-11-23 devnull s = sqrt(n) + 1;
63 bc7cb1a1 2003-11-23 devnull while(modf(n/3, ") == 0) {
64 bc7cb1a1 2003-11-23 devnull print(" 3\n");
65 bc7cb1a1 2003-11-23 devnull n = quot;
66 bc7cb1a1 2003-11-23 devnull s = sqrt(n) + 1;
68 bc7cb1a1 2003-11-23 devnull while(modf(n/5, ") == 0) {
69 bc7cb1a1 2003-11-23 devnull print(" 5\n");
70 bc7cb1a1 2003-11-23 devnull n = quot;
71 bc7cb1a1 2003-11-23 devnull s = sqrt(n) + 1;
73 bc7cb1a1 2003-11-23 devnull while(modf(n/7, ") == 0) {
74 bc7cb1a1 2003-11-23 devnull print(" 7\n");
75 bc7cb1a1 2003-11-23 devnull n = quot;
76 bc7cb1a1 2003-11-23 devnull s = sqrt(n) + 1;
79 bc7cb1a1 2003-11-23 devnull for(i=1;;) {
80 bc7cb1a1 2003-11-23 devnull d += wheel[i];
81 bc7cb1a1 2003-11-23 devnull while(modf(n/d, ") == 0) {
82 bc7cb1a1 2003-11-23 devnull print(" %.0f\n", d);
83 bc7cb1a1 2003-11-23 devnull n = quot;
84 bc7cb1a1 2003-11-23 devnull s = sqrt(n) + 1;
87 bc7cb1a1 2003-11-23 devnull if(i >= whsiz) {
89 bc7cb1a1 2003-11-23 devnull if(d > s)
93 bc7cb1a1 2003-11-23 devnull if(n > 1)
94 bc7cb1a1 2003-11-23 devnull print(" %.0f\n",n);
95 bc7cb1a1 2003-11-23 devnull print("\n");