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.

User: William P. Orrick

William P. Orrick's wiki page.

William P. Orrick has authored 20 sequences. Here are the ten most recent ones:

A368134 Characteristic numbers of Markov triples in the binary tree A368546.

Original entry on oeis.org

2, 5, 12, 13, 75, 179, 70, 34, 507, 2923, 1120, 2673, 15571, 6089, 408, 89, 3468, 51709, 19760, 113922, 1701181, 651838, 16725, 39916, 3472225, 20226717, 1354498, 529673, 3087111, 206855, 2378, 233, 23763, 925943, 353702, 5273811, 205543262, 78545995, 770133
Offset: 0

Author

William P. Orrick, Jan 11 2024

Keywords

Comments

The characteristic number u of a Markov triple (r, m, s) is the solution in (0, m) of r * x == s (mod m). It satisfies u^2 == -1 (mod m), so that v = (u^2 + 1) / m is also an integer. The other solution in (0,m) of u^2 == -1 (mod m), namely m - u, is always greater than u, so u < m / 2.
The Markov tree may be formulated in terms of a set of Cohn matrices. There is a one-parameter family of such sets, parametrized by an integer c. Given a vertex of the Markov tree with Farey triple (x, y, z) and Markov triple (r, m, s), producing characteristic number u and v = (u^2 + 1) / m, the Cohn matrix C_y(c) with parameter c is
[ c * m + u m ]
[(3 * c - c^2) * m - (2 * c - 3) * u - v (3 - c) * m - u].
Then the vertex is associated with a triple of Cohn matrices, (R, R S, S), where R = C_x(c), RS = C_y(c), and S = C_z(c). See A368546 for a description of Farey and Markov triples. The left child of the vertex is associated with the triple (R, R^2 S, RS) and the right child with (RS, R S^2, S).

Examples

			The initial rows of the binary tree are
                                  2
             5                                       12
   13                 75                    179                70
34    507        2923   1120            2673  15571       6089   408
		

References

  • Martin Aigner, Markov's theorem and 100 years of the uniqueness conjecture. A mathematical journey from irrational numbers to perfect matchings. Springer, 2013. x+257 pp. ISBN: 978-3-319-00887-5; 978-3-319-00888-2 MR3098784

Crossrefs

Cf. A368546.

Programs

  • SageMath
    rowM = [[1,5,2]]
    rowU = [[0,2,1]]
    a368134 = [2]
    for rw in range(1,6):
        prevRowM = rowM
        prevRowU = rowU
        rowM = []
        rowU = []
        for i in range(len(prevRowM)):
            [r,m,s] = prevRowM[i]
            [t,u,v] = prevRowU[i]
            ltM = [r,3*r*m - s,m]
            rtM = [m,3*m*s - r,s]
            ltU = [t,3*r*u - v,u]
            rtU = [u,3*u*s - t,v]
            rowM = rowM + [ltM,rtM]
            rowU = rowU + [ltU,rtU]
            a368134 = a368134 + [ltU[1],rtU[1]]
    a368134

Formula

Recurrence: The left child of the Markov triple (r, m, s) is (r, 3rm - s, m); the right child is (m, 3ms - r, s). The corresponding triple of characteristic numbers (t, u, v) has left child (t, 3ru - v, u) and right child (u, 3us - t, v). Initial Markov triple: (1, 5, 2), initial characteristic number triple: (0, 2, 1).

A370852 Irregular triangle read by rows: row n is the list of residues mod n that occur among the Markov numbers.

Original entry on oeis.org

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

Author

Keywords

Comments

Length of row n is A370164(n).

Examples

			The first rows are:
n
1:   0
2:   0  1
3:   1  2
4:   1  2
5:   0  1  2  3  4
6:   1  2  4  5
7:   1  2  5  6
8:   1  2  5
9:   1  2  4  5  7  8
10:  0  1  2  3  4  5  6  7  8  9
11:  1  2  4  5  6  7  9 10
12:  1  2  5 10
13:  0  1  2  3  4  5  6  7  8  9 10 11 12
14:  1  2  5  6  8  9 12 13
15:  1  2  4  5  7  8 10 11 13 14
16:  1  2  5  9 13
17:  0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16
18:  1  2  4  5  7  8 10 11 13 14 16 17
19:  1  2  3  4  5  6  8  9 10 11 13 14 15 16 17 18
20:  1  2  5  6  9 10 13 14 17 18
For n = 14 residues congruent to 0, 3, or 4 mod 7 are forbidden. (See comments to A370164 for explanation.) All other residues occur. For example, the Markov numbers 1, 2, 5, 34, 610, 1325, 194, and 13 produce the residues shown in row 14 of the triangle (mod 14).
		

