Primes in a prescribed arithmetic progression dividing the sequence a^k+b^k
Given positive integers a,b,c and d such that c and d are coprime we show that the primes p=c(mod d)dividing a^k+b^k for some k>=1 have a natural density and explicitly compute this density. We demonstrate our results by considering some claims of Fermat that he made in a 1641 letter to Mersenne.