%I#19 2014年10月16日21:10:38

%S 15,35,51,65,77,91115123141161185187201209219221235259，

%电话26730130532133934135536537381403413427437451453481，

%电话：4854975015375455896496676816896997177217237377457497637898358438498901905

%具有两个因子的半素数平均值的N个半素数。

%C该定义类似于A115585，但考虑的是算术平均值，而不是系数的总和。

%C即使是半素数A100484也不在序列中，因为（除4外）它们因子的平均值不是整数。

%H Alois P.Heinz，n的表格，n=1..1000时的a（n）</a>

%e半素数187=11*17在序列中，因为平均值（11+17）/2=14=2*7是半素数。

%e半素数267=3*89在序列中，因为平均数（3+89）/2=46=2*23是半素数。

%t半素数Q[n_]：=总数[FactorInteger[n]][2]]==2；Reap[Do[{p，e}=Transpose[FactorInteger[n]]；如果[Total[e]==2&&semiPrimeQ[Total[p]/2]，Sow[n]]，{n，1000}]][[2，1]]

%o（PARI）sopf（n）={局部（f，s=0）；f=因子（n）；对于（i=1，矩阵大小（f）[1]，s+=f[i，1]）；返回}

%o平均值（n）=｛局部（s）；s=sopf（n）/ω（n）；return（s）｝

%o{对于（n=4，10^3，m=averg（n）；如果（bigomega（n）==2，如果（m==楼层（m）&&bigomegam（m）==2，打印1（n，“，”））}

%o//安东尼奥·罗尔丹，2012年10月15日

%o（PARI）列表（lim）=我的（v=列表（））；对于素数（p=3，lim/3，对于素数，q=3，min（p-2，lim\p），if（bigomega（（p+q）/2）==2，listput（v，p*q））；集（v）\\_Charles R Greathouse IV_，2014年10月16日

%Y参考A001358。

%K nonn公司

%O 1,1号机组

%A _Antonio Roldán，2011年3月9日