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

A002982 Numbers k such that k! - 1 is prime.

Original entry on oeis.org

3, 4, 6, 7, 12, 14, 30, 32, 33, 38, 94, 166, 324, 379, 469, 546, 974, 1963, 3507, 3610, 6917, 21480, 34790, 94550, 103040, 147855, 208003
Offset: 1

Views

Author

Keywords

Comments

The corresponding primes n!-1 are often called factorial primes.

Examples

			From _Gus Wiseman_, Jul 04 2019: (Start)
The sequence of numbers n! - 1 together with their prime indices begins:
                    1: {}
                    5: {3}
                   23: {9}
                  119: {4,7}
                  719: {128}
                 5039: {675}
                40319: {9,273}
               362879: {5,5,430}
              3628799: {10,11746}
             39916799: {6,7,9,992}
            479001599: {25306287}
           6227020799: {270,256263}
          87178291199: {3610490805}
        1307674367999: {7,11,11,16,114905}
       20922789887999: {436,318519035}
      355687428095999: {8,21,10165484947}
     6402373705727999: {17,20157,25293727}
   121645100408831999: {119,175195,4567455}
  2432902008176639999: {11715,659539127675}
(End)
		

References

  • J.-M. De Koninck, Ces nombres qui nous fascinent, Entry 166, p. 53, Ellipses, Paris 2008.
  • R. K. Guy, Unsolved Problems in Number Theory, Section A2.
  • 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 118.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See entry 719 at p. 160.

Crossrefs

Cf. A002981 (numbers n such that n!+1 is prime).
Cf. A055490 (primes of form n!-1).
Cf. A088332 (primes of form n!+1).

Programs

Extensions

21480 sent in by Ken Davis (ken.davis(AT)softwareag.com), Oct 29 2001
Updated Feb 26 2007 by Max Alekseyev, based on progress reported in the Carmody web site.
Inserted missing 21480 and 34790 (see Caldwell). Added 94550, discovered Oct 05 2010. Eric W. Weisstein, Oct 06 2010
103040 was discovered by James Winskill, Dec 14 2010. It has 471794 digits. Corrected by Jens Kruse Andersen, Mar 22 2011
a(26) = 147855 from Felix Fröhlich, Sep 02 2013
a(27) = 208003 from Sou Fukui, Jul 27 2016

A225485 Number of partitions of n that have frequency depth k, an array read by rows.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 2, 3, 1, 3, 4, 3, 1, 1, 4, 8, 1, 1, 3, 6, 9, 3, 1, 2, 8, 12, 7, 1, 3, 11, 17, 10, 1, 1, 11, 26, 17, 1, 5, 19, 25, 27, 1, 1, 17, 44, 38, 1, 3, 25, 53, 52, 1, 1, 3, 29, 63, 76, 4
Offset: 1

Views

Author

Clark Kimberling, May 08 2013

Keywords

Comments

Let S = {x(1),...,x(k)} be a multiset whose distinct elements are y(1),...,y(h). Let f(i) be the frequency of y(i) in S. Define F(S) = {f(1),..,f(h)}, F(1,S) = F(S), and F(m,S) = F(F(m-1),S) for m>1. Then lim(F(m,S)) = {1} for every S, so that there is a least positive integer i for which F(i,S) = {1}, which we call the frequency depth of S.
Equivalently, the frequency depth of an integer partition is the number of times one must take the multiset of multiplicities to reach (1). For example, the partition (32211) has frequency depth 5 because we have: (32211) -> (221) -> (21) -> (11) -> (2) -> (1). - Gus Wiseman, Apr 19 2019
From Clark Kimberling, Sep 26 2023: (Start)
Below, m^n abbreviates the sum m+...+m of n terms. In the following list, the numbers p_1,...,p_k are distinct, m >= 1, and k >= 1. The forms of the partitions being counted are as follows:
column 1: [n],
column 2: [m^k],
column 3: [p_1^m,...,p_k^m],
column 4: [(p_1^m_1)^m,..., (p_k^m_k)^m], distinct numbers m_i.
Column 3 is of special interest. Assume first that m = 1, so that the form of partition being counted is p = [p_1,...,p_k], with conjugate given by [q_1,...,q_m] where q_i is the number of parts of p that are >= i. Since the p_i are distinct, the distinct parts of q are the integers 1,2,...,k. For the general case that m >= 1, the distinct parts of q are the integers m,...,km. Let S(n) denote the set of partitions of n counted by column 3. Then if a and b are in the set S*(n) of conjugates of partitions in S(n), and if a > b, then a - b is also in S*(n). Call this the subtraction property. Conversely, if a partition q has the subtraction property, then q must consist of a set of numbers m,..,km for some m. Thus, column 3 counts the partitions of n that have the subtraction property. (End)

