login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A185718
For positive n with prime decomposition n = Product_{j=1..m} (p_j^k_j) define A_n = Sum_{j=1..m} (p_j*k_j) and B_n = Sum_{j=1..m} (p_j^k_j). This sequence gives those n for which A_n and B_n are both prime and B_n = A_n + 2 (i.e., form a twin prime pair).
1
40, 88, 184, 424, 808, 1048, 1384, 1528, 1864, 2104, 2184, 3080, 4504, 4744, 5224, 5928, 6440, 6568, 7224, 8104, 8360, 8840, 9784, 10264, 10472, 11480, 11544, 11848, 12808, 12904, 14136, 14840, 14968, 16280, 16648, 18664, 19608, 20344, 21080, 22040, 23240, 23704, 24440, 24648, 24920, 26008, 26584, 27384, 27608, 27688, 28264, 28952, 29240
OFFSET
1,1
COMMENTS
Assuming the twin prime conjecture, my advisor and I are able to prove there are infinitely many of these pairs. In other words, there are infinitely many n such that A_n and B_n are prime and B_n = A_n + 2.
From Bobby Browning, Feb 14 2011: (Start)
8*A046138 is a subsequence of A185718 for the following reasons:
i) the n in A185718 for which A_n and B_n form a twin prime pair are of the form n=2^3*p_1*p_2*...p_k.
ii) the A046138 sequence consists of primes p such that p+6 and p+8 form a twin prime pair.
iii) so if p is a prime such that p+6 and p+8 form a twin prime pair and n = 2^3*p then A_n = p+6 and B_n = p+8. Thus, the integers such that n = 2^3*p are a subsequence of A185718. (End)
EXAMPLE
a(1) = 40 = 2^3*5^1, with a = 11 and b = 13.
a(2) = 88 = 2^3*11^1 with a = 17 and b = 19.
a(3) = 184 = 2^3*23^1 with a = 29 and b = 31.
MATHEMATICA
okQ[n_] := Module[{p, e, a, b}, {p, e} = Transpose[FactorInteger[n]]; a = Plus @@ (p*e); b = Plus @@ (p^e); b == a + 2 && PrimeQ[a] && PrimeQ[b]]; Select[Range[30000], okQ]
PROG
(PARI) sopfr(n)=my(f=factor(n)); sum(i=1, #f[, 1], f[i, 1]*f[i, 2]);
forstep(n=8, 1e5, 16, if(issquarefree(n/8)&&isprime(k=sopfr(n))&isprime(k+2), print1(n", ")))
CROSSREFS
Cf. A001414 (A_n), A008475 (B_n).
Sequence in context: A243803 A203855 A188335 * A043414 A044178 A044559
KEYWORD
nonn
AUTHOR
Bobby Browning, Feb 10 2011
STATUS
approved