#从百科全书的问候序列行!本次搜索:id:a087694 展示1-1的1个1 ;%I a087694;%S a087694%S a087694 1,1,1,3,4,4,1,3,4,1,3,13,4,1,3,13,4,1,3,12,1,1,12,25,1,12,25,1,13,3,16,1,1,9,37,4,39,39,1,4,39,1,1,1,39,1,1,1,12,25,25,27,27,52,;%T a087694,3,3,61,16,16,3,1,1,13,36,73,37,37,37,75,4,4,,9,1,1,48133,25,3100,1,27, %U a087694 1,52111,1,1,12121,61117,64,25,3133,4,3,13 %N a087694 x^2+xy+y^2==0(mod N)的解数。 %HA087694安德鲁·豪罗伊德,n=1..10000的n,a(n)表%F A087694与a(3^e)=3^e相乘,a(p^e)=((p-1)*e+p)*p^(e-1)如果p mod 3=1,a(p^e)=p^(2*地板(e/2)),如果p mod 3=2。-2003年9月27日,;%p A087694 A087694 A087694 A087694:=proc(n)选项记住;局部pf,p,f,e;如果n=1则1;else pf:=ifactors(n)[2];如果nops(pf)=1则f:=op(1,pf);p:=op(1,f);e:=op(2,f);p=3则n;elif p p mod 3=1然后(p-1)*e+p)*p^(e-1);其他p^(2*地板(e/2))(2*地板(e/2)))p(2*地板(e/2));p:=p=op(2,p(1,f(1,f)p:=op(2,f;end if;else mul(procname(op(1,p)^op(2,p)),p=pf);结束end if;end if;end proc: %p A087694 seq(A087694(n,n=1..70);\35; U R.J.Mathar_,2011年1月7日;%t A087694 a[n U]:=if[n==1,1,1,产品[{p,e}=pe;其中[p==3,3^e,Mod[p,3]==2,(p^2)^商[e,2],真的,((p-1)e+p)p^(e-1)],{pe,pe,1,产品[{pe,1,e{p==1,e{p=p=p=p=p=p保理商[n]}]]; %t A087694 a/@范围[1100](*u Jean-François Alcover,2019年9月20日,来自巴黎*);%o A087694(p=3,3^e,if(p%3==2,(p^2)^(e\2),((p p-3=3,3^e,if(p%p%3==2,(p^2)^(e\2),((p p-1)*e+p*p^(e-1)))))))}安德鲁·霍罗伊德,2018年7月9日,2018年7月9日,;%Y A087694 Cf A000000086。;%K A087694 mult,NONON A087694 mult,NONONON A087694 mult,NONONONONON A087694 mult,NONONONON(e-1)*(e-1)*(e-1)*(e-1)*(e-1)p(e-n %O A087694 1,3 %a A087694 Yuval Dekel(dekelyuval(AT)hotmail.com),2003年9月27日 #内容根据OEIS最终用户许可协议提供:http://OEIS.org/License