快速幂算法

Contents
说明

算法说明详见 TAOCP 4.6.3 Evaluation of Powers

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45

package main

import "fmt"

func power(x, n uint64) uint64 {
	var total uint64 = 1

	if n == 0 {
		return 1
	} else if n == 1 {
		return x
	} else {
		for n > 0 {
			if n & 1 == 1 {
				total *= x
			}
			x *= x
			n = n >> 1
		}
	}

	return total
}


func main() {
	var x, n, result uint64

	fmt.Print("please input x and n:")

	_, err := fmt.Scanf("%d %d", &x, &n)
	if err != nil {
		fmt.Println("error: ", err)
	}

	result = power(x, n)

	if result == 0 {
		fmt.Println("result is to large!")
	} else {
		fmt.Printf("The result is : %d\n", power(x, n))
	}

}
0%