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

A365314 Number of unordered pairs of distinct positive integers <= n that can be linearly combined using nonnegative coefficients to obtain n.

Original entry on oeis.org

0, 0, 1, 3, 6, 8, 14, 14, 23, 24, 33, 28, 52, 36, 55, 58, 73, 53, 95, 62, 110, 94, 105, 81, 165, 105, 133, 132, 176, 112, 225, 123, 210, 174, 192, 186, 306, 157, 223, 218, 328, 180, 354, 192, 324, 315, 288, 216, 474, 260, 383, 311, 404, 254, 491, 338, 511, 360
Offset: 0

Views

Author

Gus Wiseman, Sep 05 2023

Keywords

Comments

Is there only one case of nonzero adjacent equal parts, at position n = 6?

Examples

			We have 19 = 4*3 + 1*7, so the pair (3,7) is counted under a(19).
The a(2) = 1 through a(7) = 14 pairs:
  (1,2)  (1,2)  (1,2)  (1,2)  (1,2)  (1,2)
         (1,3)  (1,3)  (1,3)  (1,3)  (1,3)
         (2,3)  (1,4)  (1,4)  (1,4)  (1,4)
                (2,3)  (1,5)  (1,5)  (1,5)
                (2,4)  (2,3)  (1,6)  (1,6)
                (3,4)  (2,5)  (2,3)  (1,7)
                       (3,5)  (2,4)  (2,3)
                       (4,5)  (2,5)  (2,5)
                              (2,6)  (2,7)
                              (3,4)  (3,4)
                              (3,5)  (3,7)
                              (3,6)  (4,7)
                              (4,6)  (5,7)
                              (5,6)  (6,7)
		

Crossrefs

