This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A316525 #13 Sep 28 2018 15:31:01 %S A316525 2,3,4,5,7,8,9,11,13,16,17,19,20,21,23,25,27,29,31,32,33,37,41,43,44, %T A316525 47,49,53,57,59,60,61,64,67,68,69,71,73,79,81,83,85,89,93,97,101,103, %U A316525 105,107,109,112,113,116,121,125,127,128,129,131,133,137,139 %N A316525 Numbers whose average of prime factors is prime. %C A316525 Prime factors counted with multiplicity. - _Harvey P. Dale_, Sep 28 2018 %H A316525 Harvey P. Dale, <a href="/A316525/b316525.txt">Table of n, a(n) for n = 1..1000</a> %e A316525 60 = 2*2*3*5 has average of prime factors (2+2+3+5)/4 = 3, which is prime, so 60 belongs to the sequence. %t A316525 Select[Range[100],PrimeQ[Mean[If[#==1,{},Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[p,{k}]]]]]]&] %t A316525 Select[Range[200],PrimeQ[Mean[Flatten[Table[#[[1]],#[[2]]]&/@ FactorInteger[ #]]]]&] (* _Harvey P. Dale_, Sep 28 2018 *) %o A316525 (PARI) isok(n) = {my(f=factor(n)); iferr(isprime(sum(k=1, #f~, f[k,1]*f[k,2])/sum(k=1, #f~, f[k,2])), E, 0);} \\ _Michel Marcus_, Jul 06 2018 %Y A316525 Subsequence of A078175. %Y A316525 Cf. A000040, A000607, A056768, A071321, A100118, A316219, A316313. %K A316525 nonn %O A316525 1,1 %A A316525 _Gus Wiseman_, Jul 05 2018