login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A251964 For a prime p, denote by s(p,k) the odd part of the digital sum of p^k. Let k_1 be the smallest k such that s(p,k) is divisible by 5. Sequence lists primes p for which s(p,k_1)=5. 11
2, 5, 7, 11, 19, 23, 37, 41, 61, 71, 73, 101, 109, 113, 127, 131, 163, 179, 181, 211, 229, 241, 251, 271, 307, 311, 313, 383, 389, 401, 421, 433, 449, 479, 521, 523, 541, 557, 569, 571, 587, 601, 613, 631, 659, 677, 751, 811, 827, 839, 857, 929, 947, 971, 977 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Let p be a prime other than 3. If p is not in the sequence, then either s(p,k_1) >= 25 or k_1 does not exist. We conjecture that k_1=k_1(p) exists for every prime p.

LINKS

Table of n, a(n) for n=1..55.

EXAMPLE

For p=7, s(p,1) = 7, s(p,2) = 4+9 = 13, s(p,3) = (3+4+3)/2 = 5. So 7 is a term.

For p=13, s(p,1) = 1, s(p,2) = 1, s(p,3) = 19, s(p,4) = 11, s(p,5) = 25. So 13 is not in the sequence.

MATHEMATICA

s[p_, k_] := Module[{s = Total[IntegerDigits[p^k]]}, s/2^IntegerExponent[s, 2]]; f5[p_] := Module[{k = 1}, While[! Divisible[s[p, k], 5], k++]; k]; ok5Q[p_] := s[p, f5[p]] == 5; Select[Range[1000], PrimeQ[#] && ok5Q[#] &] (* Amiram Eldar, Dec 07 2018 *)

PROG

(PARI) s(p, k) = my(s=sumdigits(p^k)); s >> valuation(s, 2);

f5(p) = my(k=1); while(s(p, k) % 5, k++); k;

isok5(p) = s(p, f5(p)) == 5;

lista5(nn) = forprime(p=2, nn, if (isok5(p), print1(p, ", "))); \\ Michel Marcus, Dec 07 2018

CROSSREFS

Cf. A221858, A225039, A225093.

Sequence in context: A152216 A045350 A179273 * A045351 A051645 A116635

Adjacent sequences:  A251961 A251962 A251963 * A251965 A251966 A251967

KEYWORD

nonn,base

AUTHOR

Vladimir Shevelev, Dec 11 2014

EXTENSIONS

More terms from Peter J. C. Moses, Dec 11 2014

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 17 14:09 EDT 2019. Contains 324185 sequences. (Running on oeis4.)