login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A136474 Primes that divide 2^(3^n)+1 for some n. 3
3, 19, 163, 1459, 17497, 52489, 87211, 135433, 139483, 1220347, 5419387, 6049243, 28934011, 86093443, 227862073, 272010961 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

This sequence is a subsequence of A057719.

272010961 is the last term less than 3*10^9. The n for each prime is 0, 2, 4, 5, 7, 8, 3, 4, 5, 9, 7, 7, 8, 16, 6, 4. Some terms from A111974 are in this sequence also: 411782264189299, 116299474006080119380780339, and 84782316550432407028588866403. If p=2*3^k+1 is prime for an even k, then p is in this sequence.

LINKS

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

EXAMPLE

1220347 belongs to the sequence as it is a factor of 2^(3^9)+1 (This is the largest member of the sequence less than 5000000)

MAPLE

with(numtheory):L:=3; for p from 5 to 5000000 do if isprime(p) then q:=op(2, ifactors(order(2, p))); if nops(q)=2 then if op(1, op(1, q))=2 and op(2, op(1, q))=1 and op(1, op(2, q))=3 then L:=L, p; fi; fi; fi; od; L;

MATHEMATICA

Reap[Do[p=Prime[n]; mo=MultiplicativeOrder[2, p]; If[EvenQ[mo] && IntegerQ[Log[3, mo/2]], Sow[p]], {n, PrimePi[10^7]}]][[2, 1]]

CROSSREFS

Cf. A057719, A136475, A005621, A136473.

Sequence in context: A057719 A289258 A199559 * A337818 A183069 A256710

Adjacent sequences:  A136471 A136472 A136473 * A136475 A136476 A136477

KEYWORD

more,nonn

AUTHOR

Christopher J. Smyth, Feb 16 2008

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified September 24 08:31 EDT 2021. Contains 347623 sequences. (Running on oeis4.)