A365303 a(n) is the smallest nonnegative integer such that the sum of any seven ordered terms a(k), k<=n (repetitions allowed), is unique.

0, 1, 8, 57, 256, 1153, 4181, 14180, 47381, 115267, 307214, 737909, 1682367, 3850940, 8557010, 18311575, 37925058, 61662056

Offset

1,3

Comments

This is the greedy B_7 sequence.

Links

Table of n, a(n) for n=1..18.

J. Cilleruelo and J Jimenez-Urroz, B_h[g] sequences, Mathematika (47) 2000, pp. 109-115.

Melvyn B. Nathanson, The third positive element in the greedy B_h-set, arXiv:2310.14426 [math.NT], 2023.

Melvyn B. Nathanson and Kevin O'Bryant, The fourth positive element in the greedy B_h-set, arXiv:2311.14021 [math.NT], 2023.

Kevin O'Bryant, A complete annotated bibliography of work related to Sidon sequences, Electron. J. Combin., DS11, Dynamic Surveys (2004), 39 pp.

Example

a(3) != 7 because 7+0+0+0+0+0+0 = 1+1+1+1+1+1+1.

Prog

(Python)
def GreedyBh(h, seed, stopat):
    A = [set() for _ in range(h+1)]
    A[1] = set(seed)    # A[i] will hold the i-fold sumset
    for j in range(2, h+1): # {2, ..., h}
        for x in A[1]:
            A[j].update([x+y for y in A[j-1]])
    w = max(A[1])+1
    while w <= stopat:
        wgood = True
        for k in range(1, h):
            if wgood:
                for j in range(k+1, h+1):
                    if wgood and (A[j].intersection([(j-k)*w + x for x in A[k]]) != set()):
                        wgood = False
        if wgood:
            A[1].add(w)
            for k in range(2, h+1): # update A[k]
                for j in range(1, k):
                    A[k].update([(k-j)*w + x for x in A[j]])
        w += 1
        return A[1]
GreedyBh(7, [0], 10000)
(Python)
from itertools import count, islice, combinations_with_replacement
def A365303_gen(): # generator of terms
    aset, alist = set(), []
    for k in count(0):
        bset = set()
        for d in combinations_with_replacement(alist+[k], 6):
            if (m:=sum(d)+k) in aset:
                break
            bset.add(m)
        else:
            yield k
            alist.append(k)
            aset |= bset
A365303_list = list(islice(A365303_gen(), 10)) # Chai Wah Wu, Sep 01 2023

Crossrefs

Row 7 of A365515.

Cf. A025582, A051912, A365300, A365301, A365302, A365304, A365305.

Sequence in context: A317040 A317694 A143570 * A296721 A283693 A297453

Adjacent sequences: A365300 A365301 A365302 * A365304 A365305 A365306

Keyword

nonn,more

Author

Kevin O'Bryant, Aug 31 2023

Extensions

a(13)-a(18) from Chai Wah Wu, Sep 13 2023

Status

approved