References

  • Martin Aigner, Markov's theorem and 100 years of the uniqueness conjecture. A mathematical journey from irrational numbers to perfect matchings. Springer, 2013. x+257 pp. ISBN: 978-3-319-00887-5; 978-3-319-00888-2 MR3098784.

Crossrefs

Markov numbers: A002559.
Markov tree: A327345, A368546.
Cf. A370164.

Programs

  • SageMath
    def generateAllMarkovTreeResidues(n):
        row = [[1 % n,5 % n,2 % n]]
        residuesFound = []
        triplesFound = []
        while row != []:
            newRow = []
            for trpl in row:
                if trpl[1] not in residuesFound:
                    residuesFound.append(trpl[1])
                if trpl[2] < trpl[0]:
                    trpl.reverse()
                if trpl not in triplesFound:
                    triplesFound.append(trpl)
                    newRow.append([trpl[0],(3*trpl[0]*trpl[1]-trpl[2]) % n,trpl[1]])
                    newRow.append([trpl[1],(3*trpl[1]*trpl[2]-trpl[0]) % n,trpl[2]])
            row = newRow
        residuesFound.sort()
        return(residuesFound)
    [r for n in range(1,16) for r in generateAllMarkovTreeResidues(n)]

A370164 The number of residues mod n that occur among the Markov numbers.

Original entry on oeis.org

1, 2, 2, 2, 5, 4, 4, 3, 6, 10, 8, 4, 13, 8, 10, 5, 17, 12, 16, 10, 8, 16, 20, 6, 25, 26, 18, 8, 29, 20, 28, 9, 16, 34, 20, 12, 37, 32, 26, 15, 41, 16, 40, 16, 30, 40, 44, 10, 28, 50, 34, 26, 53, 36, 40, 12, 32, 58, 56, 20, 61, 56, 24, 18, 65, 32, 64, 34, 40
Offset: 1

Author

Keywords

Comments

