%I #17 Aug 31 2020 17:16:12
%S 3,5,7,20,31,53,86,97,108,110,121,132,143,154,165,176,187,198,209,211,
%T 222,233,244,266,277,288,299,310,323,334,345,356,367,389,400,411,413,
%U 424,435,446,457,468,479,490,501,512,514,536,547,558,569,580,591,602,613
%N Colombian numbers that are not Bogotá numbers.
%C 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 As repunits are trivially Bogotá numbers, there are not repunits in the data.
%C A336983, A336984, A336986 and this sequence form a partition of the set of positive integers N*
%H Giovanni Resta, <a href="http://www.numbersaplenty.com/set/self_number/">Self or Colombian number</a>, Numbers Aplenty.
%H <a href="https://oeis.org/wiki/Index_to_OEIS:_Section_Coi#Colombian">Index to sequences related to Colombian or self numbers</a>.
%e 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 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 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 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 (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 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 nonn,base
%O 1,1
%A _Bernard Schott_, Aug 26 2020