OFFSET
1,1
COMMENTS
Primes of the forms 2^r * b^s - 1 where b = 1, 5, 7, 11, 13 are A000668 (Mersenne prime exponents), A077313, A077314, A077315 and A173062. When b = 3 we get A005105 with initial term 2.
For n > 1, all terms are congruent to 1 (mod 3).
Also, these are prime numbers p for which (34^p)/(p+1) is an integer.
EXAMPLE
With n = 1, a(1) = 2^2 * 17^0 - 1 = 3.
With n = 4, a(4) = 2^2 * 17^1 - 1 = 67.
list of (r, s): (2, 0), (3, 0), (5, 0), (2, 1), (3, 1), (7, 0), (4, 1), (1, 2), (6, 1), (3, 2), (13, 0), (4, 3), (17, 0), (19, 0), (16, 1), (13, 2), (5, 4), (20, 1), (9, 4), (6, 5).
PROG
(GAP)
K:=10^7+1;; # to get all terms <= K.
A:=Filtered(Filtered([1..K], i->i mod 3=1), IsPrime);; I:=[17];;
B:=List(A, i->Elements(Factors(i+1)));;
C:=List([0..Length(I)], j->List(Combinations(I, j), i->Concatenation([2], i)));;
A292890:=Concatenation([3], List(Set(Flat(List([1..Length(C)], i->List([1..Length(C[i])], j->Positions(B, C[i][j]))))), i->A[i]));
(PARI) isok(p) = isprime(p) && (denominator((34^p)/(p+1)) == 1); \\ Michel Marcus, Sep 27 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Muniru A Asiru, Sep 26 2017
EXTENSIONS
More terms from Jinyuan Wang, Feb 23 2020
STATUS
approved