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

A276157 a(n) = {greatest primorial less than or equal to n} divided by {greatest primorial number which divides n} = A260188(n)/A053589(n).

Original entry on oeis.org

1, 1, 2, 1, 2, 1, 6, 3, 6, 3, 6, 1, 6, 3, 6, 3, 6, 1, 6, 3, 6, 3, 6, 1, 6, 3, 6, 3, 6, 1, 30, 15, 30, 15, 30, 5, 30, 15, 30, 15, 30, 5, 30, 15, 30, 15, 30, 5, 30, 15, 30, 15, 30, 5, 30, 15, 30, 15, 30, 1, 30, 15, 30, 15, 30, 5, 30, 15, 30, 15, 30, 5, 30, 15, 30, 15, 30, 5, 30, 15, 30
Offset: 1

Views

Author

Antti Karttunen, Aug 30 2016

Keywords

Crossrefs

Cf. A060735 (positions of ones).

Programs

  • Mathematica
    Table[(j = #; While[! Divisible[n, Times @@ Prime@ Range@ j], j--]; (Times @@ Prime@ Range@ #)/(Times @@ Prime@ Range@ j)) &@ (k = 0; While[Times @@ Prime@ Range[k + 1] <= n, k++]; k), {n, 120}] (* Michael De Vlieger, Aug 30 2016 *)
  • Scheme
    (define (A276157 n) (/ (A260188 n) (A053589 n)))

Formula

a(n) = A260188(n)/A053589(n).

A002110 Primorial numbers (first definition): product of first n primes. Sometimes written prime(n)#.

Original entry on oeis.org

1, 2, 6, 30, 210, 2310, 30030, 510510, 9699690, 223092870, 6469693230, 200560490130, 7420738134810, 304250263527210, 13082761331670030, 614889782588491410, 32589158477190044730, 1922760350154212639070, 117288381359406970983270, 7858321551080267055879090
Offset: 0

Views

Author

Keywords

Comments

See A034386 for the second definition of primorial numbers: product of primes in the range 2 to n.
a(n) is the least number N with n distinct prime factors (i.e., omega(N) = n, cf. A001221). - Lekraj Beedassy, Feb 15 2002
Phi(n)/n is a new minimum for each primorial. - Robert G. Wilson v, Jan 10 2004
Smallest number stroked off n times after the n-th sifting process in an Eratosthenes sieve. - Lekraj Beedassy, Mar 31 2005
Apparently each term is a new minimum for phi(x)*sigma(x)/x^2. 6/Pi^2 < sigma(x)*phi(x)/x^2 < 1 for n > 1. - Jud McCranie, Jun 11 2005
Let f be a multiplicative function with f(p) > f(p^k) > 1 (p prime, k > 1), f(p) > f(q) > 1 (p, q prime, p < q). Then the record maxima of f occur at n# for n >= 1. Similarly, if 0 < f(p) < f(p^k) < 1 (p prime, k > 1), 0 < f(p) < f(q) < 1 (p, q prime, p < q), then the record minima of f occur at n# for n >= 1. - David W. Wilson, Oct 23 2006
Wolfe and Hirshberg give ?, ?, ?, ?, ?, 30030, ?, ... as a puzzle.
Records in number of distinct prime divisors. - Artur Jasinski, Apr 06 2008
For n >= 2, the digital roots of a(n) are multiples of 3. - Parthasarathy Nambi, Aug 19 2009 [with corrections by Zak Seidov, Aug 30 2015]
Denominators of the sum of the ratios of consecutive primes (see A094661). - Vladimir Joseph Stephan Orlovsky, Oct 24 2009
Where record values occur in A001221. - Melinda Trang (mewithlinda(AT)yahoo.com), Apr 15 2010
It can be proved that there are at least T prime numbers less than N, where the recursive function T is: T = N - N*Sum_{i = 0..T(sqrt(N))} A005867(i)/A002110(i). This can show for example that at least 0.16*N numbers are primes less than N for 29^2 > N > 23^2. - Ben Paul Thurston, Aug 23 2010
The above comment from Parthasarathy Nambi follows from the observation that digit summing produces a congruent number mod 9, so the digital root of any multiple of 3 is a multiple of 3. prime(n)# is divisible by 3 for n >= 2. - Christian Schulz, Oct 30 2013
The peaks (i.e., local maximums) in a graph of the number of repetitions (i.e., the tally of values) vs. value, as generated by taking the differences of all distinct pairs of odd prime numbers within a contiguous range occur at regular periodic intervals given by the primorial numbers 6 and greater. Larger primorials yield larger (relative) peaks, however the range must be >50% larger than the primorial to be easily observed. Secondary peaks occur at intervals of those "near-primorials" divisible by 6 (e.g., 42). See A259629. Also, periodicity at intervals of 6 and 30 can be observed in the local peaks of all possible sums of two, three or more distinct odd primes within modest contiguous ranges starting from p(2) = 3. - Richard R. Forberg, Jul 01 2015
If a number k and a(n) are coprime and k < (prime(n+1))^b < a(n), where b is an integer, then k has fewer than b prime factors, counting multiplicity (i.e., bigomega(k) < b, cf. A001222). - Isaac Saffold, Dec 03 2017
If n > 0, then a(n) has 2^n unitary divisors (A034444), and a(n) is a record; i.e., if k < a(n) then k has fewer unitary divisors than a(n) has. - Clark Kimberling, Jun 26 2018
Unitary superabundant numbers: numbers k with a record value of the unitary abundancy index, A034448(k)/k > A034448(m)/m for all m < k. - Amiram Eldar, Apr 20 2019
Psi(n)/n is a new maximum for each primorial (psi = A001615) [proof in link: Patrick Sole and Michel Planat, proposition 1 page 2]; compare with comment 2004: Phi(n)/n is a new minimum for each primorial. - Bernard Schott, May 21 2020
The term "primorial" was coined by Harvey Dubner (1987). - Amiram Eldar, Apr 16 2021
a(n)^(1/n) is approximately (n log n)/e. - Charles R Greathouse IV, Jan 03 2023
Subsequence of A267124. - Frank M Jackson, Apr 14 2023

Examples

			a(9) = 23# = 2*3*5*7*11*13*17*19*23 = 223092870 divides the difference 5283234035979900 in the arithmetic progression of 26 primes A204189. - _Jonathan Sondow_, Jan 15 2012
		

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 50.
  • G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, Cambridge, University Press, 1940, p. 49.
  • P. Ribenboim, The Book of Prime Number Records. Springer-Verlag, NY, 2nd ed., 1989, p. 4.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 114.
  • D. Wolfe and S. Hirshberg, Underspecified puzzles, in Tribute to A Mathemagician, Peters, 2005, pp. 73-74.

Crossrefs

A034386 gives the second version of the primorial numbers.
Subsequence of A005117 and of A064807. Apart from the first term, a subsequence of A083207.
Cf. A001615, A002182, A002201, A003418, A005235, A006862, A034444 (unitary divisors), A034448, A034387, A033188, A035345, A035346, A036691 (compositorial numbers), A049345 (primorial base representation), A057588, A060735 (and integer multiples), A061742 (squares), A072938, A079266, A087315, A094348, A106037, A121572, A053589, A064648, A132120, A260188.
Cf. A061720 (first differences), A143293 (partial sums).
Cf. also A276085, A276086.
The following fractions are all related to each other: Sum 1/n: A001008/A002805, Sum 1/prime(n): A024451/A002110 and A106830/A034386, Sum 1/nonprime(n): A282511/A282512, Sum 1/composite(n): A250133/A296358.

Programs

  • Haskell
    a002110 n = product $ take n a000040_list
    a002110_list = scanl (*) 1 a000040_list
    -- Reinhard Zumkeller, Feb 19 2012, May 03 2011
    
  • Magma
    [1] cat [&*[NthPrime(i): i in [1..n]]: n in [1..20]]; // Bruno Berselli, Oct 24 2012
    
  • Magma
    [1] cat [&*PrimesUpTo(p): p in PrimesUpTo(60)]; // Bruno Berselli, Feb 08 2015
    
  • Maple
    A002110 := n -> mul(ithprime(i),i=1..n);
  • Mathematica
    FoldList[Times, 1, Prime[Range[20]]]
    primorial[n_] := Product[Prime[i], {i, n}]; Array[primorial,20] (* José María Grau Ribas, Feb 15 2010 *)
    Join[{1}, Denominator[Accumulate[1/Prime[Range[20]]]]] (* Harvey P. Dale, Apr 11 2012 *)
  • PARI
    a(n)=prod(i=1,n, prime(i)) \\ Washington Bomfim, Sep 23 2008
    
  • PARI
    p=1; for (n=0, 100, if (n, p*=prime(n)); write("b002110.txt", n, " ", p) )  \\ Harry J. Smith, Nov 13 2009
    
  • PARI
    a(n) = factorback(primes(n)) \\ David A. Corneth, May 06 2018
    
  • Python
    from sympy import primorial
    def a(n): return 1 if n < 1 else primorial(n)
    [a(n) for n in range(51)]  # Indranil Ghosh, Mar 29 2017
    
  • Sage
    [sloane.A002110(n) for n in (1..20)] # Giuseppe Coppoletta, Dec 05 2014
    
  • Scheme
    ; with memoization-macro definec
    (definec (A002110 n) (if (zero? n) 1 (* (A000040 n) (A002110 (- n 1))))) ;; Antti Karttunen, Aug 30 2016

Formula

Asymptotic expression for a(n): exp((1 + o(1)) * n * log(n)) where o(1) is the "little o" notation. - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 08 2001
a(n) = A054842(A002275(n)).
Binomial transform = A136104: (1, 3, 11, 55, 375, 3731, ...). Equals binomial transform of A121572: (1, 1, 3, 17, 119, 1509, ...). - Gary W. Adamson, Dec 14 2007
a(0) = 1, a(n+1) = prime(n)*a(n). - Juri-Stepan Gerasimov, Oct 15 2010
a(n) = Product_{i=1..n} A000040(i). - Jonathan Vos Post, Jul 17 2008
a(A051838(n)) = A116536(n) * A007504(A051838(n)). - Reinhard Zumkeller, Oct 03 2011
A000005(a(n)) = 2^n. - Carlos Eduardo Olivieri, Jun 16 2015
a(n) = A035345(n) - A005235(n) for n > 0. - Jonathan Sondow, Dec 02 2015
For all n >= 0, a(n) = A276085(A000040(n+1)), a(n+1) = A276086(A143293(n)). - Antti Karttunen, Aug 30 2016
A054841(a(n)) = A002275(n). - Michael De Vlieger, Aug 31 2016
a(n) = A270592(2*n+2) - A270592(2*n+1) if 0 <= n <= 4 (conjectured for all n by Alon Kellner). - Jonathan Sondow, Mar 25 2018
Sum_{n>=1} 1/a(n) = A064648. - Amiram Eldar, Oct 16 2020
Sum_{n>=1} (-1)^(n+1)/a(n) = A132120. - Amiram Eldar, Apr 12 2021
Theta being Chebyshev's theta function, a(0) = exp(theta(1)), and for n > 0, a(n) = exp(theta(m)) for A000040(n) <= m < A000040(n+1) where m is an integer. - Miles Englezou, Nov 26 2024

A276150 Sum of digits when n is written in primorial base (A049345); minimal number of primorials (A002110) that add to n.

Original entry on oeis.org

0, 1, 1, 2, 2, 3, 1, 2, 2, 3, 3, 4, 2, 3, 3, 4, 4, 5, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 1, 2, 2, 3, 3, 4, 2, 3, 3, 4, 4, 5, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 2, 3, 3, 4, 4, 5, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 6, 7, 7, 8, 8, 9, 3, 4, 4, 5, 5, 6, 4, 5, 5, 6, 6, 7, 5, 6, 6, 7, 7, 8, 6, 7, 7, 8, 8, 9, 7, 8, 8, 9, 9, 10, 4
Offset: 0

Views

Author

Antti Karttunen, Aug 22 2016

Keywords

Comments

The sum of digits of n in primorial base is odd if n is 1 or 2 (mod 4) and even if n is 0 or 3 (mod 4). Proof: primorials are 1 or 2 (mod 4) and a(n) can be constructed via the greedy algorithm. So if n = 4k + r where 0 <= r < 4, 4k needs an even number of primorials and r needs hammingweight(r) = A000120(r) primorials. Q.E.D. - David A. Corneth, Feb 27 2019

Examples

			For n=24, which is "400" in primorial base (as 24 = 4*(3*2*1) + 0*(2*1) + 0*1, see A049345), the sum of digits is 4, thus a(24) = 4.
		

Crossrefs

Cf. A333426 [k such that a(k)|k], A339215 [numbers not of the form x+a(x) for any x], A358977 [k such that gcd(k, a(k)) = 1].
Cf. A014601, A042963 (positions of even and odd terms), A343048 (positions of records).
Differs from analogous A034968 for the first time at n=24.

Programs

  • Mathematica
    nn = 120; b = MixedRadix[Reverse@ Prime@ NestWhileList[# + 1 &, 1, Times @@ Prime@ Range[# + 1] <= nn &]]; Table[Total@ IntegerDigits[n, b], {n, 0, nn}] (* Version 10.2, or *)
    nn = 120; f[n_] := Block[{a = {{0, n}}}, Do[AppendTo[a, {First@ #, Last@ #} &@ QuotientRemainder[a[[-1, -1]], Times @@ Prime@ Range[# - i]]], {i, 0, #}] &@ NestWhile[# + 1 &, 0, Times @@ Prime@ Range[# + 1] <= n &]; Rest[a][[All, 1]]]; Table[Total@ f@ n, {n, 0, 120}] (* Michael De Vlieger, Aug 26 2016 *)
  • PARI
    A276150(n) = { my(s=0, p=2, d); while(n, d = (n%p); s += d; n = (n-d)/p; p = nextprime(1+p)); (s); }; \\ Antti Karttunen, Feb 27 2019
  • Python
    from sympy import prime, primefactors
    def Omega(n): return 0 if n==1 else Omega(n//primefactors(n)[0]) + 1
    def a276086(n):
        i=0
        m=pr=1
        while n>0:
            i+=1
            N=prime(i)*pr
            if n%N!=0:
                m*=(prime(i)**((n%N)/pr))
                n-=n%N
            pr=N
        return m
    def a(n): return Omega(a276086(n))
    print([a(n) for n in range(201)]) # Indranil Ghosh, Jun 23 2017
    

Formula

a(n) = 1 + a(A276151(n)) = 1 + a(n-A002110(A276084(n))), a(0) = 0.
or for n >= 1: a(n) = 1 + a(n-A260188(n)).
Other identities and observations. For all n >= 0:
a(n) = A001222(A276086(n)) = A001222(A278226(n)).
a(n) >= A371091(n) >= A267263(n).
From Antti Karttunen, Feb 27 2019: (Start)
a(n) = A000120(A277022(n)).
a(A283477(n)) = A324342(n).
(End)
a(n) = A373606(n) + A373607(n). - Antti Karttunen, Jun 19 2024

A060735 a(1)=1, a(2)=2; thereafter, a(n) is the smallest number m not yet in the sequence such that every prime that divides a(n-1) also divides m.

Original entry on oeis.org

1, 2, 4, 6, 12, 18, 24, 30, 60, 90, 120, 150, 180, 210, 420, 630, 840, 1050, 1260, 1470, 1680, 1890, 2100, 2310, 4620, 6930, 9240, 11550, 13860, 16170, 18480, 20790, 23100, 25410, 27720, 30030, 60060, 90090, 120120, 150150, 180180, 210210
Offset: 1

Views

Author

Robert G. Wilson v, Apr 23 2001

Keywords

Comments

Also, numbers k at which k / (phi(k) + 1) increases.
Except for the initial 1, this sequence is a primorial (A002110) followed by its multiples until the next primorial, then the multiples of that primorial and so on. - Wilfredo Lopez (chakotay147138274(AT)yahoo.com), Dec 28 2006
a(1)=1, a(2)=2. For n >= 3, a(n) is the smallest integer > a(n-1) that is divisible by every prime which divides lcm(a(1), a(2), a(3), ..., a(n)). - Leroy Quet, Feb 23 2010
Numbers n for which A053589(n) = A260188(n), thus numbers with only one nonzero digit when written in primorial base A049345. - Antti Karttunen, Aug 30 2016
Lexicographically earliest infinite sequence of distinct positive numbers with property that every prime that divides a(n-1) also divides a(n). - N. J. A. Sloane, Apr 08 2022

Examples

			After a(2)=2 the next term must be even, so a(3)=4.
Then a(4) must be even so a(4) = 6.
Now a(5) must be a multiple of 2*3=6, so a(5)=12.
Then a(6)=18, a(7)=24, a(8)=30.
Now a(9) must be a multiple of 2*3*5 = 30, so a(9)=60. And so on.
		

Crossrefs

Indices of ones in A276157 and A267263.
One more than A343048.

Programs

  • Maple
    seq(seq(k*mul(ithprime(i),i=1..n-1),k=1..ithprime(n)-1),n=1..10); # Vladeta Jovovic, Apr 08 2004
    a := proc(n) option remember; if n=1 then return 1 fi; a(n-1);
    % + convert(numtheory:-factorset(%), `*`) end:
    seq(a(n), n=1..42); # after Zumkeller, Peter Luschny, Aug 30 2016
  • Mathematica
    a = 0; Do[ b = n/(EulerPhi[ n ] + 1); If[ b > a, a = b; Print[ n ] ], {n, 1, 10^6} ]
    f[n_] := Range[Prime[n + 1] - 1] Times @@ Prime@ Range@ n;  Array[f, 7, 0] // Flatten (* Robert G. Wilson v, Jul 22 2015 *)
  • PARI
    first(n)=my(v=vector(n),k=1,p=1,P=1); v[1]=1; for(i=2,n, v[i]=P*k++; if(k>p && isprime(k), p=k; P=v[i]; k=1)); v \\ Charles R Greathouse IV, Jul 22 2015
    
  • PARI
    is_A060735(n,P=1)={forprime(p=2,,n>(P*=p)||return(1);n%P&&return)} \\ M. F. Hasler, Mar 14 2017
    
  • Python
    from functools import cache;
    from sympy import primefactors, prod
    @cache
    def a(n): return 1 if n == 0 else a(n-1) + prod(primefactors(a(n-1)))
    print([a(n) for n in range(42)]) # Trey Deitch, Jun 08 2024

Formula

a(1) = 1, a(n) = a(n-1) + rad(a(n-1)) with rad=A007947, squarefree kernel. - Reinhard Zumkeller, Apr 10 2006
a(A101301(n)+1) = A002110(n). - Enrique Pérez Herrero, Jun 10 2012
a(n) = 1 + A343048(n). - Antti Karttunen, Nov 14 2024

Extensions

Definition corrected by Franklin T. Adams-Watters, Apr 16 2009
Simpler definition, comments, examples from N. J. A. Sloane, Apr 08 2022

A053589 Greatest primorial number (A002110) which divides n.

Original entry on oeis.org

1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 30, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 30, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 30, 1, 2, 1, 2, 1, 6, 1, 2, 1, 2, 1, 6
Offset: 1

Views

Author

Frederick Magata (frederick.magata(AT)uni-muenster.de), Jan 19 2000

Keywords

Examples

			a(30) = 30 because 30=2*3*5, a(15) = 1 because 15=3*5.
		

Crossrefs

Programs

  • Maple
    N:= 1000: # to get a(1)..a(N)
    P:= 1: p:= 1:
    A:= Vector(N,1):
    do
      p:= nextprime(p);
      P:= P*p;
      if P > N then break fi;
      A[[seq(i,i=P..N,P)]]:= P;
    od:
    convert(A,list); # Robert Israel, Aug 30 2016
  • Mathematica
    Table[k = 1; While[Divisible[n, Times @@ Prime@ Range@ k], k++]; Times @@ Prime@ Range[k - 1], {n, 120}] (* Michael De Vlieger, Aug 30 2016 *)
  • PARI
    a(n)=my(f=factor(n), r = 1, k = 1, p); while(k<=matsize(f)[1], p=prime(k); if(f[k,1]!=p,return(r));r*=p; k++) ; r
    a(n) = my(r = 1, p = 2); while(n/p==n\p, r*=p; p=nextprime(p+1));r
    \\ list of all terms up to n#.
    lista(n) = my(l = List([1]),k,s=1); forprime(i=2,n, for(j=1,i-1, for(k=1,s, listput(l,l[k]))); l[#l]*=i; s=#l); l \\ David A. Corneth, Aug 30 2016
    
  • PARI
    a(n)=my(s=1); forprime(p=2,, if(n%p, return(s), s *= p)) \\ Charles R Greathouse IV, Sep 07 2016
  • Scheme
    (define (A053589 n) (A002110 (A276084 n))) ;; Antti Karttunen, Aug 30 2016
    

Formula

From Antti Karttunen, Aug 30 2016: (Start)
a(n) = A002110(A276084(n)).
a(n) = n/A111701(n).
A276157(n) = A260188(n)/a(n).
(End)

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Oct 02 2000

A260187 a(n) = n modulo the greatest primorial <= n.

Original entry on oeis.org

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

Views

Author

Jean-Marc Rebert, Jul 18 2015

Keywords

Comments

I think this can help in finding prime numbers.
If n>2 and a(n)=0 then n is not prime.
If n>2 and 2 <= a(n) <= p with p = A006530(A260188(n)) (i.e., greatest prime factor of A260188(n)) then n is not prime.
If n>2 and (a(n) mod k) = 0 and 2 <= k <= p with p = A006530(A260188(n)) (i.e., greatest prime factor of A260188(n)) then n is not prime.
Alternative definition: count up from 0 to primorial(n)-1, prime(n+1)-1 times, where primorial(n) is A002110(n). - Franklin T. Adams-Watters, Jul 20 2015

Examples

			a(5) = 1 because 5 modulo 2# = 1 and 2# = 2 is the greatest primorial <= 5. (3# = 2*3 = 6 > 5)
		

Crossrefs

Cf. A034386 and A002110 (primorials), A260188, A257687.

Programs

  • Maple
    N:= 100: # to get a(1) to a(N)
    P:= 1: p:= 2: R:= 2:
    for n from 1 to N do
      if n >= R then
        P:= R; p:= nextprime(p); R:= P*p;
      fi;
      A[n]:= n mod P;
    od:
    seq(A[i],i=1..N); # Robert Israel, Jul 20 2015
  • Mathematica
    s = Product[Prime@ n, {n, #}] & /@ Range[0, 6]; Table[Mod[n, Last@ Select[s, # <= n &]], {n, 120}] (* Michael De Vlieger, Jul 20 2015 *)
    f[n_] := Block[{m = p = 1}, While[p*(m + 1) <= n, p = p*m; m = NextPrime@ m]; Mod[n, p]]; Array[f, 101] (* Robert G. Wilson v, Jul 21 2015 *)
  • PARI
    a(n)=my(t=1, k); forprime(p=2, , k=t*p; if(k>n, return(n%t), t=k)) \\ Charles R Greathouse IV, Jul 20 2015

Formula

a(n) = n mod A260188(n).
a(n) <= (n+1)/2. - Charles R Greathouse IV, Jul 20 2015

A258937 Define f_i as the i-th iterate of A260187. a(n) is the least prime for which f_i(a(n)) is prime for all i such that 0 <= i < n and f_n(a(n)) is not prime.

Original entry on oeis.org

2, 11, 41, 251, 2579, 32609, 543131, 10243031, 233336819, 6703033091, 207263540933, 7628002016027, 311878266460847, 13394639616667427, 628284422215925129, 33217442899664876729, 1955977793054900415107, 119244359152469819863541
Offset: 1

Views

Author

Jean-Marc Rebert, Nov 07 2015

Keywords

Comments

If p is prime, we replace p with A260187(p), until A260187(p) is not prime.
a(n) is the least prime for which the number of steps is n.

Examples

			a(3) = 41, because A260187(41) = 11, A260187(11) = 5 are prime and A260187(5) = 1 is not prime.
Trajectories of the first few terms:
2->0
11->5->1
41->11->5->1
251->41->11->5->1
2579->269->59->29->5->1
32609->2579->269->59->29->5->1
543131->32621->2591->281->71->11->5->1
10243031->543341->32831->2801->491->71->11->5->1
233336819->10243949->544259->33749->3719->1409->149->29->5->1
6703033091->233339861->10246991->547301->36791->6761->2141->41->11->5->1.
		

Crossrefs

Programs

  • PARI
    a260187(n)=my(t=1, k); forprime(p=2, , k=t*p; if(k>n, return(n%t), t=k));
    isok(k,n) = {for (j=1, n-1, nk = 260187(k); if (! isprime(nk), return (0)); k = nk;); ! isprime(a260187(k));}
    a(n) = {my(k = 2); while(! isok(k, n), k = nextprime(k+1)); k;} \\ Michel Marcus, Nov 16 2015

A377884 Composite numbers k without prime factors that are divisors of the greatest primorial less than k.

Original entry on oeis.org

25, 49, 77, 91, 119, 121, 133, 143, 161, 169, 187, 203, 209, 221, 247, 253, 289, 299, 319, 323, 341, 361, 377, 391, 403, 407, 437, 451, 473, 481, 493, 517, 527, 529, 533, 551, 559, 583, 589, 611, 629, 649, 667, 671, 689, 697, 703, 713, 731, 737, 767, 779, 781, 793, 799
Offset: 1

Views

Author

Daniel D Gibson, Nov 10 2024

Keywords

Examples

			The greatest primorial below 25 is 6, the factors of which are 2 and 3, neither of which are factors of 25=5^2.
		

Crossrefs

Intersection of A002808 and A335284.

Programs

  • Mathematica
    q[k_] := Module[{p = 2, r = 2}, If[CompositeQ[k], While[r < k && ! Divisible[k, p], p = NextPrime[p]; r *= p]; r >= k, False]]; Select[Range[800], q] (* Amiram Eldar, Nov 14 2024 *)
  • PARI
    P(n)=my(t=1, k); forprime(p=2, , k=t*p; if(k>n, return(t), t=k)); \\ A260188
    isok(k) = (k>1) && !isprime(k) && (gcd(k, P(k-1)) == 1); \\ Michel Marcus, Nov 12 2024
Showing 1-8 of 8 results.