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.

A070229 Next m>n such that m is divisible by lpf(n), lpf=A006530 largest prime factor.

Original entry on oeis.org

2, 4, 6, 6, 10, 9, 14, 10, 12, 15, 22, 15, 26, 21, 20, 18, 34, 21, 38, 25, 28, 33, 46, 27, 30, 39, 30, 35, 58, 35, 62, 34, 44, 51, 42, 39, 74, 57, 52, 45, 82, 49, 86, 55, 50, 69, 94, 51, 56, 55, 68, 65, 106, 57, 66, 63, 76, 87, 118, 65, 122
Offset: 1

Views

Author

Reinhard Zumkeller, May 07 2002

Keywords

Comments

Or, for n>1, n + (largest prime divisor of n). [Anne Robinson, daughter of Herman P. Robinson, Oct 08 1981]

Crossrefs

Cf. A036441.
Iterations: A076271 (start=1), A036441 (start=2), A123581 (start=3).

Programs

Formula

a(n) = (n/lpf(n)+1)*lpf(n).
a(n) = n+lpf(n) where lpf=A006530 largest prime factor. For example, a(14)=14+7=21 (instead of ((14/7)+1)*7). - Philippe Lallouet (philip.lallouet(AT)wanadoo.fr), Jun 14 2007

A076271 a(1) = 1, a(2) = 2, and for n > 2, a(n) = a(n-1) + gpf(a(n-1)), where gpf = greatest prime factor = A006530.

Original entry on oeis.org

1, 2, 4, 6, 9, 12, 15, 20, 25, 30, 35, 42, 49, 56, 63, 70, 77, 88, 99, 110, 121, 132, 143, 156, 169, 182, 195, 208, 221, 238, 255, 272, 289, 306, 323, 342, 361, 380, 399, 418, 437, 460, 483, 506, 529, 552, 575, 598, 621, 644, 667, 696, 725, 754, 783, 812, 841
Offset: 1

Views

Author

Reinhard Zumkeller, Oct 04 2002

Keywords

Comments

a(n+1) is the smallest number such that the largest prime divisor of a(n) is the highest common factor of a(n) and a(n+1). - Amarnath Murthy, Oct 17 2002
Essentially the same as A036441(n) = a(n+1) and A180107(n) = a(n-1) (n > 1).
The equivalent sequence with A020639 = spf instead of A006530 = gpf begins a(1) = 1, a(2) = 2, and from then on we get all even numbers: a(n) = a(2) + 2*(n-2), n > 1. - M. F. Hasler, Apr 08 2015
From David James Sycamore, Apr 27 2017: (Start)
The sequence contains only one prime; a(2)=2, all other terms (excluding a(1)=1) being composite, since if a(n) for some n > 2 is assumed to be the first prime after 2, then a(n) = a(n-1) + gpf(a(n-1))= m*q+q = q*(m+1) for some integer m > 1 and some prime q. This number is composite; contradiction. Terms after a(3)=4 alternate between even and odd values since each is created by addition of a prime (odd term).
All terms a(n) arise as consecutive multiples of consecutive primes occurring in their natural ascending order, 2,3,5,7.... (A000040). The number of (consecutive) terms which arise as multiples of p(n)= A000040(n) is 1 + p(n+1)- p(n-1), namely n-th term of the sequence: 2,4,5,7,7,7,7,7,11, etc. Example: Number of multiples of 17, the 7th prime, is 1+p(8)-p(6) = 1+19-13 = 7.
For any pair of consecutive primes, p,q (p < q) a(p+q-1) = p*q, the (semiprime) term where multiples of p end and multiples of q start. Example a(7+11-1) = a(17) = 77 = 11*7, the last multiple of 7 and first multiple of 11. Every string of multiples of prime p contains the term p^2, located at a(2*p-1). E.g.: a(3)=4, a(5)=9, a(9)=25. (End)

Crossrefs

Cf. A036441, A076272(n) = a(n+1) - a(n).
See also A180107.
Cf. A070229.

