OFFSET
1,1
COMMENTS
A number k has a coprime divisor shift s if GCD(d + s, k) = 1 for all divisors d of k.
A number k has a coprime divisor shift iff it is not divisible by any number in the sequence.
If k has no coprime divisor shift, then so is any multiple of k.
Theorem: if p is an odd prime of the form p = 3*q - 1, then 3*p is a term. Proof: p <> 3 and 3*p has 4 divisors: 1, 3, p, 3*p. A proper divisor of 3*p is either 1 or an odd prime and thus has a coprime divisor shift and is not a term. There are 3 distinct residues of the divisors of 3*p modulo 3, i.e. 0, 1, 2, which is a complete set of residues and thus 3*p has no coprime divisor shift (see A366251). The converse is not true, for instance 819 = 3^2*7*13 is a term. - Chai Wah Wu, Oct 19 2025
REFERENCES
a(1) = 2 for GCD(2 + 0, 2) > 1 and GCD(1 + 1, 2) > 1.
a(2) = 15 for GCD(3 + 0, 15) > 1, GCD(5 + 1, 15) > 1, GCD(1 + 2, 15) > 1, and any odd number between 2 and 15 has a coprime divisor shift.
LINKS
M. Farrokhi D. G., Table of n, a(n) for n = 1..10000
PROG
(Python)
from itertools import count, islice
from sympy import divisors, primerange
def A366330_gen(): # generator of terms
yield 2
a = []
for k in count(3, 2):
ds = divisors(k)
if all(k%b for b in a) and any(len({d%p for d in ds})==p for p in primerange(3, len(ds)+1)):
yield k
a.append(k)
CROSSREFS
KEYWORD
nonn
AUTHOR
M. Farrokhi D. G., Oct 07 2023
STATUS
approved
