OFFSET
1,1
COMMENTS
From Antti Karttunen, Nov 29 & 30 2013: (Start)
Proof:
A000069(n) reduces according to the bit parity of n-1 as follows:
A000069(n) = 2n - 2 when n-1 is odious.
A000069(n) = 2n - 1 when n-1 is evil.
which means that no prime in this sequence can be evil, as then p-1 would be an odious number (true for all odd primes) and A000069(p) would be 2(p-1) which obviously cannot be a prime, contradicting the requirement. Thus all primes present must belong to the set of odious primes, A027697.
As each prime p here is thus odious, it means that each p-1 is an evil number (A001969), and thus A000069(p) = 2p-1. And the stipulation that it also must be prime, is just what is required from the terms of A005382. Thus this sequence contains exactly those primes that occur in both A005382 and A027697.
Equally: this is the intersection of A000069 and A005382, thus prime p occurs here iff A000120(p) is odd and 2p-1 is prime also.
Also, apart from the first term (2), all the primes (2*a(n))-1 are also odious. This follows because for any odd number k, A000120(2k-1) = A000120(k).
(End)
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..225
EXAMPLE
7 is a prime and A000069(7) = 13, a prime also, thus 7 is in this sequence.
19 is a prime and A000069(19) = 37, a prime also, thus 19 is in this sequence.
Alternatively:
7 is a prime, 2*7-1 = 13 is also prime, and when written in binary, 7 = '111', with an odd number of 1-bits. Thus 7 is included in this sequence.
The next time this happens, is for 19, as it is a prime, 2*19-1 = 37 is also prime, and when written in binary, 19 = '10011', also has on odd number of 1-bits.
PROG
(Scheme, with Antti Karttunen's IntSeq-library and Aubrey Jaffer's SLIB-library)
(require 'factor) ;; Includes predicate prime? from SLIB-library.
;; Implementation based on the original definition:
;; Alternative implementation based on the other definition:
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Nov 27 2013
EXTENSIONS
Edited and erroneous terms removed by Antti Karttunen, Nov 29-30 2013
STATUS
approved