cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A097501 p^q + q^p for twin primes p and q.

This page as a plain text file.
%I A097501 #16 Sep 07 2019 09:05:08
%S A097501 368,94932,36314872537968,244552822542936127033092,
%T A097501 2177185942561672462146321298650240665136431700,
%U A097501 2246585380039521951243337580678537047744572047581514711375688196554564
%N A097501 p^q + q^p for twin primes p and q.
%C A097501 Except for the first term, 6 divides a(n). Let p = 3k+2 for odd k since k even implies p even, a contradiction. Then p = 6m + 5 and q = 6m+7 = 6m1 + 1. So p^q+q^p = (6m+5)^(6m1+1) + (6m1+1)^(6m+5) = 6H + 5^odd + 1^odd. Now 5 = (6-1) and (6-1)^odd + 1 = 6G -1 + 1 = 6G as stated. Are 3 and 17 the only primes in A051442(n)?
%e A097501 Consider the second twin prime pair (5,7). 5^7 + 7^5 = 94932, the 2nd entry.
%t A097501 lst={}; Do[p=Prime[n]; If[PrimeQ[q=p+2],a=(p^q+q^p); AppendTo[lst,a]],{n,2*4!}]; lst (* _Vladimir Joseph Stephan Orlovsky_, Aug 16 2009 *)
%t A097501 #[[1]]^#[[2]]+#[[2]]^#[[1]]&/@Select[Partition[Prime[Range[20]],2,1],#[[2]] - #[[1]]==2&] (* _Harvey P. Dale_, Sep 07 2019 *)
%o A097501 (PARI) f(n) = for(x=1,n,p=prime(x);q=prime(x+1);if(q-p==2,v=p^q+q^p;print1(v",")))
%Y A097501 Cf. A051442.
%K A097501 nonn
%O A097501 2,1
%A A097501 _Cino Hilliard_, Aug 25 2004