



UpperCamelCase_      			=  {
    "0": 1,
    "1": 1,
    "2": 2,
    "3": 6,
    "4": 24,
    "5": 120,
    "6": 720,
    "7": 5_040,
    "8": 40_320,
    "9": 362_880,
}

UpperCamelCase_      			=  {145: 145}

UpperCamelCase_      			=  {
    145: 0,
    169: 3,
    36_301: 3,
    1_454: 3,
    871: 2,
    45_361: 2,
    872: 2,
}







def 							A     (						__UpperCAmelCase      )			->			int:



						'''simple docstring'''






						if n in CACHE_SUM_DIGIT_FACTORIALS:
												return CACHE_SUM_DIGIT_FACTORIALS[n]
						UpperCAmelCase_					=	sum(DIGIT_FACTORIALS[let] for let in str(__UpperCAmelCase      )      )
						UpperCAmelCase_					=	ret
						return ret







def 							A     (						__UpperCAmelCase				,   __UpperCAmelCase = None      )			->			int:



						'''simple docstring'''






						UpperCAmelCase_					=	previous or set()
						if n in CHAIN_LENGTH_CACHE:
												return CHAIN_LENGTH_CACHE[n]
						UpperCAmelCase_					=	sum_digit_factorials(__UpperCAmelCase      )
						if next_number in previous:
												UpperCAmelCase_					=	0
												return 0
						else:
												previous.add(__UpperCAmelCase      )
												UpperCAmelCase_					=	1 + chain_length(__UpperCAmelCase				,   __UpperCAmelCase      )
												UpperCAmelCase_					=	ret
												return ret







def 							A     (						__UpperCAmelCase = 60				,   __UpperCAmelCase = 100_0000      )			->			int:



						'''simple docstring'''






						return sum(1 for i in range(1				,   __UpperCAmelCase      ) if chain_length(__UpperCAmelCase      ) == num_terms      )


if __name__ == "__main__":
							print(f"{solution() = }")










import numpy as np
import qiskit







def 							A     (						__UpperCAmelCase				,   __UpperCAmelCase      )			->			qiskit.QuantumCircuit:



						'''simple docstring'''






						UpperCAmelCase_					=	qiskit.QuantumCircuit(num_qubits + 1      )

						# First, let's deal with the case in which oracle is balanced
						if case == "balanced":
												# First generate a random number that tells us which CNOTs to
												# wrap in X-gates:
												UpperCAmelCase_					=	np.random.randint(1				,   2**num_qubits      )
												# Next, format 'b' as a binary string of length 'n', padded with zeros:
												UpperCAmelCase_					=	format(__UpperCAmelCase				,   f"0{num_qubits}b"      )
												# Next, we place the first X-gates. Each digit in our binary string
												# corresponds to a qubit, if the digit is 0, we do nothing, if it's 1
												# we apply an X-gate to that qubit:
												for index, bit in enumerate(__UpperCAmelCase      ):
																		if bit == "1":
																								oracle_qc.x(__UpperCAmelCase      )
        # Do the controlled-NOT gates for each qubit, using the output qubit
        # as the target:
												for index in range(__UpperCAmelCase      ):
																		oracle_qc.cx(__UpperCAmelCase				,   __UpperCAmelCase      )
												# Next, place the final X-gates
												for index, bit in enumerate(__UpperCAmelCase      ):
																		if bit == "1":
																								oracle_qc.x(__UpperCAmelCase      )

    # Case in which oracle is constant
						if case == "constant":
												# First decide what the fixed output of the oracle will be
												# (either always 0 or always 1)
												UpperCAmelCase_					=	np.random.randint(2      )
												if output == 1:
																		oracle_qc.x(__UpperCAmelCase      )

						UpperCAmelCase_					=	oracle_qc.to_gate()
						UpperCAmelCase_					=	'''Oracle'''  # To show when we display the circuit
						return oracle_gate







def 							A     (						__UpperCAmelCase				,   __UpperCAmelCase      )			->			qiskit.QuantumCircuit:



						'''simple docstring'''






						UpperCAmelCase_					=	qiskit.QuantumCircuit(num_qubits + 1				,   __UpperCAmelCase      )
						# Set up the output qubit:
						dj_circuit.x(__UpperCAmelCase      )
						dj_circuit.h(__UpperCAmelCase      )
						# And set up the input register:
						for qubit in range(__UpperCAmelCase      ):
												dj_circuit.h(__UpperCAmelCase      )
						# Let's append the oracle gate to our circuit:
						dj_circuit.append(__UpperCAmelCase				,   range(num_qubits + 1      )      )
						# Finally, perform the H-gates again and measure:
						for qubit in range(__UpperCAmelCase      ):
												dj_circuit.h(__UpperCAmelCase      )

						for i in range(__UpperCAmelCase      ):
												dj_circuit.measure(__UpperCAmelCase				,   __UpperCAmelCase      )

						return dj_circuit







def 							A     (						__UpperCAmelCase				,   __UpperCAmelCase      )			->			qiskit.result.counts.Counts:



						'''simple docstring'''






						UpperCAmelCase_					=	qiskit.Aer.get_backend('''aer_simulator'''      )

						UpperCAmelCase_					=	dj_oracle(__UpperCAmelCase				,   __UpperCAmelCase      )
						UpperCAmelCase_					=	dj_algorithm(__UpperCAmelCase				,   __UpperCAmelCase      )

						# Execute the circuit on the simulator
						UpperCAmelCase_					=	qiskit.execute(__UpperCAmelCase				,   __UpperCAmelCase				,   shots=1000      )

						# Return the histogram data of the results of the experiment.
						return job.result().get_counts(__UpperCAmelCase      )


if __name__ == "__main__":
							print(f"Deutsch Jozsa - Constant Oracle: {deutsch_jozsa('constant', 3)}")
							print(f"Deutsch Jozsa - Balanced Oracle: {deutsch_jozsa('balanced', 3)}")





