cp's OEIS Frontend

A067863 Numbers k such that k divides the sum of digits of 7^k.

%I A067863 #17 Sep 22 2024 18:08:56
%S A067863 1,13,67,94,139,220,805
%N A067863 Numbers k such that k divides the sum of digits of 7^k.
%C A067863 There are no other terms < 3000. - _Stefan Steinerberger_, Mar 14 2006
%C A067863 No more terms < 50000. - _David Wasserman_, May 30 2008
%C A067863 From _Jon E. Schoenfield_, May 29 2010: (Start)
%C A067863 No more terms < 100000. It is nearly certain that no terms exist beyond 805.
%C A067863 Let f(k) be the sum of digits of 7^k. Let d be the number of digits, i.e., d=ceiling(log_10(7^k)).
%C A067863 Let s(m) be the sum of m random digits (each drawn independently from a uniform distribution over the integers 0 through 9).
%C A067863 As k increases, the behavior of f(k)/k becomes increasingly similar to that of s(d)/k.
%C A067863 The mean and variance of s(d)/k are 4.5*d/k and 28.5*d/k^2, respectively.
%C A067863 For large values of k, the distribution of s(d)/k approaches a standard normal distribution with mean 4.5*log_10(7) (approximately 3.80294) and variance 28.5*log_10(7)/k.
%C A067863 The probability P(k) that s(d)/k departs from the mean by an amount at least sufficient to reach the nearest higher or lower integer (so that k divides the sum of digits) becomes vanishingly small (e.g., P(50000) < 10^-18, P(100000) < 10^-36, P(150000) < 10^-54), and the same is true of the sum of P(i) for all i >= k (this sum is less than 10^-33 at k=100000). (End)
%e A067863 13 divides the sum of digits of 7^13 (i.e., 9 + 6 + 8 + 8 + 9 + 0 + 1 + 0 + 4 + 0 + 7 = 52), so 13 is in the sequence.
%t A067863 For[n = 1, n < 2000, n++, a := DigitCount[7^n]; If[IntegerQ[Sum[a[[i]]*i, {i, 1, 9}]/n], Print[n]]] (* _Stefan Steinerberger_, Mar 14 2006 *)
%t A067863 Select[Range[1000],Mod[Total[IntegerDigits[7^#]],#]==0&] (* _Harvey P. Dale_, Sep 22 2024 *)
%Y A067863 Cf. A062927, A067862, A067864.
%K A067863 hard,more,nonn,base
%O A067863 1,2
%A A067863 _Shyam Sunder Gupta_ and _Amarnath Murthy_, Feb 16 2002
%E A067863 Edited by _Jon E. Schoenfield_, May 29 2010