@staticmethod
def decomposition(wires):
    decomp_ops = [qml.Toffoli(wires=wires[1:]), qml.CNOT(wires=[wires[1], wires[2]]), qml.Toffoli(wires=[wires[0], wires[2], wires[3]])]
    return decomp_ops
@staticmethod
def decomposition(wires):
    decomp_ops = [qml.Toffoli(wires=wires[1:]), qml.CNOT(wires=[wires[1], wires[2]]), qml.Toffoli(wires=[wires[0], wires[2], wires[3]])]
    return decomp_ops