|
| |
|
|
A037074
|
|
Numbers that are the product of a pair of twin primes.
|
|
75
|
|
|
|
15, 35, 143, 323, 899, 1763, 3599, 5183, 10403, 11663, 19043, 22499, 32399, 36863, 39203, 51983, 57599, 72899, 79523, 97343, 121103, 176399, 186623, 213443, 272483, 324899, 359999, 381923, 412163, 435599, 656099, 675683, 685583, 736163
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,1
|
|
|
COMMENTS
|
Each entry is the product of p and p+2 where both p and p+2 are prime, i.e., the product of the lesser and greater of a twin prime pair.
Except for the first term, all entries have digital root 8. - Lekraj Beedassy, Jun 11 2004
The above statement follows from p > 3 => (p,p+2) = (6k-1,6k+1) => p*(p+2) = 36k^2 - 1 == 8 (mod 9), and A010888 === A010878 (mod 9). - M. F. Hasler, Jan 11 2013
Albert A. Mullin states that m is a product of twin primes iff phi(m)*sigma(m) = (m-3)*(m+1), where phi(m) = A000010(m) and sigma(m) = A000203(m). Of course, for a product of distinct primes p*q we know sigma(p*q) = (p+1)*(q+1) and if p, q, are twin primes, say q = p + 2, then sigma(p*q) = (p+1)*(q+1) = (p+1)*(p+3). - Jonathan Vos Post, Feb 21 2006
Also the area of twin prime rectangles. A twin prime rectangle is a rectangle whose sides are components of twin prime pairs. E.g., the twin prime pair (3,5) produces a 3 X 5 unit rectangle which has area 15 square units. - Cino Hilliard, Jul 28 2006
If m is the product of twin primes, then sigma(m) = m + 1 + 2*sqrt(m + 1), phi(m) = m + 1 - 2*sqrt(m + 1). pmin(m) = sqrt(m + 1) - 1, pmax(m) = sqrt(m + 1) + 1. - Wesley Ivan Hurt, Jan 06 2013
|
|
|
REFERENCES
|
Albert A. Mullin, "Bicomposites, twin primes and arithmetic progression", Abstract 04T-11-48, Abstracts of AMS, Vol. 25, No. 4, 2004, p. 795.
|
|
|
LINKS
|
|
|
|
FORMULA
|
|
|
|
EXAMPLE
|
a(2)=35 because 5*7=35, that is (5,7) is the 2nd pair of twin primes.
|
|
|
MAPLE
|
ZL:=[]: for p from 1 to 863 do if (isprime(p) and isprime(p+2) ) then ZL:=[op(ZL), (p*(p+2))]; fi; od; print(ZL); # Zerinvary Lajos, Mar 07 2007
for i from 1 to 150 do if ithprime(i+1) = ithprime(i) + 2 then print({ithprime(i)*ithprime(i+1)}); fi; od; # Zerinvary Lajos, Mar 19 2007
|
|
|
MATHEMATICA
|
s = Select[ Prime@ Range@170, PrimeQ[ # + 2] &]; s(s + 2) (* Robert G. Wilson v, Feb 21 2006 *)
(* For checking large numbers, the following code is better. For instance, we could use the fQ function to determine that 229031718473564142083 is in this sequence. *) fQ[n_] := Block[{fi = FactorInteger[n]}, Last@# & /@ fi == {1, 1} && Differences[ First@# & /@ fi] == {2}]; Select[ Range[750000], fQ] (* Robert G. Wilson v, Feb 08 2012 *)
Times@@@Select[Partition[Prime[Range[500]], 2, 1], Last[#]-First[#]==2&] (* Harvey P. Dale, Oct 16 2012 *)
|
|
|
PROG
|
(PARI) g(n) = for(x=1, n, if(prime(x+1)-prime(x)==2, print1(prime(x)*prime(x+1)", "))) \\ Cino Hilliard, Jul 28 2006
(Magma) [p*(p+2): p in PrimesUpTo(1000) | IsPrime(p+2)]; // Bruno Berselli, Jul 08 2011
(Magma) IsSemiprime:=func<n | &+[d[2]: d in Factorization(n)] eq 2>; [s: n in [1..500] | IsSemiprime(s) where s is 4*n^2-1]; // Vincenzo Librandi, Apr 13 2013
(Haskell)
|
|
|
CROSSREFS
|
|
|
|
KEYWORD
|
nice,nonn
|
|
|
AUTHOR
|
|
|
|
EXTENSIONS
|
|
|
|
STATUS
|
approved
|
| |
|
|
