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 21-30 of 33 results. Next

A191741 Dispersion of A047217, (numbers >1 and congruent to 0 or 1 or 2 mod 5), by antidiagonals.

Original entry on oeis.org

1, 2, 3, 5, 6, 4, 10, 11, 7, 8, 17, 20, 12, 15, 9, 30, 35, 21, 26, 16, 13, 51, 60, 36, 45, 27, 22, 14, 86, 101, 61, 76, 46, 37, 25, 18, 145, 170, 102, 127, 77, 62, 42, 31, 19, 242, 285, 171, 212, 130, 105, 71, 52, 32, 23, 405, 476, 286, 355, 217, 176, 120
Offset: 1

Views

Author

Clark Kimberling, Jun 14 2011

Keywords

Comments

For a background discussion of dispersions and their fractal sequences, see A191426. For dispersions of congruence sequences mod 3, mod 4, or mod 5, see A191655, A191663, A191667, A191702.
...
Suppose that {2,3,4,5,6} is partitioned as {x1, x2} and {x3,x4,x5}. Let S be the increasing sequence of numbers >1 and congruent to x1 or x2 mod 5, and let T be the increasing sequence of numbers >1 and congruent to x3 or x4 or x5 mod 5. There are 10 sequences in S, each matched by a (nearly) complementary sequence in T. Each of the 20 sequences generates a dispersion, as listed here:
...
A191722=dispersion of A008851 (0, 1 mod 5 and >1)
A191723=dispersion of A047215 (0, 2 mod 5 and >1)
A191724=dispersion of A047218 (0, 3 mod 5 and >1)
A191725=dispersion of A047208 (0, 4 mod 5 and >1)
A191726=dispersion of A047216 (1, 2 mod 5 and >1)
A191727=dispersion of A047219 (1, 3 mod 5 and >1)
A191728=dispersion of A047209 (1, 4 mod 5 and >1)
A191729=dispersion of A047221 (2, 3 mod 5 and >1)
A191730=dispersion of A047211 (2, 4 mod 5 and >1)
A191731=dispersion of A047204 (3, 4 mod 5 and >1)
...
A191732=dispersion of A047202 (2,3,4 mod 5 and >1)
A191733=dispersion of A047206 (1,3,4 mod 5 and >1)
A191734=dispersion of A032793 (1,2,4 mod 5 and >1)
A191735=dispersion of A047223 (1,2,3 mod 5 and >1)
A191736=dispersion of A047205 (0,3,4 mod 5 and >1)
A191737=dispersion of A047212 (0,2,4 mod 5 and >1)
A191738=dispersion of A047222 (0,2,3 mod 5 and >1)
A191739=dispersion of A008854 (0,1,4 mod 5 and >1)
A191740=dispersion of A047220 (0,1,3 mod 5 and >1)
A191741=dispersion of A047217 (0,1,2 mod 5 and >1)
...
For further information about these 20 dispersions, see A191722.
...
Regarding the dispersions A191722-A191741, there are general formulas for sequences of the type "(a or b mod m)" and "(a or b or c mod m)" used in the relevant Mathematica programs.

Examples

			Northwest corner:
1....2....5....10...17
3....6....11...20...35
4....7....12...21...36
8....15...26...45...76
9....16...27...46...77
13...22...37...62...105
		

Crossrefs

Programs

  • Mathematica
    (* Program generates the dispersion array t of the increasing sequence f[n] *)
    r = 40; r1 = 12;  c = 40; c1 = 12;
    a=2; b=5; c2=6; m[n_]:=If[Mod[n,3]==0,1,0];
    f[n_]:=a*m[n+2]+b*m[n+1]+c2*m[n]+5*Floor[(n-1)/3]
    Table[f[n], {n, 1, 30}]  (* A047217 *)
    mex[list_] := NestWhile[#1 + 1 &, 1, Union[list][[#1]] <= #1 &, 1, Length[Union[list]]]
    rows = {NestList[f, 1, c]};
    Do[rows = Append[rows, NestList[f, mex[Flatten[rows]], r]], {r}];
    t[i_, j_] := rows[[i, j]];
    TableForm[Table[t[i, j], {i, 1, 10}, {j, 1, 10}]] (* A191741 *)
    Flatten[Table[t[k, n - k + 1], {n, 1, c1}, {k, 1, n}]] (* A191741  *)

A134153 a(n) = 15*n^2 + 9*n + 1.

Original entry on oeis.org

1, 25, 79, 163, 277, 421, 595, 799, 1033, 1297, 1591, 1915, 2269, 2653, 3067, 3511, 3985, 4489, 5023, 5587, 6181, 6805, 7459, 8143, 8857, 9601, 10375, 11179, 12013, 12877, 13771, 14695, 15649, 16633, 17647, 18691, 19765, 20869, 22003, 23167
Offset: 0

Views

Author

Artur Jasinski, Oct 10 2007

Keywords

Comments

A119617 is union of A134153 and A134154. A000538(n) is divisible by A000330(n) if and only n is congruent to {1, 3} mod 5 (see A047219). A134154(n) is case when n is congruent to 1 mod 5 see cases 2)