No Markov number is divisible by any prime congruent to 3 (mod 4). Proof: Let m be a Markov number and let p be an odd prime dividing m. Then there exist positive integers b and c such that m^2 + b^2 + c^2 = 3 * m * b * c (Markov's equation). Furthermore it is known that b and c must be coprime to m. Evaluating Markov's equation (mod p) gives b^2 = -c^2 (mod p), which implies -1 is a square (mod p). This is only possible if p = 1 (mod 4).
Given a prime p congruent to 3 (mod 4) and greater than 3, no Markov number is congruent to 1/3 * (p - 2 * sigma(p)) or 2/3 * (p + sigma(p)) (mod p), where sigma(p) = -1 if p is congruent to 7 (mod 12) and sigma(p) = 1 if p is congruent to 11 (mod 12). Proof: Let m be a Markov number congruent to one of the forbidden residues (mod p). Then evaluating Markov's equation (mod p) using the forms of the forbidden residues above and doing a few simplifications implies that -4 is a square (mod p). This contradicts that p = 3 (mod 4).
The first comment above implies that an odd Markov number is congruent to 1 (mod 4). It has also been shown that any positive even integer satisfying all of the constraints of the first two comments above is congruent to 2 (mod 32).
It is conjectured that, modulo n, all residues not forbidden by the constraints in the three comments above are actually realized by some Markov number. Specifically, mod n the forbidden residues include the following: (1) if p is a prime congruent to 3 (mod 4) that divides n, then any residue congruent to 0 (mod p) or congruent to either of the forbidden residues listed in the second comment (mod p), (2) if 4 divides n, then any odd residue congruent to 3 (mod 4), (3) if 2^r is the highest power of 2 dividing n and 2 <= r <= 5, then any even residue not congruent to 2 (mod 2^r), (4) if 2^r is the highest power of 2 dividing n and r > 5, then any even residue not congruent to 2 (mod 32). It is conjectured that all other residues occur. This has been verified for all n <= 38000.
If the conjecture in the previous comment is correct, then it follows from the Chinese remainder theorem that a(n) may be found by writing n = 2^s * 3^t * u * p_1^r_1 * p_2^r_2 * ... * p_k^r_k, where p_1, ..., p_k are distinct primes greater than 3 and congruent to 3 (mod 4) and r_1, ..., r_k are positive and where the prime divisors of u are all congruent to 1 (mod 4). Then a(n) = u * C_s * D_t * Product_{j=1..k} (p_j - 3) * p_j^(r_j - 1), where C_s = 2^s if s < 2, 1 + 2^(s-2) if 2 <= s <= 5, and 2^(s - 5) + 2^(s - 2) if s > 5, and where D_t = 1 if t = 0 and 2 * 3^(t-1) if t > 0. If the conjecture holds then a(n) is a multiplicative function.

Examples

			If n = 56 = 7 * 8 then, since only the residues 1, 2, 5, 6 are allowed (mod 7) and only the residues 1, 2, 5 are allowed (mod 8), the number of potential residues (mod 56) is 4 * 3 = 12, and these residues are 1, 2, 5, 9, 13, 26, 29, 33, 34, 37, 41, 50. That these residues are realized by Markov numbers is witnessed by 1, 2, 5, 233, 13, 194, 29, 89, 34, 1325, 433, 610.
		

References

  • Martin Aigner, Markov's theorem and 100 years of the uniqueness conjecture. A mathematical journey from irrational numbers to perfect matchings. Springer, 2013. x+257 pp. ISBN: 978-3-319-00887-5; 978-3-319-00888-2 MR3098784.

Crossrefs

Markov numbers: A002559.
Markov tree: A327345, A368546.
Triangle giving list of residues mod n: A370852.

Formula

Conjectured: write n = 2^s * 3^t * u * p_1^r_1 * p_2^r_2 * ... * p_k^r_k, where p_1, ..., p_k are distinct primes greater than 3 and congruent to 3 (mod 4) and r_1, ..., r_k are positive and where the prime divisors of u are all congruent to 1 (mod 4). Then a(n) = u * C_s * D_t * Product_{j=1..k} (p_j - 3) * p_j^(r_j - 1), where C_s = 2^s if s < 2, 1 + 2^(s-2) if 2 <= s <= 5, and 2^(s - 5) + 2^(s - 2) if s > 5, and where D_t = 1 if t = 0 and 2 * 3^(t-1) if t > 0.

A368546 Alternative version of the Markov tree A327345.

Original entry on oeis.org

5, 13, 29, 34, 194, 433, 169, 89, 1325, 7561, 2897, 6466, 37666, 14701, 985, 233, 9077, 135137, 51641, 294685, 4400489, 1686049, 43261, 96557, 8399329, 48928105, 3276509, 1278818, 7453378, 499393, 5741, 610, 62210, 2423525, 925765, 13782649, 537169541
Offset: 0

Author

William P. Orrick, Jan 04 2024

Keywords

Comments

The Markov tree is a complete, infinite binary tree. Vertices are labeled by triples. The root vertex is (1, 5, 2). The left child of (a, b, c) is (a, 3*a*b - c, b); its right child is (b, 3*b*c - a, c). The sequence is a triangle read by rows consisting of the middle element of each triple, which is always the largest element of the triple. Row r contains 2^r elements.
The tree contains contains exactly one representative of each class of permutation equivalent nonsingular solutions to Markov's equation, a^2 + b^2 + c^2 = 3 * a * b * c. Nonsingular solutions are those in which a, b, and c are three distinct numbers. The two singular triples (1, 1, 1) and (1, 2, 1) are omitted in this sequence.
A consequence of Markov's equation is that the recurrence for the tree may be reformulated as follows: the left child of (a, b, c) is (a, (a^2 + b^2) / c, b); its right child is (b, (b^2 + c^2) / a, c).
An open problem is to prove the uniqueness conjecture, which asserts that the largest element of a triple determines the other two.
Frobenius proposed assigning a rational number index in (0,1) to each vertex of the tree, and hence to each term in this sequence. This is the Farey index, obtained by assigning the triple (0/1, 1/2, 1/1) to the root vertex and using the following rules to assign triples to the rest of the tree: the vertex labeled (u/v, w/x, y/z) with w = u + u and x = v + z has left child (u/v, (u+w)/(v+x), w/x) and right child (w/x, (w+y)/(x+z), y/z). The Farey index is the center element of the triple. Each rational number in (0, 1) appears as the Farey index of exactly one vertex of the tree. The index of a(n) is A007305(n+2) / A007306(n+2).
A sequence of leftward steps in the tree produces odd-indexed Fibonacci numbers, A001519, which have Farey indices of the form 1 / n. A sequence of rightward steps in the tree produces odd-indexed Pell numbers, A001653, which have Farey indices of the form (n - 1) / n. A sequence of leftward steps followed by a single rightward step produces A350922, corresponding to Farey indices of the form 2 / (2 * n + 1). Alternating steps right, left, right, left, right, ... produces A064098, which corresponds to Farey indices of the form F(n) / F(n + 1), where F(n) is the n-th Fibonacci number.

Examples

			The initial levels of the tree are as follows. (See p. 47 of Aigner's book.)
                               (1,5,2)
             (1,13,5)                              (5,29,2)
   (1,34,13)         (13,194,5)         (5,433,29)             (29,169,2)
(1,        (34,     (13,     (194,    (5,       (433,       (29,       (169,
 89,        1325,    7561,    2897,    6466,     37666,      14701,     985,
 ,34)        13)      194)     5)       433)      29)         169)       2)
		

References

  • Martin Aigner, Markov's theorem and 100 years of the uniqueness conjecture. A mathematical journey from irrational numbers to perfect matchings. Springer, 2013. x+257 pp. ISBN: 978-3-319-00887-5; 978-3-319-00888-2 MR3098784.

Crossrefs

Other presentations of the Markov numbers, Markov triples, or the Markov tree: A002559, A253809, A291694, A305313, A305314, A327345.
Subsequences in the Markov tree: A001519, A001653, A350922, A064098.
Farey tree: A007305, A007306.

Programs

  • Python
    def Mtree(x): return(x[0],(3*x[0]*x[1])-x[2],x[1]), (x[1],(3*x[1]*x[2])-x[0],x[2])
    def A368546_rowlist(maxrow):
        A,B = [[(1,5,2)]],[]
        for n in range(maxrow+1):
            A.append([])
            for j in A[n]:
                B.append(max(j))
                for k in Mtree(j):
                    A[n+1].append(k)
        return(B) # John Tyler Rascoe, Feb 09 2024
  • SageMath
    def stripUpToFirst1(w):
        x = w
        while x % 2 == 0:
            x = x // 2
        return(x // 2)
    def stripUpToFirst0(w):
        x = w
        while x % 2 == 1:
            x = x // 2
        if x == 0:
            return(None)
        else:
            return(x // 2)
    @CachedFunction
    def markovNumber(w):
        if w == None:
            return(2)
        elif w == 0:
            return(1)
        elif w == 1:
            return(5)
        elif w % 2 == 0:
            return(3*markovNumber(stripUpToFirst1(w))*markovNumber(w//2) - markovNumber(stripUpToFirst0(w//2)))
        else:
            return(3*markovNumber(stripUpToFirst0(w))*markovNumber(w//2) - markovNumber(stripUpToFirst1(w//2)))
    [markovNumber(w) for w in range(1,38)]
    

A365687 a(n) = number of fractions m/n, 0 <= m < n, gcd(m,n) = 1 whose partial fraction decomposition has integer part 0.

Original entry on oeis.org

1, 1, 2, 2, 4, 1, 6, 4, 6, 2, 10, 2, 12, 3, 4, 8, 16, 3, 18, 4, 6, 5, 22, 4, 20, 6, 18, 6, 28, 0, 30, 16, 10, 8, 12, 6, 36, 9, 12, 8, 40, 1, 42, 10, 12, 11, 46, 8, 42, 10, 16, 12, 52, 9, 20, 12, 18, 14, 58, 2, 60, 15, 18, 32, 24, 1, 66, 16, 22, 2, 70, 12, 72, 18, 20
Offset: 1

Author

William P. Orrick, Sep 15 2023

Keywords

Comments

If n = p_1^r_1 p_2^r_2 ... p_k^r_k where p_1, ..., p_k are distinct primes, the partial fraction decomposition of m/n has the form z + Sum_{i=1..k} a_i / p_i^r_i = z + Sum_{i=1..k} Sum_{j=1..r_i} a_ij / p_i^j where 0 < a_i < p_i^r_i, gcd(a_i,p_i) = 1, and where 0 <= a_ij < p_i (a_ij > 0 when j = r_i). z is the integer part.
If 0 < m / n < 1 and gcd(m,n) = 1 then the integer part satisfies 1 - k <= z <= 0.
If n is a nonhyperbolic number, which means that Sum_{i=1..k} 1 / p_i^r_i >= 1, then the integer part is not zero for any m/n with 0 < m / n < 1, gcd(m,n) = 1.
Assuming gcd(m,n) = 1, if m/n has integer part z then (n - m)/n has integer part 1 - k - z. This means there are equally many reduced fractions with denominator n > 1 between 0 and 1 having integer part z and integer part 1 - k - z.

Examples

			a(10) = 2 because, of the four nonnegative reduced fractions less than 1 with denominator 10, two of them (7/10 and 9/10) have integer part 0:
   1/10 = -1 + 1/2 + 3/5
   3/10 = -1 + 1/2 + 4/5
   7/10 = 1/2 + 1/5
   9/10 = 1/2 + 2/5.
		

Crossrefs

Programs

  • SageMath
    def a(n):
        b = n
        fs = factor(b)
        # bzList will hold Bezout coefficients to express 1/n as combination
        # of the reciprocals of the prime power factors of n. ppList will
        # hold the prime power factors themselves.
        bzList = []
        bz0 = 1
        ppList = []
        for f in fs:
            q = f[0]^f[1]
            ppList.append(q)
            b = b / q
            bzThis = xgcd(q,b)
            bzList.append(bz0*bzThis[2])
            bz0 = bz0 * bzThis[1]
        ct = 0
        for j in n.coprime_integers(n):
            if sum(floor(j*bzList[i]/ppList[i])\
             for i in range(len(ppList))) == 0:
                ct = ct + 1
        return(ct)

Formula

a(n) = phi(n) if n is a prime power.
a(n) = phi(n) / 2 if n is the product of powers of two distinct primes.
a(n) = 0 if n is in A181629, that is, if n is a nonhyperbolic number.
a(n) = A070306(n) if n is a prime power or a product of powers of two distinct primes.

A360441 Triangle read by rows: T(n,k) is the number of pairs (c,m), where c is a covering of the 1 X (2n) grid with 1 X 2 rectangles and equal numbers of red and blue 1 X 1 squares and m is a matching between red squares and blue squares, such that exactly k matched pairs are adjacent.

Original entry on oeis.org

1, 1, 2, 7, 8, 4, 71, 78, 36, 8, 1001, 1072, 504, 128, 16, 18089, 19090, 9080, 2480, 400, 32, 398959, 417048, 199980, 56960, 10320, 1152, 64, 10391023, 10789982, 5204556, 1523480, 295120, 38304, 3136, 128, 312129649, 322520672, 156264304, 46629632, 9436000, 1336832, 130816, 8192, 256
Offset: 0

Author

William P. Orrick, Mar 08 2023

Keywords

Comments

If row elements are divided by row sums, one obtains a probability distribution that approaches a Poisson distribution with expected value 1 as n approaches infinity.

Examples

			Triangle begins:
         1
         1        2
         7        8       4
        71       78      36       8
      1001     1072     504     128     16
     18089    19090    9080    2480    400    32
    398959   417048  199980   56960  10320  1152   64
  10391023 10789982 5204556 1523480 295120 38304 3136 128
		

Crossrefs

Column 1 is |A002119|.
T(n,k) equals 2^k times the corresponding element of the triangle of A168422.
Sum of row n is A001517(n).
Cf. A253667.

Programs

  • SageMath
    def T(n,k):
        return(2^k*sum((-1)^(j-k)*binomial(2*n-j,n)*binomial(n,j)\
         *binomial(j,k)*factorial(n-j) for j in range(k,n+1)))

Formula

T(n,k) equals 2^k times the corresponding element of the triangle of A168422.
T(n,k) = 2^k * Sum_{j=k..n} (-1)^(j-k) * C(2*n-j,n) * C(n,j) * C(j,k) * (n-j)!.
Recurrence: T(n,k) = (1/k!) * Sum_{j=0..k} T(n-j,0) * (-1)^j * C(k,j) * Sum_{t=0..min(j,k-j)} (-1)^(j-t) * C(j,t) * (k-j)! / (k-j-t)!
= (1/k!) * Sum_{j=0..k} T(n-j,0) * (-1)^j * C(k,j) * R(k,j) where R(k,j) is an element of the triangle of A253667.

A360440 Square array read by antidiagonals upwards: T(n,k), n>=0, k>=0, is the number of ways of choosing nonnegative numbers for k indistinguishable A063008(n)-sided dice so that it is possible to roll every number from 0 to (A063008(n))^k-1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 7, 15, 1, 1, 1, 1, 10, 71, 105, 1, 1, 1, 1, 42, 280, 1001, 945, 1, 1, 1, 1, 115, 3660, 15400, 18089, 10395, 1, 1, 1, 1, 35, 20365, 614040, 1401400, 398959, 135135, 1, 1
Offset: 0

Author

William P. Orrick, Feb 19 2023

Keywords

Comments

The number of configurations depends on the number of sides on the dice only through its prime signature. A063008 provides a canonical representative of each prime signature.
Also the number of Krasner factorizations of (x^(A063008(n))^k)-1) / (x-1) into k polynomials each having A063008(n) nonzero terms all with coefficient +1. (Krasner and Ranulac, 1937)

Examples

			A063008(2) = 4. There are 3 ways to assign numbers to two 4-sided dice:
 {{0, 1, 2, 3}, {0, 4, 8, 12}},
 {{0, 1, 8, 9}, {0, 2, 4,  6}},
 {{0, 1, 4, 5}, {0, 2, 8, 10}}.
The table begins:
  1  1    1       1          1              1                 1  ...
  1  1    1       1          1              1                 1  ...
  1  1    3      15        105            945             10395  ...
  1  1    7      71       1001          18089            398959  ...
  1  1   10     280      15400        1401400         190590400  ...
  1  1   42    3660     614040      169200360       69444920160  ...
  1  1  115   20365    6891361     3815893741     3141782433931  ...
  1  1   35    5775    2627625     2546168625     4509264634875  ...
  1  1  230  160440  299145000  1175153779800  8396156461492800  ...
  ...
The rows shown enumerate configurations for dice of 1, 2, 4, 6, 8, 12, 30, 16, and 24 sides, which represent the prime signatures {}, {1}, {2}, {1,1}, {3}, {2,1}, {1,1,1}, {4}, and {3,1}.
		

Crossrefs

The concatenation of all prime signatures, listed in the order that corresponds to the rows of T(n,k), is A080577.
T(3,k) is |A002119(k)|. Starting with k = 1, T(1,k), T(2,k), T(4,k), and T(7,k) are given by columns 1-4 of A060540.
Row n is row A063008(n) of A360098.

Programs

  • SageMath
    @cached_function
    def r(i,M):
        kminus1 = len(M)
        u = tuple([1 for j in range(kminus1)])
        if i > 1 and M == u:
            return(1)
        elif M != u:
            divList = divisors(i)[:-1]
            return(sum(r(M[j],tuple(sorted(M[:j]+tuple([d])+M[j+1:])))\
             for d in divList for j in range(kminus1)))
    def f(n,k):
        if n == 1 or k == 0:
            return(1)
        else:
            return(r(n,tuple([n for j in range(k-1)]))) / factorial(k-1)
    # The following function produces the top left corner of the table:
    def TArray(maxn,maxk):
        retArray = []
        primesList = []
        ptnSum = 0
        ptnItr = Partitions(ptnSum)
        ptn = ptnItr.first()
        n = 0
        while n <= maxn:
            if ptn == None:
                primesList.append(Primes()[ptnSum])
                ptnSum = ptnSum + 1
                ptnItr = Partitions(ptnSum)
                ptn = ptnItr.first()
            prdct = prod(primesList[j]^ptn[j] for j in range(len(ptn)))
            retArray.append([f(prdct,k) for k in range(maxk+1)])
            n = n + 1
            ptn = ptnItr.next(ptn)
        return(retArray)

Formula

T(n,k) = f(A063008(n),k), where f(n,k) is the table given by A360098.

A360439 Square array read by antidiagonals upwards: T(n,k), n>=0, k>=0, is the number of ways of choosing nonnegative numbers for k indistinguishable (p^n*q)-sided dice so that it is possible to roll every number from 0 to (p^n*q)^k-1, where p and q are distinct primes.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 42, 71, 1, 1, 1, 230, 3660, 1001, 1, 1, 1, 1190, 160440, 614040, 18089, 1, 1, 1, 5922, 6387150, 299145000, 169200360, 398959, 1, 1, 1, 28644, 238504266, 127534407000, 1175153779800, 69444920160, 10391023, 1
Offset: 0

Author

William P. Orrick, Feb 18 2023

Keywords

Comments

Also the number of Krasner factorizations of (x^((p^n*q)^k)-1) / (x-1) into k polynomials each having p^n*q nonzero terms all with coefficient +1. (Krasner and Ranulac, 1937)

Examples

			For two ten-sided dice we have k = 2 and n = 1 since 10 = 2^1*5. The seven configurations are
  {{0,1,2,3,4,5,6,7,8,9}, {0,10,20,30,40,50,60,70,80,90}},
  {{0,1,2,3,4,50,51,52,53,54}, {0,5,10,15,20,25,30,35,40,45}},
  {{0,1,2,3,4,25,26,27,28,29}, {0,5,10,15,20,50,55,60,65,70}},
  {{0,1,10,11,20,21,30,31,40,41}, {0,2,4,6,8,50,52,54,56,58}},
  {{0,1,20,21,40,41,60,61,80,81}, {0,2,4,6,8,10,12,14,16,18}},
  {{0,1,2,3,4,10,11,12,13,14}, {0,5,20,25,40,45,60,65,80,85}},
  {{0,1,4,5,8,9,12,13,16,17}, {0,2,20,22,40,42,60,62,80,82}}.
Array begins:
  1  1      1           1                  1                         1  ...
  1  1      7          71               1001                     18089  ...
  1  1     42        3660             614040                 169200360  ...
  1  1    230      160440          299145000             1175153779800  ...
  1  1   1190     6387150       127534407000          6888547183518000  ...
  1  1   5922   238504266     49829456981304      36179571823974699120  ...
  1  1  28644  8507955456  18306027156441024  175934152220744900062080  ...
  ...
		

Crossrefs

For a table with the number of sides not restricted to the form p^n*q see A360098.
T(n,2) = A349427(n+1).
T(1,k) = |A002119(k)|.

Programs

  • SageMath
    def T(n,k):
        return(factorial(k*n)/factorial(n)^k/factorial(k)\
         * sum((-n)^(k-j)*binomial(n*k+j,j)*falling_factorial(k,j)\
         for j in range(k+1)))

Formula

T(n,k) = (n*k)!/((n!)^k*k!) * Sum_{j=0}^k (-n)^(k-j)*binomial(n*k+j,j)*k!/(k-j)!.
T(n,k) = A060540(k,n) * Sum_{j=0}^k (-n)^(k-j)*binomial(n*k+j,j)*k!/(k-j)! for n>=1, k>=1.

A360864 Number of unlabeled connected multigraphs with circuit rank n and degree >= 3 at each node, loops allowed.

Original entry on oeis.org

0, 3, 15, 111, 1076, 13870, 220520, 4185406, 92235118, 2314204852, 65129484278, 2032179006943, 69640160993587, 2600585852722150, 105127528809344785, 4574251821427917425, 213171992131468465801, 10593983324971249199532, 559293301762878627195807, 31259896932477899016109585, 1844062168535890557437809526
Offset: 1

Author

Keywords

Comments

The terms up to a(21) were computed by Michael Borinsky and Jos Vermaseren using a program written in FORM. The graphs enumerated by a(n) are called admissible (Borinsky and Vogtmann, 2023).

Crossrefs

Diagonal sums of A360862.

Formula

a(n) = Sum_{k>=1} A360862(n + k - 1, k).

A360098 Square array read by antidiagonals upwards: T(n,k) is the number of ways of choosing nonnegative numbers for k n-sided dice, k >= 0, n >= 1, so that summing the faces can give any integer from 0 to n^k - 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 15, 1, 1, 1, 1, 1, 7, 1, 105, 1, 1, 1, 1, 1, 1, 71, 1, 945, 1, 1, 1, 1, 1, 10, 1, 1001, 1, 10395, 1, 1, 1, 1, 1, 3, 280, 1, 18089, 1, 135135, 1, 1, 1, 1, 1, 7, 15, 15400, 1, 398959, 1
Offset: 1

Author

William P. Orrick, Jan 25 2023

Keywords

Comments

T(n,k) depends on n only through its prime signature. (See A118914.) For example, for fixed k, T(p*q^2,k) will be the same for any pair of distinct primes p and q. Hence we may define T(n,k) = R(s(n),k), where s(n) is the prime signature of n.
Also the number of Krasner factorizations of (x^(n^k)-1) / (x-1) into k polynomials each having n nonzero terms all with coefficient +1. (Krasner and Ranulac, 1937)

Examples

			There are 3 ways to assign numbers to two 4-sided dice:
 {{0, 1, 2, 3}, {0, 4, 8, 12}},
 {{0, 1, 8, 9}, {0, 2, 4,  6}},
 {{0, 1, 4, 5}, {0, 2, 8, 10}}.
The northwest corner of T(n,k) begins:
  1   1    1    1       1          1             1 ...     (s(1)  = {})
  1   1    1    1       1          1             1 ...     (s(2)  = {1})
  1   1    1    1       1          1             1 ...     (s(3)  = {1})
  1   1    3   15     105        945         10395 ...     (s(4)  = {2})
  1   1    1    1       1          1             1 ...     (s(5)  = {1})
  1   1    7   71    1001      18089        398959 ...     (s(6)  = {1,1})
  1   1    1    1       1          1             1 ...     (s(7)  = {1})
  1   1   10  280   15400    1401400     190590400 ...     (s(8)  = {3})
  1   1    3   15     105        945         10395 ...     (s(9)  = {2})
  1   1    7   71    1001      18089        398959 ...     (s(10) = {1,1})
  1   1    1    1       1          1             1 ...     (s(11) = {1})
  1   1   42  3660 614040  169200360   69444920160 ...     (s(12) = {1,2})
  1   1    1    1       1          1             1 ...     (s(13) = {1})
  1   1    7   71    1001      18089        398959 ...     (s(14) = {1,1})
  1   1    7   71    1001      18089        398959 ...     (s(15) = {1,1})
  1   1   35 5775 2627625 2546168625 4509264634875 ...     (s(16) = {4})
  ...
		

Crossrefs

For rows of index n = p^j, p prime, or equivalently, for rows of signature {j} we have T(p^2,k) = R({2},k) = A001147(k), T(p^3,k) = R({3},k) = A025035(k), T(p^4,k) = R({4},k) = A025036(k), and, generally, T(p^j,k) = R({j},k) = the k-th element of the j-th column of the square array A060540.
For n = p * q, p and q distinct primes, we have T(p*q,k) = R({1,1},k) = |A002119(k)|.
Column correspondences are T(n,2) = A273013(n) and T(n,3) = A131514(n).
Cf. A118914.

Programs

  • SageMath
    @cached_function
    def r(i,M):
        kminus1 = len(M)
        u = tuple([1 for j in range(kminus1)])
        if i > 1 and M == u:
            return(1)
        elif M != u:
            divList = divisors(i)[:-1]
            return(sum(r(M[j],tuple(sorted(M[:j]+tuple([d])+M[j+1:]))) for d in divList for j in range(kminus1)))
    def T(n,k):
        if n == 1 or k == 0:
            return(1)
        else:
            return(r(n,tuple([n for j in range(k-1)]))) / factorial(k-1)

Formula

Use M to denote a (k-1)-element multiset of positive integers. Let U denote the (k-1)-element multiset whose elements all equal 1 and let N denote the (k-1)-element multiset whose elements all equal n. For i in M, let M_{i,j} denote the result of replacing i with j in M. Then T(1,k) = T(n,0) = 1, while for n > 1 and k > 0 we have T(n,k) = r(n,N) / (k-1)! where r(i,M) is given by the recurrence
r(i,U) = 1 for i > 1,
r(i,M) = Sum_{m in M} Sum_{d|i,d