利用awk实现整数的因数分解---找问题？

UltimateSniper

#计算整数的因数分解，一行列出一个
#语法：awk -f factorize.awk Numbers
#使用举例：awk -f factorize.awk 10000 100001

#factorize.awk内容如下
{
        n = int($1)
        m = n = (n >= 2) ? n : 2
        factors = ""
        for ( k=2; (m>1)&&(k^2 <= m); )
        {
                if ( int(m%n) != 0 )
                {
                        k++
                        continue
                }
                m /= k
                factors = (factors=="") ? ("" k) : (factors " * " k)
        }
        if ((1<m) && (m<n))
                factors = factors " * " m
        print n, (factors == "") ? "is prime" : ("= " factors)
}

#但是运行总会出问题，提示为：10000无法打开。为什么将10000当做文件来处理了呢？希望帮忙解释一下。

sk1418

这样呢？

1. awk -f factorize.awk <<<"10000 100001"

复制代码

或者

1. echo "10000 100001"|awk -f factorize.awk

复制代码

seesea2517

如果你那个是按行处理的话，楼上的办法就行了。如果你那个是想用参数传递，可以试试两个方案（运行结果好像不对，那是代码逻辑的事啦）：

1. ##########################################
   
2. ## 用传参数的方式
   
3. 
   
4. [seesea@UC ~]$ awk -f factorize.awk 10000
   
5. 10000 = 2 * 5000
   
6. [seesea@UC ~]$ awk -f factorize.awk
   
7. No argument.
   
8. [seesea@UC ~]$ cat factorize.awk
   
9. BEGIN {
   
10.     if (ARGC < 2)
    
11.     {
    
12.         print "No argument."
    
13.         exit
    
14.     }
    
15. 
    
16.     fac(ARGV[1])
    
17. }
    
18. 
    
19. function fac(num) {
    
20.         n = int(num)
    
21.         m = n = (n >= 2) ? n : 2
    
22.         factors = ""
    
23.         for ( k=2; (m>1)&&(k^2 <= m); )
    
24.         {
    
25.                 if ( int(m%n) != 0 )
    
26.                 {
    
27.                         k++
    
28.                         continue
    
29.                 }
    
30.                 m /= k
    
31.                 factors = (factors=="") ? ("" k) : (factors " * " k)
    
32.         }
    
33.         if ((1<m) && (m<n))
    
34.                 factors = factors " * " m
    
35.         print n, (factors == "") ? "is prime" : ("= " factors)
    
36. }
    
37. 
    
38. 
    
39. #########################################
    
40. ## 或者用传变量的方式
    
41. 
    
42. [seesea@UC ~]$ awk -v num=1000 -f factorize.awk
    
43. 1000 = 2 * 500
    
44. [seesea@UC ~]$ awk  -f factorize.awk            
    
45. No argument.
    
46. [seesea@UC ~]$ cat factorize.awk
    
47. BEGIN {
    
48.     if (num=="")
    
49.     {
    
50.         print "No argument."
    
51.         exit
    
52.     }
    
53.     fac(num)
    
54. }
    
55. 
    
56. function fac(num) {
    
57.         n = int(num)
    
58.         m = n = (n >= 2) ? n : 2
    
59.         factors = ""
    
60.         for ( k=2; (m>1)&&(k^2 <= m); )
    
61.         {
    
62.                 if ( int(m%n) != 0 )
    
63.                 {
    
64.                         k++
    
65.                         continue
    
66.                 }
    
67.                 m /= k
    
68.                 factors = (factors=="") ? ("" k) : (factors " * " k)
    
69.         }
    
70.         if ((1<m) && (m<n))
    
71.                 factors = factors " * " m
    
72.         print n, (factors == "") ? "is prime" : ("= " factors)
    
73. }

复制代码

UltimateSniper

回复 2# sk1418


    想问一下，为什么是3个<呢？ 能不能解释一下，这点不太清楚。

   

awk -f factorize.awk <<<"10000 100001"

UltimateSniper

回复 3# seesea2517


    逻辑错误找到了：

1. for ( k=2; (m>1)&&(k^2 <= m); )
   
2.         {
   
3.                 if ( int(m%k) != 0 )
   
4.                 {
   
5.                         k++
   
6.                         continue
   
7.                 }
   
8.                 m /= k
   
9.                 factors = (factors=="") ? ("" k) : (factors " * " k)
   
10.         }

复制代码

应该是： m%k 而不是 m%n 我的失误。

PS：程序现在还只能一次读取一个数，我得再去修改修改，增强一下功能。不过还是谢谢了！

seesea2517

回复 4# UltimateSniper


    搜索 here string 看看有没有收获。

UltimateSniper

回复 6# seesea2517


    十分感谢，让我有机会再次看到 “Here Document” ，重新复习了一下。

    here string 应该算是 here document 的特例了吧？
   

bash，ksh或zsh中也可以用here-字串：
$ tr a-z A-Z <<<"Yes it is a string"
YES IT IS A STRING

seesea2517

回复 7# UltimateSniper


    是不是特例这个我还真没想过，只是觉得看上去长得像，用处也像

UltimateSniper

回复 8# seesea2517


    最后的那个例子我是在 Here文档的中文wiki中看到的：

http://zh.wikipedia.org/wiki/Here%E6%96%87%E6%A1%A3

seesea2517

回复 9# UltimateSniper


    学习了，感谢分享。这样看来，应该是特例。