    c_square = c ** c

    linear = {}
    quadratic = {}
    offset = 0.0
    vartype = dimod.BINARY
    for index, value in enumerate(numbers):
        linear[index + 1] = value * (value - c)
    for index1, value1 in enumerate(numbers[:-1]):
        for index2 in range(index1 + 1, len(numbers)):
            value = value1 * numbers[index2]
            idx = (index1 + 1, index2 + 1)
            quadratic[idx] = quadratic[tuple(reversed(idx))] = value
