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

A244725 a(n) = 5*n^3.

Original entry on oeis.org

0, 5, 40, 135, 320, 625, 1080, 1715, 2560, 3645, 5000, 6655, 8640, 10985, 13720, 16875, 20480, 24565, 29160, 34295, 40000, 46305, 53240, 60835, 69120, 78125, 87880, 98415, 109760, 121945, 135000, 148955, 163840, 179685, 196520, 214375, 233280, 253265
Offset: 0

Views

Author

Vincenzo Librandi, Jul 05 2014

Keywords

Crossrefs

Cf. similar sequences of the type k*n^3: A000578 (k=1), A033431 (k=2), A117642 (k=3), A033430 (k=4), this sequence (k=5), A244726 (k=6), A244727 (k=7), A016743 (k=8), A244728 (k=9), A244729 (k=10), A016767 (k=27), A016803 (k=64), A016851 (k=125), A016911 (k=216), A016983 (k=343), A017067 (k=512), A017163 (k=729), A017271 (k=1000), A017391 (k=1331), A017523 (k=1728).

Programs

  • Magma
    [5*n^3: n in [0..40]];
    
  • Magma
    I:=[0,5,40,135]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..40]];
    
  • Mathematica
    Table[5 n^3, {n, 0, 40}] (* or *) CoefficientList[Series[5 x (1 + 4 x + x^2)/(1 - x)^4, {x, 0, 40}], x]
  • PARI
    a(n)=5*n^3 \\ Charles R Greathouse IV, Oct 07 2015

Formula

G.f.: 5*x*(1 + 4*x + x^2)/(1 - x)^4.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n>3.

A332821 One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 1 (mod 3).

Original entry on oeis.org

2, 5, 9, 11, 12, 16, 17, 21, 23, 28, 30, 31, 39, 40, 41, 47, 49, 52, 54, 57, 59, 66, 67, 70, 72, 73, 75, 76, 83, 87, 88, 91, 96, 97, 100, 102, 103, 109, 111, 116, 126, 127, 128, 129, 130, 133, 135, 136, 137, 138, 148, 149, 154, 157, 159, 165, 167, 168, 169, 172, 175, 179, 180, 183, 184, 186, 190, 191, 197, 203, 211, 212
Offset: 1

Views

Author

Antti Karttunen and Peter Munn, Feb 25 2020

Keywords

Comments

The positive integers are partitioned between A332820, this sequence and A332822.
For each prime p, the terms include exactly one of p and p^2. The primes alternate between this sequence and A332822. This sequence has the primes with odd indexes, those in A031368.
The terms are the even numbers in A332822 halved. The terms are also the numbers m such that 5m is in A332822, and so on for alternate primes: 11, 17, 23 etc. Likewise, the terms are the numbers m such that 3m is in A332820, and so on for alternate primes: 7, 13, 19 etc.
The numbers that are half of the even terms of this sequence are in A332820, which consists exactly of those numbers. The numbers that are one third of the terms that are multiples of 3 are in A332822, which consists exactly of those numbers. For larger primes, an alternating pattern applies as described in the previous paragraph.
If we take each odd term of this sequence and replace each prime in its factorization by the next smaller prime, the resulting number is in A332822, which consists entirely of those numbers.
The product of any 2 terms of this sequence is in A332822, the product of any 3 terms is in A332820, and the product of a term of A332820 and a term of this sequence is in this sequence. So if a number k is present, k^2 is in A332822, k^3 is in A332820, and k^4 is in this sequence.
If k is an even number, exactly one of {k/2, k, 2k} is in the sequence (cf. A191257 / A067368 / A213258); and generally if k is a multiple of a prime p, exactly one of {k/p, k, k*p} is in the sequence.

Crossrefs

Positions of ones in A332823; equivalently, numbers in row 3k+1 of A277905 for some k >= 0.
Subsequences: intersection of A026478 and A066208, A031368 (prime terms), A033431\{0}, A052934\{1}, A069486, A099800, A167747\{1}, A244725\{0}, A244728\{0}, A338911 (semiprime terms).

