結果

問題 No.657 テトラナッチ数列 Easy
コンテスト
ユーザー aru aru
提出日時 2020-10-11 21:34:19
言語 Go
(1.23.4)
結果
AC  
実行時間 184 ms / 2,000 ms
コード長 2,354 bytes
コンパイル時間 11,956 ms
コンパイル使用メモリ 239,212 KB
実行使用メモリ 6,272 KB
最終ジャッジ日時 2024-07-20 17:38:33
合計ジャッジ時間 14,080 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 13
権限があれば一括ダウンロードができます

ソースコード

diff # 

package main

import (
	"bufio"
	"fmt"
	"os"
	"sort"
	"strconv"
)

var sc = bufio.NewScanner(os.Stdin)
var wr = bufio.NewWriter(os.Stdout)

func out(x ...interface{}) {
	fmt.Fprintln(wr, x...)
}

func getI() int {
	sc.Scan()
	i, e := strconv.Atoi(sc.Text())
	if e != nil {
		panic(e)
	}
	return i
}

func getF() float64 {
	sc.Scan()
	i, e := strconv.ParseFloat(sc.Text(), 64)
	if e != nil {
		panic(e)
	}
	return i
}

func getInts(N int) []int {
	ret := make([]int, N)
	for i := 0; i < N; i++ {
		ret[i] = getI()
	}
	return ret
}

func getS() string {
	sc.Scan()
	return sc.Text()
}

// min, max, asub, absなど基本関数
func max(a, b int) int {
	if a > b {
		return a
	}
	return b
}

func min(a, b int) int {
	if a < b {
		return a
	}
	return b
}

func asub(a, b int) int {
	if a > b {
		return a - b
	}
	return b - a
}

func abs(a int) int {
	if a >= 0 {
		return a
	}
	return -a
}

func lowerBound(a []int, x int) int {
	idx := sort.Search(len(a), func(i int) bool {
		return a[i] >= x
	})
	return idx
}

func upperBound(a []int, x int) int {
	idx := sort.Search(len(a), func(i int) bool {
		return a[i] > x
	})
	return idx
}

func main() {
	defer wr.Flush()
	sc.Split(bufio.ScanWords)
	sc.Buffer([]byte{}, 1000000)
	// this template is new version.
	// use getI(), getS(), getInts(), getF()
	A := make([][]int, 4)
	for i := 0; i < 4; i++ {
		A[i] = make([]int, 4)
	}
	A[0][0], A[0][1], A[0][2], A[0][3] = 1, 1, 1, 1
	A[1][0] = 1
	A[2][1] = 1
	A[3][2] = 1

	Q := getI()
	for i := 0; i < Q; i++ {
		n := getI()
		// out(n)
		if n <= 3 {
			out(0)
			continue
		}
		a := powModMatrix(A, n-4)
		// for i := 0; i < 4; i++ {
		// 	out(a[i])
		// }
		t1, t2, t3, t4 := 0, 0, 0, 1
		t := a[0][0]*t4 + a[0][1]*t3 + a[0][2]*t2 + a[0][3]*t1
		out(t % mod)

	}

}

const mod = 17

func mulMod(A, B [][]int) [][]int {
	H := len(A)
	W := len(B[0])
	K := len(A[0])
	C := make([][]int, W)
	for i := 0; i < W; i++ {
		C[i] = make([]int, W)
	}

	for i := 0; i < H; i++ {
		for j := 0; j < W; j++ {
			for k := 0; k < K; k++ {
				C[i][j] += A[i][k] * B[k][j]
				C[i][j] %= mod
			}
		}
	}

	return C
}

func powModMatrix(A [][]int, p int) [][]int {
	N := len(A)
	ret := make([][]int, N)
	for i := 0; i < N; i++ {
		ret[i] = make([]int, N)
		ret[i][i] = 1
	}

	for p > 0 {
		if p&1 == 1 {
			ret = mulMod(ret, A)
		}
		A = mulMod(A, A)
		p >>= 1
	}

	return ret
}
0