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-10 of 11 results. Next

A379730 a(1)=1. For n>1 if a(n-1) is a novel term a(n)=A083720(a(n-1)), else a(n) = k*a(n-1), where k(>1) is the number of times a(n-1) has occurred.

Original entry on oeis.org

1, 1, 2, 1, 3, 2, 4, 1, 4, 8, 1, 5, 6, 1, 6, 12, 1, 7, 30, 1, 8, 16, 1, 9, 2, 6, 18, 1, 10, 3, 6, 24, 1, 11, 210, 1, 12, 24, 48, 1, 13, 2310, 1, 14, 15, 2, 8, 24, 72, 1, 15, 30, 60, 1, 16, 32, 1, 17, 30030, 1, 18, 36, 1, 19, 510510, 1, 20, 3, 9, 18, 54, 1, 21
Offset: 1

Views

Author

David James Sycamore, Dec 31 2024

Keywords

Comments

In other words if a(n-1) is a novel term, a(n) is the product of all primes < Gpf(a(n-1)) which do not divide a(n-1), else if a(n-1) seen k times already, up to and including itself then a(n)=k*a(n-1). All positive integers appear eventually in the sequence, and the first occurrences of primes appear in order.
Every squarefree number appears infinitely many times (consequent to both conditions of the definition), whereas numbers m which are not squarefree can appear only from the second condition, and therefore appear finitely many (at most A000005(m)-1) times; see Example.

Examples

			a(1)=1 is a primorial term so a(2)=1. 1 has now occurred twice so a(3)=2*1=2.
a(3)=2 is a novel primorial term so a(4)=1, and then a(5)=3. Since 3 is novel a(6)=A083720(3)=2.
4 appears twice = A5(4)-1 (as 2*2=a(7), and 4*1=a(9)). It never appears again since A083720(m)!=4 for any m, so 4 can only arise as consequence of the second condition of the definition. Similarly 8 appears only 3 (A5(8)-1) times: 2*4=a(10), 8*1=a(21), and 4*2 a(47). However 16 appears not 5 (A5(16)-1) times but 4 times because 4*4 is prevented since 4 appears only twice.
		

Crossrefs

