%I#22 2023年11月21日05:23:06

%S 1,1,3,4,1,3,13,4,9,1,12,25,13,3,16,1,9,37,4,39,11,12,25,25,27,52，

%电话：1,3,61,16,3,1,13,36,73,37,75,4,1,39,85,4,9,1,1,48133,25,3100,1,27，

%U 1,52111,1,12121,61117,64,25,3133,4,3,13

%N x^2+xy+y^2==0（mod N）的解的数目。

%H Andrew Howroyd，n的表，n=1..10000的a（n）</a>

%F与a（3^e）=3^e相乘，如果p mod 3=1，a（p^e）=（（p-1）*e+p）*p^（e-1），如果p mod 3=2，a（p^e）=p^（2*floor（e/2））。-_Vladeta Jovovic，2003年9月27日

%F和{k=1..n}a（k）~c*n^2/2，其中c=A073010/A086724=0.77383581325017004332….-_Amiram Eldar，2023年11月21日

%p A087694:=过程（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 mod 3=1，则（（p-1）*e+p）*p^（e-1）；否则p^（2*层（e/2））；结束条件：；else mul（进程名（op（1，p）^op（2，p）），p=pf）；结束条件：；结束条件：；结束进程：

%p序列（A087694（n），n=1..70）；#_R.J.Mathar，2011年1月7日

%t a[n]：=如果[n=1，1，乘积[｛p，e｝=pe；哪个[p=3，3^e，Mod[p，3]==2，（p^2）^商[e，2]，True，（（p-1）e+p）p^（e-1）]，｛pe，FactorInteger[n]｝]]；

%t a/@Range[1100]（*_Jean-François Alcover_，2019年9月20日，来自PARI*）

%o（PARI）a（n）={my（f=系数（n））；prod（i=1，#f~，my（p=f[i，1]，e=f[i,2]）；if（p=3，3^e，if（p%3==2，（p^2）^（e\2），（（p-1）*e+p）*p^（e-1）））}\\ Andrew Howroyd_，2018年7月9日

%Y参考A000086、A073010、A086724。

%K mult、nonn、easy

%氧1,3

%A Yuval Dekel（dekelyuval（AT）hotmail.com），2003年9月27日