%I#20 2020年2月15日23:22:15

%第6971241151316491921232933281364939613977409743694633页，

%电话：4777532156175729600164976817695370817361763377697913，

%电话：8249863388579553967398099811005710081102171048910537

%N正好两个不同素数的乘积，与1模8（A007519）一致。

%C半素子集A001358。{d=p_1*p_2*…*p_m的子集，其中p_i==1（mod 8），1<=i<=m是不同的素数}，如Wei，p.2所示。

%H Vincenzo Librandi，n的表，n=1..1000的a（n）</a>

%H Dasheng Wei，<a href=“http://arxiv.org/abs/102.3811“>关于方程x^2-Dy^2=n</a>，2011年2月18日。

%对于i<j}，F{A007519（i）*A007518（j）。

%对于i<j，F{A000040（i）*A000040。

%e 10001位于该序列中，因为10001=73*137=A007519（3）*A007518（7）。

%t p=选择[Prime[Range[200]]，Mod[#，8]==1&]；排序[Reap[Do[n=p[i]]p[[j]]；如果[n<=p[[1]]p[[-1]]，母猪[n]]，{i，2，长度[p]}，{j，i-1}]][[2,1]]]

%o（PARI）列表（lim）=我的（v=列表（），P=列表（（），t）；对于素数（p=2，lim\17，if（p%8==1，listput（p，p）））；对于（i=2，#P，my（P=P[i]）；对于（j=1，i-1，t=p*p[j]；如果（t>lim，break）；列表（v，t））；集（v）\\_Charles R Greathouse IV_，2016年7月3日

%Y参见A000040、A001358、A007519、A017077。

%K nonn，简单

%O 1,1号机组

%A _Jonathan Vos Post，2011年2月20日