Examples

			The first 9 rows:
  n = 1 .... 0
  n = 2 .... 1..1
  n = 3 .... 1..1..1
  n = 4 .... 1..2..1..1
  n = 5 .... 1..1..2..3
  n = 6 .... 1..3..4..3
  n = 7 .... 1..1..4..8..1
  n = 8 .... 1..3..6..9..3
  n = 9 .... 1..2..8.12..7
For the 7 partitions of 5, successive frequencies are shown here:
  5 -> 1 (depth 1)
  41 -> 11 -> 2 -> 1 (depth 3)
  32 -> 11 -> 2 -> 1 (depth 3)
  311 -> 12 -> 11 -> 2 -> 1 (depth 4)
  221 -> 12 -> 11 -> 2 -> 1 (depth 4)
  2111 -> 13 -> 11 -> 2 -> 1 (depth 4)
  11111 -> 5 -> 1 (depth 2)
Summary: 1 partition has depth 1; 1 has depth 2; 2 have 3; and 3 have 4, so that the row for n = 5 is 1..1..2..3 .
		

Crossrefs

Row sums are A000041.
Column k = 2 is A032741.
Column k = 3 is A325245.
a(n!) = A325272(n).
Integer partition triangles: A008284 (first omega), A116608 (second omega), A325242 (third omega), A325268 (second-to-last omega), A225485 or A325280 (length/frequency depth).

Programs

  • Mathematica
    c[s_] := c[s] = Select[Table[Count[s, i], {i, 1, Max[s]}], # > 0 &]
    f[s_] := f[s] = Drop[FixedPointList[c, s], -2]
    t[s_] := t[s] = Length[f[s]]
    u[n_] := u[n] = Table[t[Part[IntegerPartitions[n], i]],
      {i, 1, Length[IntegerPartitions[n]]}];
    Flatten[Table[Count[u[n], k], {n, 2, 25}, {k, 1, Max[u[n]]}]]

A071626 Number of distinct exponents in the prime factorization of n!.

Original entry on oeis.org

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

Views

Author

Labos Elemer, May 29 2002

Keywords

Comments

Erdős proved that there exist two constants c1, c2 > 0 such that c1 (n / log(n))^(1/2) < a(n) < c2 (n / log(n))^(1/2). - Carlo Sanna, May 28 2019
R. Heyman and R. Miraj proved that the cardinality of the set { floor(n/p) : p <= n, p prime } is same as the number of distinct exponents in the prime factorization of n!. - Md Rahil Miraj, Apr 05 2024

Examples

			n=7: 7! = 5040 = 2*2*2*2*3*3*5*7; three different exponents arise: 4, 2 and 1; a(7)=3.
n=7: { floor(7/p) : p <= 7, p prime } = {3,2,1}. So, its cardinality is 3. - _Md Rahil Miraj_, Apr 05 2024
		

Crossrefs

Programs

  • Mathematica
    ffi[x_] := Flatten[FactorInteger[x]] lf[x_] := Length[FactorInteger[x]] ep[x_] := Table[Part[ffi[x], 2*w], {w, 1, lf[x]}] Table[Length[Union[ep[w! ]]], {w, 1, 100}]
    Table[Length[Union[Last/@If[n==1,{},FactorInteger[n!]]]],{n,30}] (* Gus Wiseman, May 15 2019 *)
  • PARI
    a(n) = #Set(factor(n!)[, 2]); \\ Michel Marcus, Sep 05 2017

Formula

a(n) = A071625(n!) = A323023(n!,3). - Gus Wiseman, May 15 2019

A325273 Prime omicron of n!.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 18 2019

Keywords

Comments

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1).
The prime omicron of n (A304465) is 0 if n is 1, 1 if n is prime, and otherwise the second-to-last part of the omega-sequence of n. For example, the prime omicron of 180 is 2.
Conjecture: all terms after a(10) = 4 are less than 4.
From James Rayman, Apr 17 2021: (Start)
The conjecture is false. a(3804) = 4. In fact, there are 91 values of n < 10000 such that a(n) = 4.
The first value of n such that a(n) = 5 is 37934. For any other n < 5*10^5, a(n) < 5. (End)

