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

A273673 Square array A(n,k) = (n / prime(1+A084558(k))^e) * prime(1+A084558(k)-A099563(k))^e, where e = A249344((1+A084558(k)), n) = the exponent of the largest power of prime(1+A084558(k)) which divides n. Array is read by descending antidiagonals as A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), etc.

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 1, 2, 3, 4, 1, 2, 3, 4, 5, 1, 2, 3, 4, 3, 4, 1, 2, 3, 4, 3, 6, 7, 1, 2, 3, 4, 2, 6, 7, 8, 1, 2, 3, 4, 2, 6, 7, 8, 4, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 6, 11, 1, 2, 3, 4, 5, 6, 5, 8, 9, 6, 11, 8, 1, 2, 3, 4, 5, 6, 5, 8, 9, 4, 11, 12, 13, 1, 2, 3, 4, 5, 6, 5, 8, 9, 4, 11, 12, 13, 14, 1, 2, 3, 4, 5, 6, 5, 8, 9, 10, 11, 12, 13, 14, 10
Offset: 1

Views

Author

Antti Karttunen, Aug 09 2016

Keywords

Comments

Informally: "clear" the exponent of prime(1+A084558(k)) and add it (the old value of exponent) to the exponent of prime(1+A084558(k)-A099563(k)) in the prime factorization of n.
Auxiliary function for computing array A275723.

Examples

			The top left 6 x 15 corner of the array:
   1,  1,  1,  1,  1,  1
   2,  2,  2,  2,  2,  2
   2,  3,  3,  3,  3,  3
   4,  4,  4,  4,  4,  4
   5,  3,  3,  2,  2,  5
   4,  6,  6,  6,  6,  6
   7,  7,  7,  7,  7,  5
   8,  8,  8,  8,  8,  8
   4,  9,  9,  9,  9,  9
  10,  6,  6,  4,  4, 10
  11, 11, 11, 11, 11, 11
   8, 12, 12, 12, 12, 12
  13, 13, 13, 13, 13, 13
  14, 14, 14, 14, 14, 10
  10,  9,  9,  6,  6, 15
		

Crossrefs

Programs

Formula

A(n,k) = (n / prime(1+A084558(k))^e) * prime(1+A084558(k)-A099563(k))^e, where e = A249344((1+A084558(k)), n), the exponent of the largest power prime(1+A084558(k)) which divides n.

A104244 Suppose m = Product_{i=1..k} p_i^e_i, where p_i is the i-th prime number and each e_i is a nonnegative integer. Then we can define P_m(x) = Sum_{i=1..k} e_i*x^(i-1). The sequence is the square array A(n,m) = P_m(n) read by descending antidiagonals.

Original entry on oeis.org

0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 1, 2, 3, 1, 0, 2, 4, 2, 4, 1, 0, 1, 3, 9, 2, 5, 1, 0, 3, 8, 4, 16, 2, 6, 1, 0, 2, 3, 27, 5, 25, 2, 7, 1, 0, 2, 4, 3, 64, 6, 36, 2, 8, 1, 0, 1, 5, 6, 3, 125, 7, 49, 2, 9, 1, 0, 3, 16, 10, 8, 3, 216, 8, 64, 2, 10, 1, 0, 1, 4, 81, 17, 10, 3, 343, 9, 81, 2, 11, 1, 0, 2, 32, 5
Offset: 1

Views

Author

Olaf Voß, Feb 26 2005

Keywords

Comments

From Antti Karttunen, Jul 29 2015: (Start)
The square array A(row,col) is read by downwards antidiagonals as: A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), etc.
A(n,m) (entry at row=n, column=m) gives the evaluation at x=n of the polynomial (with nonnegative integer coefficients) bijectively encoded in the prime factorization of m. See A206284, A206296 for the details of that encoding. (The roles of variables n and m were accidentally swapped in this description, corrected by Antti Karttunen, Oct 30 2016)
(End)
Each row is a completely additive sequence, row n mapping prime(m) to n^(m-1). - Peter Munn, Apr 22 2022

Examples

			a(13) = 3 because 3 = p_1^0 * p_2^1 * p_3^0 * ..., so P_3(x) = 0*x^(1-1) + 1*x^(2-1) + 0*x^(3-1) + ... = x. Hence a(13) = A(3,3) = P_3(3) = 3. [Elaborated by _Peter Munn_, Aug 13 2022]
