A227203 Prime(n)^2 mod (prime(n) + prime(n+1)).

4, 1, 1, 13, 1, 19, 1, 25, 9, 1, 9, 43, 1, 49, 9, 9, 1, 9, 73, 1, 9, 85, 9, 109, 103, 1, 109, 1, 115, 49, 133, 9, 1, 25, 1, 9, 9, 169, 9, 9, 1, 25, 1, 199, 1, 241, 253, 229, 1, 235, 9, 1, 25, 9, 9, 9, 1, 9, 283, 1, 25, 49, 313, 1, 319, 49, 9, 25, 1, 355, 9, 379

Offset

1,1

Comments

a(n) = 1 if and only if prime(n) and prime(n+1) are twin prime pair (that is, prime(n+1) = 2 + prime(n)).

Links

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

Mathematica

Table[Mod[Prime[n]^2, Prime[n] + Prime[n + 1]], {n, 200}]
PowerMod[#[[1]], 2, Total[#]]&/@Partition[Prime[Range[80]], 2, 1] (* Harvey P. Dale, Dec 26 2017 *)

Prog

(PARI) for(n=1, 200, print1(prime(n)^2%(prime(n)+prime(n+1))", "))

Crossrefs

Cf. A029707, A167770.

Sequence in context: A051433 A163366 A181145 * A140070 A158815 A101275

Adjacent sequences: A227200 A227201 A227202 * A227204 A227205 A227206

Keyword

nonn

Author

Zak Seidov, Sep 18 2013

Status

approved