from qiskit import QuantumRegister, ClassicalRegister, QuantumCircuit, execute
from qiskit import Aer

qreg = QuantumRegister(5, 'q')      # REPLACE 5 BY THE NUMBER OF QUBITS IN THE CIRCUIT
creg = ClassicalRegister(2, 'c')    # REPLACE 1 BY THE NUMBER OF OUTPUT QUBITS
circuit = QuantumCircuit(qreg, creg)
ONE = [0, 1]
ZERO = [1, 0]

#QUBITS INITIALIZATION
#Use some many lines as qubits you need
circuit.initialize(ZERO, 0)          # Set to ZERO or ONE depending on the desired input
circuit.initialize(ZERO, 1)
circuit.initialize(ONE, 2)

#Set all output qubits to ZERO
circuit.initialize(ZERO, 3)
circuit.initialize(ZERO, 4)         

#HERE YOU HAVE TO ADD YOUR QUANTUM OPERATIONS


#UPDATE TO MEASURE DE QUBITS
circuit.measure(qreg[3], creg[1])
circuit.measure(qreg[4], creg[0])


backend = Aer.get_backend('statevector_simulator')
job = backend.run(circuit, shots=100) 
result = job.result()
histogram = result.get_counts()
print(histogram)