The top left corner of the array:
0, 1,  1, 2,   1,  2,   1,  3,  2,   2,     1,  3,      1,    2,   2, 4
0, 1,  2, 2,   4,  3,   8,  3,  4,   5,    16,  4,     32,    9,   6, 4
0, 1,  3, 2,   9,  4,  27,  3,  6,  10,    81,  5,    243,   28,  12, 4
0, 1,  4, 2,  16,  5,  64,  3,  8,  17,   256,  6,   1024,   65,  20, 4
0, 1,  5, 2,  25,  6,  125, 3, 10,  26,   625,  7,   3125,  126,  30, 4
0, 1,  6, 2,  36,  7,  216, 3, 12,  37,  1296,  8,   7776,  217,  42, 4
0, 1,  7, 2,  49,  8,  343, 3, 14,  50,  2401,  9,  16807,  344,  56, 4
0, 1,  8, 2,  64,  9,  512, 3, 16,  65,  4096, 10,  32768,  513,  72, 4
0, 1,  9, 2,  81, 10,  729, 3, 18,  82,  6561, 11,  59049,  730,  90, 4
0, 1, 10, 2, 100, 11, 1000, 3, 20, 101, 10000, 12, 100000, 1001, 110, 4
...
		

Crossrefs

Cf. A000720.
Transpose: A104245.
Main diagonal: A090883.
Row 1: A001222, row 2: A048675, row 3: A090880, row 4: A090881, row 5: A090882, row 10: A054841; and, in the extrapolated table, row 0: A007814, row -1: A195017.
Other completely additive sequences with prime(k) mapped to a function of k include k: A056239, k-1: A318995, k+1: A318994, k^2: A289506, 2^k-1: A293447, k!: A276075, F(k-1): A265753, F(k-2): A265752.
For completely additive sequences with primes p mapped to a function of p, see A001414.
For completely additive sequences where some primes are mapped to 1, the rest to 0 (notably, some ruler functions) see the cross-references in A249344.
For completely additive sequences, s, with primes p mapped to a function of s(p-1) and maybe s(p+1), see A352957.
See the formula section for the relationship to A073133, A206296.
See the comments for the relevance of A206284.
A297845 represents multiplication of the relevant polynomials.
Cf. A090884, A248663, A265398, A265399 for other related sequences.
A167219 lists columns that contain their own column number.

Formula

A(n,A206296(k)) = A073133(n,k). [This formula demonstrates how this array can be used with appropriately encoded polynomials. Note that A073133 reads its antidiagonals by ascending order, while here the order is opposite.] - Antti Karttunen, Oct 30 2016
From Peter Munn, Apr 05 2021: (Start)
The sequence is defined by the following identities:
A(n, 3) = n;
A(n, m*k) = A(n, m) + A(n, k);
A(n, A297845(m, k)) = A(n, m) * A(n, k).
(End)

Extensions

Starting offset changed from 0 to 1 by Antti Karttunen, Jul 29 2015
Name edited (and aligned with rest of sequence) by Peter Munn, Apr 23 2022

A060175 Square array A(n,k) = exponent of the largest power of k-th prime which divides n, read by falling antidiagonals.

Original entry on oeis.org

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

Views

Author

Henry Bottomley, Mar 14 2001

Keywords

Examples

			a(12,1) = 2 since 4 = 2^2 = p_1^2 divides 12 but 8 = 2^3 does not.
a(12,2) = 1 since 3 = p_2 divides 12 but 9 = 3^2 does not.
See also examples in A249344, which is transpose of this array.
The top-left corner of the array:
n\k | 1  2  3  4  5  6  7  8
----+------------------------
1   | 0, 0, 0, 0, 0, 0, 0, 0,
2   | 1, 0, 0, 0, 0, 0, 0, 0,
3   | 0, 1, 0, 0, 0, 0, 0, 0,
4   | 2, 0, 0, 0, 0, 0, 0, 0,
5   | 0, 0, 1, 0, 0, 0, 0, 0,
6   | 1, 1, 0, 0, 0, 0, 0, 0,
7   | 0, 0, 0, 1, 0, 0, 0, 0,
8   | 3, 0, 0, 0, 0, 0, 0, 0,
9   | 0, 2, 0, 0, 0, 0, 0, 0,
10  | 1, 0, 1, 0, 0, 0, 0, 0,
11  | 0, 0, 0, 0, 1, 0, 0, 0,
12  | 2, 1, 0, 0, 0, 0, 0, 0,
...
		

Crossrefs

Transpose: A249344.
Column 1: A007814.
Column 2: A007949.
Column 3: A112765.
Column 4: A214411.
Row sums: A001222.