Programs

  • Haskell
    a076271 n = a076271_list !! (n-1)
    a076271_list = iterate a070229 1  -- Reinhard Zumkeller, Nov 07 2015
  • Mathematica
    NestList[#+FactorInteger[#][[-1,1]]&,1,60] (* Harvey P. Dale, May 11 2015 *)
  • PARI
    print1(n=1);for(i=1,199,print1(","n+=A006530(n))) \\ M. F. Hasler, Apr 08 2015
    

Formula

a(A076274(n)) = A008578(n)^2 for all n.
a(n+1) = A070229(a(n)). - Reinhard Zumkeller, Nov 07 2015

Extensions

Edited by M. F. Hasler, Apr 08 2015

A180107 Partial sums of terms in A180101.

Original entry on oeis.org

0, 1, 1, 2, 4, 6, 9, 12, 15, 20, 25, 30, 35, 42, 49, 56, 63, 70, 77, 88, 99, 110, 121, 132, 143, 156, 169, 182, 195, 208, 221, 238, 255, 272, 289, 306, 323, 342, 361, 380, 399, 418, 437, 460, 483, 506, 529, 552, 575, 598, 621, 644, 667, 696, 725, 754, 783, 812, 841, 870, 899, 930, 961, 992, 1023, 1054, 1085, 1116, 1147, 1184, 1221, 1258, 1295, 1332, 1369, 1406, 1443, 1480, 1517, 1558, 1599, 1640
Offset: 0

Views

Author

N. J. A. Sloane, Jan 16 2011

Keywords

Comments

Apart from initial terms, same as A036441 and A076271, but the simple definition warrants an independent entry.

Crossrefs

A123581 a(1) = 3, a(n) = a(n-1) + greatest prime factor of a(n-1).

Original entry on oeis.org

3, 6, 9, 12, 15, 20, 25, 30, 35, 42, 49, 56, 63, 70, 77, 88, 99, 110, 121, 132, 143, 156, 169, 182, 195, 208, 221, 238, 255, 272, 289, 306, 323, 342, 361, 380, 399, 418, 437, 460, 483, 506, 529, 552, 575, 598, 621, 644, 667, 696, 725, 754, 783, 812, 841, 870
Offset: 1

Views

Author

Ben Paul Thurston, Nov 12 2006

Keywords

Examples

			a(16) = 88 because a(15) is 77 whose largest prime factor is 11 so 77 + 11 = 88.
		

Crossrefs

Essentially the same as A036441 and A076271.
Cf. A070229.

Programs

  • Haskell
    a123581 n = a123581_list !! (n-1)
    a123581_list = iterate a070229 3  -- Reinhard Zumkeller, Nov 07 2015
  • Maple
    A123581:= proc(n) option remember;
    local t;
    t:= procname(n-1);
    t + max(numtheory[factorset](t));
    end proc;
    A123581(1):= 3;
    seq(A123581(n),n=1..100); # Robert Israel, May 18 2014
  • Mathematica
    a[1] = 3; a[n_] := a[n] = a[n - 1] + FactorInteger[a[n - 1]][[ -1, 1]]; Array[a, 56] (* Robert G. Wilson v *)
  • PARI
    {print1(a=3,",");for(n=2,57,print1(a=a+vecmax(factor(a)[,1]),","))} \\ Klaus Brockhaus, Nov 19 2006
    

Formula

a(n+1) = A070229(a(n)). - Reinhard Zumkeller, Nov 07 2015

Extensions

More terms from Robert G. Wilson v and Klaus Brockhaus, Nov 18 2006

A367504 a(1) = 2; for n > 1, a(n) = a(n-1) + 2*gpf(a(n-1)), where gpf(k) = A006530(k) = greatest prime dividing k.

Original entry on oeis.org

2, 6, 12, 18, 24, 30, 40, 50, 60, 70, 84, 98, 112, 126, 140, 154, 176, 198, 220, 242, 264, 286, 312, 338, 364, 390, 416, 442, 476, 510, 544, 578, 612, 646, 684, 722, 760, 798, 836, 874, 920, 966, 1012, 1058, 1104, 1150, 1196, 1242, 1288, 1334, 1392, 1450, 1508, 1566, 1624, 1682, 1740, 1798
Offset: 1

Views

Author

Scott R. Shannon, Nov 21 2023

Keywords

Comments

Conjecture: with the requirement that the prime factorization of each term is written so that the primes are ordered from smallest to largest, the sequence is the lexicographically earliest infinite sequence of distinct positive numbers such that gpf(a(n-1)) * lpf(a(n)) = |a(n) - a(n-1)|, where gpf(k) = A006530(k) = greatest prime factor of k and lpf(k) = A020639(k) = least prime factor of k. In this way the sequence is the ordered prime factorization version of the 'Commas sequence', A121805. One can show that, for such a sequence to be infinite, no odd number can appear. Although for many terms a lower even number can be chosen for the following term, which can lead to even lower numbers for further terms, it is conjectured all such choices will ultimately halt the sequence as a number is eventually reached for which no unused next number exists which follows the required rule for the difference between the terms. Therefore all terms must be larger than the previous, and the earliest such infinite sequence is the given sequence.
See A367465 for the sequence when the requirement that the primes in the factorization of each term must be in order is removed.

Examples

			a(7) = 40 as a(6) = 30 = 2*3*5, thus A006530(30) = 5 and a(7) = a(6) + 2*5 = 30 + 2*5 = 40.
		

Crossrefs

Programs

  • Mathematica
    NestList[#+2FactorInteger[#][[-1,1]]&,2,100] (* Paolo Xausa, Dec 31 2023 *)

Formula

a(n) = 2*A123581(n). The exponent of 2 in a(2n) is 1+A367624(n). - N. J. A. Sloane, Dec 06 2023
Showing 1-5 of 5 results.