Update cmd

Author: itsubaki

cmd/counting/main.go

@@ -76,17 +76,13 @@ func main() {
 	qsim.Swap(c...)
 	qsim.InvQFT(c...)
 
-	// results
-	N := 1 << len(r)
+	// estimate
+	N, Q := 1<<len(r), 1<<t
 	for _, state := range q.Top(qsim.State(c, r, s, a), top) {
-		phi := ldexp(state.Int()[0], -t)                 // phi = k / 2**t
+		phi := float64(state.Int()[0]) / float64(Q)      // phi = k / 2**t
 		theta := 2 * math.Pi * phi                       // theta = 2*pi*phi
 		M := float64(N) * math.Pow(math.Sin(theta/2), 2) // M = N*(sin(theta/2))**2
 
 		fmt.Printf("%v; phi=%.4f, theta=%.4f, M=%.4f, M'=%.4f\n", state, phi, theta, M, float64(N)-M)
 	}
 }
-
-func ldexp(a, b int) float64 {
-	return math.Ldexp(float64(a), b)
-}