Crossrefs

a(n) = A055396(A325275(n)/2).
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number), A325249 (sum).

Programs

  • Mathematica
    omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
    omicron[n_]:=Switch[n,1,0,?PrimeQ,1,,omseq[n][[-2]]];
    Table[omicron[n!],{n,0,100}]
  • Python
    from sympy.ntheory import *
    def red(v):
        r = {}
        for i in v: r[i] = r.get(i, 0) + 1
        return r
    def omicron(v):
        if len(v) == 0: return 0
        if len(v) == 1: return v[0]
        else: return omicron(list(red(v).values()))
    f, a_list = {}, []
    for i in range(101):
        a_list.append(omicron(list(f.values())))
        g = factorint(i+1)
        for k in g: f[k] = f.get(k, 0) + g[k]
    print(a_list) # James Rayman, Apr 17 2021

Extensions

More terms from James Rayman, Apr 17 2021

A115627 Irregular triangle read by rows: T(n,k) = multiplicity of prime(k) as a divisor of n!.

Original entry on oeis.org

1, 1, 1, 3, 1, 3, 1, 1, 4, 2, 1, 4, 2, 1, 1, 7, 2, 1, 1, 7, 4, 1, 1, 8, 4, 2, 1, 8, 4, 2, 1, 1, 10, 5, 2, 1, 1, 10, 5, 2, 1, 1, 1, 11, 5, 2, 2, 1, 1, 11, 6, 3, 2, 1, 1, 15, 6, 3, 2, 1, 1, 15, 6, 3, 2, 1, 1, 1, 16, 8, 3, 2, 1, 1, 1, 16, 8, 3, 2, 1, 1, 1, 1
Offset: 2

Views

Author

Keywords

Comments

The factorization of n! is n! = 2^T(n,1)*3^T(n,2)*...*p_(pi(n))^T(n,pi(n)) where p_k = k-th prime, pi(n) = A000720(n).
Nonzero terms of A085604; T(n,k) = A085604(n,k), k = 1..A000720(n). - Reinhard Zumkeller, Nov 01 2013
For n=2, 3, 4 and 5, all terms of the n-th row are odd. Are there other such rows? - Michel Marcus, Nov 11 2018
From Gus Wiseman, May 15 2019: (Start)
Differences between successive rows are A067255, so row n is the sum of the first n row-vectors of A067255 (padded with zeros on the right so that all n row-vectors have length A000720(n)). For example, the first 10 rows of A067255 are
{}
1
0 1
2 0
0 0 1
1 1 0
0 0 0 1
3 0 0 0
0 2 0 0
1 0 1 0
with column sums (8,4,2,1), which is row 10.
(End)
For all prime p > 7, 3*p > 2*nextprime(p), so for any n > 21 there will always be a prime p dividing n! with exponent 2 and there are no further rows with all entries odd. - Charlie Neder, Jun 03 2019

Examples

			From _Gus Wiseman_, May 09 2019: (Start)
Triangle begins:
   1
   1  1
   3  1
   3  1  1
   4  2  1
   4  2  1  1
   7  2  1  1
   7  4  1  1
   8  4  2  1
   8  4  2  1  1
  10  5  2  1  1
  10  5  2  1  1  1
  11  5  2  2  1  1
  11  6  3  2  1  1
  15  6  3  2  1  1
  15  6  3  2  1  1  1
  16  8  3  2  1  1  1
  16  8  3  2  1  1  1  1
  18  8  4  2  1  1  1  1
