题目37：找出所有11个可以双向裁剪的质数的和

渡风

此代码使用matlab编程
Problem37所用时间为： 14.9352秒
Problem37的答案为： 748317

1. %% Problem37.m
   
2. % 最后编辑时间：17-06-14 22:34
   
3. % 找出11个双向裁剪质数的和
   
4. % Problem37所用时间为： 14.9352秒
   
5. % Problem37的答案为： 748317
   
6. function Output = Problem37()
   
7. tic
   
8. Limit = 1000000;
   
9. Vector = GetPrime(Limit);
   
10. Vector = Vector(5:end);
    
11. N = length(Vector);
    
12. 
    
13. Output = 0;
    
14. 
    
15. Time = 0;
    
16. 
    
17. for ii = 1:N
    
18.     Judge = 1;
    
19.     L = length(num2str(Vector(ii)));
    
20.     for jj = 1:L-1
    
21.         if isprime(floor(Vector(ii)/(10^jj))) == 0
    
22.             Judge = 0;
    
23.             break
    
24.         end
    
25.     end
    
26.    
    
27.     if Judge == 1
    
28.         Str = num2str(Vector(ii));
    
29.         for jj = 2:L
    
30.             if isprime(str2double(Str(jj:L))) == 0
    
31.                Judge = 0;
    
32.                break
    
33.             end
    
34.         end
    
35.     end
    
36.    
    
37.     if Judge == 1
    
38.         Output = Output + Vector(ii);
    
39.         Time = Time + 1;
    
40.         disp(Vector(ii))
    
41.         if Time == 11
    
42.             break
    
43.         end        
    
44.     end
    
45.          
    
46. end
    
47. 
    
48. toc
    
49. disp('此代码使用matlab编程')
    
50. disp(['Problem37所用时间为： ',num2str(toc),'秒'])
    
51. disp(['Problem37的答案为： ',num2str(Output)])
    
52. end
    
53. %% 得到n以下的所有质数
    
54. function Vector = GetPrime(n)
    
55. if nargin == 0
    
56. n = 10000;
    
57. end
    
58. 
    
59. Judge = ones(1,n-1);
    
60. Judge(1) = 0;
    
61. for ii = 1:n-1
    
62.     if Judge(ii) == 1
    
63.         for jj = 2*ii : ii : n-1
    
64.             Judge(jj) = 0;
    
65.         end
    
66.     end
    
67. end
    
68. Vector = find(Judge == 1);
    
69. end

复制代码