%I #31 Sep 08 2022 08:46:17
%S 1,7,11,17,71,77,111,117,171,177,711,717,771,777,1111,1117,1171,1177,
%T 1711,1717,1771,1777,7111,7117,7171,7177,7711,7717,7771,7777,11111,
%U 11117,11171,11177,11711,11717,11771,11777,17111,17117,17171,17177,17711,17717,17771,17777
%N Numbers using only digits 1 and 7.
%C Numbers k such that the product of digits of k is a power of 7.
%C There are no prime terms whose number of digits is divisible by 3: for every d that is a multiple of 3, every d-digit number j consisting of no digits other than 1's and 7's will have a digit sum divisible by 3, so j will also be divisible by 3. - _Mikk Heidemaa_, Mar 27 2021
%H David A. Corneth, <a href="/A276039/b276039.txt">Table of n, a(n) for n = 1..10000</a>
%e 7717 is in the sequence because 7*7*1*7 = 343 = 7^3.
%t Select[Range[20000], IntegerQ[Log[7, Times@@(IntegerDigits[#])]] &] (* or *) Flatten[Table[FromDigits/@Tuples[{1, 7}, n], {n, 6}]]
%o (Magma) [n: n in [1..24000] | Set(Intseq(n)) subset {1,7}];
%o (PARI) is(n) = my(d=digits(n), e=[0, 2, 3, 4, 5, 6, 8, 9]); if(#setintersect(Set(d), Set(e))==0, return(1), return(0)) \\ _Felix Fröhlich_, Aug 19 2016
%o (PARI) a(n) = { my(b = binary(n + 1)); b = b[^1]; b = apply(x -> 6*x + 1, b); fromdigits(b) } \\ _David A. Corneth_, Mar 27 2021
%o (Python)
%o def a(n):
%o b = bin(n+1)[3:]
%o return int("".join(b.replace("1", "7").replace("0", "1")))
%o print([a(n) for n in range(1, 47)]) # _Michael S. Branicky_, Mar 27 2021
%Y Cf. similar sequences listed in A276037.
%Y Cf. A020455, A002275, A002281, A034054.
%K nonn,easy,base
%O 1,2
%A _Vincenzo Librandi_, Aug 19 2016
|