A356627 Primes whose powers appear in A332979.

2, 3, 5, 7, 11, 17, 29, 37, 41, 59, 67, 71, 97, 127, 149, 191, 223, 269, 307, 347, 419, 431, 557, 563, 569, 587, 593, 599, 641, 727, 809, 937, 967, 1009, 1213, 1277, 1423, 1861, 1973, 2237, 2267, 2657, 3163, 3299, 3449, 3457, 3527, 3907, 4001, 4211, 4441, 4637

Offset: 1

Author

Michael De Vlieger, Sep 27 2022

Keywords

• nonn

Comments

Maxima of row n > 0 of A005940, A182944, and A182945 are powers of these primes.

Indices k of primes, A000040(k), listed here show an interesting correlation with the function f(k) = A000040(k) - A302334(k). - Peter Munn, Sep 29 2022

Examples

			5 | A332979(5..7), thus 5 is in the sequence.
7 | A332979(8), thus 7 is in the sequence.
13 does not divide any term in A332979, so it is not a term in this sequence.
		

Links

• Michael De Vlieger, Table of n, a(n) for n = 1..100
• Michael De Vlieger, Plot A332979(n) = p^e at (e, pi(p)) for n = 1..360. The primes p are labeled in bold and appear in this sequence. The least and greatest exponents of p^e are labeled in italic.

Crossrefs

Cf. A000040, A000961, A005940, A180944, A180945, A302334, A332979.

Programs

Mathematica

Prime@ Union@ Table[MaximalBy[Table[{k, n - k}, {k, n}], Prime[#1]^#2 & @@ # &][[1, 1]], {n, 2^10}]
(* or use concise file in A332979 *)
Prime /@ Union@ Rest@ Map[ToExpression@ StringTrim[#, "p"] & @@ StringSplit[#, "^"] &, Import["https://oeis.org/A332979/a332979.txt", "Data"][[All, -1]]]