Crossrefs

Programs

  • Mathematica
    Table[1 + 9 n + 15 n^2, {n, 0, 50}]
    Table[Sum[k^4, {k, 1, 5m + 1}]/Sum[k^2, {k, 1, 5m + 1}], {m, 0, 10}]
  • PARI
    a(n)=15*n^2+9*n+1 \\ Charles R Greathouse IV, Jun 17 2017

Formula

a(n) = 15*n^2 + 9*n + 1.
a(n) = (3*(5*n + 1)^2 + 3*(5*n + 1) - 1)/5.
a(n) = (Sum_{k=1..5*n+1} k^4) / (Sum_{k=1..5*n+1} k^2).
G.f.: -(1+22*x+7*x^2)/(-1+x)^3. - R. J. Mathar, Nov 14 2007

Extensions

Offset corrected and some punctuation added by R. J. Mathar, Jul 09 2009

A134154 a(n) = 15*n^2 - 9*n + 1.

Original entry on oeis.org

1, 7, 43, 109, 205, 331, 487, 673, 889, 1135, 1411, 1717, 2053, 2419, 2815, 3241, 3697, 4183, 4699, 5245, 5821, 6427, 7063, 7729, 8425, 9151, 9907, 10693, 11509, 12355, 13231, 14137, 15073, 16039, 17035, 18061, 19117, 20203, 21319, 22465, 23641
Offset: 0

Views

Author

Artur Jasinski, Oct 10 2007

Keywords

Comments

A119617 is union of A134153 and A134154 A000538(n) is divisible by A000330(n) if and only n is congruent to {1, 3} mod 5 (see A047219) A134154(n) is case when n is congruent to 3 mod 5 see cases 2)

Crossrefs

Programs

  • Mathematica
    Table[1 - 9 n + 15 n^2, {n, 0, 50}]
    Table[Sum[k^4, {k, 1, 5m + 3}]/Sum[k^2, {k, 1, 5m + 3}], {m, 0, 30}]
  • PARI
    a(n)=15*n^2-9*n+1 \\ Charles R Greathouse IV, Jun 17 2017

Formula

a(n) = 15*n^2 - 9*n + 1.
a(n+1) = (3*(5*n + 3)^2 + 3*(5*n + 3) - 1)/5.
a(n+1) = (Sum_{k=1..5*n+3} k^4) / (Sum_{k=1..5*n+3} k^2).
G.f.: -(1+4*x+25*x^2)/(-1+x)^3. - R. J. Mathar, Nov 14 2007

A012132 Numbers z such that x*(x+1) + y*(y+1) = z*(z+1) is solvable in positive integers x,y.

Original entry on oeis.org

3, 6, 8, 10, 11, 13, 15, 16, 18, 20, 21, 23, 26, 27, 28, 31, 33, 36, 37, 38, 40, 41, 43, 44, 45, 46, 48, 49, 51, 52, 53, 54, 55, 56, 57, 58, 59, 61, 62, 63, 64, 66, 67, 68, 71, 73, 74, 75, 76, 77, 78, 80, 81, 83, 86, 88, 89, 91, 92, 93
Offset: 1

Views

Author

Sander van Rijnswou (sander(AT)win.tue.nl)

Keywords

Comments

