This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A336985 #17 Aug 31 2020 17:16:12 %S A336985 3,5,7,20,31,53,86,97,108,110,121,132,143,154,165,176,187,198,209,211, %T A336985 222,233,244,266,277,288,299,310,323,334,345,356,367,389,400,411,413, %U A336985 424,435,446,457,468,479,490,501,512,514,536,547,558,569,580,591,602,613 %N A336985 Colombian numbers that are not Bogotá numbers. %C A336985 Equivalently, numbers m that are not of the form k + sum of digits of k for any k (A003052), and that are not of the form q * product of digits of q for any q (complement of A336826). %C A336985 As repunits are trivially Bogotá numbers, there are not repunits in the data. %C A336985 A336983, A336984, A336986 and this sequence form a partition of the set of positive integers N* %H A336985 Giovanni Resta, <a href="http://www.numbersaplenty.com/set/self_number/">Self or Colombian number</a>, Numbers Aplenty. %H A336985 <a href="https://oeis.org/wiki/Index_to_OEIS:_Section_Coi#Colombian">Index to sequences related to Colombian or self numbers</a>. %e A336985 7 is a term because there are not k < 7 such that 7 = k + sum of digits of k, and that are not q such that 7 = q * product of digits of q. %e A336985 13 is not of the form q * product of digits of q for any q <= 13, so 13 is not a Bogotá number, but 13 = 11 + (1+1) is not Colombian, hence 13 is not a term. %e A336985 42 is Colombian because there does not exist m < 42 such that 42 = m + sum of digits of m, but as 42 = 21 * (2*1) is a Bogota number, 42 is not a term. %t A336985 m = 600; Intersection[Complement[Range[m], Select[Union[Table[n + Plus @@ IntegerDigits[n], {n, 1, m}]], # <= m &]], Complement[Range[m], Select[Union[Table[n * Times @@ IntegerDigits[n], {n, 1, m}]], # <= m &]]] (* _Amiram Eldar_, Aug 26 2020 *) %o A336985 (PARI) lista(nn) = Vec(setintersect(setminus([1..nn], Set(vector(nn, k, k+sumdigits(k)))), setminus([1..nn], Set(vector(nn, k, k*vecprod(digits(k))))))); \\ _Michel Marcus_, Aug 26 2020 %Y A336985 Cf. A003052 (Colombian), A176995 (not Colombian), A336826 (Bogotá numbers), A336983 (Bogotá not Colombian), A336984 (Bogotá and Colombian), this sequence (Colombian not Bogotá), A336986 (not Colombian and not Bogotá). %K A336985 nonn,base %O A336985 1,1 %A A336985 _Bernard Schott_, Aug 26 2020