(End)
m such that 5^m||101!: floor(log(101)/log(5)) = 2 terms. floor(101/5) = 20. floor(20/5) = 4. So m = u_1 + u_2 = 20 + 4 = 24. - _David A. Corneth_, Jun 22 2014
		

Crossrefs

Row lengths are A000720.
Row-sums are A022559.
Row-products are A135291.
Row maxima are A011371.

Programs

  • Haskell
    a115627 n k = a115627_tabf !! (n-2) !! (k-1)
    a115627_row = map a100995 . a141809_row . a000142
    a115627_tabf = map a115627_row [2..]
    -- Reinhard Zumkeller, Nov 01 2013
    
  • Maple
    A115627 := proc(n,k) local d,p; p := ithprime(k) ; n-add(d,d=convert(n,base,p)) ; %/(p-1) ; end proc: # R. J. Mathar, Oct 29 2010
  • Mathematica
    Flatten[Table[Transpose[FactorInteger[n!]][[2]], {n, 2, 20}]] (* T. D. Noe, Apr 10 2012 *)
    T[n_, k_] := Module[{p, jm}, p = Prime[k]; jm = Floor[Log[p, n]]; Sum[Floor[n/p^j], {j, 1, jm}]]; Table[Table[T[n, k], {k, 1, PrimePi[n]}], {n, 2, 20}] // Flatten (* Jean-François Alcover, Feb 23 2015 *)
  • PARI
    a(n)=my(i=2);while(n-primepi(i)>1,n-=primepi(i);i++);p=prime(n-1);sum(j=1,log(i)\log(p),i\=p) \\ David A. Corneth, Jun 21 2014

Formula

T(n,k) = Sum_{i=1..inf} floor(n/(p_k)^i). (Although stated as an infinite sum, only finitely many terms are nonzero.)
T(n,k) = Sum_{i=1..floor(log(n)/log(p_k))} floor(u_i) where u_0 = n and u_(i+1) = floor((u_i)/p_k). - David A. Corneth, Jun 22 2014

A336416 Number of perfect-power divisors of n!.

Original entry on oeis.org

1, 1, 1, 1, 3, 3, 7, 7, 11, 18, 36, 36, 47, 47, 84, 122, 166, 166, 221, 221, 346, 416, 717, 717, 1001, 1360, 2513, 2942, 4652, 4652, 5675, 5675, 6507, 6980, 13892, 17212, 20408, 20408, 39869, 45329, 51018, 51018, 68758, 68758, 105573, 138617, 284718, 284718, 338126, 421126
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

A number is a perfect power iff it is 1 or its prime exponents (signature) are not relatively prime.

Examples

			The a(1) = 0 through a(9) = 18 divisors:
       1: 1
       2: 1
       6: 1
      24: 1,4,8
     120: 1,4,8
     720: 1,4,8,9,16,36,144
    5040: 1,4,8,9,16,36,144
   40320: 1,4,8,9,16,32,36,64,128,144,576
  362880: 1,4,8,9,16,27,32,36,64,81,128,144,216,324,576,1296,1728,5184
		

Crossrefs

The maximum among these divisors is A090630, with quotient A251753.
The version for distinct prime exponents is A336414.
The uniform version is A336415.
Replacing factorials with Chernoff numbers (A006939) gives A336417.
Prime powers are A000961.
Perfect powers are A001597, with complement A007916.
Prime power divisors are counted by A022559.