Theorem (Sierpinski, 1963): n is a term iff n^2+(n+1)^2 is a composite number. - N. J. A. Sloane, Feb 29 2020
For n > 1, A047219 is a subset of this sequence. This is because n^2 + (n+1)^2 is divisible by 5 if n is (1 or 3) mod 5 (also see A027861). - Dmitry Kamenetsky, Sep 02 2008
From Hermann Stamm-Wilbrandt, Sep 10 2014: (Start)
For n > 0, A212160 is a subset of this sequence (n^2 + (n+1)^2 is divisible by 13 if n == (2 or 10) (mod 13)).
For n >= 0, A212161 is a subset of this sequence (n^2 + (n+1)^2 is divisible by 17 if n == (6 or 10) (mod 17)).
The above are for divisibility by 5, 13, 17; notation (1,3,5), (2,10,13), (6,10,17). Divisibility by p for a and p-a-1; notation (a,p-a-1,p). These are the next tuples: (8,20,29), (15,21,37), (4,36,41), (11,41,53), ... . The corresponding sequences are a subset of this sequence (8,20,37,49,66,78,... for (8,20,29)). These sequences have no entries in the OEIS yet. For any prime of the form 4*k+1 there is exactly one of these tuples/sequences.
For n > 1, A000217 (triangular numbers) is a subset of this sequence (3,6,10,15,...); z=A000217(n), y=z-1, x=n.
For n > 0, A001652 is a subset of this sequence; z=A001652(n), x=y=A053141(n).
For n > 1, A001108(=A115598) is a subset of this sequence; z=A001108(n), x=A076708(n), y=x+1.
For n > 0, A124124(2*n+1)(=A098790(2*n)) is a subset of this sequence (6,37,218,...); z=A124124(2*n+1), x=a(n)-1, y=a(n)+1, a(m) = 6*a(m-1) - a(m-2) + 2, a(0)=0, a(1)=4.
(End)

References

  • Aviezri S. Fraenkel, Diophantine equations involving generalized triangular and tetrahedral numbers, pp. 99-114 of A. O. L. Atkin and B. J. Birch, editors, Computers in Number Theory. Academic Press, NY, 1971.

Crossrefs

Complement of A027861. - Michael Somos, Jun 08 2000

Programs

Extensions

More terms and references from Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Feb 09 2000

A013642 Numbers k such that the continued fraction for sqrt(k) has period 2.

Original entry on oeis.org

3, 6, 8, 11, 12, 15, 18, 20, 24, 27, 30, 35, 38, 39, 40, 42, 48, 51, 56, 63, 66, 68, 72, 80, 83, 84, 87, 90, 99, 102, 104, 105, 110, 120, 123, 132, 143, 146, 147, 148, 150, 152, 156, 168, 171, 182, 195, 198, 200, 203, 210, 224, 227, 228, 230, 231, 235, 240, 255, 258, 260, 264
Offset: 1

Views

Author

Keywords

Comments

This sequence is identical to the sequence of numbers of the form k = a^2 + b, where a and b are positive integers and b is a factor of 2a greater than 1, in which case the continued fraction expansion of sqrt(k) is [a; [2a/b, 2a]]. - David Terr, Jun 11 2004

References

  • Kenneth H. Rosen, Elementary Number Theory and Its Applications, Addison-Wesley, 1984, page 426 (but beware of errors!).

Crossrefs

Programs

  • Mathematica
    cf2Q[n_]:=Module[{s=Sqrt[n]},If[IntegerQ[s],1,Length[ ContinuedFraction[ s][[2]]]]==2]; Select[Range[300],cf2Q] (* Harvey P. Dale, Jun 21 2017 *)

A301563 Expansion of Product_{k>=0} (1 + x^(5*k+1))*(1 + x^(5*k+3)).

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 1, 2, 1, 2, 2, 1, 3, 2, 2, 4, 3, 4, 4, 4, 6, 4, 6, 7, 5, 9, 8, 8, 11, 9, 12, 12, 12, 16, 13, 17, 19, 17, 23, 21, 24, 27, 24, 32, 30, 32, 40, 35, 43, 45, 44, 53, 50, 59, 62, 61, 75, 70, 78, 87, 83, 99, 97, 105, 118, 112, 133, 134, 138, 159, 153
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 23 2018

Keywords

Comments

Number of partitions of n into distinct parts congruent to 1 or 3 mod 5.

