login
Largest prime power that divides the average of twin-prime pairs.
1

%I #11 May 25 2024 04:23:21

%S 4,3,4,9,5,7,5,9,17,27,23,25,9,64,11,19,16,27,47,13,29,7,27,11,29,19,

%T 25,103,107,11,81,137,23,13,49,17,43,25,59,13,128,41,71,43,31,11,17,

%U 121,19,31,67,81,139,283,41,149,16,313,23,25,37,169,347,29,64

%N Largest prime power that divides the average of twin-prime pairs.

%H Amiram Eldar, <a href="/A118701/b118701.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = A034699(A014574(n)). - _Amiram Eldar_, May 25 2024

%e a(10) = 27 since the 10th twin-prime pair is (107, 109) and 108 = 4*27.

%t Max[Power @@@ FactorInteger[#]] & /@ (Select[Prime[Range[350]], PrimeQ[# + 2] &] + 1) (* _Amiram Eldar_, May 25 2024 *)

%o (PARI) maxpow(n) = {my(f = factor(n)); if(n == 1, 1, vecmax(vector(#f~, i, f[i,1]^f[i,2])));}

%o lista(pmax) = {my(prev = 2); forprime(p = 3, pmax, if(p == prev + 2, print1(maxpow(p - 1), ", ")); prev = p);} \\ _Amiram Eldar_, May 25 2024

%Y Cf. A014574, A034699, A060210, A246655.

%K nonn

%O 1,1

%A _Greg Huber_, Jul 13 2006

%E More terms from _Amiram Eldar_, May 25 2024