求1000以內(nèi)所有的完數(shù)（例如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

復(fù)制代碼

我寫(xiě)的，運(yùn)行起來(lái)根本就沒(méi)有輸出

elu_ligao

回復(fù) 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個(gè)問(wèn)題，一個(gè)是沒(méi)有重新初始化sum的值，另一個(gè)是沒(méi)有把最后判斷sum和i的if語(yǔ)句放在子循環(huán)外。

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 ~]#

復(fù)制代碼

我把它改成awk，運(yùn)行速度快多了~

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 ~]#

復(fù)制代碼

關(guān)陰月飛

回復(fù) 1# hmchzb19

如果要效率的話，試試這個(gè)：

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
   

復(fù)制代碼

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

復(fù)制代碼

elu_ligao

學(xué)習(xí)了， 強(qiáng)大

reyleon

回復(fù) 4# 關(guān)陰月飛


    不明覺(jué)厲

hmchzb19

回復(fù) 4# 關(guān)陰月飛


   
這個(gè)太高效了

dn833

目測(cè)4樓又開(kāi)掛了

關(guān)陰月飛

@dn833@hmchzb19@reyleon

這里用到了一個(gè)求完全數(shù)的公式(百度上找的)：
設(shè)p為一個(gè)整數(shù)：
如果p是質(zhì)數(shù)，且2^p-1也是質(zhì)數(shù)，那么（2^p-1）X 2^(p-1)便是一個(gè)完全數(shù)。
例如p=2，是一個(gè)質(zhì)數(shù)，2^p-1=3也是質(zhì)數(shù)，（2^p-1）X2^（p-1）=3X2=6，是完全數(shù)。
例如p=3，是一個(gè)質(zhì)數(shù)，2^p-1=7也是質(zhì)數(shù)，（2^p-1）X2^（p-1）=7X4=28，是完全數(shù)。
例如p=5，是一個(gè)質(zhì)數(shù)，2^p-1=31也是質(zhì)數(shù)，（2^p-1）X2^（p-1）=31X16=496是完全數(shù)。

yestreenstars

回復(fù) 8# dn833
不愧是好基友，一眼就看穿了~