单位分数难题续

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#一致了

1. 
   
2. N()=
   
3. {
   
4.   7
   
5. }
   
6. 
   
7. K()=
   
8. {
   
9.     N()-2
   
10. }
    
11. 
    
12. BEGIN()=
    
13. {
    
14.     0
    
15. }
    
16. 
    
17. END()=
    
18. {
    
19.     -1
    
20. }
    
21. 
    
22. OUTFILENAME()=
    
23. {
    
24.     "/home/zdu/N7d.txt"
    
25. }
    
26. 
    
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. }
    
53. 
    
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. }
    
75. 
    
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楼所言，不输出，只是统计，效果会怎样呢