%I #23 Mar 01 2018 18:36:10
%S 25,26,30,32,47,58,79,81,87,89,102,123,141,144,145,151,164,176,178,
%T 193,201,227,239,242,257,264,282,289,300,306,319,324,329,335,336,338,
%U 348,351,358,365,395,403,437,441,450,460,468,484,489,492,495,517,518,541,542,544,554,555,563,565,570
%N Numbers n such that the sum of digits of 3^n (A004166) is divisible by 7.
%H Chai Wah Wu, <a href="/A292931/b292931.txt">Table of n, a(n) for n = 1..10000</a>
%H Mathematics StackExchange, <a href="https://math.stackexchange.com/questions/2433244">Is sum of digits of 3^1000 divisible by 7?</a>
%e a(3) = 30 is in the sequence because 3^30 = 205891132094649 has sum of digits 63, which is divisible by 7.
%p select(n -> convert(convert(3^n,base,10),`+`) mod 7 = 0, [$1..1000]);
%t Select[Range[600],Divisible[Total[IntegerDigits[3^#]],7]&] (* _Harvey P. Dale_, Mar 01 2018 *)
%o (PARI) isok(n) = !(sumdigits(3^n) % 7); \\ _Michel Marcus_, Sep 27 2017
%o (Python)
%o from __future__ import division
%o A292931_list = [n for n in range(1000) if not sum(int(d) for d in str(3**n)) % 7] # _Chai Wah Wu_, Sep 28 2017
%Y Cf. A004166.
%K nonn,base
%O 1,1
%A _Robert Israel_, Sep 27 2017
|