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-5 of 5 results.

A159008 Positive numbers k such that k^2 == 2 (mod 89).

Original entry on oeis.org

25, 64, 114, 153, 203, 242, 292, 331, 381, 420, 470, 509, 559, 598, 648, 687, 737, 776, 826, 865, 915, 954, 1004, 1043, 1093, 1132, 1182, 1221, 1271, 1310, 1360, 1399, 1449, 1488, 1538, 1577, 1627, 1666, 1716, 1755, 1805, 1844, 1894, 1933, 1983, 2022, 2072
Offset: 1

Views

Author

Vincenzo Librandi, Jun 30 2009

Keywords

Comments

Numbers congruent to {25, 64} mod 89. - Amiram Eldar, Feb 26 2023

Links

Crossrefs

Cf. A159007.

Programs

  • Magma
    I:=[25, 64, 114]; [n le 3 select I[n] else Self(n-1)+Self(n-2)-Self(n-3): n in [1..60]]; // Vincenzo Librandi, Mar 02 2012
  • Mathematica
    LinearRecurrence[{1, 1, -1}, {25, 64, 114}, 50] (* Vincenzo Librandi, Mar 02 2012 *)
Select[Range[2100],PowerMod[#,2,89]==2&] (* Harvey P. Dale, May 09 2019 *)
  • PARI
    for(n=1, 60, print1((89+11*(-1)^(n-1)+178*(n-1))/4", ")); \\ Vincenzo Librandi, Mar 02 2012

Formula

a(n) = a(n-1) + a(n-2) - a(n-3).
G.f.: x*(25 + 39*x + 25*x^2)/((1+x)*(x-1)^2).
a(n) = (89 + 11*(-1)^(n-1) + 178*(n-1))/4.
Sum_{n>=1} (-1)^(n+1)/a(n) = tan(39*Pi/178)*Pi/89. - Amiram Eldar, Feb 26 2023

Extensions

Slightly edited by R. J. Mathar, Jul 26 2009

A204769 a(n) = 151*(n-1) - a(n-1) with n>1, a(1)=46.

Original entry on oeis.org

46, 105, 197, 256, 348, 407, 499, 558, 650, 709, 801, 860, 952, 1011, 1103, 1162, 1254, 1313, 1405, 1464, 1556, 1615, 1707, 1766, 1858, 1917, 2009, 2068, 2160, 2219, 2311, 2370, 2462, 2521, 2613, 2672, 2764, 2823
Offset: 1

Views

Author

Vincenzo Librandi, Mar 08 2012

Keywords

Comments

Positive numbers k such that k^2 == 2 (mod 151), where the prime 151 == -1 (mod 8).
Equivalently, numbers k such that k == 46 or 105 (mod 151). - Bruno Berselli, Mar 08 2012

Links

Crossrefs

Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k-1: A047341, A155450, A164131, A164135, A167533, A167534, A177044, A177046, A204766.
Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k+1: A155449, A158803, A159007, A159008, A176010, A206525.

Programs

  • Magma
    [(-151-33*(-1)^n+302*n)/4: n in [1..60]];
  • Mathematica
    LinearRecurrence[{1,1,-1}, {46,105,197}, 40] (* or *) CoefficientList[Series[x*(46+59*x+46*x^2)/((1+x)*(x-1)^2),{x,0,33}],x] (* or *) a[1] = 46; a[n_] := a[n] = 151*(n-1) - a[n-1]; Table[a[n], {n, 1, 40}]
  • PARI
    a(n)=(-151-33*(-1)^n+302*n)/4 \\ Charles R Greathouse IV, Oct 16 2015

Formula

G.f.: x*(46+59*x+46*x^2)/((1+x)*(x-1)^2).
a(n) = (-151-33*(-1)^n+302*n)/4.
a(n) = a(n-1) +a(n-2) -a(n-3).
Sum_{n>=1} (-1)^(n+1)/a(n) = tan(59*Pi/302)*Pi/151. - Amiram Eldar, Feb 28 2023

A206525 a(n) = 113*(n-1) - a(n-1) with n>1, a(1)=51.

Original entry on oeis.org

51, 62, 164, 175, 277, 288, 390, 401, 503, 514, 616, 627, 729, 740, 842, 853, 955, 966, 1068, 1079, 1181, 1192, 1294, 1305, 1407, 1418, 1520, 1531, 1633, 1644, 1746, 1757, 1859, 1870, 1972, 1983, 2085, 2096, 2198, 2209
Offset: 1

Views

Author

Vincenzo Librandi, Mar 09 2012

Keywords

Comments

Positive numbers k such that k^2 == 2 (mod 113), where the prime 113 == 1 (mod 8).
Equivalently, numbers k such that k == 51 or 62 (mod 113).

Links

Crossrefs

Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k+1: A155449, A158803, A159007, A159008, A176010.
Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k-1: A047341, A155450, A164131, A164135, A167533, A167534, A177044, A177046, A204769.

Programs

  • Magma
    [(-113-91*(-1)^n+226*n)/4: n in [1..60]];
  • Mathematica
    LinearRecurrence[{1, 1, -1}, {51, 62, 164}, 40] (* or *) CoefficientList[Series[x*(51+11*x+51*x^2)/((1+x)*(x-1)^2), {x, 0, 40}], x] (* or *) a[1] = 51; a[n_] := a[n] = 113*(n-1) - a[n-1]; Table[a[n], {n, 1, 40}]

Formula

a(n) = a(n-2) + 113.
G.f.: x*(51+11*x+51*x^2)/((1+x)*(x-1)^2).
a(n) = (-113-91*(-1)^n+226*n)/4.
a(n) = a(n-1)+a(n-2)-a(n-3).
Sum_{n>=1} (-1)^(n+1)/a(n) = tan(11*Pi/226)*Pi/113. - Amiram Eldar, Feb 28 2023

A204766 a(n) = 167*(n-1)-a(n-1) with n>1, a(1)=13.

Original entry on oeis.org

13, 154, 180, 321, 347, 488, 514, 655, 681, 822, 848, 989, 1015, 1156, 1182, 1323, 1349, 1490, 1516, 1657, 1683, 1824, 1850, 1991, 2017, 2158, 2184, 2325, 2351, 2492, 2518, 2659, 2685, 2826, 2852, 2993, 3019, 3160
Offset: 1

Views

Author

Vincenzo Librandi, Mar 09 2012

Keywords

Comments

Positive numbers k such that k^2 == 2 (mod 167), where the prime 167 == -1 (mod 8).
Equivalently, numbers k such that k == 13 or 154 (mod 167).

Links

Crossrefs

Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k-1: A047341, A155450, A164131, A164135, A167533, A167534, A177044, A177046, A204769.
Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k+1: A155449, A158803, A159007, A159008, A176010, A206525, A206526.

Programs

  • Magma
    [(-167+115*(-1)^n+334*n)/4: n in [1..60]];
  • Mathematica
    CoefficientList[Series[x*(13+141*x+13*x^2)/((1+x)*(x-1)^2), {x, 0, 40}], x] (* or *) LinearRecurrence[{1, 1, -1}, {13, 154, 180}, 40]

Formula

G.f.: x*(13+141*x+13*x^2)/((1+x)*(x-1)^2).
a(n) = (-167+115*(-1)^n+334*n)/4.
a(n) = a(n-1)+a(n-2)-a(n-3).
Sum_{n>=1} (-1)^(n+1)/a(n) = cot(13*Pi/167)*Pi/167. - Amiram Eldar, Feb 28 2023

A206526 a(n) = 137*(n-1) - a(n-1) with n>1, a(1)=31.

Original entry on oeis.org

31, 106, 168, 243, 305, 380, 442, 517, 579, 654, 716, 791, 853, 928, 990, 1065, 1127, 1202, 1264, 1339, 1401, 1476, 1538, 1613, 1675, 1750, 1812, 1887, 1949, 2024, 2086, 2161, 2223, 2298, 2360, 2435, 2497, 2572, 2634, 2709, 2771, 2846, 2908, 2983
Offset: 1

Views

Author

Vincenzo Librandi, Mar 09 2012

Keywords

Comments

Positive numbers k such that k^2 == 2 (mod 137), where the prime 137 == 1 (mod 8).
Equivalently, numbers k such that k == 31 or 106 (mod 137).
The subsequence of primes begins: 31, 853, 1613, 1949, 2161. - Jonathan Vos Post, Mar 09 2012

Links

Crossrefs

Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k+1: A155449, A158803, A159007, A159008, A176010, A206525.
Sequences of the type n^2 == 2 (mod p), where p is a prime of the form 8k-1: A047341, A155450, A164131, A164135, A167533, A167534, A177044, A177046, A204769.

Programs

  • Magma
    [(-137+13*(-1)^n+274*n)/4: n in [1..60]];
  • Magma
    [n: n in [1..3000] | n^2 mod 137 eq 2]; // Vincenzo Librandi, Mar 31 2016
  • Mathematica
    LinearRecurrence[{1, 1, -1}, {31, 106, 168}, 40] (* or *) CoefficientList[Series[x*(31+75*x+31*x^2)/((1+x)*(x-1)^2), {x, 0, 50}], x] (* or *) a[1] = 31; a[n_] := a[n] = 137*(n-1) - a[n-1]; Table[a[n], {n, 1, 40}]

Formula

a(n) = a(n-2) + 137.
G.f.: x*(31+75*x+31*x^2)/((1+x)*(x-1)^2).
a(n) = (-137+13*(-1)^n+274*n)/4.
a(n) = a(n-1)+a(n-2)-a(n-3).
Sum_{n>=1} (-1)^(n+1)/a(n) = cot(31*Pi/137)*Pi/137. - Amiram Eldar, Feb 28 2023
Showing 1-5 of 5 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.