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”).

A179993
Numbers m with the property that, when a and b are positive integers such that a*b = m, |a-b| is prime.
3
3, 8, 14, 18, 38, 62, 98, 138, 230, 258, 278, 318, 338, 390, 398, 402, 458, 542, 678, 710, 770, 798, 822, 938, 1022, 1118, 1202, 1238, 1298, 1322, 1490, 1622, 1658, 2030, 2222, 2238, 2378, 2438, 2522, 2558, 2618, 2858, 2910, 3002, 3218, 3230, 3698, 4058, 4178
OFFSET
1,1
COMMENTS
All numbers in this sequence are congruent to 0 or 2 mod 3.
It is not known if this sequence is infinite. For n > 1 all terms are even.
The intersection with A080715 seems to be empty. Is this provable ?
From Amiram Eldar, Nov 15 2021: (Start)
The nonsquarefree terms of this sequence, 8, 18, 98, 338, ..., are numbers of the form 2*p^2, where p is in A349327.
The least terms with 1, 2, 3, 4 and 5 distinct prime divisors are 3, 14, 138, 390 and 13576178, respectively. Are there terms with more than 5 distinct prime divisors? (End)
All terms have either 6 (for a(n) = 2*A349327^2) or 2^k (for a(n) in A005117) divisors. - Samuel Harkness, Mar 02 2023
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Alois P. Heinz)
EXAMPLE
Example : For n = 5, the possible values of |a-b| are 17 = 19-2 and 37 = 38-1.
MATHEMATICA
m=1; While[m < 10000, m++; If[Mod[m, 3] == 1, , V = Divisors[m]; L = Length[V]; j = 0; While[j < L/2, j++; x = (m/V[[j]]) - V[[j]]; If[PrimeQ[x], , j = L]]; If[j == L/2, X = Append[X, m], ]]]; X
q[n_] := AllTrue[Divisors[n], #^2 > n || PrimeQ[Abs[# - n/#]] &]; Select[Range[4000], q] (* Amiram Eldar, Nov 15 2021 *)
PROG
(Python)
from itertools import islice, takewhile, count
from sympy import isprime, divisors
def A179993(): # generator of terms
for m in count(1):
if all(isprime(m//a-a) for a in takewhile(lambda x: x*x <= m, divisors(m))):
yield m
A179993_list = list(islice(A179993(), 20)) # Chai Wah Wu, Nov 15 2021
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Emmanuel Vantieghem, Aug 05 2010, Aug 06 2010
STATUS
approved