A085971 Union of primes and numbers that are not prime powers (A000040, A024619).

2, 3, 5, 6, 7, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78

Offset

1,1

Comments

Complement of A025475;

A085972(n) = Max{k: a(k)<=n};

different from A007916 and A052485, as a(28)=36;

A085818(a(n)) = A000040(n).

Links

Table of n, a(n) for n=1..68.

Formula

a(n) = n + o(sqrt n). - Charles R Greathouse IV, Oct 19 2015

Mathematica

With[{nn=100}, Union[Join[Prime[Range[PrimePi[nn]]], DeleteCases[Range[2, 80], _?(PrimePowerQ[#]&)]]]] (* Harvey P. Dale, May 15 2019 *)

Prog

(PARI) is(n)=isprimepower(n)<2 && n>1 \\ Charles R Greathouse IV, Oct 19 2015

Crossrefs

Sequence in context: A089105 A028769 A094784 * A175082 A007916 A052485

Adjacent sequences: A085968 A085969 A085970 * A085972 A085973 A085974

Keyword

nonn,easy

Author

Reinhard Zumkeller, Jul 06 2003

Status

approved