|
|
A226927
|
|
Values of n such that L(7) and N(7) are both prime, where L(k) = (n^2+n+1)*2^(2*k) + (2*n+1)*2^k + 1, N(k) = (n^2+n+1)*2^k + n.
|
|
1
|
|
|
9, -15, -95, 109, -243, -297, 297, 457, 459, -477, 583, 723, -761, -771, -983, 1045, 1047, 1077, 1287, 1305, -1373, 1389, -1481, -2063, -2073, 2085, 2173, -2223, 2263, -2367, 2503, -2591, -2615, -2643, 2643, 2755, -2955, -2957, 2995, 3105, -3153, 3153, 3237, 3243, -3267, -3491, 3667, 3699, 3847, 3919, -4023, -4217, 4255, -4275, 4377, 4483, 4549, 4743, 4827, 4855, -4973
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
Computed with PARI using commands similar to those used to compute A226921.
|
|
LINKS
|
|
|
MATHEMATICA
|
k = 7; (* adjust for related sequences *) fL[n_] := (n^2 + n + 1)*2^(2*k) + (2*n + 1)*2^k + 1; fN[n_] := (n^2 + n + 1)*2^k + n; nn = 5000; A = {}; For[n = -nn, n <= nn, n++, If[PrimeQ[fL[n]] && PrimeQ[fN[n]], AppendTo[A, n]]]; cmpfunc[x_, y_] := If[x == y, Return[True], ax = Abs[x]; ay = Abs[y]; If[ax == ay, Return[x < y], Return[ ax < ay]]]; Sort[A, cmpfunc] (* Jean-François Alcover, Jul 17 2013, translated and adapted from Joerg Arndt's Pari program in A226921 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|