OFFSET
1,1
COMMENTS
Subsequence of A309769 and A309780. a(1) is the only nonsquarefree known term. Conjecture: For n > 1 the minimal condition requires m to be squarefree and for every odd prime divisor p of m to be such that m/p - 1 is composite with least prime divisor q < p (k=p-q). No term is divisible by 3.
Not all terms are coprime to 7 and the terms aren't all 97-smooth. For n = 16..18, there are upper bounds on a(n): 16808251841257353520347590, 1627869069994521415245268370, 202089221222977079276742661490. - David A. Corneth, Sep 26 2019
EXAMPLE
a(2) = 110 = (2*5)*11; q = 3 < 11; also 110=(2*11)*5; q = 3 < 5.
MATHEMATICA
kQ[n_, p_] := Module[{ans = False}, Do[If[Divisible[n - k, p - k], ans = True; Break[]], {k, 1, p - 2}]; ans]; aQ[n_] := EvenQ[n] && Length[(p = FactorInteger[ n][[2 ;; -1, 1]])] > 0 && AllTrue[p, kQ[n, #] &]; oddomega[n_] := PrimeNu[n / 2^IntegerExponent[n, 2]]; s = {}; om = 1; Do[If[oddomega[n] == om && aQ[n], AppendTo[s, n]; om++], {n, 2, 10^16, 2}]; s (* Amiram Eldar, Aug 17 2019 *)
PROG
(PARI) getk(p, m) = {for (k=1, p-2, if (((m-k) % (p-k)) == 0, return(k)); ); }
isok1(m) = {if ((m % 2) == 0, my(f = factor(m)[, 1]~); if (#f == 1, return (0)); for (i=2, #f, if (!getk(f[i], m), return(0)); ); return (1); ); }
isok(k, n) = (omega(k/(2^valuation(k, 2))) == n) && isok1(k);
a(n) = {my(k=2*prod(k=2, n+1, prime(k))); while (!isok(k, n), k+=2); k; } \\ Michel Marcus, Aug 27 2019
CROSSREFS
KEYWORD
nonn,more
AUTHOR
David James Sycamore, Aug 17 2019
EXTENSIONS
a(8) from Michel Marcus, Sep 25 2019
a(9)-a(15) from David A. Corneth, Sep 26 2019
STATUS
approved