OFFSET
1,2
COMMENTS
Subsequence of A246445.
Generated by k = 1, 2, 3, 4, 5, 6, 7, 11, 12, 13, 15, 17, 19, 23, 28, 29, 31,. ..
This set of k contains all terms of A122781 and all primes. [It contains the primes because j^p == j (mod p) for every integer j if p is prime; see e.g. the corollary 4.4 to the Lagrange theorem in Jones et al.]
LINKS
G. A. Jones and J. M. Jones, Congruences with a prime-power modulus, p 65-82 in "Elementary Number Theory", Springer Undergraduate Mathematics Series, (1988).
EXAMPLE
a(9) = 1398101 because (4^12 - 4)/12 = 1398101 for k = 12.
PROG
(PARI) lista(nn) = {for (k=1, nn, va = (4^k - 4)/k; if (type(va) == "t_INT", print1(va, ", ")); ); } \\ Michel Marcus, Sep 12 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Sep 11 2014
STATUS
approved