# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a136333 Showing 1-1 of 1 %I A136333 #26 Feb 15 2024 01:56:01 %S A136333 1,3,7,9,11,13,17,19,31,33,37,39,71,73,77,79,91,93,97,99,111,113,117, %T A136333 119,131,133,137,139,171,173,177,179,191,193,197,199,311,313,317,319, %U A136333 331,333,337,339,371,373,377,379,391,393,397,399,711,713,717,719,731 %N A136333 Numbers containing only digits coprime to 10 in their decimal representation. %C A136333 Numbers containing digits 1,3,7,9 only, or, numbers written in base 4 (cf. A007090) with digits mapped by: 0->1, 1->3, 2->7 and 3->9. - _Reinhard Zumkeller_, Jul 17 2014 %H A136333 Reinhard Zumkeller, Table of n, a(n) for n = 1..10000 %H A136333 Robert Baillie and Thomas Schmelzer, Summing Kempner's Curious (Slowly-Convergent) Series, Mathematica Notebook kempnerSums.nb, Wolfram Library Archive, 2008. %H A136333 Index entries for 10-automatic sequences. %F A136333 Sum_{n>=1} 1/a(n) = 2.395867871130444522329053889312125689319669370758630349552737883715872077555... (calculated using Baillie and Schmelzer's kempnerSums.nb, see Links). - _Amiram Eldar_, Feb 15 2024 %t A136333 fQ[n_] := Block[{s = {1, 3, 7, 9}}, Union[Join[s, IntegerDigits@ n]] == s]; Select[ Range@ 1000, fQ] (* or *) %t A136333 depth = 3; FromDigits@# & /@ FlattenAt[ Table[ Tuples[{1, 3, 7, 9}, n], {n, depth}], {#} & /@ Range[depth]] (* _Robert G. Wilson v_, Jul 02 2014 *) %o A136333 (Haskell) %o A136333 import Data.List (intersect) %o A136333 a136333 n = a136333_list !! (n-1) %o A136333 a136333_list = filter (null . intersect "024568" . show) [1..] %o A136333 -- _Reinhard Zumkeller_, Jul 17 2014 %o A136333 (PARI) isok(m) = my(d=digits(m)); apply(x->gcd(x, 10), d) == vector(#d, k, 1); \\ _Michel Marcus_, Feb 25 2022 %Y A136333 Cf. A007090, A091633 (primes), A245193. %K A136333 nonn,base,easy %O A136333 1,2 %A A136333 _Reinhard Zumkeller_, Mar 26 2008 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE