OFFSET
1,1
COMMENTS
Both divisors d and d/k are composite, since primes p either divide or are coprime to another number, and all numbers smaller than p are coprime to p.
Consider composite k, m, k != m. Define a "neutral" relation to be such that 1 < gcd(k,m) and not equal to either k or m. Then neither k nor m divides the other, and k and m are not coprime. If k is neutral to m, then m is neutral to k, since order does not matter. Then either the squarefree kernel of one divides the other or it does not. Thus, there are 3 kinds of neutral relation:
Type A: Though gcd(k,m) > 1, k has a factor P that does not divide m, and m has a factor Q that does not divide k.
Type B: rad(k) = rad(m), yet neither k divides m nor m divides k, where rad = A007947 is the squarefree kernel.
Type C: Squarefree kernel of one number divides the other, while the other has a factor that does not divide the former.
A378769, subset of this sequence, contains numbers k that have all 3 types of neutral relation between at least 1 divisor pair (d, k/d) for each.
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Chart of (d, a(n)/d) for a(n) = 1..144, showing only the smallest d for each type of neutral relation, where type A is shown in gray, type B in black, and type C in either blue or gold.
FORMULA
EXAMPLE
a(1) = 24 = 2^3 * 3 = 4*6, both composite; gcd(4,6) = 2, 4 does not divide 6 (type C).
a(2) = 40 = 2^3 * 5 = 4*10, gcd(4,10) = 2 (type C).
a(3) = 48 = 2^4 * 3 = 6*8, gcd(6,8) = 2 (type C).
a(6) = 60 = 2^2 * 3 * 5 = 6*10, gcd(6,10) = 2 (type A).
a(12) = 96 = 2^5 * 3 = 6*16 = 8*12, both type C.
a(38) = 216 = 2^3 * 3^3 = 4*54 (type C) = 9*24 (type C) = 12*18 (type B)
a(1605) = 5400 = 2^3 * 3^3 * 5^2 = 4*1350 (type C) = 24*225 (type A) = 60*90 (type B) = A378769(1).
MATHEMATICA
nn = 300; mm = Floor@ Sqrt[nn]; p = 2; q = 3;
Complement[
Select[Range[nn], And[#2 > #1 > 1, #2 > 3] & @@ {PrimeNu[#], PrimeOmega[#]} &],
Union[Reap[
While[p <= mm, q = NextPrime[p];
While[p*q <= mm, If[p != q, Sow[p*q]]; q = NextPrime[q]];
p = NextPrime[p] ] ][[-1, 1]] ]^2 ]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Dec 24 2024
STATUS
approved