from qiskit.opflow import Zero, One

state = One ^ Zero ^ One ^ Zero ^ Zero
print(state)

from qiskit.opflow import X, Z

pauli_piXZ = -(X ^ Z)
print(pauli_piXZ)

print(pauli_piXZ.to_matrix())

from qiskit.opflow import I, X, H, CX

op = (CX ^ I) @ (I ^ CX) @ (I ^ I ^ H)
print(op)

from qiskit.quantum_info import Statevector

qc = op.to_circuit()
sv = Statevector(qc)
print(sv.sample_counts(1000))