Programs

  • Mathematica
    T[n_, k_] := IntegerExponent[n, Prime[k]];
    Table[T[n-k+1, k], {n, 1, 15}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Nov 18 2019 *)
  • PARI
    a(n, k) = valuation(n, prime(k)); \\ Michel Marcus, Jun 24 2017
  • Python
    from sympy import prime
    def a(n, k):
        p=prime(n)
        i=z=0
        while p**i<=k:
            if k%(p**i)==0: z=i
            i+=1
        return z
    for n in range(1, 10): print([a(n - k + 1, k) for k in range(1, n + 1)]) # Indranil Ghosh, Jun 24 2017
    
  • Scheme
    (define (A060175 n) (A249344bi (A004736 n) (A002260 n)))
    (define (A249344bi row col) (let ((p (A000040 row))) (let loop ((n col) (i 0)) (cond ((not (zero? (modulo n p))) i) (else (loop (/ n p) (+ i 1)))))))
    ;; Antti Karttunen, Oct 28 2014
    

Formula

A(n, k) = log(A060176(n, k))/log(A000040(k)) = k-th digit from right of A054841(n).

Extensions

Erroneous example corrected and more terms computed by Antti Karttunen, Oct 28 2014
Name clarified by Antti Karttunen, Jan 16 2025

A249421 A(n,k) = exponent of the largest power of n-th prime which divides the product of the elements on row (k-1) of Pascal's triangle; a square array read by antidiagonals.

Original entry on oeis.org

0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 5, 2, 0, 0, 0, 2, 1, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 4, 4, 0, 0, 0, 0, 0, 17, 2, 3, 0, 0, 0, 0, 0, 0, 10, 0, 2, 0, 0, 0, 0, 0, 0, 0, 12, 14, 1, 6, 0, 0, 0, 0, 0, 0, 0, 4, 10, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 18, 6, 8, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 8, 13, 6, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 8, 4, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Antti Karttunen, Oct 28 2014

Keywords

Comments

Square array A(n,k), where n = row, k = column, read by antidiagonals: A(1,1), A(1,2), A(2,1), A(1,3), A(2,2), A(3,1), ....
A(n,k) is A000040(n)-adic valuation of A001142(k-1).

Examples

			The top left corner of the array:
0, 0, 1, 0, 5, 2, 4, 0, 17, 10, 12,  4, 18,  8, 11,  0, 49, 34, 36, 20, 42,
0, 0, 0, 2, 1, 0, 4, 2,  0, 14, 10,  6, 13,  8,  3, 12,  6,  0, 28, 20, 12,
0, 0, 0, 0, 0, 4, 3, 2,  1,  0,  8,  6,  4,  2,  0, 12,  9,  6,  3,  0, 16,
0, 0, 0, 0, 0, 0, 0, 6,  5,  4,  3,  2,  1,  0, 12, 10,  8,  6,  4,  2,  0,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0, 10,  9,  8,  7,  6,  5,  4,  3,  2,  1,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0,  0,  0, 12, 11, 10,  9,  8,  7,  6,  5,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 16, 15, 14, 13,
0, 0, 0, 0, 0, 0, 0, 0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0,  0, 18, 17,
...
		

Crossrefs

Transpose: A249422.
Row 1: A187059, Row 2: A249343, Row 3: A249345, Row 4 A249347. (Cf. also A249346).

Programs

Formula

A(n, k) = A249344(n, A001142(k-1)).

A248601 Square array T(n, k) read by antidiagonals, n > 0 and k > 0: for any number n > 0, let f(n) be the function that associates k to the prime(k)-adic valuation of n (where prime(k) denotes the k-th prime); f establishes a bijection between the positive numbers and the arithmetic functions with nonnegative integer values and a finite number of nonzero values; let g be the inverse of f; T(n, k) = g(f(n) * f(k)) (where i * j denotes the Dirichlet convolution of i and j).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 7, 4, 1, 1, 5, 9, 9, 5, 1, 1, 6, 13, 16, 13, 6, 1, 1, 7, 21, 25, 25, 21, 7, 1, 1, 8, 19, 36, 23, 36, 19, 8, 1, 1, 9, 27, 49, 65, 65, 49, 27, 9, 1, 1, 10, 49, 64, 37, 126, 37, 64, 49, 10, 1, 1, 11, 39, 81, 125, 133, 133, 125
Offset: 1

Views

Author

Rémy Sigrist, Dec 21 2017

Keywords

Comments