Programs

  • Mathematica
    Select[Range@ 212, Mod[Total@ #, 3] == 1 &@ Map[#[[-1]]*2^(PrimePi@ #[[1]] - 1) &, FactorInteger[#]] &] (* Michael De Vlieger, Mar 15 2020 *)
  • PARI
    isA332821(n) =  { my(f = factor(n)); (1==((sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2)%3)); };

Formula

{a(n) : n >= 1} = {2 * A332820(k) : k >= 1} U {A003961(A332822(k)) : k >= 1}.
{a(n) : n >= 1} = {A332822(k)^2 : k >= 1} U {A331590(2, A332820(k)) : k >= 1}.

A287335 Nonnegative numbers k such that 3*k + 2 is a cube.

Original entry on oeis.org

2, 41, 170, 443, 914, 1637, 2666, 4055, 5858, 8129, 10922, 14291, 18290, 22973, 28394, 34607, 41666, 49625, 58538, 68459, 79442, 91541, 104810, 119303, 135074, 152177, 170666, 190595, 212018, 234989, 259562, 285791, 313730, 343433, 374954, 408347, 443666, 480965
Offset: 1

Views

Author

Bruno Berselli, May 23 2017

Keywords

Comments

Corresponding cubes are listed in A016791.
Primes in the sequence: 2, 41, 443, 1637, 22973, 34607, 91541, 234989, ...

Crossrefs

Subsequence of A047292.
Cf. A244728: nonnegative k such that 3*k is a cube.
Cf. A121628: nonnegative k such that 3*k + 1 is a cube.

Programs

  • Magma
    [9*n^3-9*n^2+3*n-1: n in [1..40]];
  • Mathematica
    Table[9 n^3 - 9 n^2 + 3 n - 1, {n, 0, 40}]
    LinearRecurrence[{4,-6,4,-1},{2,41,170,443},40] (* Harvey P. Dale, Aug 28 2021 *)
  • Maxima
    makelist(9*n^3-9*n^2+3*n-1, n, 1, 40);
    
  • Python
    [9*n**3-9*n**2+3*n-1 for n in range(1,40)]
    
  • Sage
    [9*n^3-9*n^2+3*n-1 for n in (1..40)]
    

Formula

O.g.f.: x*(2 + 33*x + 18*x^2 + x^3)/(1 - x)^4.
E.g.f.: 1 - (1 - 3*x - 18*x^2 - 9*x^3)*exp(x).
a(n) = 9*n^3 - 9*n^2 + 3*n - 1.
a(n) = A131476(3*n-1) = A212069(3*n-1).

A286930 Integers whose double is a square and whose triple is a cube.

Original entry on oeis.org

0, 72, 4608, 52488, 294912, 1125000, 3359232, 8470728, 18874368, 38263752, 72000000, 127552392, 214990848, 347530248, 542126592, 820125000, 1207959552, 1737904968, 2448880128, 3387303432, 4608000000, 6175160712, 8163353088, 10658584008, 13759414272, 17578125000
Offset: 1

Views

Author

Michel Marcus, May 16 2017

Keywords

Examples

			From _Michael De Vlieger_, May 16 2017: (Start)
72 is a term because 2*72 = 144 = 12^2 and 3*72 = 216 = 6^3.4608 is a term because 2*4608 = 96^2 and 3*4608 = 24^3. (End)
		

Crossrefs

Cf. A001014.
Intersection of A001105 and A244728.

Programs

  • Mathematica
    Array[72 (# - 1)^6 &, 26] (* Michael De Vlieger, May 16 2017 *)
    LinearRecurrence[{7,-21,35,-35,21,-7,1},{0,72,4608,52488,294912,1125000,3359232},30] (* Harvey P. Dale, May 07 2022 *)
  • PARI
    isok(x) = issquare(2*x) && ispower(3*x, 3);
    
  • PARI
    concat(0, Vec(72*x^2*(1 + x)*(1 + 56*x + 246*x^2 + 56*x^3 + x^4) / (1 - x)^7 + O(x^30))) \\ Colin Barker, May 17 2017

Formula

a(n) = 72*(n-1)^6. - David A. Corneth, May 16 2017
O.g.f.: 72*x^2*(1 + x)*(1 + 56*x + 246*x^2 + 56*x^3 + x^4) / (1 - x)^7. - Colin Barker, May 17 2017
E.g.f.: 72*(-1 + (1 - x + x^2 + 10*x^3 + 20*x^4 + 9*x^5 + x^6)*exp(x)). - Bruno Berselli, May 17 2017

Extensions

More terms from Michael De Vlieger, May 16 2017
Showing 1-4 of 4 results.