a prime question

gxqcn

原帖由 mathe 于 2009-3-9 21:12 发表
(%i6) primep(8*5+1);(%o6) true(%i7) primep(2*3*7*11*13*29*31*37+1);(%o7) true
(%i9) primep(4*7+1);(%o9) true(%i10) primep(4*3*5*11*13*29*31*37+1);(%o10) true
(%i29) primep(4*37+1);(%o29) true(%i30)  ...

不懂这种语法，
请以第一行为例给我们讲解一下，扫扫盲。。。

mathe

(%i6)是maxima的提示符
然后我输入primep(8*5+1)用来判断8*5+1是否素数。
而maxima下的换行符贴过来好像没了，所以没有显示换行。
接下去(%o6)true是maxima的输出，表示上面的计算结果是true,也就是8*5+1是素数。
当然中间还有i8等结果是false的判断，我没有贴出来

kofeffect

1. 
   
2. {7,11,59}
   
3. {3,31,59}
   
4. {2,61,59}
   
5. 
   
6. {7,13,71}
   
7. {5,19,71}
   
8. {3,37,71}
   
9. {2,73,71}
   
10. 
    
11. {7,41,239}
    
12. {5,61,239}
    
13. {2,241,239}
    
14. 
    
15. {7,53,311}
    
16. {5,79,311}
    
17. {3,157,311}
    
18. {2,313,311}
    
19. 
    
20. {11,37,359}
    
21. {7,61,359}
    
22. {3,181,359}
    
23. 
    
24. {7,71,419}
    
25. {3,211,419}
    
26. {2,421,419}
    
27. 
    
28. {13,37,431}
    
29. {7,73,431}
    
30. {5,109,431}
    
31. {2,433,431}
    
32. 
    
33. {17,31,479}
    
34. {13,41,479}
    
35. {3,241,479}
    
36. 
    
37. {19,29,503}
    
38. {13,43,503}
    
39. {5,127,503}
    
40. 
    
41. {11,61,599}
    
42. {7,101,599}
    
43. {5,151,599}
    
44. 
    
45. {19,37,647}
    
46. {7,109,647}
    
47. {5,163,647}
    
48. 
    
49. {23,31,659}
    
50. {11,67,659}
    
51. {3,331,659}
    
52. 
    
53. {19,41,719}
    
54. {13,61,719}
    
55. {11,73,719}
    
56. {5,181,719}
    
57. 
    
58. {13,71,839}
    
59. {5,211,839}
    
60. {3,421,839}
    
61. 
    
62. {17,73,1151}
    
63. {13,97,1151}
    
64. {7,193,1151}
    
65. 
    
66. {31,43,1259}
    
67. {19,71,1259}
    
68. {11,127,1259}
    
69. {7,211,1259}
    
70. 
    
71. {23,73,1583}
    
72. {19,89,1583}
    
73. {5,397,1583}
    
74. 
    
75. {31,67,1979}
    
76. {11,199,1979}
    
77. {7,331,1979}
    
78. 
    
79. {37,59,2087}
    
80. {7,349,2087}
    
81. {5,523,2087}
    
82. 
    
83. {41,61,2399}
    
84. {17,151,2399}
    
85. {11,241,2399}
    
86. {7,401,2399}
    
87. 
    
88. {37,73,2591}
    
89. {17,163,2591}
    
90. {7,433,2591}
    
91. 
    
92. {41,73,2879}
    
93. {31,97,2879}
    
94. {17,181,2879}
    
95. {13,241,2879}
    
96. 
    
97. {41,79,3119}
    
98. {11,313,3119}
    
99. {7,521,3119}
    
100. 
     
101. {17,211,3359}
     
102. {13,281,3359}
     
103. {11,337,3359}
     
104. 
     
105. {31,127,3779}
     
106. {19,211,3779}
     
107. {11,379,3779}
     
108. 
     
109. {41,103,4079}
     
110. {31,137,4079}
     
111. {11,409,4079}
     
112. 
     
113. {61,79,4679}
     
114. {37,131,4679}
     
115. {31,157,4679}
     
116. 
     
117. {47,109,4967}
     
118. {37,139,4967}
     
119. {19,277,4967}
     
120. 
     
121. {71,73,5039}
     
122. {41,127,5039}
     
123. {29,181,5039}
     
124. {19,281,5039}
     
125. {13,421,5039}
     
126. 
     
127. {61,89,5279}
     
128. {23,241,5279}
     
129. {17,331,5279}
     
130. 
     
131. {31,197,5879}
     
132. {29,211,5879}
     
133. {13,491,5879}
     
134. 
     
135. {31,199,5939}
     
136. {23,271,5939}
     
137. {19,331,5939}
     
138. 
     
139. {59,109,6263}
     
140. {19,349,6263}
     
141. {13,523,6263}
     
142. 
     
143. {71,97,6719}
     
144. {61,113,6719}
     
145. {29,241,6719}
     
146. {17,421,6719}
     
147. 
     
148. {73,97,6911}
     
149. {37,193,6911}
     
150. {17,433,6911}
     
151. 
     
152. {67,109,7127}
     
153. {37,199,7127}
     
154. {19,397,7127}
     
155. 
     
156. {71,109,7559}
     
157. {61,127,7559}
     
158. {43,181,7559}
     
159. {37,211,7559}
     
160. {29,271,7559}
     
161. {19,421,7559}
     
162. 
     
163. {71,127,8819}
     
