小练习：20160620 算出所有不能写成两个过剩数之和的正整数之和

小剑剑

1. import time
   
2. begin=time.clock()
   
3. a=list(filter(lambda x: sum([i  for  i in range(1,x) if not x%i])>x , range(28111)))
   
4. b=list(set([i+ii for i in a for ii in a if i+ii<28123]))
   
5. c=[i for i in range(1,28123) if i not in b]
   
6. d=sum(c)
   
7. end=time.clock()
   
8. print(end-begin,d)
   

复制代码

凯布里克恩

1. import math
   
2. import time
   
3. start = time.perf_counter()
   
4. def truefactor(n):
   
5.             fac_list = [1]
   
6.             sum1 = 0
   
7.             halfnum = n // 2 + 1
   
8.             for i in range(2, halfnum):
   
9.                         if n % i == 0:
   
10.                                     fac_list.append(i)
    
11.             sum1 = sum(fac_list)
    
12.             if sum1 > n:
    
13.                         return True
    
14.             else:
    
15.                         return False
    
16. 
    
17. def getfactor():
    
18.             factorlist = []
    
19.             factorlist1 = []
    
20.             sum_allintegers = 0
    
21.             sum_allfactors = 0
    
22. 
    
23.             #获得所有的过剩数
    
24.             for n in range(1, 28124):
    
25.                         sum_allintegers += n
    
26.                         if truefactor(n):
    
27.                                     factorlist.append(n)
    
28. 
    
29.             #获得过剩数相加的和，然后得到过剩数之和，总和减去过去过剩数之后为所有不是过剩数的正整数之和
    
30.             length = len(factorlist)
    
31.             for i in range(length):
    
32.                         for j in range(length):
    
33.                                     a = factorlist[i] + factorlist[j]
    
34.                                     if  a > 28123:
    
35.                                                 break
    
36.                                     else:
    
37.                                                 factorlist1.append(a)                        
    
38.             factorlist1 = sorted(list(set(factorlist1)))
    
39.             sum_allfactors = sum(factorlist1)
    
40.             return sum_allintegers - sum_allfactors
    
41.                         
    
42. print(getfactor())
    
43. print('Running time is %s Seconds' %(time.perf_counter() - start))
    

复制代码

bacon6581

1. from time import time
   
2. start=time()
   
3. 
   
4. lst=[]
   
5. for i in range(12,28124):
   
6.     l=[1]
   
7.     j=i**0.5
   
8.     if j==int(j):
   
9.         l.append(j)
   
10.     else:
    
11.         j=int(j)+1
    
12. 
    
13.     for x in range(2,int(j)):
    
14.         if i%x==0:
    
15.             l.append(x)
    
16.             l.append(i/x)
    
17. 
    
18.     sm=0
    
19.     for each in l:
    
20.         sm+=each
    
21. 
    
22.     if sm>i:
    
23.         lst.append(i)
    
24. 
    
25. #print(lst)
    
26. 
    
27. lst1=list(range(1,28124))
    
28. for m in lst:
    
29.     for n in range(lst.index(m),len(lst)):
    
30.         if m+lst[n]>28124:
    
31.             break
    
32.         if m + lst[n] in lst1:
    
33.             lst1.remove(m + lst[n])
    
34. 
    
35. result=0
    
36. for each in lst1:
    
37.     result+=each
    
38. print(result)
    
39. print(time()-start)
    

复制代码

小剑剑

挥舞乾坤 发表于 2016-6-25 18:07

为什么是20612不是28123啊

挥舞乾坤

小剑剑 发表于 2016-6-29 15:47
为什么是20612不是28123啊

百度一下盈数，有介绍

yellowfgglk

我仿佛看到一尊大神。。。是不是数学不好的基本和编程无缘的