an={};
he[a_,b_,e_]:=1/(2e)平方[(-((a-b-e)(a+b-e)
paX[e_]:={e,0}(*顶点A坐标*)
pbX[a_,b_,e_]:={(-a^2+b^2+e^2)/(2e),he[a,b,e]}(*顶点b坐标*)
pc={0,0};(*顶点C坐标*)
pdX[c_,d_,e_]:={(c^2-d^2+e^2)/(2e),-he[c,d,e]}(*顶点d坐标*)
concaveQ[{bx_,by_},{dx_,dy_}
gQ[x_,y_]:=模块[{z=x-y,res=False},Do[If[z[[i]]>0,res=True;中断[],
如果[z[[i]]<0,中断[]]],{i,1,4}];资源]
规范Q[{a_,b_,c_,d_}]:=模块[{m={a,b,c,d}},If[(gQ[{b,a,d,c},m]| | gQ[[{d,c,b,a},m]| | gQ[{c,d,a,b},m]),False,True]]
做[cnt=0;
Do[pa=paX[e];pb=pbX[a,b,e];pd=pdX[c,d,e];
如果[(f=Sqrt[(pb-pd).(pb-pd)];整数Q[f])&&concaveQ[pb,pd,e]&&canonicalQ[{a,b,c,d}],cnt++
(*;打印[{{a,b,c,d,e,f},图形[Line[{pa,pb,pc,pd,pa}]}]*)],
{b,1,a},{e,a-b+1,a-1},{c,1,a},}d,Abs[e-c]+1,Min[a,e+c-1]}];
附加到[an,cnt],
{a,1,75}
]
一个
|