Programs

  • Mathematica
    perpouQ[n_]:=Or[n==1,GCD@@FactorInteger[n][[All,2]]>1];
    Table[Length[Select[Divisors[n!],perpouQ]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, (d==1) || ispower(d)); \\ Michel Marcus, Aug 19 2020
    
  • PARI
    addhelp(val, "exponent of prime p in n!")
    val(n, p) = my(r=0); while(n, r+=n\=p);r
    a(n) = {if(n<=3, return(1)); my(pr = primes(primepi(n\2)), v = vector(#pr, i, val(n, pr[i])), res = 1, cv); for(i = 2, v[1], if(issquarefree(i), cv = v\i; res-=(prod(i = 1, #cv, cv[i]+1)-1)*(-1)^omega(i) ) ); res } \\ David A. Corneth, Aug 19 2020

Formula

a(p) = a(p-1) for prime p. - David A. Corneth, Aug 19 2020

Extensions

a(26)-a(34) from Jinyuan Wang, Aug 19 2020
a(35)-a(49) from David A. Corneth, Aug 19 2020

A325617 Multinomial coefficient of the prime signature of n!.

Original entry on oeis.org

1, 1, 1, 2, 4, 20, 105, 840, 3960, 51480, 675675, 10810800, 139675536, 2793510720, 58663725120, 1799020903680, 26985313555200, 782574093100800, 25992639520848000, 857757104187984000, 30021498646579440000, 1563341744336692320000, 64179292662243158400000
Offset: 0

Views

Author

Gus Wiseman, May 12 2019

Keywords

Comments

Number of permutations of the multiset of prime factors of n!.

Examples

			The a(5) = 20 permutations of {2,2,2,3,5}:
  (22235)  (32225)  (52223)
  (22253)  (32252)  (52232)
  (22325)  (32522)  (52322)
  (22352)  (35222)  (53222)
  (22523)
  (22532)
  (23225)
  (23252)
  (23522)
  (25223)
  (25232)
  (25322)
		

Crossrefs

Programs

  • Mathematica
    Table[Multinomial@@Last/@FactorInteger[n!],{n,0,15}]

Formula

a(n) = A318762(A181819(n!)).

A336414 Number of divisors of n! with distinct prime multiplicities.

Original entry on oeis.org

1, 1, 2, 3, 7, 10, 20, 27, 48, 86, 147, 195, 311, 390, 595, 1031, 1459, 1791, 2637, 3134, 4747, 7312, 10766, 12633, 16785, 26377, 36142, 48931, 71144, 82591, 112308, 128023, 155523, 231049, 304326, 459203, 568095, 642446, 812245, 1137063, 1441067, 1612998, 2193307, 2429362
Offset: 0

Views

Author

Gus Wiseman, Jul 22 2020

Keywords

Comments

A number has distinct prime multiplicities iff its prime signature is strict.

Examples

			The first and second columns below are the a(6) = 20 counted divisors of 6! together with their prime signatures. The third column shows the A000005(6!) - a(6) = 10 remaining divisors.
      1: ()      20: (2,1)    |    6: (1,1)
      2: (1)     24: (3,1)    |   10: (1,1)
      3: (1)     40: (3,1)    |   15: (1,1)
      4: (2)     45: (2,1)    |   30: (1,1,1)
      5: (1)     48: (4,1)    |   36: (2,2)
      8: (3)     72: (3,2)    |   60: (2,1,1)
      9: (2)     80: (4,1)    |   90: (1,2,1)
     12: (2,1)  144: (4,2)    |  120: (3,1,1)
     16: (4)    360: (3,2,1)  |  180: (2,2,1)
     18: (1,2)  720: (4,2,1)  |  240: (4,1,1)
		

Crossrefs

Perfect-powers are A001597, with complement A007916.
Numbers with distinct prime multiplicities are A130091.
Divisors with distinct prime multiplicities are counted by A181796.
The maximum divisor with distinct prime multiplicities is A327498.
Divisors of n! with equal prime multiplicities are counted by A336415.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n!],UnsameQ@@Last/@FactorInteger[#]&]],{n,0,15}]
  • PARI
    a(n) = sumdiv(n!, d, my(ex=factor(d)[,2]); #vecsort(ex,,8) == #ex); \\ Michel Marcus, Jul 24 2020

Formula

a(n) = A181796(n!).

Extensions

a(21)-a(41) from Alois P. Heinz, Jul 24 2020

A325508 Product of primes indexed by the prime exponents of n!.

Original entry on oeis.org

1, 1, 2, 4, 10, 20, 42, 84, 204, 476, 798, 1596, 3828, 7656, 12276, 24180, 36660, 73320, 120840, 241680, 389424, 785680, 1294440, 2588880, 3848880, 7147920, 11264760, 15926040, 26057304, 52114608, 74421648, 148843296, 187159392, 340949280, 527531760, 926505360
Offset: 0

Views

Author

Gus Wiseman, May 08 2019

Keywords

Comments

The prime indices of a(n) are the signature of n!, which is row n of A115627.

Examples

			We have 7! = 2^4 * 3^2 * 5^1 * 7^1, so a(7) = prime(4)*prime(2)*prime(1)*prime(1) = 84.
The sequence of terms together with their prime indices begins:
          1: {}
          1: {}
          2: {1}
          4: {1,1}
         10: {1,3}
         20: {1,1,3}
         42: {1,2,4}
         84: {1,1,2,4}
        204: {1,1,2,7}
        476: {1,1,4,7}
        798: {1,2,4,8}
       1596: {1,1,2,4,8}
       3828: {1,1,2,5,10}
       7656: {1,1,1,2,5,10}
      12276: {1,1,2,2,5,11}
      24180: {1,1,2,3,6,11}
      36660: {1,1,2,3,6,15}
      73320: {1,1,1,2,3,6,15}
     120840: {1,1,1,2,3,8,16}
     241680: {1,1,1,1,2,3,8,16}
		

Crossrefs

Programs

  • Mathematica
    Table[Times@@Prime/@Last/@If[(n!)==1,{},FactorInteger[n!]],{n,0,30}]

Formula

a(n) = A181819(n!).
A001221(a(n)) = A071626(n).
A001222(a(n)) = A000720(n).
A056239(a(n)) = A022559(n).
A003963(a(n)) = A135291(n).
A061395(a(n)) = A011371(n).
A007814(a(n)) = A056171(n).
a(n) = A122111(A307035(n)). - Antti Karttunen, Nov 19 2019

A076934 Smallest integer of the form n/k!.

Original entry on oeis.org

1, 1, 3, 2, 5, 1, 7, 4, 9, 5, 11, 2, 13, 7, 15, 8, 17, 3, 19, 10, 21, 11, 23, 1, 25, 13, 27, 14, 29, 5, 31, 16, 33, 17, 35, 6, 37, 19, 39, 20, 41, 7, 43, 22, 45, 23, 47, 2, 49, 25, 51, 26, 53, 9, 55, 28, 57, 29, 59, 10, 61, 31, 63, 32, 65, 11, 67, 34, 69, 35, 71
Offset: 1

Views

Author

Amarnath Murthy, Oct 19 2002

Keywords

Comments

Equivalently, n divided by the largest factorial divisor of n.
Also, the smallest r such that n/r is a factorial number.
Positions of 1's are the factorial numbers A000142. Is every positive integer in this sequence? - Gus Wiseman, May 15 2019
Let m = A055874(n), the largest integer such that 1,2,...,m divides n. Then a(n*m!) = n since m+1 does not divide n, showing that every integer is part of the sequence. - Etienne Dupuis, Sep 19 2020

Crossrefs

Programs

  • Mathematica
    Table[n/Max@@Intersection[Divisors[n],Array[Factorial,n]],{n,100}] (* Gus Wiseman, May 15 2019 *)
    a[n_] := Module[{k=1}, While[Divisible[n, k!], k++]; n/(k-1)!]; Array[a, 100] (* Amiram Eldar, Dec 25 2023 *)
  • PARI
    first(n) = {my(res = [1..n]); for(i = 2, oo, k = i!; if(k <= n, for(j = 1, n\k, res[j*k] = j ) , return(res) ) ) } \\ David A. Corneth, Sep 19 2020

Formula

From Amiram Eldar, Dec 25 2023: (Start)
a(n) = n/A055881(n)!.
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = BesselI(2, 2) = 0.688948... (A229020). (End)

Extensions

More terms from David A. Corneth, Sep 19 2020
Showing 1-10 of 41 results. Next