For any n > 0, f(n) corresponds to the function k -> A249344(k, n).
For any n > 0 and m > 0, f(n * m) = f(n) + f(m).
Also, f(1) = A000004 and f(2) corresponds to k -> A000007(k-1).
The function f can be naturally extended to the set of positive rational numbers: if r = u/v (not necessarily in reduced form), then f(r) = f(u) - f(v); as such, f is an homomorphism from the multiplicative group of positive rational numbers to the additive group of arithmetic functions with integer values and a finite number of nonzero values.
For any arithmetic function with integer values and a finite number of nonzero values j, g(j) = Product_{k > 0} A000040(k)^j(k).
See A296857 for the main diagonal of T.

Examples

			Array T(n, k) begins:
  n\k|  1   2   3    4    5    6    7     8     9    10
  ---+-------------------------------------------------
    1|  1   1   1    1    1    1    1     1     1     1  -> A000012
    2|  1   2   3    4    5    6    7     8     9    10  -> A000027
    3|  1   3   7    9   13   21   19    27    49    39  -> A297002
    4|  1   4   9   16   25   36   49    64    81   100  -> A000290
    5|  1   5  13   25   23   65   37   125   169   115
    6|  1   6  21   36   65  126  133   216   441   390
    7|  1   7  19   49   37  133   53   343   361   259
    8|  1   8  27   64  125  216  343   512   729  1000  -> A000578
    9|  1   9  49   81  169  441  361   729  2401  1521
   10|  1  10  39  100  115  390  259  1000  1521  1150
		

Crossrefs

Cf. A000004, A000007, A000012 (first row/column), A000027 (second row/column), A000040, A000290 (fourth row/column), A000578 (eighth row/column), A001221, A001222, A055396, A061395, A249344, A289508, A296857 (main diagonal), A297002 (third row/column).

