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.

Previous Showing 31-34 of 34 results.

A185438 a(n) = 8*n^2 - 2*n + 1.

Original entry on oeis.org

1, 7, 29, 67, 121, 191, 277, 379, 497, 631, 781, 947, 1129, 1327, 1541, 1771, 2017, 2279, 2557, 2851, 3161, 3487, 3829, 4187, 4561, 4951, 5357, 5779, 6217, 6671, 7141, 7627, 8129, 8647, 9181, 9731, 10297, 10879, 11477, 12091, 12721, 13367, 14029, 14707, 15401, 16111, 16837, 17579
Offset: 0

Views

Author

Paul Curtz, Feb 03 2011

Keywords

Comments

Odd numbers (A005408) written clockwise as a square spiral:
.
41--43--45--47--49--51
| |
39 13--15--17--19 53
| | | |
37 11 1---3 21 55
| | | | |
35 9---7---5 23 57
| | |
33--31--29--27--25 59
|
71--69--67--65--63--61
.
Walking in straight lines away from the center:
1, 17, 49, ... = A069129(n+1) = 1 - 8*n + 8*n^2,
1, 3, 21, ... = A033567(n) = 1 - 6*n + 8*n^2,
1, 15, 45, ... = A014634(n) = 1 + 6*n + 8*n^2,
1, 5, 25, ... = A080856(n) = 1 - 4*n + 8*n^2,
1, 13, 41, ... = A102083(n) = 1 + 4*n + 8*n^2,
1, 7, 29, ... = a(n) = 1 - 2*n + 8*n^2,
1, 11, 37, ... = A188135(n) = 1 + 2*n + 8*n^2,
1, 9, 33, ... = A081585(n) = 1 + 8*n^2,
5, 29, 69, ... = A108928(n+1) = -3 + 8*n^2,
7, 31, 71, ... = A157914(n+1) = -1 + 8*n^2,
9, 35, 77, ... = A033566(n+1) = -1 + 2*n + 8*n^2.
All are quadrisections of sequences in A181407(n) (example: A014634(n) and A033567(n) in A064038(n+1)) or of this family (?): a(n) is a quadrisection of f(n) = 1,1,1,1,2,7,11,8,11,29,37,23,28,67,79,46,... f(n) is just before A064038(n+1) (fifth vertical) in A181407(n). The companion to a(n) is A188135(n), another quadrisection of f(n). Two last quadrisections of f(n) are A054552(n) and A033951(n).
For n >= 1, bisection of A193867. - Omar E. Pol, Aug 16 2011
Also the sequence may be obtained by starting with the segment (1, 7) followed by the line from 7 in the direction 7, 29, ... in the square spiral whose vertices are the generalized hexagonal numbers (A000217). - Omar E. Pol, Aug 01 2016

Crossrefs

Programs

Formula

a(n) = a(n-1) + 16*n - 10 (n > 0).
a(n) = 2*a(n-1) - a(n-2) + 16 (n > 1).
a(n) = 3*(n-1) - 3*a(n-2) + a(n-3) (n > 2).
G.f.: (-1 - 4*x - 11*x^2)/(x-1)^3. - R. J. Mathar, Feb 03 2011
a(n) = A014635(n) + 1. - Bruno Berselli, Apr 09 2011
E.g.f.: exp(x)*(1 + 6*x + 8*x^2). - Elmo R. Oliveira, Nov 17 2024

A195241 Expansion of (1-x+19*x^3-3*x^4)/(1-x)^3.

Original entry on oeis.org

1, 2, 3, 23, 59, 111, 179, 263, 363, 479, 611, 759, 923, 1103, 1299, 1511, 1739, 1983, 2243, 2519, 2811, 3119, 3443, 3783, 4139, 4511, 4899, 5303, 5723, 6159, 6611, 7079, 7563, 8063, 8579, 9111, 9659, 10223, 10803, 11399, 12011, 12639, 13283, 13943
Offset: 0

Views

Author