Programs

  • Mathematica
    nn = 120; c[_] := 0; j = 1;
    {j}~Join~Reap[Do[
      If[c[j] == 0,
        k = Product[Prime[i], {i, PrimePi[#[[-1]] ] } ]/Apply[Times, #] &[
          FactorInteger[j][[All, 1]] ]; c[j]++,
        k = ++c[j]*j ];
      j = Sow[k], {nn}] ][[-1, 1]]

Formula

If a(n-1) is a novel term and rad(a(n-1)) is primorial then a(n)=1.

A079067 Number of primes less than greatest prime factor of n but not dividing n.

Original entry on oeis.org

0, 0, 1, 0, 2, 0, 3, 0, 1, 1, 4, 0, 5, 2, 1, 0, 6, 0, 7, 1, 2, 3, 8, 0, 2, 4, 1, 2, 9, 0, 10, 0, 3, 5, 2, 0, 11, 6, 4, 1, 12, 1, 13, 3, 1, 7, 14, 0, 3, 1, 5, 4, 15, 0, 3, 2, 6, 8, 16, 0, 17, 9, 2, 0, 4, 2, 18, 5, 7, 1, 19, 0, 20, 10, 1, 6, 3, 3, 21, 1, 1, 11, 22, 1, 5, 12, 8, 3, 23, 0, 4, 7, 9, 13, 6, 0
Offset: 1

Views

Author

Reinhard Zumkeller, Dec 20 2002

Keywords

Comments

For n >= 2, a(n) is the largest part minus the number of distinct parts of the partition having Heinz number n. The Heinz number of a partition [i_1, i_2, ..., i_r] is defined as Product_{j=1..r} (i_j-th prime) (concept used by Alois P. Heinz in A215366 as an encoding of a partition). For example, for the partition [1, 1, 1, 4] we get 2*2*2*7 = 56; a(56) = 4 - #{1,4} = 2. - Emeric Deutsch, Jun 09 2015 [edited by Peter Munn, Apr 09 2024]

Crossrefs

See the formula section for the relationships with A000040, A001221, A002110, A006530, A049084, A061395, A083720.

Programs

  • Maple
    with(numtheory): a := proc (n) local B: B := proc (n) local nn, j, m: nn := op(2, ifactors(n)): for j to nops(nn) do m[j] := op(j, nn) end do: [seq(seq(pi(op(1, m[i])), q = 1 .. op(2, m[i])), i = 1 .. nops(nn))] end proc: max(B(n))-nops(convert(B(n), set)) end proc: 0, seq(a(n), n = 2 .. 96); # The subprogram B yields the partition having Heinz number n. # Emeric Deutsch, Jun 09 2015
    # second Maple program:
    with(numtheory):
    a:= n-> (s-> pi(max(0, s))-nops(s))(factorset(n)):
    seq(a(n), n=1..100);  # Alois P. Heinz, Sep 03 2019
  • Mathematica
    a[1] = 0; a[n_] := With[{fi = FactorInteger[n]}, PrimePi[fi][[-1, 1]] - Length[fi]]; Array[a, 100] (* Jean-François Alcover, Jan 08 2016 *)
  • PARI
    a(n) = if (n==1, 0, my(pf=factor(n)[,1]); primepi(vecmax(pf)) - #pf); \\ Michel Marcus, May 05 2017

Formula

a(n) = A049084(A006530(n)) - A001221(n) = A061395(n) - A001221(n).
a(n) = 0 iff n = m*prime(k)#, where prime(k)# is the k-th primorial (A002110(k)) and A006530(m) <= A000040(k).
a(A000040(k)) = k - 1.
a(n) = A001221(A083720(n)). - Peter Munn, Apr 09 2024

A083722 Product of primes greater than the greatest prime factor of n but not greater than n.

Original entry on oeis.org

1, 1, 1, 3, 1, 5, 1, 105, 35, 7, 1, 385, 1, 143, 1001, 15015, 1, 85085, 1, 323323, 46189, 4199, 1, 37182145, 7436429, 7429, 37182145, 1062347, 1, 215656441, 1, 100280245065, 86822723, 392863, 955049953, 33426748355, 1, 765049, 247110827, 247357937827, 1, 1448810778701, 1
Offset: 1

Views

Author

Reinhard Zumkeller, May 04 2003

Keywords

Comments

a(n) = 1 iff n is prime or n is 1.
Apart from 1-terms, the other duplicated terms are: a(24) = a(27), a(120) = a(125), a(140) = a(147), a(528) = a(539), etc, whose positions are listed by A293893 and A293894. - Antti Karttunen, Nov 01 2017

Crossrefs

Cf. A293892 (restricted growth sequence transform), A293893, A293894.

Programs

  • Mathematica
    Array[Times @@ Select[Prime@ Range[#1, #1 + #2], Function[p, p <= #3]] & @@ {PrimePi@ NextPrime[FactorInteger[#][[-1, 1]]], PrimePi@ #, #} &, 43] (* Michael De Vlieger, Nov 01 2017 *)
  • PARI
    a(n) = {if (n==1, return (1)); my(gpf = vecmax(factor(n)[,1])); my(pp = 1); forprime(p=gpf+1, n, pp *= p;); pp;} \\ Michel Marcus, Jun 26 2016

Formula

Extensions

More terms from Michel Marcus, Jun 26 2016

A137795 Smallest positive m such that m*n is free of prime gaps in canonical factorization.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 15, 1, 1, 1, 1, 1, 3, 5, 105, 1, 1, 1, 1155, 1, 15, 1, 1, 1, 1, 35, 15015, 1, 1, 1, 255255, 385, 3, 1, 5, 1, 105, 1, 4849845, 1, 1, 1, 3, 5005, 1155, 1, 1, 7, 15, 85085, 111546435, 1, 1, 1, 3234846615, 5, 1, 77, 35, 1, 15015, 1616615, 3, 1, 1
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 11 2008

Keywords

Examples

			n=42: A073490(42) = A073490([2*3]*[7]) = 1,
the gap is filled by a(42) = 5: A073490(42*5) = 0.
		

Crossrefs

Programs

  • PARI
    A137795(n) = if(1==n,1, my(f = factor(n), p = f[1, 1], gpf = f[#f~, 1], m = 1); while(pAntti Karttunen, Sep 06 2018

Formula

A073490(n*a(n)) = 0; A137794(n*a(n)) = 1.
For m < a(n), A073490(n*m) > 0 and A137794(n*m) = 0.
a(A073491(n)) = 1; a(A073492(n)) > 1.
a(n) = A083720(n) / A034386(A020639(n)-1). - Peter Munn, Feb 24 2024

A364154 Lexicographically earliest sequence of distinct positive integers such that a(n) is least novel multiple m of the product of all primes less than the greatest prime factor of a(n-1) which do not divide a(n-1); a(1) = 1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 30, 8, 9, 10, 12, 11, 210, 13, 2310, 14, 15, 16, 17, 30030, 18, 19, 510510, 20, 21, 40, 24, 22, 105, 26, 1155, 28, 45, 32, 23, 9699690, 25, 36, 27, 34, 15015, 38, 255255, 42, 35, 48, 29, 223092870, 31, 6469693230, 33, 70, 39, 770, 51, 10010
Offset: 1

Views

Author

David James Sycamore, Jul 11 2023

Keywords

Comments

It follows from the definition that the sequence is infinite, and that the records (outside of the first 7 terms) are all primorial numbers, meaning that it grows very quickly.
When there are no primes less than the greatest prime factor of a(n-1) which do not divide a(n-1) then m is the least novel multiple of 1, the empty product, and therefore a(n) = u, the least unused number in the sequence so far. The only way a prime can enter the sequence is as u. When a(n-1) = prime(k), a(n) is A002110(k-1), and any primorial term is followed by u. Thus: prime —> primorial —> u.
Sequence is a permutation of the positive integers since by the definition no number appears more than once and m = 1 eventually introduces any number not already placed by the first part of the definition (m > 1).

Examples

			a(1) = 1 and there are no primes < 1 which divide 1 therefore m = 1 so a(2) = 2, the least unused number. Likewise a(3) = 3.
a(4) = 2*2 = 4 since 2 is the only prime < 3 which does not divide 3 and 2 has already occurred.
Since a(7) = 7, a(8) = 2*3*5 = 30.
		

Crossrefs

Programs

  • Mathematica
    nn = 120; c[] := False; m[] := 1; a[1] = j = 1; c[1] = True;
    Do[k = Times @@ Complement[Prime@ Range[PrimePi@ Last[#] - 1], #] &[
       FactorInteger[j][[All, 1]] ];
     While[c[k m[k]], m[k]++]; k *= m[k];
     Set[{a[n], c[k], j}, {k, True, k}], {n, 2, nn}];
    Array[a, nn]
  • PARI
    lista(nn) = my(c, m, v=List([1, 2])); for(k=3, nn, c=m=1; forprime(p=2, vecmax(factor(v[k-1])[, 1]), if(v[k-1]%p, m*=p)); while(setsearch(Set(v), c*m), c++); listput(v, c*m)); Vec(v) \\ Jinyuan Wang, Jul 11 2023

Extensions

More terms from Jinyuan Wang, Jul 11 2023

A375579 Lexicographically earliest sequence of distinct positive integers such that the products of two consecutive terms belong to A055932.

Original entry on oeis.org

1, 2, 3, 4, 6, 5, 12, 8, 9, 10, 15, 14, 30, 7, 60, 16, 18, 20, 21, 40, 24, 25, 36, 27, 32, 45, 28, 75, 42, 35, 48, 50, 54, 64, 72, 70, 33, 140, 63, 80, 81, 90, 49, 120, 56, 105, 22, 210, 11, 420, 44, 315, 66, 175, 84, 55, 126, 100, 96, 108, 125, 144, 128, 135
Offset: 1

Views

Author

Rémy Sigrist, Aug 19 2024

Keywords

Comments

a(n+1) is a multiple of A083720(a(n)).
This sequence has similarities with A175343; here we consider prime factors of consecutive terms, there ones in binary expansions of consecutive terms.

Examples

			The first terms, alongside their prime factors, are:
  n   a(n)  2  3  5  7
  --  ----  -  -  -  -
   1     1
   2     2  X
   3     3     X
   4     4  X
   5     6  X  X
   6     5        X
   7    12  X  X
   8     8  X
   9     9     X
  10    10  X     X
  11    15     X  X
  12    14  X        X
		

Crossrefs

Programs

  • PARI
    \\ See Links section.

A375615 Lexicographically earliest sequence of distinct positive integers such that for any n > 0, n * a(n) belongs to A055932.

Original entry on oeis.org

1, 2, 4, 3, 6, 5, 30, 8, 10, 9, 210, 12, 2310, 15, 14, 16, 30030, 18, 510510, 21, 20, 105, 9699690, 24, 36, 1155, 32, 45, 223092870, 7, 6469693230, 27, 70, 15015, 42, 25, 200560490130, 255255, 770, 48, 7420738134810, 35, 304250263527210, 315, 28, 4849845
Offset: 1

Views

Author

Rémy Sigrist, Aug 21 2024

Keywords

Comments

This sequence is a self-inverse permutation of the positive integers.
The n-th term is a multiple of A083720(n).

Examples

			The first terms, alongside the prime factors of n*a(n), are:
  n   a(n)  Prime factors
  --  ----  --------------------
   1     1  {}
   2     2  {2}
   3     4  {2, 3}
   4     3  {2, 3}
   5     6  {2, 3, 5}
   6     5  {2, 3, 5}
   7    30  {2, 3, 5, 7}
   8     8  {2}
   9    10  {2, 3, 5}
  10     9  {2, 3, 5}
  11   210  {2, 3, 5, 7, 11}
  12    12  {2, 3}
  13  2310  {2, 3, 5, 7, 11, 13}
  14    15  {2, 3, 5, 7}
  15    14  {2, 3, 5, 7}
		

Crossrefs

Programs

  • PARI
    \\ See Links section.

A372699 a(1) = 1, j = a(n-1). For n > 1, a(n) is the least novel m such that rad(m*j) is the smallest primorial number divisible by A007947(j).

Original entry on oeis.org

1, 2, 3, 4, 6, 5, 12, 10, 9, 8, 18, 15, 16, 24, 20, 27, 32, 36, 25, 30, 7, 60, 14, 45, 40, 48, 50, 54, 75, 64, 72, 80, 81, 96, 90, 21, 70, 42, 35, 84, 100, 108, 120, 28, 105, 56, 135, 128, 144, 125, 150, 49, 180, 63, 140, 126, 160, 162, 200, 192, 225, 216, 240
Offset: 1

Views

Author

Keywords

Comments

a(n) depends upon if rad(j) = A002110(k) for some k (equivalently A083720(j) = 1), or not. If so a(n) is least novel m such that rad(m*j) = A002110(k+1). Otherwise a(n) = least novel m such that rad(m*j) = A002110(A000720(q)), where q = gpf(j).
Put otherwise, if p = nextprime(q), and A = A083720, then for n > 1 if A(j) = 1, a(n) is the least novel p-smooth number divisible by p, and if A(j) = w > 1, a(n) is the least novel q-smooth number divisible by w.
If j is a term in A002110, a(n) = smallest prime which has not yet appeared in the sequence (e.g., 1-->2, 2-->3, 6-->5, 30-->7, 210-->11, and so on).
Primes are in order and if p is prime and p|a(n) there is an i <= n such that a(i) = p (no multiple of p appears prior to p). Sequence is conjectured to have "Property S" of A368900. Also, for integers x, y with x < y and rad(x) = rad(y), x appears in the sequence before y. Conjecture: Sequence is a permutation of the positive integers which preserves the above mentioned properties of A000027.

Examples

			a(1) = 1 = A002110(0) so a(2) = 2 since rad(2*1) = 2 = A002110(1).
a(2) = 2 = A002110(1) so a(3) = 3 since then rad(3*2) = 6 = A002110(2).
a(6) = 5 is not primorial so a(7) = 12, since rad(12*5) = 30 = A002110(3) and no smaller novel number has this property.
		

Crossrefs

Programs

  • Mathematica
    nn = 120;
    c[] := False; m[] := 1;
    Array[Set[{a[#], c[#], m[#]}, {#, True, 2}] &, 2]; j = a[2];
    f[x_] := f[x] = FactorInteger[x][[All, 1]];
    Do[(If[Divisible[j, Set[q, Times @@ Prime@ Range[#2]]],
           q = Times @@ Prime@ Range[#2 + 1]];
           k = q/#3; mm = m[k];
           While[Or[c[k*mm], ! Divisible[q, Times @@ f[mm]]], mm++];
           k *= mm; While[c[k*m[k]], m[k]++]) & @@
         {#, PrimePi@ Last[#], Times @@ #} &[f[j]];
      Set[{a[n], c[k], j}, {k, True, k}], {n, 3, nn}];
    Array[a, nn]

Formula

a(n) = A002110(k) --> a(n+1) = prime(k+1); a(n+2) = 2*a(n), a(n+3) = 2*prime(k+1) (n >= 1).

Extensions

More terms from Michael De Vlieger, May 18 2024

A377193 Lexicographically earliest infinite sequence of distinct positive integers such that any term j = a(n-1) with primorial kernel is followed by a prime, whereas any other term is followed by a number with prime factors p < q = Gpf(j) which do not divide j.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 11, 9, 16, 13, 10, 27, 32, 17, 12, 19, 14, 15, 64, 23, 18, 29, 20, 81, 128, 31, 21, 25, 24, 37, 22, 35, 36, 41, 26, 33, 28, 45, 256, 43, 30, 47, 34, 39, 40, 243, 512, 53, 38, 49, 48, 59, 42, 125, 54, 61, 44, 63, 50, 729, 1024, 67, 46, 51
Offset: 1

Views

Author

Keywords

Comments

Following j = a(n-1), a term in A005932, a(n) is the smallest prime not already listed. Otherwise a(n) = smallest novel product of powers of non divisor primes of j; a number of the form: Product_{i = 0..k} p_i^e_i; p_i a prime < q = Gpf(j) which does not divide j, e_i >= 0, k = the number of primes p_i < q which do not divide j.
Adjacent terms are coprime and the greedy algorithm implied by the definition forces naked prime p to appear in advance of any multiple m*p of p; m >1.
Prime powers enter the sequence early, consequent to j having a single non divisor prime. A power of 3 is always followed by a power of 2.
Conjectures:
(i) A permutation of the positive integers in which the primes appear in order.
(ii)The sequence obeys Selcoe's theorem (see A280864) regarding numbers that have the same squarefree kernel, namely: Construct a sequence S_r = { m*r : rad(m) | r } = { k : rad(k) = r }, squarefree r. Terms w in S_r appear in this sequence in order. This is to say, for example, that for r = 6, terms in A033845 = {6, 12, 18, 24, 36, 48, 54, ...} appear in order.

Examples

			a(1) = 1 implies a(2) = 2 since A007947(1) = A002110(1) = 1, and 2 is the earliest unrecorded prime so far, and likewise a(3) = 3. Since rad(3) = 3 is not a primorial number a(4) = 2^2 = 4, the smallest novel number derived from 2, the only non divisor prime of 3 and < 3.
a(8) = 8 implies a(9) = 11 because 8 is a term in A055932. The non divisor primes of 11 and < 11 are 2,3,5,7  and the smallest number which can be composed using some or all of these primes is a(10) = 3^2 = 9 (since 2,3,4,5,6,7,8 have all occurred previously). Consequently a(11) = 2^4 = 16, the smallest novel power of 2.
a(195) = 154 = 2*7*11, the non divisor primes < 11 are 3 and 5, so a(196) = 405 = 3^4*5 since all smaller candidates (3,5,9,15,25,45,75,81,125,135,243,375) have already appeared.
		

Crossrefs

A115333 Sum of primes that do not divide n and are less than the largest prime dividing n.

Original entry on oeis.org

0, 0, 2, 0, 5, 0, 10, 0, 2, 3, 17, 0, 28, 8, 2, 0, 41, 0, 58, 3, 7, 15, 77, 0, 5, 26, 2, 8, 100, 0, 129, 0, 14, 39, 5, 0, 160, 56, 25, 3, 197, 5, 238, 15, 2, 75, 281, 0, 10, 3, 38, 26, 328, 0, 12, 8, 55, 98, 381, 0, 440, 127, 7, 0, 23, 12, 501, 39, 74, 3, 568, 0, 639, 158, 2, 56, 10
Offset: 1

Views

Author

Leroy Quet, Mar 05 2006

Keywords

Comments

When n is prime, n = largest prime dividing n; hence a(n) is the sum of all primes less than n = A034387(n)-n. a(n) = SUM{p such that p is in A000040 AND NOT(p|n) AND p < A006530(n)}. - Jonathan Vos Post, Mar 08 2006
The zeros give A055932: All prime divisors are consecutive primes starting at 2. - Robert G. Wilson v, May 01 2006

Examples

			The primes < 7 and coprime to 7 are 2, 3 and 5. So a(7) = 2+3+5 = 10.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Plus @@ Complement[Prime@ Range@ PrimePi[ Max[First /@ FactorInteger@n] - 1], First /@ FactorInteger@n]; Array[f, 77] (* Hans Havermann, Mar 06 2006 *)

Extensions

More terms from Hans Havermann, Mar 06 2006
Showing 1-10 of 11 results. Next