def rps(n):
resources = create_empty_resources()
# \epsilon = \delta^2 / 4
l = 2 + (2*n)
# Run PPE to get a random k
ppe_res = ppe(l = l, s = s, n = n)
# Solve for p given s and l
div_res = divmod(n = n)
resources = combine_resources(ppe_res, div_res)
# Apply A^p to enact desired phase shift
add_res = dkrs_add(n=n)
resources = combine_resources(resources, add_res)
# Add up the initial Hadamards to create \ket{\nu}
resources['single'] += 2
resources['depth'] += 1
return resources
from ppe import * s = 517 n = 24 print "s= " + str(s) print "n= " + str(n) for i in range(6,16): resources = ppe(s=s, n=n, l=i) print "l= " + str(i) print "resources= " + str(resources)
