求1000以内所有的完数（例如6=1+2+3，1*2*3=6）

hmchzb19

1. #! /bin/bash
   
2. 
   
3. 
   
4. 
   
5. for ((  i=1;i<=1000;i++ ))
   
6. do
   
7.     for (( j=i;j>1;j-- ))
   
8.     do
   
9.         let temp=$i%$j
   
10.         let temp1=$i/$j
    
11.         if [ $temp -eq 0 ];then
    
12.             k=$temp1
    
13.             let sum=$sum+$k
    
14.             if [ $sum -eq $i ];then
    
15.                 echo $i
    
16.             fi
    
17.         fi
    
18.     done
    
19. done

复制代码

我写的，运行起来根本就没有输出

elu_ligao

回复 1# hmchzb19


  
#! /bin/sh -


for (( i=1; i<=1000; ++i ))
do
    s=0
    for (( j=1; 2*j<=i; ++j ))
    do
        let tmp=$i%$j
        if [[ $tmp -eq 0 ]]; then
            let s=$s+$j
        fi
    done
    if [[ $s -eq $i ]]; then
        echo $i
    fi
done

yestreenstars

你的代码有2个问题，一个是没有重新初始化sum的值，另一个是没有把最后判断sum和i的if语句放在子循环外。

1. [root@localhost ~]# cat t
   
2. #! /bin/bash
   
3. 
   
4. for((i=1;i<=1000;i++))
   
5. do
   
6.         sum=0
   
7.         for((j=i;j>1;j--))
   
8.         do
   
9.                 let temp=i%j
   
10.                 let temp1=i/j
    
11.                 if [ $temp -eq 0 ];then
    
12.                         let sum+=temp1
    
13.                 fi
    
14.         done
    
15.         if [ $sum -eq $i ];then
    
16.                 echo $i
    
17.         fi
    
18. done
    
19. [root@localhost ~]# ./t
    
20. 6
    
21. 28
    
22. 496
    
23. [root@localhost ~]#

复制代码

我把它改成awk，运行速度快多了~

1. [root@localhost ~]# awk 'BEGIN{for(i=1;i<=1000;i++){s=0;for(j=i;j>1;j--){t1=i%j;t2=i/j;if(t1==0)s+=t2}if(s==i)print i}}'
   
2. 6
   
3. 28
   
4. 496
   
5. [root@localhost ~]#

复制代码

关阴月飞

回复 1# hmchzb19

如果要效率的话，试试这个：

1. 
   
2. #!/bin/bash
   
3. for (( i=1; $[(2**i-1)*(2**(i-1))]<=$1;i++ ))
   
4. do
   
5. if [ `factor $i $[2**i-1] |awk 'NF==2' |wc -l` -eq 2 ];then
   
6.         echo $[(2**i-1)*(2**(i-1))]
   
7. fi
   
8. done
   

复制代码

1. [root@master awktest]# cat a.sh
   
2. #!/bin/bash
   
3. for (( i=1; $[(2**i-1)*(2**(i-1))]<=$1;i++ ))
   
4. do
   
5. if [ `factor $i $[2**i-1] |awk 'NF==2' |wc -l` -eq 2 ];then
   
6.         echo $[(2**i-1)*(2**(i-1))]
   
7. fi
   
8. done
   
9. [root@master awktest]# time ./a.sh 1000
   
10. 6
    
11. 28
    
12. 496
    
13. 
    
14. real    0m0.034s
    
15. user    0m0.002s
    
16. sys     0m0.023s
    
17. [root@master awktest]# time ./a.sh 10000
    
18. 6
    
19. 28
    
20. 496
    
21. 8128
    
22. 
    
23. real    0m0.048s
    
24. user    0m0.007s
    
25. sys     0m0.027s
    
26. [root@master awktest]# time ./a.sh 10000000000000000
    
27. 6
    
28. 28
    
29. 496
    
30. 8128
    
31. 33550336
    
32. 8589869056
    
33. 137438691328
    
34. 
    
35. real    0m0.160s
    
36. user    0m0.059s
    
37. sys     0m0.031s

复制代码

elu_ligao

学习了， 强大

reyleon

回复 4# 关阴月飞


    不明觉厉

hmchzb19

回复 4# 关阴月飞


   
这个太高效了

dn833

目测4楼又开挂了

关阴月飞

@dn833@hmchzb19@reyleon

这里用到了一个求完全数的公式(百度上找的)：
设p为一个整数：
如果p是质数，且2^p-1也是质数，那么（2^p-1）X 2^(p-1)便是一个完全数。
例如p=2，是一个质数，2^p-1=3也是质数，（2^p-1）X2^（p-1）=3X2=6，是完全数。
例如p=3，是一个质数，2^p-1=7也是质数，（2^p-1）X2^（p-1）=7X4=28，是完全数。
例如p=5，是一个质数，2^p-1=31也是质数，（2^p-1）X2^（p-1）=31X16=496是完全数。

yestreenstars

回复 8# dn833
不愧是好基友，一眼就看穿了~