%I#37 2022年8月15日02:07:44

%序号1,4,7,23,41,79,89,97127167191223239241257281439577593

%N形式为x^2+xy+y^2的Pell数指数。

%这个序列的项大多数时候是质数。

%e4是一个术语，因为A000129（4）=12和12=2^2+2*2+2^2。

%o（PARI）a000129（n）=（[2,1；1,0]^n）[2,1]；

%o是A003136（n）=#bnfisint范数（bnfinit（z^2+z+1），n）；

%o表示（n=1，1e3，if（isA0003136（a000129（n）），print1（n，“，”））；

%Y参考A000129、A003136。

%K nonn，更多

%O 1,2号机组

%A Altug Alkan，2016年3月18日

%E a（17）-a（19），来自王金源，2022年8月14日