Examples

			a(14) = 3 because we have [13, 1], [11, 3] and [8, 6].
		

Crossrefs

Programs

  • Mathematica
    nmax = 72; CoefficientList[Series[Product[(1 + x^(5 k + 1)) (1 + x^(5 k + 3)), {k, 0, nmax}], {x, 0, nmax}], x]
    nmax = 72; CoefficientList[Series[QPochhammer[-x, x^5] QPochhammer[-x^3, x^5], {x, 0, nmax}], x]
    nmax = 72; CoefficientList[Series[Product[(1 + Boole[MemberQ[{1, 3}, Mod[k, 5]]] x^k), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

G.f.: Product_{k>=1} (1 + x^A047219(k)).
a(n) ~ exp(Pi*sqrt(2*n/15)) / (2^(21/20) * 15^(1/4) * n^(3/4)). - Vaclav Kotesovec, Mar 24 2018

A212160 Numbers that are congruent to {2, 10} mod 13.

Original entry on oeis.org

2, 10, 15, 23, 28, 36, 41, 49, 54, 62, 67, 75, 80, 88, 93, 101, 106, 114, 119, 127, 132, 140, 145, 153, 158, 166, 171, 179, 184, 192, 197, 205, 210, 218, 223, 231, 236, 244, 249, 257, 262, 270, 275, 283, 288, 296, 301, 309, 314, 322, 327, 335, 340, 348
Offset: 0

Views

Author

Wolfdieter Lang, May 09 2012

Keywords

Comments

A001844(N) = N^2 + (N+1)^2 = 4*A000217(N) + 1 is divisible by 13 if and only if N=a(n), n>=0. For the proof it suffices to show that only N=2 and N=10 from {0,1,..,12} satisfy A001844(N)== 0 (mod 13). Note that only primes of the form p= 4*k+1 (A002144) can be divisors of A001844 (see a Wolfdieter Lang comment there giving the reference).
Partial sums of the sequence [2,5,8,5,8,5,8,5,8,...] (see the o.g.f., and subtract 2 to see the 5,8 periodicity).

Examples

			Divisibility of A001844 by 13:
n=0: A001844(2) = 13 == 0 (mod 13).
n=3: A001844(23) = 1105 = 85*13 == 0 (mod 13).
However, 8^2 + 9^2 = 145 == 2 (mod 13) is not divisible by 13 because 8 is not a member of the present sequence.
		

Crossrefs

Cf. A047219 (case p=5).

Programs

Formula

Bisection: a(2*n) = 13*n + 2, a(2*n+1) = 13*n + 10, n>=0.
O.g.f.: (2 + 8*x + 3*x^2)/((1-x)*(1-x^2)).
E.g.f.: ((26*x + 11)*exp(x) - 3*exp(-x))/4. - David Lovler, Aug 09 2022

A212161 Numbers congruent to 6 or 10 mod 17.

Original entry on oeis.org

6, 10, 23, 27, 40, 44, 57, 61, 74, 78, 91, 95, 108, 112, 125, 129, 142, 146, 159, 163, 176, 180, 193, 197, 210, 214, 227, 231, 244, 248, 261, 265, 278, 282, 295, 299, 312, 316, 329, 333, 346, 350
Offset: 0

Views

Author

Wolfdieter Lang, May 09 2012

Keywords

Comments

A001844(N) = N^2 + (N+1)^2 = 4*A000217(N) + 1 is divisible by 17 if and only if N = a(n), n >= 0. For the proof it suffices to show that only N=6 and N=10 from {0,1,...,16} satisfy A001844(N) == 0 (mod 17). Note that only primes of the form p = 4*k+1 (A002144) can be divisors of A001844 (see a Wolfdieter Lang comment there giving the reference). Note also that if N^2 + (N+1)^2 == 0 (mod p), with any prime p (necessarily from A002144), then also p-1-N satisfies this congruence. This explains why 10 = 17-1-6 is the (incongruent) companion of 6.
Partial sums of the sequence 6,4,13,4,13,4,13,4,13,4,13,... (see the o.g.f., and subtract 6 to see the remaining 4, 13=17-4 periodicity).

Examples

			Divisibility of A001844 by 17:
n=0: A001844(6) = 85 = 5*17 == 0 (mod 17).
n=2: A001844(23) = 1105 = 5*13*17 == 0 (mod 17).
However, 8^2 + 9^2 = 145 == 9 (mod 17) is not divisible by 17 because 8 is not a term of the present sequence.
		

Crossrefs

Cf. A047219 (p=5), A212160 (p=13).

Programs

  • Magma
    [1/4*(34*n+9*(-1)^n+15): n in [0..60]]; // Vincenzo Librandi, May 24 2012
    
  • Mathematica
    Table[1/4*(34*n+9*(-1)^n+15),{n,0,60}] (* Vincenzo Librandi, May 24 2012 *)
  • PARI
    a(n) = (34*n + 9*(-1)^n + 15)/4 \\ David Lovler, Aug 09 2022

Formula

Bisection: a(2*n) = 17*n + 6, a(2*n+1) = 17*n + 10, n >= 0.
O.g.f.: (6 + 4*x + 7*x^2)/((1-x)*(1-x^2)).
E.g.f.: ((34*x + 15)*exp(x) + 9*exp(-x))/4. - David Lovler, Aug 09 2022

A174814 a(n) = n*(n+1)*(5*n+1)/3.

Original entry on oeis.org

0, 4, 22, 64, 140, 260, 434, 672, 984, 1380, 1870, 2464, 3172, 4004, 4970, 6080, 7344, 8772, 10374, 12160, 14140, 16324, 18722, 21344, 24200, 27300, 30654, 34272, 38164, 42340, 46810, 51584, 56672, 62084, 67830, 73920, 80364, 87172, 94354, 101920, 109880
Offset: 0

Views

Author

Bruno Berselli, Dec 01 2010 - Dec 02 2010

Keywords

Comments

Also zero followed by bisection (even part) of A088003.
Numbers ending in 0, 2 or 4 (cf. 2*A053796(n)). Therefore we can easily see that a(m)^(2*k+1)==-1 (mod 5) only for m in A047219, while a(m)^(2*k)==-1 (mod 5) only for m in A016873 and k odd.

Programs

Formula

G.f.: 2*x*(2+3*x)/(1-x)^4.
a(n) = 2*A033994(n) for n>0.
a(n) = n*A147875(n+1)-sum(k=1..n, A147875(k)) for n>0.
a(-n) = -A144945(n).

A144652 Triangle, read by rows, where T(m,n) = floor((2mn+m+n)/2) with m >= n >= 1.

Original entry on oeis.org

2, 3, 6, 5, 8, 12, 6, 11, 15, 20, 8, 13, 19, 24, 30, 9, 16, 22, 29, 35, 42, 11, 18, 26, 33, 41, 48, 56, 12, 21, 29, 38, 46, 55, 63, 72, 14, 23, 33, 42, 52, 61, 71, 80, 90, 15, 26, 36, 47, 57, 68, 78, 89, 99, 110, 17, 28, 40, 51, 63, 74, 86, 97, 109, 120, 132, 18, 31, 43, 56, 68
Offset: 1

Views

Author

Vincenzo Librandi, Jan 27 2009

Keywords

Comments

From Vincenzo Librandi, Nov 16 2012: (Start)
First column: A007494(n+1);
second column: A047219(n+2);
third column: A047383(n+3);
fourth column: A193910(n+4).
Conjecture: If h does not belong to the sequence, then 4*h+1 is prime. (End)

Examples

			Triangle begins:
2;
3,  6;
5,  8,  12;
6,  11, 15, 20;
8,  13, 19, 24, 30;
9,  16, 22, 29, 35, 42;
11, 18, 26, 33, 41, 48, 56; etc.
		

Crossrefs

Programs

  • Magma
    [Floor((2*n*k+n+k)/2): k in [1..n], n in [1..11]]; // Vincenzo Librandi, Nov 16 2012
  • Mathematica
    Flatten[Table[Floor[(2*n*m + m + n)/2], {n, 1, 20}, {m, n}]] (* Vincenzo Librandi, Nov 16 2012 *)

Extensions

Definition edited (specifying m >= n >= 1), and terms recomputed to match definition, as was done with the similar sequence A140869, by Jon E. Schoenfield, Jun 24 2010
Previous Showing 21-30 of 33 results. Next