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

A002808 The composite numbers: numbers n of the form x*y for x > 1 and y > 1.

Original entry on oeis.org

4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 25, 26, 27, 28, 30, 32, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 49, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 81, 82, 84, 85, 86, 87, 88
Offset: 1

Views

Author

Keywords

Comments

The natural numbers 1,2,... are divided into three sets: 1 (the unit), the primes (A000040) and the composite numbers (A002808).
The number of composite numbers <= n (A065855) = n - pi(n) (A000720) - 1.
n is composite iff sigma(n) + phi(n) > 2n. This is a nice result of the well known theorem: For all positive integers n, n = Sum_{d|n} phi(d). For the proof see my contribution to puzzle 76 of Carlos Rivera's Primepuzzles. - Farideh Firoozbakht, Jan 27 2005, Jan 18 2015
The composite numbers have the semiprimes A001358 as primitive elements.
A211110(a(n)) > 1. - Reinhard Zumkeller, Apr 02 2012
A060448(a(n)) > 1. - Reinhard Zumkeller, Apr 05 2012
A086971(a(n)) > 0. - Reinhard Zumkeller, Dec 14 2012
Composite numbers n which are the product of r=A001222(n) prime numbers are sometimes called r-almost primes. Sequences listing r-almost primes are: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), A069275 (r = 14), A069276 (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), A069280 (r = 19), A069281 (r = 20). - Jason Kimberley, Oct 02 2011
a(n) = A056608(n) * A160180(n). - Reinhard Zumkeller, Mar 29 2014
Degrees for which there are irreducible polynomials which are reducible mod p for all primes p, see Brandl. - Charles R Greathouse IV, Sep 04 2014
An integer is composite if and only if it is the sum of strictly positive integers in arithmetic progression with common difference 2: 4 = 1 + 3, 6 = 2 + 4, 8 = 3 + 5, 9 = 1 + 3 + 5, etc. - Jean-Christophe Hervé, Oct 02 2014
This statement holds since k+(k+2)+...+k+2(n-1) = n*(n+k-1) = a*b with arbitrary a,b (taking n=a and k=b-a+1 if b>=a). - M. F. Hasler, Oct 04 2014
For n > 4, these are numbers n such that n!/n^2 = (n-1)!/n is an integer (see A056653). - Derek Orr, Apr 16 2015
Let f(x) = Sum_{i=1..x} Sum_{j=2..i-1} cos((2*Pi*x*j)/i). It is known that the zeros of f(x) are the prime numbers. So these are the numbers n such that f(n) > 0. - Michel Lagneau, Oct 13 2015
Numbers n that can be written as solutions of the Diophantine equation n = (x+2)(y+2) where {x,y} in N^2, pairs of natural numbers including zero (cf. Mathematica code and Davis). - Ron R Spencer and Bradley Klee, Aug 15 2016
Numbers n with a partition (containing at least two summands) so that its summands also multiply to n. If n is prime, there is no way to find those two (or more) summands. If n is composite, simply take a factor or several, write those divisors and fill with enough 1's so that they add up to n. For example: 4 = 2*2 = 2+2, 6 = 1*2*3 = 1+2+3, 8 = 1*1*2*4 = 1+1+2+4, 9 = 1*1*1*3*3 = 1+1+1+3+3. - Juhani Heino, Aug 02 2017

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, page 2.
  • A. E. Bojarincev, Asymptotic expressions for the n-th composite number, Univ. Mat. Zap. 6:21-43 (1967). - In Russian.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See p. 127.
  • Martin Davis, "Algorithms, Equations, and Logic", pp. 4-15 of S. Barry Cooper and Andrew Hodges, Eds., "The Once and Future Turing: Computing the World", Cambridge 2016.
  • R. K. Guy, Unsolved Problems Number Theory, Section A.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers. 3rd ed., Oxford Univ. Press, 1954, p. 2.
  • D. R. Hofstadter, Goedel, Escher, Bach: an Eternal Golden Braid, Random House, 1980, p. 66.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 51.
  • 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).

Crossrefs