164. {43,211,8819}
     
165. {19,491,8819}
     
166. 
     
167. {101,109,10799}
     
168. {73,151,10799}
     
169. {61,181,10799}
     
170. {41,271,10799}
     
171. 
     
172. {97,131,12479}
     
173. {53,241,12479}
     
174. {41,313,12479}
     
175. 
     
176. {67,211,13859}
     
177. {43,331,13859}
     
178. {31,463,13859}
     
179. 
     
180. {79,241,18719}
     
181. {61,313,18719}
     
182. {37,521,18719}
     
183. 
     
184. {71,277,19319}
     
185. {47,421,19319}
     
186. {43,461,19319}
     
187. 
     
188. {79,281,21839}
     
189. {71,313,21839}
     
190. {53,421,21839}
     
191. {43,521,21839}
     
192. 
     
193. {127,181,22679}
     
194. {109,211,22679}
     
195. {61,379,22679}
     
196. 
     
197. {113,241,26879}
     
198. {97,281,26879}
     
199. {61,449,26879}
     
200. 
     
201. {137,199,26927}
     
202. {89,307,26927}
     
203. {67,409,26927}
     
204. 
     
205. {137,211,28559}
     
206. {103,281,28559}
     
207. {71,409,28559}
     
208. 
     
209. {151,197,29399}
     
210. {71,421,29399}
     
211. {61,491,29399}
     
212. 
     
213. {163,211,34019}
     
214. {127,271,34019}
     
215. {71,487,34019}
     
216. 
     
217. {181,197,35279}
     
218. {127,281,35279}
     
219. {73,491,35279}
     
220. 
     
221. {127,353,44351}
     
222. {113,397,44351}
     
223. {97,463,44351}
     
224. 
     
225. {197,271,52919}
     
226. {127,421,52919}
     
227. {109,491,52919}
     
228. {211,313,65519}
     
229. {157,421,65519}
     
230. {127,521,65519}
     
231. 
     
232. {281,337,94079}
     
233. {211,449,94079}
     
234. {193,491,94079}
     
235. 
     
236. {331,337,110879}
     
237. {281,397,110879}
     
238. {241,463,110879}
     

复制代码

mathe

Kofeffect  是用穷举的方法？

无心人

2的23
3的59
4的71
5的5039
6的7559

medie2005

呵呵，那就再求７到１０的吧．

kofeffect

原帖由 mathe 于 2009-3-10 08:25 发表
Kofeffect  是用穷举的方法？

看出来了啊，呵呵。

无心人

似乎kofeffect求的不对

看我的结果

1. 
   
2. 
   
3.   import Primes
   
4.   import List
   
5.   
   
6.   main = do
   
7.         let a = sort [(r, p, q) | p <- primesTo10000, q <- primesTo10000, q > p, let r = p * q - (p + q), isPrime r]
   
8.         let b = map (\l -> (head l, length l))  \$ group \$ map (\(r, _, _) -> r) a
   
9.         let c = nub \$ sort \$ map (\(_, l) -> l) b
   
10.         let maxc = maximum c
    
11.         print maxc
    
12.         let d = [(n, l) | n <- [1..maxc], let hds = filter (\(_, l) -> l == n) b, not (null hds), let hd = fst \$ head hds,let l = filter (\(r, _, _) -> r == hd) a]
    
13.         print \$ show d
    

复制代码

=============================================================================
输出

"[(1,[(3,2,5)]),(2,[(11,2,13),(11,3,7)]),(3,[(59,2,61),(59,3,31),(59,7,11)]),(4,[(71,2,73),(71,3,37),(71,5,19),(71,7,13)]),(5,[
1151,2,1153),(1151,3,577),(1151,7,193),(1151,13,97),(1151,17,73)]),(6,[(2399,3,1201),(2399,5,601),(2399,7,401),(2399,11,241),(2
99,17,151),(2399,41,61)]),(7,[(9719,2,9721),(9719,3,4861),(9719,7,1621),(9719,13,811),(9719,19,541),(9719,37,271),(9719,61,163)
),(8,[(19319,3,9661),(19319,5,4831),(19319,7,3221),(19319,11,1933),(19319,29,691),(19319,43,461),(19319,47,421),(19319,71,277)]
,(9,[(7559,2,7561),(7559,11,757),(7559,13,631),(7559,19,421),(7559,29,271),(7559,37,211),(7559,43,181),(7559,61,127),(7559,71,1
9)]),(10,[(166319,19,9241),(166319,23,7561),(166319,37,4621),(166319,41,4159),(166319,67,2521),(166319,71,2377),(166319,73,2311
,(166319,127,1321),(166319,271,617),(166319,397,421)]),(11,[(347759,37,9661),(347759,47,7561),(347759,71,4969),(347759,73,4831)
(347759,109,3221),(347759,139,2521),(347759,181,1933),(347759,211,1657),(347759,271,1289),(347759,421,829),(347759,461,757)]),(
2,[(453599,61,7561),(453599,71,6481),(453599,73,6301),(453599,109,4201),(453599,113,4051),(453599,163,2801),(453599,181,2521),(
53599,211,2161),(453599,281,1621),(453599,379,1201),(453599,433,1051),(453599,601,757)])]"

无心人

2，3，4的对
5的最小是1151
6的最小是2399
7的最小是9719
这个确定是最小的了

无心人

8的是19319
9的最小是7559（已确定）
10的166319
11的347759
12的453599

8，10-12的无法确定最小