%I#22 2020年11月9日21:33:48
%S 2,14,33,951612794736589441139149117381826288429964318,
%电话:458550046191668378498413103141074611394131571339316013,
%电话:16566189361978320376239462705727872043088335541352323638439832456715085851363500595509756040
%N a(N)是A002145(N)的二次无余量之和(第N素数==3 mod 4)。
%H Robert Israel,n的表,n的a(n)=1..10000</a>
%H Christian Aebi,Grant Cairns,<a href=“http://arxiv.org/abs/1512.00896“>二次剩余和非剩余之和,arXiv:1512.00896[math.NT],2015。
%F a(n)=和{k=1..(A002145(n)-1)/2}(-k^2)mod A002145(n).-_J.M.Bergot和_Robert Israel,2020年11月9日
%p(数字理论):
%p a:=[];m: =[];d: =[];
%i1的p从1到200 do
%p p:=i素数(i1);
%p如果(p mod 4)=3,则
%p sp:=0;sm:=0;
%从1到p-1的j的p do
%p如果legendre(j,p)=1,则sp:=sp+j;否则sm:=sm+j;fi;od;
%p a:=[op(a),sp];m: =[op(m),sm];d: =[op(d),sm-sp];
%p fi;
%日期:
%p a;米;d、 #个A282035、A282036和A282037
%p#备选方案:
%pf:=p->添加(-k^2 mod p,k=1..(p-1)/2):
%p映射(f,选择(isprime,[seq(p,p=3..1000.4)]);#_罗伯特·伊斯雷尔(Robert Israel),2020年11月9日
%t f[p_]:=总计[范围[p-1]~补码~表[Mod[k^2,p],{k,(p-1)/2}]];f/@Select[范围[3,1000,4],PrimeQ](*_Jean-François Alcover_,2018年2月16日,在_Robert Israel_*之后)
%o(PARI)列表a(nn)=素数(p=2,nn,如果(p%4==3,打印1(总和(k=1,p-1,如果(!issquare(Mod(k,p),k)),“,”));\\_Michel Marcus,2020年11月9日
%Y残留物、非残留物及其差值之和,对于p==1 mod 4、p==3 mod 5和所有p:A171555;A282035、A282036和A282037;A076409、A125615、A282038。
%Y参考A002145。
%K非n
%O 1,1号机组
%A _N.J.A.Sloane,2017年2月20日
|