The unrestricted version is A000217, ranks A001358.
For all subsets instead of just pairs we have A365073, complement A365380.
For strict partitions we have A365311, complement A365312.
The case of positive coefficients is A365315, for all subsets A088314.
The binary complement is A365320, positive A365321.
For partitions we have A365379, complement A365378.
A004526 counts partitions of length 2, shift right for strict.
A007865 counts sum-free subsets, complement A093971.
A179822 and A326080 count sum-closed subsets.
A364350 counts combination-free strict partitions.
A364914/A365046 count combination-full subsets, complement A326083/A124506.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,0,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]];
    Table[Length[Select[Subsets[Range[n],{2}], combs[n,#]!={}&]],{n,0,30}]
  • Python
    from itertools import count
    from sympy import divisors
    def A365314(n):
        a = set()
        for i in range(1,n+1):
            if not n%i:
                a.update(tuple(sorted((i,j))) for j in range(1,n+1) if j!=i)
            else:
                for j in count(0,i):
                    if j > n:
                        break
                    k = n-j
                    for d in divisors(k):
                        if d>=i:
                            break
                        a.add((d,i))
        return len(a) # Chai Wah Wu, Sep 12 2023

A365315 Number of unordered pairs of distinct positive integers <= n that can be linearly combined using positive coefficients to obtain n.

Original entry on oeis.org

0, 0, 0, 1, 2, 4, 5, 8, 10, 12, 15, 18, 20, 24, 28, 28, 35, 37, 42, 44, 49, 49, 60, 59, 66, 65, 79, 74, 85, 84, 93, 93, 107, 100, 120, 104, 126, 121, 142, 129, 145, 140, 160, 150, 173, 154, 189, 170, 196, 176, 208, 193, 223, 202, 238, 203, 241, 227, 267, 235
Offset: 0

Views

Author

Gus Wiseman, Sep 06 2023

Keywords

Comments

We consider (for example) that 2x + y + 3z is a positive linear combination of (x,y,z), but 2x + y is not, as the coefficient of z is 0.

Examples

			We have 19 = 4*3 + 1*7, so the pair (3,7) is counted under a(19).
For the pair p = (2,3), we have 4 = 2*2 + 0*3, so p is counted under A365314(4), but it is not possible to write 4 as a positive linear combination of 2 and 3, so p is not counted under a(4).
The a(3) = 1 through a(10) = 15 pairs:
  (1,2)  (1,2)  (1,2)  (1,2)  (1,2)  (1,2)  (1,2)  (1,2)
         (1,3)  (1,3)  (1,3)  (1,3)  (1,3)  (1,3)  (1,3)
                (1,4)  (1,4)  (1,4)  (1,4)  (1,4)  (1,4)
                (2,3)  (1,5)  (1,5)  (1,5)  (1,5)  (1,5)
                       (2,4)  (1,6)  (1,6)  (1,6)  (1,6)
                              (2,3)  (1,7)  (1,7)  (1,7)
                              (2,5)  (2,3)  (1,8)  (1,8)
                              (3,4)  (2,4)  (2,3)  (1,9)
                                     (2,6)  (2,5)  (2,3)
                                     (3,5)  (2,7)  (2,4)
                                            (3,6)  (2,6)
                                            (4,5)  (2,8)
                                                   (3,4)
                                                   (3,7)
                                                   (4,6)
		

Crossrefs

The unrestricted version is A000217, ranks A001358.
For all subsets instead of just pairs we have A088314, complement A365322.
For strict partitions we have A088571, complement A088528.
The case of nonnegative coefficients is A365314, for all subsets A365073.
The (binary) complement is A365321, nonnegative A365320.
A004526 counts partitions of length 2, shift right for strict.
A007865 counts sum-free subsets, complement A093971.
A179822 and A326080 count sum-closed subsets.
A326083 and A124506 appear to count combination-free subsets.
A364350 counts combination-free strict partitions.
A364914 and A365046 count combination-full subsets.

Programs

  • Mathematica
    combp[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,1,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]];
    Table[Length[Select[Subsets[Range[n],{2}],combp[n,#]!={}&]],{n,0,30}]
  • Python
    from itertools import count
    from sympy import divisors
    def A365315(n):
        a = set()
        for i in range(1,n+1):
            for j in count(i,i):
                if j >= n:
                    break
                for d in divisors(n-j):
                    if d>=i:
                        break
                    a.add((d,i))
        return len(a) # Chai Wah Wu, Sep 13 2023

A365382 Number of relatively prime integer partitions with sum < n that cannot be linearly combined using nonnegative coefficients to obtain n.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 4, 4, 2, 4, 12, 8, 20, 11, 14, 26, 43, 19, 38, 53, 51, 48, 101, 48, 124, 96, 121, 159, 134, 103, 241, 261, 244, 175, 401, 229, 488, 358, 328
Offset: 0

Views

Author

Gus Wiseman, Sep 08 2023

Keywords

Examples

			The a(11) = 2 through a(18) = 8 partitions:
  (5,4)  .  (6,5)  (6,5)   (7,6)  (7,5)   (7,4)     (7,5)
  (7,3)     (7,4)  (8,5)   (9,4)  (7,6)   (7,6)     (8,7)
            (7,5)  (9,4)          (9,5)   (8,5)     (10,7)
            (8,3)  (10,3)         (11,3)  (8,7)     (11,4)
                                          (9,5)     (11,5)
                                          (9,7)     (12,5)
                                          (10,3)    (13,4)
                                          (11,4)    (7,5,5)
                                          (11,5)
                                          (13,3)
                                          (7,4,4)
                                          (10,3,3)
		

Crossrefs

Relatively prime partitions are counted by A000837, ranks A289509.
This is the relatively prime case of A365378.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A116861 and A364916 count linear combinations of strict partitions.
A364350 counts combination-free strict partitions, non-strict A364915.
A364839 counts combination-full strict partitions, non-strict A364913.

Programs

  • Mathematica
    combsu[n_,y_]:=With[{s=Table[{k,i},{k,Union[y]},{i,0,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]];
    Table[Length[Select[Join@@IntegerPartitions/@Range[n-1],GCD@@#==1&&combsu[n,#]=={}&]],{n,0,20}]
  • Python
    from math import gcd
    from sympy.utilities.iterables import partitions
    def A365382(n):
        a = {tuple(sorted(set(p))) for p in partitions(n)}
        return sum(1 for m in range(1,n) for b in partitions(m) if gcd(*b.keys()) == 1 and not any(set(d).issubset(set(b)) for d in a)) # Chai Wah Wu, Sep 13 2023

Extensions

a(21)-a(45) from Chai Wah Wu, Sep 13 2023
Showing 1-3 of 3 results.