a最大值=150;(*警告:O(n^3)*)
做[
如果[x\[NotElement]基本原理,请继续[]];
{d1,d2}=1/(-(1/a)+1/b+1/c+{1,-1}2x);
如果[GCD[a,b,c,d1]==1,{a,b、c,d1}//Sow];
如果[d2>c||d2==d1,则继续[]];
如果[GCD[a,b,c,d2]==1,{a,b、c,d2}//Sow];
,{a,aMax},{b,2,a-2},}c,Min[b,a-b]}
,{x,{Sqrt[(-a+b+c)/(-a b c)]}}]//Reap//Last//Last//TableForm
|