This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A271101 Squarefree semiprimes (A006881(n)) whose average prime factor is prime. 0

%I

%S 21,33,57,69,85,93,129,133,145,177,205,213,217,237,249,253,265,309,

%T 393,417,445,469,489,493,505,517,553,565,573,597,633,669,685,697,753,

%U 781,793,813,817,865,889,913,933,949,973,985,993,1057,1077,1137,1149,1177,1257,1273,1285,1329

%N Squarefree semiprimes (A006881(n)) whose average prime factor is prime.

%C Sum of factors of a(n) if semiprime (product 2*p, with p prime).

%C This sequence is subsequence of A006881, A089765, A187073, A108633 and A213015.

%C This sequence is also subsequence of A045835, because sopfr(omega(a(n))) = omega(sopfr(a(n))): sopfr(omega(a(n)))=sopfr(2)=2, and omega(sopfr(a(n)))=omega(2*p)=2 (p prime, p>2, average prime factor).

%e 133 is in the sequence because 133 is squarefree semiprime: 133=7*19, and (7+19)/2=13, prime number.

%t Select[Select[Range@ 1330, SquareFreeQ@ # && PrimeOmega@ # == 2 &], PrimeQ@ Mean[First /@ FactorInteger@ #] &] (* _Michael De Vlieger_, Mar 30 2016 *)

%o (PARI)

%o sopf(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]); return(s) }

%o {for (n=6, 2*10^3, if(bigomega(n)==2&&omega(n)==2, m=sopf(n)/2;if(m==truncate(m),if(isprime(m), print1(n, ", ")))))}

%Y Cf. A006881, A046315, A046388, A115585, A187073, A089765, A108633, A213015.

%K nonn

%O 1,1

%A _Antonio RoldÃ¡n_, Mar 30 2016

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified October 23 03:21 EDT 2019. Contains 328335 sequences. (Running on oeis4.)