%I#18 2023年3月22日12:02:02
%S 8,12,5,24,28,40,44,13,56,60,17,76,21,88,92104,29120124,33136140,
%电话:37152156,41168172184188204,53220,57232236,61248,65264,
%电话:268,69280284,73296,77312316328332,85344348,89364,93376380,97
%实二次数域的判别式。
%D David A.Cox,“形式x^2+ny^2的素数”,威利,1989年,第103页。
%F对于无平方n>=2,如果n=1,则列出n,否则列出4n。
%p with(numtheory):a:=proc(n)如果issqrfree(n),则返回(逐段(n mod 4=1,n,4*n)),否则返回(NULL)fi:结束:seq(a(n)、n=2..150);#C.罗纳尔多(aga_new_ac(AT)hotmail.com),2005年1月3日
%t收获[For[n=2,n<=100,n++,If[SquareFreeQ[n],Sow[If[Mod[n,4]==1,n,4n]]][2,1]](*Jean-François Alcover_,2023年3月22日*)
%Y参考A144338。
%K nonn,简单
%O 1,1号机组
%A _N.J.A.斯隆_
%E更多来自C.Ronaldo(aga_new_ac(AT)hotmail.com)的条款,2005年1月3日
|