QUBO 多項式の展開

quboの計算で、以下のように多項式を展開する必要があります。かつ、$q_n$ は、0または1なので、$q_n^2=q_n$と置き換えないといけません。もちろん、筆算でもいいのですが、このご時世ですから、プログラムで処理したいですよね。

$$(q_0 + q_1 – 1)^2 + (q_2 + q_3 – 1)^2= 2q_0q_1 -q_0-q_1+1 + 2q_2q_3 -q_2-q_3+1$$

そこで、sympy を使って、コーディングしてみました。