```
// test project main.go
package main
import (
"fmt"
"math/big"
"os"
"runtime"
"sync"
"time"
)
var file, _ = os.Create("20-1.txt")
var round = 1
var input = 3
var size = 1048576
var block = 8192
var s_box = [16][16]int{
{16, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, //0
{0, 0, 0, 4, 0, 0, 0, 4, 0, 4, 0, 0, 0, 4, 0, 0}, //1
{0, 0, 0, 2, 0, 4, 2, 0, 0, 0, 2, 0, 2, 2, 2, 0}, //2
{0, 2, 0, 2, 2, 0, 4, 2, 0, 0, 2, 2, 0, 0, 0, 0}, //3
{0, 0, 0, 0, 0, 4, 2, 2, 0, 2, 2, 0, 2, 0, 2, 0}, //4
{0, 2, 0, 0, 2, 0, 0, 0, 0, 2, 2, 2, 4, 2, 0, 0}, //5
{0, 0, 2, 0, 0, 0, 2, 0, 2, 0, 0, 4, 2, 0, 0, 4}, //6
{0, 4, 2, 0, 0, 0, 2, 0, 2, 0, 0, 0, 2, 0, 0, 4}, //7
{0, 0, 0, 2, 0, 0, 0, 2, 0, 2, 0, 4, 0, 2, 0, 4}, //8
{0, 0, 2, 0, 4, 0, 2, 0, 2, 0, 0, 0, 2, 0, 4, 0}, //9
{0, 0, 2, 2, 0, 4, 0, 0, 2, 0, 2, 0, 0, 2, 2, 0}, //A
{0, 2, 0, 0, 2, 0, 0, 0, 4, 2, 2, 2, 0, 2, 0, 0}, //B
{0, 0, 2, 0, 0, 4, 0, 2, 2, 2, 2, 0, 0, 0, 2, 0}, //C
{0, 2, 4, 2, 2, 0, 0, 2, 0, 0, 2, 2, 0, 0, 0, 0}, //D
{0, 0, 2, 2, 0, 0, 2, 2, 2, 2, 0, 0, 2, 2, 0, 0}, //E
{0, 4, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 4}} //F
func SPN(a []big.Rat, b [][]big.Rat, base int) {
for k := base * block; k < (base+1)*block; k++ {
//fmt.Println(k)
if a[k].Cmp(big.NewRat(0, 1)) > 0 {
first := k / 65536
second := k / 4096 % 16
third := k / 256 % 16
fourth := k / 16 % 16
fifth := k % 16
for m := 0; m < 16; m++ {
for n := 0; n < 16; n++ {
for x := 0; x < 16; x++ {
for y := 0; y < 16; y++ {
for z := 0; z < 16; z++ {
per := int64(s_box[first][m] * s_box[second][n] * s_box[third][x] * s_box[fourth][y] * s_box[fifth][z])
if per > 0 {
temp := big.NewRat(per, 1048576)
temp.Mul(temp, &a[k])
tt := m*65536 + n*4096 + x*256 + y*16 + z
var pos [20]int
for i := 0; i < 20; i++ {
pos[i] = tt % 2
tt >>= 1
tt = 524288*pos[19] + 16384*pos[18] + 512*pos[17] + 16*pos[16] + 262144*pos[15] + 8192*pos[14] + 256*pos[13] + 8*pos[12] + 131072*pos[11] + 4096*pos[10] +
128*pos[9] + 4*pos[8] + 65536*pos[7] + 2048*pos[6] + 64*pos[5] + 2*pos[4] + 32768*pos[3] + 1024*pos[2] + 32*pos[1] + 1*pos[0]
b[base][tt].Add(temp, &b[base][tt])
}
}
}
}
}
}
}
}
}
}
func main() {
start := time.Now()
runtime.GOMAXPROCS(128)
var wg sync.WaitGroup
var a = make([]big.Rat, size)
var b = make([][]big.Rat, 128)
var c = make([]big.Rat, size)
var dot string
for i := 0; i < 128; i++ {
b[i] = make([]big.Rat, size)
}
a[input].SetFrac(big.NewInt(1), big.NewInt(1))
change := a[:]
fmt.Printf("%v", a[input].String())
for i := 0; i < round; i++ {
for j := 0; j < 128; j++ {
//fmt.Printf("%d", j)
wg.Add(1)
go func(wg *sync.WaitGroup, j int) {
//SPN(a, b[j], j)
SPN(a, b, j)
for i := block * j; i < block*(j+1); i++ {
a[i].SetFrac(big.NewInt(0), big.NewInt(1))
}
for i := 0; i < size; i++ {
c[i].Add(&c[i], &b[j][i])
b[j][i].SetFrac(big.NewInt(0), big.NewInt(1))
}
defer wg.Done()
}(&wg, j)
}
wg.Wait()
change = a
a = c
c = change
/*for i := 0; i < 128; i++ {
for j := 0; j < size; j++ {
a[j].Add(&a[j], &b[i][j])
}
}
for i := 0; i < 128; i++ {
for j := 0; j < size; j++ {
b[i][j].SetFrac(big.NewInt(0), big.NewInt(1))
}
}*/
}
for i := 1; i < size; i++ {
dot = a[i].FloatString(20)
file.WriteString(dot + "\r\n")
}
elapsed := time.Since(start)
fmt.Printf("\nTook %s\n", elapsed)
}
```
错误如下
有疑问加站长微信联系（非本文作者）
