login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A221982 Primes p == 2 (mod 5) for which 4*p+1 is also prime. 2

%I #23 Feb 12 2013 11:54:24

%S 7,37,67,97,127,277,307,487,577,727,997,1087,1297,1327,1567,1597,1777,

%T 1987,2017,2437,2647,2677,2767,2887,3037,3067,3307,3457,3637,3907,

%U 4057,4297,4447,4567,4987,5197,5527,5557,6007,6247,6337,6367,6397,6547,6577,7027,7057,7237,7417,7507,7717,7867

%N Primes p == 2 (mod 5) for which 4*p+1 is also prime.

%C The corresponding primes 4*p+1 are Chebyshev's subsequence A221981 of the primes with primitive root 10.

%D P. L. Chebyshev, Theory of congruences. Elements of number theory, Chelsea, 1972, p. 306.

%D R. K. Guy, Unsolved Problems in Number Theory, F9.

%H Paolo P. Lava, <a href="/A221982/b221982.txt">Table of n, a(n) for n = 1..10000</a>

%H P. Moree, <a href="http://arxiv.org/abs/math/0412262v2">Artin's primitive root conjecture - a survey</a>, arXiv 2004, revised 2012, p. 1.

%H <a href="/index/Pri#primes_root">Index entries for primes by primitive root</a>

%F a(n) = (A221981(n) - 1)/4.

%e 7 is a member because 7 == 2 (mod 5) and 29 = 4*7 + 1 are both prime.

%p A221982:=proc(q)

%p local n;

%p for n from 1 to q do

%p if isprime(n) and isprime(4*n+1) and (n mod 5)=2 then print(n) fi; od; end:

%p A221982 (10000); # _Paolo P. Lava_, Feb 12 2013

%t Select[ Prime[ Range[1000]], Mod[#, 5] == 2 && PrimeQ[4 # + 1] &]

%Y Cf. A001913, A006883, A045380, A106849, A221981, A222008.

%K nonn

%O 1,1

%A _Jonathan Sondow_, Feb 02 2013

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 24 22:17 EDT 2024. Contains 371964 sequences. (Running on oeis4.)