MAX_NUM_QUBITS = 10
NUM_QUBITS = 5
NUM_STATE_DIMS = 2**NUM_QUBITS
def comp_basis_states(num_qubits):
num_qb = min(num_qubits, MAX_NUM_QUBITS)
basis_states = []
for idx in range(2**num_qb):
state = format(idx, '0' + str(num_qb) + 'b')
basis_states.append(state)
return basis_states