Complement of A008578. - Omar E. Pol, Dec 16 2016
Cf. A073783 (first differences), A073445 (second differences).
Boustrophedon transforms: A230954, A230955.
Cf. A163870 (nontrivial divisors).
Related sequences:
Primes (p) and composites (c): A000040, A002808, A000720, A065855.
Primes between p(n) and 2*p(n): A063124, A070046; between c(n) and 2*c(n): A376761; between n and 2*n: A035250, A060715, A077463, A108954.
Composites between p(n) and 2*p(n): A246514; between c(n) and 2*c(n): A376760; between n and 2*n: A075084, A307912, A307989, A376759.

Programs

  • Haskell
    a002808 n = a002808_list !! (n-1)
    a002808_list = filter ((== 1) . a066247) [2..]
    -- Reinhard Zumkeller, Feb 04 2012
    
  • Magma
    [n: n in [2..250] | not IsPrime(n)]; // G. C. Greubel, Feb 24 2024
    
  • Maple
    t := []: for n from 2 to 20000 do if isprime(n) then else t := [op(t),n]; fi; od: t; remove(isprime,[$3..89]); # Zerinvary Lajos, Mar 19 2007
    A002808 := proc(n) option remember ; local a ; if n = 1 then 4; else for a from procname(n-1)+1 do if not isprime(a) then return a; end if; end do ; end if; end proc; # R. J. Mathar, Oct 27 2009
  • Mathematica
    Select[Range[2,100], !PrimeQ[#]&] (* Zak Seidov, Mar 05 2011 *)
    With[{nn=100},Complement[Range[nn],Prime[Range[PrimePi[nn]]]]] (* Harvey P. Dale, May 01 2012 *)
    Select[Range[100], CompositeQ] (* Jean-François Alcover, Nov 07 2021 *)
  • PARI
    A002808(n)=for(k=0,primepi(n),isprime(n++)&&k--);n \\ For illustration only: see below. - M. F. Hasler, Oct 31 2008
    
  • PARI
    A002808(n)= my(k=-1); while(-n + n += -k + k=primepi(n),); n \\ For n=10^4 resp. 3*10^4, this is about 100 resp. 500 times faster than the former; M. F. Hasler, Nov 11 2009
    
  • PARI
    forcomposite(n=1, 1e2, print1(n, ", ")) \\ Felix Fröhlich, Aug 03 2014
    
  • PARI
    for(n=1, 1e3, if(bigomega(n) > 1, print1(n, ", "))) \\ Altug Alkan, Oct 14 2015
    
  • Python
    from sympy import primepi
    def A002808(n):
        m, k = n, primepi(n) + 1 + n
        while m != k:
            m, k = k, primepi(k) + 1 + n
        return m # Chai Wah Wu, Jul 15 2015, updated Apr 14 2016
    
  • Python
    from sympy import isprime
    def ok(n): return n > 1 and not isprime(n)
    print([k for k in range(89) if ok(k)]) # Michael S. Branicky, Nov 07 2021
    
  • Python
    next_A002808=lambda n: next(n for n in range(n,n*5)if not isprime(n)) # next composite >= n > 0; next_A002808(n)==n <=> iscomposite(n). - M. F. Hasler, Mar 28 2025
    is_A002808=lambda n:not isprime(n) and n>1 # where isprime(n) can be replaced with: all(n%d for d in range(2, int(n**.5)+1))
    # generators of composite numbers:
    A002808_upto=lambda stop=1<<59: filter(is_A002808, range(2,stop))
    A002808_seq=lambda:(q:=2)and(n for p in primes if (o:=q)<(q:=p) for n in range(o+1,p)) # with, e.g.: primes=filter(isprime,range(2,1<<59)) # M. F. Hasler, Mar 28 2025
    
  • SageMath
    [n for n in (2..250) if not is_prime(n)] # G. C. Greubel, Feb 24 2024

Formula

a(n) = pi(a(n)) + 1 + n, where pi is the prime counting function.
a(n) = A136527(n, n).
A000005(a(n)) > 2. - Juri-Stepan Gerasimov, Oct 17 2009
A001222(a(n)) > 1. - Juri-Stepan Gerasimov, Oct 30 2009
A000203(a(n)) < A007955(a(n)). - Juri-Stepan Gerasimov, Mar 17 2011
A066247(a(n)) = 1. - Reinhard Zumkeller, Feb 05 2012
Sum_{n>=1} 1/a(n)^s = Zeta(s)-1-P(s), where P is prime zeta. - Enrique Pérez Herrero, Aug 08 2012
n + n/log n + n/log^2 n < a(n) < n + n/log n + 3n/log^2 n for n >= 4, see Panaitopol. Bojarincev gives an asymptotic version. - Charles R Greathouse IV, Oct 23 2012
a(n) > n + A000720(n) + 1. - François Huppé, Jan 08 2025

Extensions

Deleted an incomplete and broken link. - N. J. A. Sloane, Dec 16 2010

A328028 Nonprime numbers n whose proper divisors (greater than 1 and less than n) have no consecutive divisibilities.

Original entry on oeis.org

1, 4, 6, 9, 10, 12, 14, 15, 21, 22, 24, 25, 26, 30, 33, 34, 35, 36, 38, 39, 45, 46, 48, 49, 51, 55, 57, 58, 60, 62, 63, 65, 69, 70, 72, 74, 77, 82, 84, 85, 86, 87, 90, 91, 93, 94, 95, 96, 105, 106, 108, 111, 115, 118, 119, 120, 121, 122, 123, 129, 132, 133, 134
Offset: 1

Views

Author

Gus Wiseman, Oct 06 2019

Keywords

Examples

			The proper divisors of 18 are {2, 3, 6, 9}, and {3, 6} are a consecutive divisible pair, so 18 does not belong to the sequence.
The proper divisors of 60 are {2, 3, 4, 5, 6, 10, 12, 15, 20, 30}, and none of {2, 3}, {3, 4}, {4, 5}, {5, 6}, {6, 10}, {10, 12}, {12, 15}, {15, 20}, or {20, 30} are divisible pairs, so 60 belongs to the sequence.
		

Crossrefs

Positions of 0's or 2's in A328026.
1 and positions of 1's in A328194.
The version including primes is A328161.
Partitions with no consecutive divisibilities are A328171.
Numbers whose proper divisors have no consecutive successions are A088725.
Contains A001358.

Programs

  • Maple
    filter:= proc(n) local D,i;
      if isprime(n) then return false fi;
      D:= sort(convert(numtheory:-divisors(n) minus {1,n}, list));
      for i from 1 to nops(D)-1 do if (D[i+1]/D[i])::integer then return false fi od:
      true
    end proc:
    select(filter, [$1..300]); # Robert Israel, Oct 11 2019
  • Mathematica
    Select[Range[100],!PrimeQ[#]&&!MatchQ[DeleteCases[Divisors[#],1|#],{_,x_,y_,_}/;Divisible[y,x]]&]

A106708 a(n) is the concatenation of its nontrivial divisors.

Original entry on oeis.org

0, 0, 0, 2, 0, 23, 0, 24, 3, 25, 0, 2346, 0, 27, 35, 248, 0, 2369, 0, 24510, 37, 211, 0, 2346812, 5, 213, 39, 24714, 0, 23561015, 0, 24816, 311, 217, 57, 234691218, 0, 219, 313, 24581020, 0, 23671421, 0, 241122, 35915, 223, 0, 23468121624, 7, 251025, 317
Offset: 1

Views

Author

N. J. A. Sloane, Jul 20 2007

Keywords

Crossrefs

Cf. A037278, A120712, A037279, A131983 (records), A131984 (where records occur).

Programs

  • Haskell
    a106708 1           = 0
    a106708 n
       | a010051 n == 1 = 0
       | otherwise = read $ concat $ (map show) $ init $ tail $ a027750_row n
    -- Reinhard Zumkeller, May 01 2012
    
  • Maple
    A106708 := proc(n) local dvs ; if isprime(n) or n = 1 then 0; else dvs := [op(numtheory[divisors](n) minus {1,n} )] ; dvs := sort(dvs) ; cat(op(dvs)) ; fi ; end: seq(A106708(n),n=1..80) ; # R. J. Mathar, Aug 01 2007
  • Mathematica
    Table[If[CompositeQ[n],FromDigits[Flatten[IntegerDigits/@Rest[ Most[ Divisors[ n]]]]],0],{n,60}] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Jun 22 2020 *)
  • PARI
    {map(n) = local(d); d=divisors(n); if(#d<3, 0, d[1]=""; eval(concat(vecextract(d, concat("..", #d-1)))))}
    for(n=1,51,print1(map(n),",")) /* Klaus Brockhaus, Aug 05 2007 */
    
  • Python
    from sympy import divisors
    def a(n):
      nontrivial_divisors = [d for d in divisors(n)[1:-1]]
      if len(nontrivial_divisors) == 0: return 0
      else: return int("".join(str(d) for d in nontrivial_divisors))
    print([a(n) for n in range(1, 52)]) # Michael S. Branicky, Dec 31 2020

Formula

a(n) = A037279(n) * A010051(n). - R. J. Mathar, Aug 01 2007

Extensions

More terms from R. J. Mathar and Klaus Brockhaus, Aug 01 2007
Name edited by Michael S. Branicky, Dec 31 2020

A328161 Numbers n that are prime or whose proper divisors (greater than 1 and less than n) have no consecutive divisibilities.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 17, 19, 21, 22, 23, 24, 25, 26, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 43, 45, 46, 47, 48, 49, 51, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 67, 69, 70, 71, 72, 73, 74, 77, 79, 82, 83, 84, 85, 86, 87, 89, 90, 91
Offset: 1

Views

Author

Gus Wiseman, Oct 06 2019

Keywords

Examples

			The proper divisors of 18 are {2, 3, 6, 9}, and {3, 6} are a consecutive divisible pair, so 18 does not belong to the sequence.
The proper divisors of 60 are {2, 3, 4, 5, 6, 10, 12, 15, 20, 30}, and none of {2, 3}, {3, 4}, {4, 5}, {5, 6}, {6, 10}, {10, 12}, {12, 15}, {15, 20}, or {20, 30} are divisible pairs, so 60 belongs to the sequence.
		

Crossrefs

Equals the union of A328028 and A000040.
Complement of A328189.
One, primes, and positions of 1's in A328194.
Partitions with no consecutive divisibilities are A328171.

Programs

  • Maple
    filter:= proc(n) local D,i;
      if isprime(n) then return true fi;
      D:= sort(convert(numtheory:-divisors(n) minus {1,n}, list));
      for i from 1 to nops(D)-1 do if (D[i+1]/D[i])::integer then return false fi od:
      true
    end proc:
    select(filter, [$1..100]); # Robert Israel, Oct 11 2019
  • Mathematica
    Select[Range[100],!MatchQ[DeleteCases[Divisors[#],1|#],{_,x_,y_,_}/;Divisible[y,x]]&]

A328189 Numbers n with at least one pair of consecutive divisible nontrivial divisors (greater than 1 and less than n).

Original entry on oeis.org

8, 16, 18, 20, 27, 28, 32, 40, 42, 44, 50, 52, 54, 56, 64, 66, 68, 75, 76, 78, 80, 81, 88, 92, 98, 99, 100, 102, 104, 110, 112, 114, 116, 117, 124, 125, 126, 128, 130, 136, 138, 140, 147, 148, 152, 153, 156, 160, 162, 164, 170, 171, 172, 174, 176, 184, 186
Offset: 1

Views

Author

Gus Wiseman, Oct 13 2019

Keywords

Examples

			The nontrivial divisors of 42 are {2, 3, 6, 7, 14, 21}, with pairs of consecutive divisible divisors {3, 6} and {7, 14}, so 42 belongs to the sequence.
		

Crossrefs

Complement of A328161.
Positions of terms greater than 1 in A328194.
Partitions with a pair of consecutive divisible parts are A328221.

Programs

  • Mathematica
    Select[Range[200],MatchQ[DeleteCases[Divisors[#],1|#],{_,x_,y_,_}/;Divisible[y,x]]&]
    Select[Range[2,200],AnyTrue[Partition[Most[Rest[Divisors[#]]],2,1],Mod[#[[2]],#[[1]]] == 0&]&] (* Harvey P. Dale, Mar 14 2023 *)

A328194 Maximum length of a divisibility chain of consecutive nontrivial divisors of n (greater than 1 and less than n).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 14 2019

Keywords

Comments

The nontrivial divisors of n are row n of A163870.

Examples

			The nontrivial divisors of 272 are {2, 4, 8, 16, 17, 34, 68, 136} with divisibility chains {{2, 4, 8, 16}, {17, 34, 68, 136}}, so a(272) = 4.
		

Crossrefs

Positions of 1's are A328028 without 1.
The version with all divisors allowed is A328162.
Allowing n as a divisor of n gives A328195.
Indices of terms greater than 1 are A328189.
The maximum run-length of divisors of n is A055874(n).

Programs

  • Mathematica
    Table[Switch[n,1,0,?PrimeQ,0,,Max@@Length/@Split[DeleteCases[Divisors[n],1|n],Divisible[#2,#1]&]],{n,100}]
  • PARI
    A328194(n) = if(1==n || isprime(n), 0, my(divs=divisors(n), rl=0,ml=1); for(i=2,#divs-1,if(!(divs[i]%divs[i-1]), rl++, ml = max(rl,ml); rl=1)); max(ml,rl)); \\ Antti Karttunen, Dec 07 2024

Extensions

Data section extended up to a(105) by Antti Karttunen, Dec 07 2024

A144925 Number of nontrivial divisors of the n-th composite number.

Original entry on oeis.org

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

Views

Author

Huen Yeong Kong (cosmology(AT)pacific.net.sg), Sep 25 2008

Keywords

Comments

1 and the number itself are excluded as divisors.
First occurrence of k: 1, 2, 9, 6, 45, 14, 24, 32, 851, 42, 3531, 148, 109, 89, 58993, 138, ..., which corresponds to the composite number (A005179): 4, 6, 16, 12, 64, 24, 36, 48, 1024, 60, 4096, 192, 144, 120, 65536, 180, ..., . - Robert G. Wilson v, Aug 30 2009
Row lengths of table in A163870. - Reinhard Zumkeller, Mar 29 2014

Crossrefs

Programs

  • Haskell
    a144925 = length . a163870_row  -- Reinhard Zumkeller, Mar 29 2014
  • Mathematica
    Composite[n_Integer] := FixedPoint[n + PrimePi@# + 1 &, n + PrimePi@n + 1]; f[n_] := DivisorSigma[0, n] - 2; Table[f@ Composite@ n, {n, 101}] (* Robert G. Wilson v, Aug 30 2009 *)
    DivisorSigma[0,#]-2&/@Select[Range[300],CompositeQ] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Nov 15 2018 *)
  • PARI
    k=1;vector(120,n,while(isprime(k++),0);numdiv(k)-2)
    

Formula

a(n) = A070824(A002808(n)) = A000005(A002808(n)) - 2.
A144925(n) = A070824(A002808(n)) = A000005(A002808(n)) - 2. - Robert G. Wilson v, Aug 30 2009

Extensions

Sequence extended by Juri-Stepan Gerasimov, Aug 05 2009
Edited and extended by Franklin T. Adams-Watters, Aug 30 2009

A328195 Maximum length of a divisibility chain of consecutive divisors of n greater than 1.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Oct 14 2019

Keywords

Comments

Also the maximum length of a divisibility chain of consecutive divisors of n less than n.
The divisors of n (except 1) are row n of A027749.

Examples

			The divisors of 272 greater than 1 are {2, 4, 8, 16, 17, 34, 68, 136, 272}, with divisibility chains {{2, 4, 8, 16}, {17, 34, 68, 136, 272}}, so a(272) = 5.
		

Crossrefs

Allowing 1 as a divisor gives A328162.
Forbidding n as a divisor of n gives A328194.
Positions of 1's are A000040 (primes).
Indices of terms greater than 1 are A002808 (composite numbers).
The maximum run-length of divisors of n is A055874(n).

Programs

  • Mathematica
    Table[If[n==1,0,Max@@Length/@Split[DeleteCases[Divisors[n],1],Divisible[#2,#1]&]],{n,100}]
  • PARI
    A328195(n) = if(1==n, 0, my(divs=divisors(n), rl=0,ml=1); for(i=2,#divs,if(!(divs[i]%divs[i-1]), rl++, ml = max(rl,ml); rl=1)); max(ml,rl)); \\ Antti Karttunen, Dec 07 2024

Extensions

Data section extended up to a(105) by Antti Karttunen, Dec 07 2024

A160180 Largest proper divisor of the n-th composite number.

Original entry on oeis.org

2, 3, 4, 3, 5, 6, 7, 5, 8, 9, 10, 7, 11, 12, 5, 13, 9, 14, 15, 16, 11, 17, 7, 18, 19, 13, 20, 21, 22, 15, 23, 24, 7, 25, 17, 26, 27, 11, 28, 19, 29, 30, 31, 21, 32, 13, 33, 34, 23, 35, 36, 37, 25, 38, 11, 39, 40, 27, 41, 42, 17, 43, 29, 44, 45, 13, 46, 31, 47, 19, 48, 49, 33, 50, 51, 52
Offset: 1

Views

Author

Kyle Stern, May 03 2009, May 04 2009

Keywords

Comments

Old name: The n-th positive composite number divided by its lowest nontrivial factor.

Examples

			a(1) = 4/2 = 2, a(2) = 6/2 = 3, a(3) = 8/2 = 4, a(4) = 9/3 = 3, a(5) = 10/2 = 5.
		

Crossrefs

Programs

  • Haskell
    a160180 = a032742 . a002808  -- Reinhard Zumkeller, Mar 29 2014
  • MATLAB
    function [a] = A160180(k) j = 0; n = 1; while j < k if isprime(n) == 1 skip elseif isprime(n) == 0 j = j + 1; factors = factor(n); lowfactor = factors(1,1); a(j,1) = n/lowfactor; end n = n + 1; end - Kyle Stern, May 04 2009
    
  • Mathematica
    f[n_] := Block[{k = n + PrimePi@ n + 1}, While[k != n + PrimePi@ k + 1, k++ ]; k/FactorInteger[k][[1, 1]]]; Array[f, 75] (* Robert G. Wilson v, May 11 2012 *)
    Divisors[#][[-2]]&/@Select[Range[200],CompositeQ] (* Harvey P. Dale, Dec 06 2021 *)
    (# / FactorInteger[#][[1, 1]])& /@ Select[Range[300], CompositeQ] (* Amiram Eldar, Jun 18 2022 *)

Formula

a(n) = A032742(A002808(n)) = A002808(n) / A056608(n) = A163870(n,A144925(n)). - Reinhard Zumkeller, Mar 29 2014

Extensions

Indices of b-file corrected, more terms added using b-file. - N. J. A. Sloane, Aug 31 2009
New name from Reinhard Zumkeller, Mar 29 2014
Incorrect formula removed by Ridouane Oudra, Oct 15 2021

A328458 Maximum run-length of the nontrivial divisors (greater than 1 and less than n) of n.

Original entry on oeis.org

1, 0, 0, 1, 0, 2, 0, 1, 1, 1, 0, 3, 0, 1, 1, 1, 0, 2, 0, 2, 1, 1, 0, 3, 1, 1, 1, 1, 0, 2, 0, 1, 1, 1, 1, 3, 0, 1, 1, 2, 0, 2, 0, 1, 1, 1, 0, 3, 1, 1, 1, 1, 0, 2, 1, 2, 1, 1, 0, 5, 0, 1, 1, 1, 1, 2, 0, 1, 1, 1, 0, 3, 0, 1, 1, 1, 1, 2, 0, 2, 1, 1, 0, 3, 1, 1, 1, 1, 0, 2, 1, 1, 1, 1, 1, 3, 0, 1, 1, 2, 0, 2, 0, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Oct 17 2019

Keywords

Comments

By convention, a(1) = 1, and a(p) = 0 for p prime.

Examples

			The non-singleton runs of the nontrivial divisors of 1260 are: {2,3,4,5,6,7} {9,10} {14,15} {20,21} {35,36}, so a(1260) = 6.
		

Crossrefs

Positions of first appearances are A328459.
Positions of 0's and 1's are A088723.
The version that looks at all divisors is A055874.
The number of successive pairs of divisors > 1 of n is A088722(n).
The Heinz number of the multiset of run-lengths of divisors of n is A328166(n).

Programs

  • Mathematica
    Table[Switch[n,1,1,?PrimeQ,0,,Max@@Length/@Split[DeleteCases[Divisors[n],1|n],#2==#1+1&]],{n,100}]
  • PARI
    A328458(n) = if(1==n,n,my(rl=0,pd=0,m=0); fordiv(n, d, if(1(1+pd), m = max(m,rl); rl=0); pd=d; rl++)); max(m,rl)); \\ Antti Karttunen, Feb 23 2023

Extensions

Data section extended up to a(105) by Antti Karttunen, Feb 23 2023
Showing 1-10 of 15 results. Next