cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A201998 Positive numbers n such that n^2 + n + 41 is composite and there are no positive integers c such that n = c*x^2 + (c + 1)*x + c*41 for an integer x.

Original entry on oeis.org

244, 249, 251, 266, 270, 295, 301, 336, 344, 389, 399, 407, 416, 418, 445, 449, 454, 466, 489, 494, 496, 500, 506, 527, 531, 545, 547, 563, 570, 571, 582, 583, 585, 611, 612, 620, 622, 624, 628, 630, 636, 652, 661, 662, 663, 679, 693, 699
Offset: 1

Views

Author

Matt C. Anderson, Dec 07 2011

Keywords

Comments

The composition of functions k(x) factors. k(x) = (x^2 + x + 41)*(c^2*x^2 + (c^2 + 2*c)*x + c^2*41 + c + 1). So k(x) is the product of two integers greater than one and thus composite.

References

  • John Stillwell, Elements of Number Theory, Springer, 2003, page 3.

Links

Crossrefs

Cf. A007634 (n^2 + n + 41 is composite).
Cf. A235381 (similar to this sequence).

Programs

  • Maple
    maxn:=1000:
A:={}:
for n from 1 to maxn do
g:=n^2+n+41:
if isprime(g)=false then
A:=A union {n}:
end if:
end do:
# The set A contains values n such that n^2+n+41 is composite and n < maxn.
c:=1:
x:=-1:
p:=41:
q:=c*x^2-(c+1)*x+c*p:
A2:=A:
while q < maxn do
while q < maxn do
A2:=A2 minus {q}:
A2:=A2 minus {c*x^2+(c+1)*x+c*p}:
x:=x+1:
q:=c*x^2-(c+1)*x+c*p:
end do:
c:=c+1:
x:=-1:
q:=c*x^2-(c+1)*x+c*p:
end do:
A2;
  • Mathematica
    Reap[For[n=1, n<700, n++, If[!PrimeQ[n^2+n+41], If[Reduce[c>0 && n == c*x^2+(c+1)*x+41*c , {c, x}, Integers] === False, Sow[n]]]]][[2, 1]] (* Jean-François Alcover, Apr 30 2014 *)

A241529 Positive numbers k such that k^2 + k + 41 is composite and there are no positive integers a,c,d such that k = c*a*z^2 + ((((d+2)*(1/3))*c-2)*a/d+1)*z + ((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2 - (((d-1)*(1/3))*c+1)/d)/c for an integer z.

Original entry on oeis.org

2887, 2969, 3056, 3220, 3365, 3464, 3565, 3611, 3719, 3746, 3814, 3836, 3874, 3879, 3955, 4142, 4147, 4211, 4277, 4371, 4403, 4483, 4564, 4572, 4661, 4730, 4813, 4881, 4888, 4902, 4906, 4965, 4982, 5132, 5175, 5208, 5410, 5431, 5509, 5527, 5564, 5624, 5669
Offset: 1

Views

Author

Matt C. Anderson, Apr 27 2014

Keywords

Comments

This sequence has a restriction involving 4 variables. More composite cases are described with a better restrictive expression. The expression for k(a,c,d,z) will force k^2 + k + 41 to be either a fraction or a composite number.
The condition on k(a,c,d,z) was determined by quadratic curve fitting. It has been automated with the Maple Interactive() command. The ultimate motivation is to try to find a closed-from expression that generates all the composite cases of k^2 + k + 41 for integer k.
What is the smallest value of n where this sequence's a(n) < 2n? (For A194634, this value is 2358.) - J. Lowell, Feb 25 2019

References

  • John Stillwell, Elements of Number Theory, Springer, 2003, page 3.

Links

Crossrefs

Cf. A007634, A055390, A201998, and with division, A235381.

Programs

  • Maple
    # Euler considered the prime values for n^2 + n + 41;
# This is a 76 second calculation on a 2.93 GHz machine
h := n^2+n+41;
y := c*a*z^2+((((d+2)*(1/3))*c-2)*a/d+1)*z+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c;
y2 := subs(n = y, h);
y3 := factor(y2);
# note that y is an expression in 4 variables.
# After a composition of functions, an algebraic factorization
# can be observed in y3.  As long as y3 is an integer, it will
# be composite.  This is because y3 factors and both factors
# are integers bigger than one.
maxn := 6000;
A := {}:
for n to maxn do
g := n^2+n+41:
if isprime(g) = false then A := `union`(A, {n}) end if :
end do:
# now the A set contains composite values of the form
# n^2 + n + 41 less than maxn.
c := 1: a := 1: d := 1: z := -1: p := 41:
q := c*a*z^2+((((d+2)*(1/3))*c-2)*a/d+1)*z+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c:
A2 := A:
while q < maxn do
while `and`(q < maxn, d < 100) do
while q < maxn do while
q < maxn do
A2 := `minus`(A2, {q});
A2 := `minus`(A2, {c*a*z^2+((((d+2)*(1/3))*c-2)*a/d+1)*z+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c});
z := z+1;
A2 := `minus`(A2, {c*a*z^2-((((d+2)*(1/3))*c-2)*a/d+1)*(1*z)+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c}); q := c*a*z^2+((((d+2)*(1/3))*c-2)*a/d+1)*z+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c
end do;
a := a+1; z := -1;
q := c*a*z^2+((((d+2)*(1/3))*c-2)*a/d+1)*z+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c :
end do;
d := d+1: a := 1:
q := c*a*z^2+((((d+2)*(1/3))*c-2)*a/d+1)*z+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c :
end do:
c := c+1: d := 1:
q := c*a*z^2+((((d+2)*(1/3))*c-2)*a/d+1)*z+((((367*d^2+d+1)*(1/9))*c^2-((d+2)*(1/3))*c+1)*a/d^2-(((d-1)*(1/3))*c+1)/d)/c :
end do:
A2;
# Matt C. Anderson, May 13 2014
Showing 1-2 of 2 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.