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 16 results. Next

A048670 Jacobsthal function A048669 applied to the product of the first n primes (A002110).

Original entry on oeis.org

2, 4, 6, 10, 14, 22, 26, 34, 40, 46, 58, 66, 74, 90, 100, 106, 118, 132, 152, 174, 190, 200, 216, 234, 258, 264, 282, 300, 312, 330, 354, 378, 388, 414, 432, 450, 476, 492, 510, 538, 550, 574, 600, 616, 642, 660, 686, 718, 742, 762, 798, 810, 834, 858, 876, 908, 926, 954
Offset: 1

Views

Author

Keywords

Comments

Pintz shows that j(x#) >= (2*e^gamma + o(1)) x log x log log log x / (log log x)^2 and hence a(n) >= (2*e^gamma + o(1)) n log^2 n log log log n / (log log n)^2 by the Prime Number Theorem. - Charles R Greathouse IV, Sep 08 2012
Jacobsthal conjectures that a(n) >= j(k) := A048669(k) for any k with n prime factors, which would make this the RECORDS transform of A048669. Hajdu & Saradha disprove the conjecture, showing that this fails for n = 24 where j(k) = 236 > 234 = a(24) for any k divisible by 76964283982898776138308824190 and with 24 prime factors in total. - Charles R Greathouse IV, Sep 08 2012 / Edited by Jan Kristian Haugland, Feb 02 2019
Ford, Green, Konyagin, Maynard, & Tao show that j(x#) >> x log x log log log x / log log x and hence a(n) >> n log^2 n log log log n / log log n. - Charles R Greathouse IV, Mar 29 2018
Computation of a(62)-a(64) was supported by Google Cloud. - Andrzej Bozek, Mar 14 2021

References

  • L. E. Dickson, History of the Theory of Numbers, Vol. 1, p. 439, Chelsea, 1952.

Crossrefs

Programs

  • Mathematica
    (* This program is not suitable to compute more than a few terms *) primorial[n_] := Product[Prime[k], {k, 1, n}]; j[n_] := Module[{L = 1, m = 1}, For[k = 2, k <= n + 1, k++, If[GCD[k, n] == 1, If[L + m < k, m = k - L]; L = k]]; m]; a[n_] := a[n] = j[primorial[n]]; Table[Print["a(", n, ") = ", a[n]]; a[n], {n, 1, 10}] (* Jean-François Alcover, Sep 27 2013, after M. F. Hasler *)

Formula

a(n) = A058989(n) + 1.
a(n) << n^2*(log n)^2, see Iwaniec. - Charles R Greathouse IV, Sep 08 2012
a(n) >= (2*e^gamma + o(1)) n log^2 n log log log n / (log log n)^2, see Pintz.
a(n) = 2 * A072752(n) + 2. - Mario Ziller, Dec 08 2016
Maier & Pomerance conjecture that Max_{n <= x} A048669(n) = log(x)*(log log x)^(2+o(1)) which suggests a(n) = n*(log n)^(3+o(1)). - Charles R Greathouse IV, Mar 29 2018
a(n) = largest (or last) term in row n of A331118. - Michael De Vlieger, Dec 11 2020

Extensions

a(21)-a(24) from Max Alekseyev, Apr 09 2006
a(25)-a(49) from Thomas Hagedorn, Feb 21 2007
a(46) corrected (published value in Hagedorn's 2009 Mathematics of Computation article was correct) and a(50)-a(54) added by Mario Ziller, Dec 08 2016
a(55)-a(57) from Robert Gerbicz, Apr 10 2017
a(58)-a(64) from Andrzej Bozek, Mar 14 2021

A128708 Let j(n) be the Jacobsthal function, A048669. Then a(n) is the number of times that the gap j(n) appears between consecutive numbers <= n+1 and coprime to n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 4, 3, 1, 1, 2, 1, 1, 2, 8, 1, 3, 1, 2, 2, 1, 1, 4, 5, 1, 9, 2, 1, 2, 1, 16, 2, 1, 2, 6, 1, 1, 2, 4, 1, 2, 1, 2, 6, 1, 1, 8, 7, 5, 2, 2, 1, 9, 2, 4, 2, 1, 1, 4, 1, 1, 6, 32, 2, 2, 1, 2, 2, 2, 1, 12, 1, 1, 10, 2, 2, 2, 1, 8, 27, 1, 1, 4, 2, 1, 2, 4, 1, 6, 2, 2, 2, 1, 2, 16, 1, 7, 6, 10
Offset: 1

Views

Author

T. D. Noe, Mar 24 2007

Keywords

Comments

Differs from A087653 starting at n=35. For prime n and e>0, a(n^e)=n^(e-1). The closely-related sequence A128707 satisfies the inequality a(n)*A128707(n) <= n-1, with equality for prime n. If m is the squarefree kernel of n (A007947), then a(n)/a(m) = n/m.

Examples

			The numbers coprime to 15 are 1,2,4,7,8,11,13,14,16,17,19,22,... Observe that the differences are periodic: 1,2,3,1,3,2,1,2,1,2,3,... The maximum value is 3, which occurs twice in the first period. Hence a(15)=2.
		

Programs

  • Mathematica
    JacobsthalCount[n_] := Module[{g,d,mx}, g=Select[Range[n+1], GCD[n,# ]==1&]; d=Rest[g]-Most[g]; mx=Max@@d; Count[d,mx]]; Table[JacobsthalCount[n], {n,100}]

A128759 Least k such that the Jacobsthal function A048669(k) = n.

Original entry on oeis.org

1, 2, 15, 6, 105, 30, 1155, 770, 36465, 210, 15015, 6006, 255255, 2310, 8580495, 102102, 4849845, 72930, 20056049013, 74364290, 5898837945, 30030, 3234846615, 881790, 195282582495, 510510, 218257003965, 20281170, 100280245065, 17160990, 934482952262145, 6614136163635
Offset: 1

Views

Author

T. D. Noe, Mar 28 2007

Keywords

Comments

Except for a(1) and a(2), the same as A070971. It appears that a(2n)=2a(n) for odd n. Because the primorial numbers (A002110) yield record values of the Jacobsthal function, we have a(A048670(n))=A002110(n). Note that numbers in this sequence up to n=18 have the form p#, p#/2, p#/q, or p#/(2q), where p and q are primes with 2

Crossrefs

Extensions

a(19) onwards from Don Reble, Oct 17 2013

A285182 L = A132468(n) = A048669(n)-1 is the length of the longest run of consecutive numbers that have a common factor with n; a(n) = smallest k >= 0 which starts such a run.

Original entry on oeis.org

0, 0, 0, 0, 2, 0, 0, 0, 4, 0, 2, 0, 6, 5, 0, 0, 2, 0, 4, 6, 10, 0, 2, 0, 12, 0, 6, 0, 2, 0, 0, 11, 16, 14, 2, 0, 18, 12, 4, 0, 6, 0, 10, 5, 22, 0, 2, 0, 4, 17, 12, 0, 2, 10, 6, 18, 28, 0, 2, 0, 30, 6, 0, 25, 8, 0, 16, 23, 4, 0, 2, 0, 36, 5, 18, 21, 12, 0, 4, 0, 40, 0, 6, 34, 42, 29, 10, 0, 2, 13, 22, 30
Offset: 2

Author

N. J. A. Sloane, Apr 19 2017

Keywords

Comments

Related to the Jacobsthal function A048669.
a(n) depends only on the radical A007947(n).

Examples

			If n = 6, among the numbers k = 0,1,2,3,4,5,6,7,... there is a run of L = 3 consecutive numbers, 2,3,4, all with gcd(k,6)>1, starting at k=2, so a(6) = 2.
If n is a prime (or prime power), a(n)=0.
		

Crossrefs

Programs

  • Maple
    acf:=[]; wcf:=[];
    for n from 2 to 120 do
    s:=[seq(j,j=0..4*n)];
    rec:=0;
       for st from 0 to n do
       len:=0;
          for i from 1 to n while gcd(s[st+i],n)>1 do len:=len+1; od:
          if len>rec then rec:=len; w:=st; fi;
       od:
    acf:=[op(acf),rec];
    wcf:=[op(wcf),w];
    od:
    acf; # A132468
    wcf; # A285182

A040001 1 followed by {1, 2} repeated.

Original entry on oeis.org

1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2
Offset: 0

Comments

Continued fraction for sqrt(3).
Also coefficient of the highest power of q in the expansion of the polynomial nu(n) defined by: nu(0)=1, nu(1)=b and for n>=2, nu(n)=b*nu(n-1)+lambda*(n-1)_q*nu(n-2) with (b,lambda)=(1,1), where (n)_q=(1+q+...+q^(n-1)) and q is a root of unity. - Y. Kelly Itakura (yitkr(AT)mta.ca), Aug 21 2002
nu(0)=1 nu(1)=1; nu(2)=2; nu(3)=3+q; nu(4)=5+3q+2q^2; nu(5)=8+7q+6q^2+4q^3+q^4; nu(6)=13+15q+16q^2+14q^3+11q^4+5q^5+2q^6.
From Jaroslav Krizek, May 28 2010: (Start)
a(n) = denominators of arithmetic means of the first n positive integers for n >= 1.
See A026741(n+1) or A145051(n) - denominators of arithmetic means of the first n positive integers. (End)
From R. J. Mathar, Feb 16 2011: (Start)
This is a prototype of multiplicative sequences defined by a(p^e)=1 for odd primes p, and a(2^e)=c with some constant c, here c=2. They have Dirichlet generating functions (1+(c-1)/2^s)*zeta(s).
Examples are A153284, A176040 (c=3), A040005 (c=4), A021070, A176260 (c=5), A040011, A176355 (c=6), A176415 (c=7), A040019, A021059 (c=8), A040029 (c=10), A040041 (c=12). (End)
a(n) = p(-1) where p(x) is the unique degree-n polynomial such that p(k) = A000325(k) for k = 0, 1, ..., n. - Michael Somos, May 12 2012
For n > 0: denominators of row sums of the triangular enumeration of rational numbers A226314(n,k) / A054531(n,k), 1 <= k <= n; see A226555 for numerators. - Reinhard Zumkeller, Jun 10 2013
From Jianing Song, Nov 01 2022: (Start)
For n > 0, a(n) is the minimal gap of distinct numbers coprime to n. Proof: denote the minimal gap by b(n). For odd n we have A058026(n) > 0, hence b(n) = 1. For even n, since 1 and -1 are both coprime to n we have b(n) <= 2, and that b(n) >= 2 is obvious.
The maximal gap is given by A048669. (End)

Examples

			1.732050807568877293527446341... = 1 + 1/(1 + 1/(2 + 1/(1 + 1/(2 + ...))))
G.f. = 1 + x + 2*x^2 + x^3 + 2*x^4 + x^5 + 2*x^6 + x^7 + 2*x^8 + x^9 + ...
		

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 186.
  • Jan Gullberg, Mathematics from the Birth of Numbers, W. W. Norton & Co., NY & London, 1997, §4.4 Powers and Roots, p. 144.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 276.

Crossrefs

Cf. A000034, A002194, A133566, A083329 (binomial Transf).
Apart from a(0) the same as A134451.

Programs

  • Haskell
    a040001 0 = 1; a040001 n = 2 - mod n 2
    a040001_list = 1 : cycle [1, 2]  -- Reinhard Zumkeller, Apr 16 2015
  • Maple
    Digits := 100: convert(evalf(sqrt(N)),confrac,90,'cvgts'):
  • Mathematica
    ContinuedFraction[Sqrt[3],300] (* Vladimir Joseph Stephan Orlovsky, Mar 04 2011 *)
    PadRight[{1},120,{2,1}] (* Harvey P. Dale, Nov 26 2015 *)
  • PARI
    {a(n) = 2 - (n==0) - (n%2)} /* Michael Somos, Jun 11 2003 */
    
  • PARI
    { allocatemem(932245000); default(realprecision, 12000); x=contfrac(sqrt(3)); for (n=0, 20000, write("b040001.txt", n, " ", x[n+1])); } \\ Harry J. Smith, Jun 01 2009
    

Formula

Multiplicative with a(p^e) = 2 if p even; 1 if p odd. - David W. Wilson, Aug 01 2001
G.f.: (1 + x + x^2) / (1 - x^2). E.g.f.: (3*exp(x)-2*exp(0)+exp(-x))/2. - Paul Barry, Apr 27 2003
a(n) = (3-2*0^n +(-1)^n)/2. a(-n)=a(n). a(2n+1)=1, a(2n)=2, n nonzero.
a(n) = sum{k=0..n, F(n-k+1)*(-2+(1+(-1)^k)/2+C(2, k)+0^k)}. - Paul Barry, Jun 22 2007
Row sums of triangle A133566. - Gary W. Adamson, Sep 16 2007
Euler transform of length 3 sequence [ 1, 1, -1]. - Michael Somos, Aug 04 2009
Moebius transform is length 2 sequence [ 1, 1]. - Michael Somos, Aug 04 2009
a(n) = sign(n) + ((n+1) mod 2) = 1 + sign(n) - (n mod 2). - Wesley Ivan Hurt, Dec 13 2013

A134451 Ternary digital root of n.

Original entry on oeis.org

0, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2
Offset: 0

Author

Reinhard Zumkeller, Oct 27 2007

Keywords

Comments

Continued fraction expansion of sqrt(3) - 1. - N. J. A. Sloane, Dec 17 2007. Cf. A040001, A048878/A002530.
Minimum number of terms required to express n as a sum of odd numbers.
Shadow transform of even numbers A005843. - Michel Marcus, Jun 06 2013
From Jianing Song, Nov 01 2022: (Start)
For n > 0, a(n) is the minimal gap of distinct numbers coprime to n. Proof: denote the minimal gap by b(n). For odd n we have A058026(n) > 0, hence b(n) = 1. For even n, since 1 and -1 are both coprime to n we have b(n) <= 2, and that b(n) >= 2 is obvious.
The maximal gap is given by A048669. (End)

Examples

			n=42: A007089(42) = '1120', A053735(42) = 1+1+2+0 = 4,
A007089(4)='11', A053735(4)=1+1=2: therefore a(42) = 2.
0.732050807568877293527446341... = 0 + 1/(1 + 1/(2 + 1/(1 + 1/(2 + ...)))). - _Harry J. Smith_, May 31 2009
		

Crossrefs

Cf. A000010, A055034, A134452, A160390 (decimal expansion).
Apart from a(0) the same as A040001.
Related base-3 sequences: A053735, A134451, A230641, A230642, A230643, A230853, A230854, A230855, A230856, A230639, A230640, A010063 (trajectory of 1).

Programs

Formula

a(n) = n if n <= 2, otherwise a(A053735(n)).
a(A005408(n)) = 1; a(A005843(n)) = 2 for n>0;
a(n) = 0 if n=0, otherwise A000034(n-1).
a(n) = ((n+1) mod 2) + 2*sign(n) - 1. - Wesley Ivan Hurt, Dec 06 2013
Multiplicative with a(2^e) = 2, a(p^e) = 1 for odd prime p. - Andrew Howroyd, Aug 06 2018
a(0) = A055034(1) / A000010(1), a(n) = A000010(n+1) / A055034(n+1), n>1. - Torlach Rush, Oct 29 2019
Dirichlet g.f.: zeta(s)*(1+1/2^s). - Amiram Eldar, Jan 01 2023

A058989 Largest number of consecutive integers such that each is divisible by a prime <= the n-th prime.

Original entry on oeis.org

1, 3, 5, 9, 13, 21, 25, 33, 39, 45, 57, 65, 73, 89, 99, 105, 117, 131, 151, 173, 189, 199, 215, 233, 257, 263, 281, 299, 311, 329, 353, 377, 387, 413, 431, 449, 475, 491, 509, 537, 549, 573, 599, 615, 641, 659, 685, 717, 741, 761, 797, 809, 833, 857, 875, 907, 925, 953, 977, 1001, 1029, 1057, 1097, 1109
Offset: 1

Author

Jud McCranie, Jan 16 2001

Keywords

Comments

Marty Weissman conjectured that a(n)=2q-1, where q is the largest prime smaller than the n-th prime. The conjecture holds for the first few terms, but then a(n) is larger than 2q-1. Phil Carmody proved a(n)>=2q-1. Terms were calculated by Weissman, Carmody and McCranie.
A049300(n) is the smallest value of the mentioned consecutive integers. - Reinhard Zumkeller, Jun 14 2003
By Lemma 5 of Maynard, there is a constant C > 0 such that, for any n, there is a prime p <= C*exp(prime(n)) such that q - p >= a(n) where q is the smallest prime larger than p. (Probably C = 7/e^3 = 0.35... is admissible.) - Charles R Greathouse IV, Aug 01 2024

Examples

			The 4th prime is 7. Nine is the maximum number of consecutive integers such that each is divisible by 2, 3, 5 or 7. (Example: 2 through 10) So a(4)=9.
		

References

  • Dickson, L. E., History of the Theory of Numbers, Vol. 1, p. 439, Chelsea, 1952.

Crossrefs

Programs

  • Mathematica
    (* This program is not suitable to compute more than a few terms *)
    primorial[n_] := Product[Prime[k], {k, 1, n}];
    j[n_] := Module[{L = 1, m = 1}, For[k = 2, k <= n+1, k++, If[GCD[k, n] == 1, If[L+m < k, m = k-L]; L = k]]; m];
    a[1] = 1;
    a[n_] := a[n] = j[primorial[n]] - 1;
    Table[Print["a(", n, ") = ", a[n]]; a[n], {n, 1, 10}] (* Jean-François Alcover, Sep 05 2017 *)
  • PARI
    do(n,P,R)=for(i=1,#R, if(n==R[i], return(do(n+1,P,R)))); if(#P==0, return(n-1)); my(b=0); for(i=1,#P,my(t=do(n+1,setminus(P,[P[i]]), concat(R,Mod(n,P[i])))); if(t>b,b=t)); b
    a(n)=do(1,primes(n),[]) \\ Charles R Greathouse IV, Aug 08 2024

Formula

a(n) = A048670(n) - 1. See that entry for additional information.
Iwaniec proved that a(n) << n^2*(log n)^2. - Charles R Greathouse IV, Sep 08 2012
a(n) >= (2e^gamma + o(1)) n log^2 n log log log n / (log log n)^2, see A048670. - Charles R Greathouse IV, Sep 08 2012
a(n) = 2 * A072752(n) + 1. - Mario Ziller, Dec 08 2016
See A048669 for many other bounds and references. - N. J. A. Sloane, Apr 19 2017

Extensions

Laison and Schick reference from Parthasarathy Nambi, Oct 19 2007
More terms from A048670 added by Max Alekseyev, Feb 07 2008
a(46) corrected and a(50)-a(54) added by Mario Ziller, Dec 08 2016
a(55)-a(64) from A048670 added by Constantino Calancha, Aug 05 2023

A049300 Smallest number starting a longest interval of consecutive integers, each of which is divisible by at least one of the first n primes.

Original entry on oeis.org

2, 2, 2, 2, 114, 9440, 217128, 60044, 20332472, 417086648, 74959204292, 187219155594, 79622514581574, 14478292443584, 6002108856728918, 12288083384384462, 5814429911995661690, 14719192159220252523420
Offset: 1

Author

Keywords

Comments

The length of such longest interval of consecutive integers is given by A058989(n), which is the first maximal gaps A048670(n) minus 1 in the reduced residue system of consecutive primorial numbers.
Let j(m)=A048669(m) be the Jacobsthal function, i.e., the maximal distance between integers relatively prime to m. Let m=2*3*5*...*prime(n). Then a(n) is the least k>0 such that k,k+1,k+2,...,k+j(m)-2 are not coprime to m. Note that a(n) begins (or is inside) a large gap between primes. - T. D. Noe, Mar 29 2007

Examples

			Between 1 and 7, all 5 numbers (2,3,4,5,6) are divisible either by 2,3 or 5. Thus a(3)=2, the initial term. Between 113 and 127 the 13 consecutive integers are divisible by 2,5,2,3,2,7,2,11,2,3,2,5,2, each from {2,3,5,7,11}. Thus a(5)=114, the smallest with this property.
		

Crossrefs

Formula

a(n) = 1 + A128707(A002110(n)). - T. D. Noe, Mar 29 2007

Extensions

More terms from T. D. Noe, Mar 29 2007
a(11)-a(12) from Donovan Johnson, Oct 13 2009
a(13) from Donovan Johnson, Oct 20 2009
a(14) and beyond from Max Alekseyev, Nov 14 2009

A070194 List the phi(n) numbers from 1 to n-1 which are relatively prime to n; sequence gives size of maximal gap.

Original entry on oeis.org

1, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 3, 2, 1, 4, 1, 4, 3, 4, 1, 4, 2, 4, 2, 4, 1, 6, 1, 2, 3, 4, 3, 4, 1, 4, 3, 4, 1, 6, 1, 4, 3, 4, 1, 4, 2, 4, 3, 4, 1, 4, 3, 4, 3, 4, 1, 6, 1, 4, 3, 2, 3, 6, 1, 4, 3, 6, 1, 4, 1, 4, 3, 4, 3, 6, 1, 4, 2, 4, 1, 6, 3, 4, 3, 4, 1, 6, 3, 4, 3, 4, 3, 4, 1, 4, 3, 4, 1, 6, 1, 4, 5, 4, 1
Offset: 3

Author

N. J. A. Sloane, May 13 2002

Keywords

Comments

Maximal gap in reduced residue system mod n.
It is an unsolved problem to determine the rate of growth of this sequence.

Examples

			For n = 10 the reduced residues are 1, 3, 7, 9; the maximal gap is a(10) = 7-3 = 4.
		

References

  • H. L. Montgomery, Ten Lectures on the Interface Between Analytic Number Theory and Harmonic Analysis, Amer. Math. Soc., 1996, p. 200.

Crossrefs

Cf. A000010.
Cf. A038566.

Programs

  • Haskell
    a070194 n = maximum $ zipWith (-) (tail ts) ts where ts = a038566_row n
    -- Reinhard Zumkeller, Oct 01 2012
  • Mathematica
    f[n_] := Block[{a = Select[ Table[i, {i, n - 1}], GCD[ #, n] == 1 & ], b = {}, k = 1, l = EulerPhi[n]}, While[k < l, b = Append[b, Abs[a[[k]] - a[[k + 1]]]]; k++ ]; Max[b]]; Table[ f[n], {n, 3, 100}]
  • PARI
    A070194(n)={my(L=1,m=1);for(k=2,n-1,gcd(k,n)>1&next;L+mM. F. Hasler, Sep 08 2012
    

Formula

a(n) = max(A048669(n),2) for all n>2. Indeed A048669 is obtained when going up to n+1 instead of only n-1 (because after n+1, the gaps among numbers coprime to n repeat). Since n-1 and n+1 are both coprime to n, this means that A048669(n)=2 whenever a(n)=1, but in all other cases, there is equality. - M. F. Hasler, Sep 08 2012

Extensions

More terms from Robert G. Wilson v and John W. Layman, May 13 2002

A070971 a(n) is the smallest positive integer m for which A070194(m) (i.e., the maximal gap in {k|gcd(k,m) = 1, 1 <= k <= m-1}) is n.

Original entry on oeis.org

3, 4, 15, 6, 105, 30, 1155, 770, 36465, 210, 15015, 6006, 255255, 2310, 8580495, 102102, 4849845, 72930, 20056049013, 74364290, 5898837945, 30030, 3234846615, 881790, 195282582495, 510510, 218257003965, 20281170, 100280245065, 17160990, 934482952262145, 6614136163635
Offset: 1

Author

John W. Layman, May 17 2002

Keywords

Comments

a(n) is the least x such that maximal gap in RRS of x equals n: a(n) = max{x: A070194(x) = n}
For n > 2, same as A128759, which gives the least k such that the Jacobsthal function A048669(k) equals n. See A128759 for more comments. - T. D. Noe, Mar 28 2007

Examples

			A070194 begins with 1,2,1,4,... with offset 3, so a(4)=6.
		

Crossrefs

Programs

  • Mathematica
    gw[x_] := Table[GCD[w, x], {w, 1, x}] rrs[x_] := Flatten[Position[gw[x], 1]] dr[x_] := Delete[RotateLeft[rrs[x]]-rrs[x], -1] t=Table[0, {25}]; Do[s=Max[dr[n]]; If[s<26&&t[[s]]==0, t[[s]]=n], {n, 3, 10000}]; t (* Labos Elemer, Oct 09 2002 *)

Extensions

a(13)-a(18) from T. D. Noe, Mar 28 2007
a(19) onwards from Don Reble, Oct 17 2013
Showing 1-10 of 16 results. Next