login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A070006 Composite numbers that are not a prime power and whose distinct prime divisors' arithmetic mean is a prime. 5
21, 33, 57, 63, 69, 85, 93, 99, 105, 129, 133, 145, 147, 171, 177, 189, 195, 205, 207, 213, 217, 231, 237, 249, 253, 265, 279, 297, 309, 315, 363, 387, 393, 417, 425, 441, 445, 465, 469, 483, 489, 493, 505, 513, 517, 525, 531, 553, 565, 567, 573, 585, 597 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Subsequence of A070005.

LINKS

Michael De Vlieger, Table of n, a(n) for n = 1..10000

EXAMPLE

n = 993 = 3*331, mean = 334/2 = 167, a prime.

MATHEMATICA

ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ba[x_] := Table[Part[ffi[x], 2*w-1], {w, 1, lf[x]}] Do[s=Apply[Plus, ba[n]]/lf[n]; If[PrimeQ[s]&&Greater[lf[n], 1], Print[n]], {n, 2, 1000}]

(* Second program: *)

Select[Range@ 600, And[! PrimePowerQ@ #, PrimeQ@ Mean[FactorInteger[#][[All, 1]]]] &] (* Michael De Vlieger, Jul 18 2017 *)

PROG

(PARI) lista(nn) = {for (n=2, nn, f = factor(n); if ((#f~ != 1) && (type(q=sum(k=1, #f~, f[k, 1])/#f~) == "t_INT") && isprime(q), print1(n, ", ")); ); } \\ Michel Marcus, Mar 28 2015

CROSSREFS

Cf. A001414, A008472, A070005.

Sequence in context: A330441 A176945 A333857 * A189986 A190299 A280262

Adjacent sequences:  A070003 A070004 A070005 * A070007 A070008 A070009

KEYWORD

easy,nonn

AUTHOR

Labos Elemer, Apr 11 2002

STATUS

approved

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 26 10:21 EDT 2021. Contains 346294 sequences. (Running on oeis4.)