A352083 Numbers k such that (3^k - k^3)/2 is prime.

5, 13, 205, 409, 413, 545, 88649

Offset

1,1

Comments

All terms must be odd. - Michael S. Branicky, Mar 03 2022

No further terms less than 25000. - Michael S. Branicky, Mar 04 2022

Links

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

Example

5 is a term since (3^5 - 5^3)/2 = 59 is prime.

Mathematica

Select[Range[1, 600, 2], PrimeQ[(3^# - #^3)/2] &] (* Amiram Eldar, Mar 03 2022 *)

Prog

(Sage) for n in srange(1, 10000):
             if ((3^n-n^3)//2).is_prime():
                  print(n)
(PARI) isok(k) = if (denominator(x=(3^k-k^3)/2) == 1, ispseudoprime(x)); \\ Michel Marcus, Mar 03 2022
(Python)
from sympy import isprime
def afind(limit):
    for k in range(1, limit+1, 2):
        if isprime((3**k - k**3)//2):
            print(k, end=", ")
afind(1000) # Michael S. Branicky, Mar 03 2022

Crossrefs

Sequence in context: A117077 A124924 A209271 * A355763 A124878 A085554

Adjacent sequences: A352080 A352081 A352082 * A352084 A352085 A352086

Keyword

nonn,more

Author

Brennan G. Benfield, Mar 03 2022

Extensions

a(7) from Michael S. Branicky, Oct 03 2024

Status

approved