Programs

  • PARI
    T(n,k) = my(fn=factor(n), pn=apply(primepi,fn[,1]~), fk=factor(k), pk=apply(primepi,fk[,1]~)); prod(i=1, #pn, prod(j=1, #pk, prime(pn[i]*pk[j])^(fn[i,2]*fk[j,2])))

Formula

T is completely multiplicative in both parameters:
- for any n > 0
- and k > 0 with prime factorization Prod_{i > 0} prime(i)^e_i:
- T(prime(n), k) = T(k, prime(n)) = Prod_{i > 0} prime(n * i)^e_i.
For any m > 0, n > 0 and k > 0:
- T(n, k) = T(k, n) (T is commutative),
- T(m, T(n, k)) = T(T(m, n), k) (T is associative),
- T(n, 1) = 1 (1 is an absorbing element for T),
- T(n, 2) = n (2 is an identity element for T),
- T(n, 2^i) = n^i for any i >= 0,
- T(n, 4) = n^2 (A000290),
- T(n, 8) = n^3 (A000578),
- T(n, 3) = A297002(n),
- T(n, 3^i) = A297002(n)^i for any i >= 0,
- A001221(T(n, k)) <= A001221(n) * A001221(k),
- A001222(T(n, k)) = A001222(n) * A001222(k),
- A055396(T(n, k)) = A055396(n) * A055396(k),
- A061395(T(n, k)) = A061395(n) * A061395(k),
- A289508(T(n, k)) = A289508(n) * A289508(k),
- T(A000040(n), A000040(k)) = A000040(n * k),
- T(A000040(n)^i, A000040(k)^j) = A000040(n * k)^(i * j) for any i >= 0 and j >= 0.

A352957 Triangle read by rows: Row n is the lexicographically earliest strictly monotonic completely additive sequence of length n.

Original entry on oeis.org

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

Views

Author

Peter Munn, Apr 11 2022

Keywords

Comments

Each sequence consists of nonnegative integers indexed from 1.
Note in particular in the formula section, the lower bound, floor(n/k), for first differences between terms in a row. This follows (using the additive property) from the strict monotonicity of floor(n/k)+1 consecutive terms near the end of the row.
For any k, with increasing length n >= k, the first k terms of the sequences approach similarity with a real-valued logarithmic function defined on the integers. For example, the asymptote of T(n,3)/T(n,2) is log(3)/log(2), A020857.

Examples

			(For row 4.) A completely additive sequence requires T(4,1) = 0. Strict monotonicity requires T(4,4) > T(4,3) > T(4,2). So T(4,4) >= T(4,2) + 2. Using the additivity this becomes T(4,2) + T(4,2) >= T(4,2) + T(4,1) + 2. Subtracting T(4,2) and substituting 0 for T(4,1) we get T(4,2) >= 2. So from T(4,4) > T(4,3) > T(4,2), we see T(4,3) >= 3, T(4,4) >= 4. So row 4 = (0, 2, 3, 4) as it is strictly monotonic and completely additive and from the preceding arguments is seen to be the lexicographically earliest such.
Triangle starts:
0;
0, 1;
0, 1,  2;
0, 2,  3,  4;
0, 2,  3,  4,  5;
0, 3,  5,  6,  7,  8;
0, 3,  5,  6,  7,  8,  9;
0, 4,  6,  8,  9, 10, 11, 12;
0, 5,  8, 10, 11, 13, 14, 15, 16;
0, 5,  8, 10, 12, 13, 14, 15, 16, 17;
0, 5,  8, 10, 12, 13, 14, 15, 16, 17, 18;
0, 7, 11, 14, 16, 18, 19, 21, 22, 23, 24, 25;
0, 7, 11, 14, 16, 18, 19, 21, 22, 23, 24, 25, 26;
0, 7, 11, 14, 16, 18, 20, 21, 22, 23, 24, 25, 26, 27;
0, 8, 13, 16, 19, 21, 23, 24, 26, 27, 28, 29, 30, 31, 32;
0, 9, 14, 18, 21, 23, 25, 27, 28, 30, 31, 32, 33, 34, 35, 36;
		

Crossrefs

Cf. A020857.
Completely additive sequences, s, with primes p mapped to a function of s(p-1) and maybe s(p+1): A064097, A344443, A344444; and for functions of earlier terms, see A334200.
For completely additive sequences with primes p mapped to a function of p, see A001414.
For completely additive sequences with prime(k) mapped to a function of k, see A104244.
For completely additive sequences where some primes are mapped to 1, the rest to 0 (notably, some ruler functions) see the cross-references in A249344.

Formula

The definition specifies: T(n,j*k) = T(n,j) + T(n,k); for k > 1, T(n,k) > T(n,k-1).
T(n,1) = 0, otherwise T(n,k) >= T(n,k-1) + floor(n/k).
For prime p, T(p,p) = T(p-1,p-1) + 1, otherwise T(p,k) = T(p-1,k).
T(n,2) >= 2*floor(n/4) + floor(n/9).
T(n,3) >= ceiling( (3*T(n,2) + floor(n/9)) / 2).
T(11,k) = A344443(k).
For k <> 13, T(23,k) = A344444(k).

A375537 Square array A(n, k) (n, k >= 1) read by antidiagonals in ascending order: A(n, k) = Max_{i = 1..n} v_prime(i)(k), where v_p(k) is the p-adic valuation of k.

Original entry on oeis.org

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

Views

Author

Amiram Eldar, Aug 19 2024

Keywords

Comments

For a given n, A(n, k) is the sequence that gives the maximum exponent in the prime factorization of the largest prime(n)-smooth divisor of k.

Examples

			Array begins:
   n | n-th row
  ---+-----------------------------
   1 | 0, 1, 0, 2, 0, 1, 0, 3, 0, 1
   2 | 0, 1, 1, 2, 0, 1, 0, 3, 2, 1
   3 | 0, 1, 1, 2, 1, 1, 0, 3, 2, 1
   4 | 0, 1, 1, 2, 1, 1, 1, 3, 2, 1
   5 | 0, 1, 1, 2, 1, 1, 1, 3, 2, 1
   6 | 0, 1, 1, 2, 1, 1, 1, 3, 2, 1
   7 | 0, 1, 1, 2, 1, 1, 1, 3, 2, 1
   8 | 0, 1, 1, 2, 1, 1, 1, 3, 2, 1
   9 | 0, 1, 1, 2, 1, 1, 1, 3, 2, 1
  10 | 0, 1, 1, 2, 1, 1, 1, 3, 2, 1
		

Crossrefs

Programs

  • Mathematica
    A[n_, k_] := Max[IntegerExponent[k, Prime[Range[n]]]]; Table[A[n - k + 1, k], {n, 1, 14}, {k, 1 n}] // Flatten
  • PARI
    A(n, k) = vecmax(apply(x -> valuation(k, x), primes(n)));

Formula

A(n, k) = Max_{i=1..n} A249344(i, k).
A(n, k) = A051903(k) for n >= A000720(A006530(k)).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{i=1..m} A(n, i) = A375538(n)/A375539(n).
Showing 1-7 of 7 results.