%I M2884 N2039#54 2023年10月16日23:41:29

%S 3,11,23,83131179191239251359419431443491659683719743，

%电话：9111019103111031223143914511499151115591583181119312003，

%电话：20392063233923512399245925432699281929032939296330239299

%N个卢卡斯素数：p=3（mod 4），其中2*p+1个素数。

%C 2*p+1除以A000225（p），即第p个梅森数_Lekraj Beedassy，2003年6月23日

%C还素数p，使得2^（2*p+1）-1除以2^_Arkadiusz Wesolowski，2011年5月26日

%A005384（Sophie Germain素数）与A002145的交点_Jeppe Stig Nielsen，2020年8月3日

%D A.J.C.坎宁安，《关于梅森的数字》，《英国科学促进协会的报告》，1894年，第563-564页。

%D·L·E·迪克森，《数字理论史》。卡内基公共研究所。256，华盛顿特区，第1卷，1919年；第2卷，1920年；1923年第3卷，见第1卷，第27页。

%D Daniel Shanks，“数论中已解决和未解决的问题”，第四版，切尔西出版公司，纽约，1993年，第28页。

%D N.J.A.斯隆，《整数序列手册》，学术出版社，1973年（包括该序列）。

%D N.J.A.Sloane和Simon Plouffe，《整数序列百科全书》，学术出版社，1995年（包括该序列）。

%H Marius A.Burtea，n表，n（n）表示n=1..10000（术语n=1..1000来自T.D.Noe）

%t选择[Range[10^4]，Mod[#，4]==3&&PrimeQ[#]&&PrimeQ[2#+1]&]

%t选择[Prime[Range[500]]，Mod[#，4]==3&&PrimeQ[2#+1]&]（*哈维·P·戴尔，2016年3月15日*）

%o（PARI）是（n）=n%4>2&&isprime（n）&&isprime（2*n+1）\\查尔斯·格里特豪斯IV_，2013年4月1日

%o（岩浆）[p:p in PrimesUpTo（6000）|IsPrime（2*p+1）and p mod 4 in[3]；//_文森佐·利班迪（Vincenzo Librandi），2016年9月3日

%o（MATLAB）p=素数（1500）；m=1；

%o对于u=1：长度（p）

%o if and（i素数（2*p（u）+1）==1，mod（p（u，4）==3）；溶胶（m）=p（u）；m=m+1；结束；

%o端

%o溶胶%_Marius A.Burtea，2019年3月26日

%Y参考A002145、A005384。

%K nonn公司

%O 1,1号机组

%A _N.J.A.斯隆_

%E更多条款摘自Robert G.Wilson v_，2002年3月7日