# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a242720 Showing 1-1 of 1 %I A242720 #69 Nov 02 2018 09:26:14 %S A242720 12,38,80,212,224,440,440,854,1250,1460,1742,2282,2282,3434,4190,4664, %T A242720 4760,4760,6890,8054,8054,8054,12374,12830,12830,13592,13592,14282, %U A242720 17402,17402,18212,22502,22502,22502,25220,28202,28202,32234,32402,32402,38012 %N A242720 Smallest even k such that the pair {k-3,k-1} is not a twin prime pair and lpf(k-1) > lpf(k-3) >= prime(n), where lpf = least prime factor (A020639). %C A242720 The sequence is nondecreasing. See comment in A242758. %C A242720 a(n) >= prime(n)^2+3. Conjecture: a(n) <= prime(n)^4. - _Vladimir Shevelev_, Jun 01 2014 %C A242720 Conjecture. There are only a finite number of composite numbers of the form a(n)-1. _Peter J. C. Moses_ found only two: a(16)-1 = 4189 = 59*71 and a(20)-1 = 6889 = 83^2 and no others up to a(2501). Most likely, there are no others. - _Vladimir Shevelev_, Jun 09 2014 %H A242720 Peter J. C. Moses, Table of n, a(n) for n = 2..2501 %H A242720 V. Shevelev, Theorems on twin primes-dual case, arXiv:0912.4006 [math.GM], 2009-2014, (Section 10). %F A242720 Conjecturally, a(n) ~ (prime(n))^2, as n goes to infinity (cf. A246748, A246821). - _Vladimir Shevelev_, Sep 02 2014 %F A242720 For n>=3, a(n) >= (prime(n)+1)^2 + 2. Equality holds for terms of A246824. - _Vladimir Shevelev_, Sep 04 2014 %t A242720 lpf[n_] := FactorInteger[n][[1, 1]]; %t A242720 Clear[a]; a[n_] := a[n] = For[k = If[n <= 2, 2, a[n-1]], True, k = k+2, If[Not[PrimeQ[k-3] && PrimeQ[k-1]] && lpf[k-1] > lpf[k-3] >= Prime[n], Return[k]]]; %t A242720 Table[a[n], {n, 2, 50}] (* _Jean-François Alcover_, Nov 02 2018 *) %o A242720 (PARI) %o A242720 lpf(k) = factorint(k)[1,1]; %o A242720 vector(60, n, k=6; while((isprime(k-3) && isprime(k-1)) || lpf(k-1)<=lpf(k-3) || lpf(k-3)