OFFSET
1,2
COMMENTS
For any n > 0, a(n) <= n * (A053669(n) - 1).
a(n) = 1 iff omega(n) < omega(n+1).
a(p) = 1 for any prime power p not in A006549.
The scatterplot of the sequence shows segments of slope -1, corresponding to frequent values of n+a(n); these segments correspond to the strands in the plot of the ordinal transform of n+a(n) (see plots in Links section).
LINKS
Rémy Sigrist, Logarithmic scatterplot of the first 10000 terms
Rémy Sigrist, Pin plot of the first 10000 terms
Rémy Sigrist, Ordinal transform of the first 10000 terms of n+a(n)
EXAMPLE
For n=2:
- omega(2+1) = 1 = omega(2),
- omega(2+2) = 1 = omega(2),
- omega(2+3) = 1 = omega(2),
- omega(2+4) = 2 > max(omega(2), omega(4)) = 1,
- hence, a(2) = 4.
PROG
(PARI) a(n) = my (on=omega(n)); for (k=1, oo, if (omega(n+k) > max(on, omega(k)), return (k)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Rémy Sigrist, Oct 26 2017
STATUS
approved