login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A254616
Let i < k and prime(k) == x(i) mod prime(i). This sequence gives the indices k of the primes prime(k) such that Sum(i=1..k-1, x(i)) is prime.
0
3, 6, 9, 30, 33, 36, 38, 54, 66, 89, 92, 93, 98, 121, 125, 128, 154, 155, 167, 189, 198, 231, 291, 296, 300, 310, 323, 334, 345, 357, 373, 407, 436, 437, 444, 456, 481, 483, 492, 493, 515, 516, 541, 554, 581, 628, 649, 681, 713, 714, 749, 750, 809, 818, 826, 834, 864, 868, 872, 881, 888
OFFSET
1,1
EXAMPLE
prime(3) is 5, 5 == 1 mod 2 == 2 mod 3. 1 + 2 = 3, which is prime. Therefore, the index of the prime(3) is in the sequence.
MATHEMATICA
lst1={}; lst2={}; i=1; primes=Prime[Range[1000]];
Do[
While[i<k, AppendTo[lst1, Mod[primes[[k]], primes[[i]]]]; i++];
If[PrimeQ[Total[lst1]], AppendTo[lst2, k]];
i=1; lst1={},
{k, 3, Length[primes]}
];
lst2
CROSSREFS
Cf. A000040.
Sequence in context: A349597 A178467 A193802 * A195205 A045638 A305322
KEYWORD
easy,nonn
AUTHOR
Ivan N. Ianakiev, Feb 03 2015
STATUS
approved