OFFSET
1,1
MATHEMATICA
s[p_, k_] := Module[{s = Total[IntegerDigits[p^k]]}, s/2^IntegerExponent[s, 2]]; f11[p_] := Module[{k = 1}, While[! Divisible[s[p, k], 11], k++]; k]; ok11Q[p_] := s[p, f11[p]] == 11; Select[Range[1000], PrimeQ[#] && ok11Q[#] &] (* Amiram Eldar, Dec 07 2018 *)
PROG
(PARI) s(p, k) = my(s=sumdigits(p^k)); s >> valuation(s, 2);
f11(p) = my(k=1); while(s(p, k) % 11, k++); k;
isok11(p) = s(p, f11(p)) == 11;
lista11(nn) = forprime(p=2, nn, if (isok11(p), print1(p, ", "))); \\ Michel Marcus, Dec 07 2018
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Vladimir Shevelev and Peter J. C. Moses, Dec 16 2014
STATUS
approved