cp's OEIS Frontend

A225017 Odd part of digit sum of 5^n divided by maximal possible power of 5.

%I A225017 #19 Oct 27 2023 22:11:46
%S A225017 1,1,7,1,13,11,19,23,1,13,1,19,7,23,17,11,29,7,1,59,61,31,67,37,41,77,
%T A225017 79,89,17,83,91,13,53,89,103,23,109,13,31,67,13,137,29,149,151,29,7,1,
%U A225017 29,79,151,19,13,119,127,167,49,43,211,191,199,97,187,17,83
%N A225017 Odd part of digit sum of 5^n divided by maximal possible power of 5.
%C A225017 Does the sequence contain every prime greater than 5?
%H A225017 Alois P. Heinz, <a href="/A225017/b225017.txt">Table of n, a(n) for n = 0..10000</a> (terms n = 0..999 from Peter J. C. Moses)
%F A225017 a(n) = A132740(A055566(n)). - _Michel Marcus_, Dec 10 2018
%p A225017 a:= proc(n) local m, r; m, r:= 0, 5^n;
%p A225017       while r>0 do m:= m+irem(r, 10, 'r') od;
%p A225017       while irem(m, 2, 'r')=0 do m:=r od;
%p A225017       while irem(m, 5, 'r')=0 do m:=r od; m
%p A225017     end:
%p A225017 seq(a(n), n=0..80);  # _Alois P. Heinz_, Apr 24 2013
%t A225017 Map[#/(2^IntegerExponent[#,2] 5^IntegerExponent[#,5])&[Total[ IntegerDigits[5^#]]]&,Range[0,99]] (* _Peter J. C. Moses_, Apr 24 2013 *)
%o A225017 (PARI) a(n) = my(x = sumdigits(5^n)); x/5^valuation(x, 5) >> valuation(x, 2); \\ _Michel Marcus_, Dec 10 2018
%Y A225017 Cf. A055566, A132740.
%K A225017 nonn,base
%O A225017 0,3
%A A225017 _Vladimir Shevelev_, Apr 24 2013