OFFSET
1,1
COMMENTS
k = 1 is the "trivial" case: sod(prime(n)) = prime(m)
n = 2, prime(2) = 3: 3^k is for k > 1 a multiple of 3^2.
REFERENCES
M. Fujiwara, Y. Ogawa: Introduction to truly beautiful Mathematics, Chikuma Shobo, Tokyo 2005.
Theo Kempermann, Zahlentheoretische Kostproben, Harri Deutsch, 2. aktualisierte Auflage 2005.
Hans Schubart: Einfuehrung in die klassische und moderne Zahlentheorie, Vieweg, Braunschweig 1974.
EXAMPLE
sod(2^5)=5, sod(5^2)=7, sod(7^2)=13, sod(11^9)=53, sod(13^3)=19, sod(17^2)=19, sod(19^5)=37, sod(23^3)=17, sod(29^2)=13, sod(31^7)=31, sod(37^2)=19, sod(41^4)=31, sod(43^5)=31, sod(47^2)=13, sod(53^2)=19, sod(59^5)=47, sod(61^2)=13, sod(67^3)=19, sod(71^6)=37, sod(73^2)=19, sod(79^2)=13, sod(83^2)=31, sod(89^2)=19, sod(97^4)=43, sod(101^8)=67, sod(103^4)=31, sod(107^2)=19, sod(109^2)=19, sod(113^4)=31, sod(127^2)=19, sod(131^8)=61, sod(137^2)=31, sod(139^3)=37, sod(149^2)=7, sod(151^2)=13, sod(157^4)=31, sod(163^4)=37, sod(167^6)=73, sod(173^2)=31, sod(179^4)=37, sod(181^2)=19, sod(191^10)=97, sod(193^3)=37, sod(197^4)=37, sod(199^2)=19, sod(211^3)=37, sod(223^2)=31, sod(227^4)=43, sod(229^3)=37.
PROG
(Magma) S:=[ 5, 0 ]; for n in [3..103] do j:=2; while not IsPrime(&+Intseq(NthPrime(n)^j)) do j+:=1; end while; Append(~S, j); end for; S; // Klaus Brockhaus, Jan 29 2010
CROSSREFS
KEYWORD
easy,nonn,base
AUTHOR
Ulrich Krug (leuchtfeuer37(AT)gmx.de), Jan 23 2010
EXTENSIONS
More terms from Klaus Brockhaus, Jan 29 2010
Edited by Charles R Greathouse IV, Aug 02 2010
STATUS
approved