%I#8 2016年10月11日08:06:37
%S 4,0,5,4,3,2,1,0,8,8,8,8,8,1,8,4,0,8,8、8,2,1,8,6,5,4、3,2,0,8、8,
%T 8,8,8,1,8,8-8,5,8,2,1,8,4,0,8,8,8,3,2,8,88,8、8,2,8,1,8、7,6,5,4,
%U 3,2,1,0,8,8,8,8,8
%N对于无限数量的整数对(x,y),8x^2-ny^2取的最小非负值。
%t p=8;f[n_,z_]:=查找实例[x>0&y>0&p*x^2-n*y^2==z,{x,y},整数,1];a[n]:=对于[z=0,z<=p,z=z+GCD[n,p],其中[n==49,z=4;fz=$中止,n==81,z=8;fz=$中止,n==83,z=5;fz=201,z=3;fz:$中止,True,fz=TimeConstrained[f[n,z],300]];fz=如果[fz===$中止,{{x->Null,y->Null}},fz];如果[fz=!={},打印[“a(”,n,“)=”,z,“{x,y}=”,{x,y}/.fz[[1]];返回[z]]];表[a[n],{n,1105}](*_Jean-François Alcover_,2016年10月11日*)
%K nonn公司
%O 1,1号机组
%A·热心的W·威尔逊_
|