%I#8 2016年10月11日16:03:14
%S 1,2,3,4,0,1,2,3,1,5,6,3,5,10,11,3,13,14,0,1,2,2,4,5,7,5,10,11,
%电话:3,7,5,15,4,5,2,19,20,9,22,23,24,0,1,2,3,1,5,6,7,10,4,3,5,6,15,16,
%U 3,7,11,5,19,22,3,1,25,11,27,5,5,30,31,7,33,34,0,12,3,4,5,1,7,8,5,10,11
%N无穷多个整数对(x,y)的N*x^2-5*y^2所取的最小非负值。
%t p=5;f[n_,z_]:=查找实例[x>0&y>0&&n*x^2-p*y^2==z,{x,y},整数,1];a[n_]:=对于[z=0,True,z=z+GCD[n,p],fz=时间约束[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,1,91}](*Jean-François Alcover_,2016年10月11日*)
%K nonn公司
%O 1,2号机组
%A·热心的W·威尔逊_
|