A071557 Numbers k such that A065876(k) = k^2-k+1.

0, 2, 3, 4, 5, 6, 7, 9, 10, 11, 14, 15, 16, 19, 20, 24, 25, 26, 29, 35, 36, 39, 40, 41, 45, 49, 51, 54, 56, 59, 61, 65, 66, 69, 71, 74, 79, 84, 85, 90, 94, 95, 101, 110, 116, 120, 121, 124, 126, 130, 131, 134, 139, 141, 145, 146, 150, 156, 159, 160, 165, 169, 170, 171

Offset

1,2

Links

Amiram Eldar, Table of n, a(n) for n = 1..1000

Mathematica

q[n_] := Module[{m = n+1}, While[!Divisible[m^2 + 1, n^2 + 1], m++]; m == n^2 - n + 1]; Select[Range[200], q] (* Amiram Eldar, May 05 2022 *)

Prog

(PARI) for(n=0, 300, q=n+1; while((q^2+1)%(n^2+1)>0, q++); if(q==n^2-n+1, print1(n, ", ")))

Crossrefs

Cf. A065876.

Sequence in context: A174328 A272570 A123101 * A331802 A271108 A179401

Adjacent sequences: A071554 A071555 A071556 * A071558 A071559 A071560

Keyword

easy,nonn

Author

Benoit Cloitre, May 30 2002

Extensions

a(1) = 0 added by Amiram Eldar, May 05 2022

Status

approved