单位分数难题续

mathe

N=7的结果为227907个,结果文件共7.2M 看来N=8是不能保留结果文件的了

数学星空

medie2005

http://www.research.att.com/~njas/sequences/A002967

数学星空

呵呵,看来只有计算n=9才能打破记录 可能只有mathe....

无心人

mathe有计算这个的先例 不穷举也可以计算的

无心人

另外，肚子去美国了？ 怎么IP显示在美国？

mathe

IP不能说明任何问题,现在呢? 发现我还有遗漏,而前面计数过多,是我不知道gp写文件,不会覆盖原先的文件.所以如果原先文件中还有数据,会将多次数据累加

mathe

这下N=5和N=6的数目同12#一致了

2. N()=
3. {
4. 7
5. }
7. K()=
8. {
9. N()-2
10. }
12. BEGIN()=
13. {
14. 0
15. }
17. END()=
18. {
19. -1
20. }
22. OUTFILENAME()=
23. {
24. "/home/zdu/N7d.txt"
25. }
27. output()=
28. {
29. local(V,MM,RR,uu,f,v,v2,a,b,d);
30. MM=M[K()];
31. RR=R[K()];
32. V=MM;
33. c=c+1;
34. if(c>=BEGIN()&&(END()<0||c<=END()),
35. f=divisors(V);
36. for(u=1,length(f),
37. for(u2=u,length(f),
38. v=f[u];v2=f[u2];
39. if((v+v2)%RR!=0, next());
40. if(gcd(v,v2)>1,next());
41. d=(v+v2)/RR*(V/(v*v2));
42. a=d*v;b=d*v2;
43. if(a>=n[K()],
44. n[K()+1]=a;
45. n[K()+2]=b;
46. write(OUTFILENAME(),n);
47. )
48. )
49. )
50. );
51. if(c%100==0, print("process " c " lines"))
52. }
54. search(level)=
55. {
56. local(ll,uu,dd);
57. if(level>=K(),
58. output(),
59. ll=floor(M[level]/R[level])+1;
60. uu=floor((N()-level)*M[level]/R[level]);
61. if(ll<=n[level],ll=n[level]);
62. for(i=ll,uu,
63. n[level+1]=i;
64. M[level+1]=M[level]*i;
65. R[level+1]=R[level]*i-M[level];
66. dd=gcd(M[level+1],R[level+1]);
67. M[level+1]=M[level+1]/dd;
68. R[level+1]=R[level+1]/dd;
69. A[i%4+1]=A[i%4+1]+1;
70. search(level+1);
71. A[i%4+1]=A[i%4+1]-1
72. )
73. )
74. }
76. search0()=
77. {
78. n=vector(N());
79. M=vector(N());
80. R=vector(N());
81. A=vector(4);
82. c=0;
83. for(i=2,8,
84. n[1]=i;
85. R[1]=i-1;
86. M[1]=i;
87. A[i%4+1]=A[i%4+1]+1;
88. search(1);
89. A[i%4+1]=A[i%4+1]-1;
90. );
91. }

复制代码

mathe

呵呵,看来只有计算n=9才能打破记录 可能只有mathe.... 数学星空 发表于 2010-1-19 09:52

用现在穷举的算法是无法破记录的了.倒是计算n=8问题不大.设置即使将所有数据保存下来都有可能,只是文件太大了,估计要超1G

wayne

用现在穷举的算法是无法破记录的了.倒是计算n=8问题不大.设置即使将所有数据保存下来都有可能,只是文件太大了,估计要超1G mathe 发表于 2010-1-19 10:29

要是如15楼所言，不输出，只是统计，效果会怎样呢