Bruno Berselli, Sep 13 2011 - based on remarks and sequences by Omar E. Pol

Keywords

Comments

Sequence found by reading the line 1, 2, 3, 23,.. in the square spiral whose vertices are the triangular numbers (A000217) - see Pol's comments in other sequences visible in this numerical spiral.
This is a subsequence of A110326 (without signs) and A047838 (apart from the second term, 2).

Crossrefs

Programs

  • Magma
    m:=44; R:=PowerSeriesRing(Integers(), m); Coefficients(R!((1-x+19*x^3-3*x^4)/(1-x)^3));
    
  • Mathematica
    CoefficientList[Series[(1 - x + 19 x^3 - 3 x^4)/(1 - x)^3, {x, 0, 50}], x] (* Vincenzo Librandi, Mar 26 2013 *)
    LinearRecurrence[{3,-3,1},{1,2,3,23,59},50] (* Harvey P. Dale, Dec 04 2022 *)
  • Maxima
    makelist(coeff(taylor((1-x+19*x^3-3*x^4)/(1-x)^3, x, 0, n), x, n), n, 0, 43);
  • PARI
    Vec((1-x+19*x^3-3*x^4)/(1-x)^3+O(x^44))
    

Formula

G.f.: (1-x+19*x^3-3*x^4)/(1-x)^3.
a(n) = 8*n^2-20*n+11 for n>1; a(0)=1, a(1)=2.

A264823 Centered 16-gonal (or hexadecagonal) primes.

Original entry on oeis.org

17, 97, 241, 337, 449, 577, 881, 1249, 3041, 3361, 3697, 4049, 4801, 6961, 7937, 9521, 10657, 13121, 14449, 15137, 16561, 18049, 23761, 25537, 26449, 31249, 32257, 37537, 39761, 40897, 46817, 48049, 49297, 54449, 58481, 64081, 65521, 66977, 68449, 77617, 79201
Offset: 1

Views

Author

Ilya Gutkovskiy, Nov 26 2015

Keywords

Comments

Primes of the form 8*k^2 - 8*k + 1.

Crossrefs

Programs

  • Mathematica
    Select[Table[8k^2-8k+1,{k,0,100}],PrimeQ] (* Harvey P. Dale, Mar 08 2023 *)
  • PARI
    for(n=1, 1e3, if(isprime(k=8*n^2-8*n+1), print1(k,", "))) \\ Altug Alkan, Nov 26 2015

A368340 Take the solution to Pellian equation x^2 - 8*n*y^2 = 1 with smallest positive y and x >= 0; sequence gives a(n) = x, or 1 if n is twice a positive square. A368339 gives values of y.

Original entry on oeis.org

3, 1, 5, 17, 19, 7, 15, 1, 17, 9, 197, 49, 51, 127, 11, 577, 35, 1, 37, 721, 13, 199, 24335, 97, 99, 649, 485, 15, 19603, 31, 63, 1, 65, 33, 251, 17, 3699, 57799, 53, 161, 163, 55, 10405, 77617, 19, 1151, 2143295, 4801, 99, 1, 101, 5201, 32080051, 1351, 21, 127
Offset: 1

Views

Author

Arkadiusz Wesolowski, Dec 21 2023

Keywords

Examples

			For n = 1, 2, 3, 4, 5 solutions are (x,y) = (3, 1), (1, 0), (5, 1), (17, 3), (19, 3).
		

Crossrefs

Programs

  • PARI
    pellsolve(n)={if(issquare(n/2), return(1), q=bnfinit('x^2-8*n, 1); i=-1; until(y&&x==floor(x)&&y==floor(y)&&x^2-8*n*y^2==1, f=lift(q.fu[1]^i); x=abs(polcoeff(f, 0)); y=abs(polcoeff(f, 1)); i++); return(x))};

Formula

a(n) = A002350(8*n).
a(n) = sqrt(8*n*A368339(n)^2 + 1).
a(A000217(n)) = 2*n + 1, n >= 1.
Previous Showing 31-34 of 34 results.