#549984 (Go) No.194 フィボナッチ数列の理解(1)

提出ソース

結果

問題	No.194 フィボナッチ数列の理解(1)
ユーザー	aru aru
提出日時	2020-09-08 22:29:12
言語	Go (1.22.1)
結果	AC
実行時間	25 ms / 5,000 ms
コード長	2,799 bytes
コンパイル時間	10,453 ms
コンパイル使用メモリ	239,460 KB
実行使用メモリ	5,376 KB
最終ジャッジ日時	2024-05-07 16:05:32
合計ジャッジ時間	11,906 ms
ジャッジサーバーID （参考情報）	judge5 / judge4

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 5,248 KB
testcase_01	AC	1 ms 5,248 KB
testcase_02	AC	24 ms 5,248 KB
testcase_03	AC	3 ms 5,376 KB
testcase_04	AC	9 ms 5,376 KB
testcase_05	AC	8 ms 5,376 KB
testcase_06	AC	9 ms 5,376 KB
testcase_07	AC	15 ms 5,376 KB
testcase_08	AC	2 ms 5,376 KB
testcase_09	AC	11 ms 5,376 KB
testcase_10	AC	5 ms 5,376 KB
testcase_11	AC	5 ms 5,376 KB
testcase_12	AC	8 ms 5,376 KB
testcase_13	AC	4 ms 5,376 KB
testcase_14	AC	1 ms 5,376 KB
testcase_15	AC	19 ms 5,376 KB
testcase_16	AC	16 ms 5,376 KB
testcase_17	AC	5 ms 5,376 KB
testcase_18	AC	18 ms 5,376 KB
testcase_19	AC	25 ms 5,376 KB
testcase_20	AC	16 ms 5,376 KB
testcase_21	AC	15 ms 5,376 KB
testcase_22	AC	15 ms 5,376 KB
testcase_23	AC	2 ms 5,376 KB
testcase_24	AC	8 ms 5,376 KB
testcase_25	AC	8 ms 5,376 KB
testcase_26	AC	7 ms 5,376 KB
testcase_27	AC	9 ms 5,376 KB
testcase_28	AC	3 ms 5,376 KB
testcase_29	AC	14 ms 5,376 KB
testcase_30	AC	23 ms 5,376 KB
testcase_31	AC	1 ms 5,376 KB
testcase_32	AC	7 ms 5,376 KB
testcase_33	AC	11 ms 5,376 KB
testcase_34	AC	8 ms 5,376 KB
testcase_35	AC	7 ms 5,376 KB
testcase_36	AC	17 ms 5,376 KB
testcase_37	AC	2 ms 5,376 KB
testcase_38	AC	20 ms 5,376 KB
testcase_39	AC	8 ms 5,376 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

package main

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

func out(x ...interface{}) {
	fmt.Println(x...)
}

var sc = bufio.NewScanner(os.Stdin)

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

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

func getString() 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
}

const mod = int(1e9 + 7)

func solve1(N, K int, a []int) {
	K = min(1000000, K)
	sum := 0
	for i := 0; i < N; i++ {
		sum += a[i]
	}
	f := sum
	s := sum
	for i := N; i < K; i++ {
		f = sum
		s += sum
		s %= mod
		// out(f, s)
		x := a[i%N]
		a[i%N] = sum
		sum += sum - x
		sum %= mod
		if sum < 0 {
			sum += mod
		}
	}
	out(f, s)
}

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
}

func solve2(N, K int, x []int) {
	a := make([][]int, N)
	for i := 0; i < N; i++ {
		a[i] = make([]int, N)
	}
	for i := 0; i < N; i++ {
		a[0][i] = 1
	}
	for i := 1; i < N; i++ {
		a[i][i-1] = 1
	}
	a = powModMatrix(a, K-N)
	f := 0
	for i := 0; i < N; i++ {
		f += x[N-1-i] * a[0][i]
		f %= mod
	}

	a = make([][]int, N+1)
	for i := 0; i < N+1; i++ {
		a[i] = make([]int, N+1)
	}
	a[0][0] = 2
	a[0][N] = -1
	for i := 1; i < N+1; i++ {
		a[i][i-1] = 1
	}
	// for i := 0; i < N+1; i++ {
	// 	out(a[i])
	// }
	a = powModMatrix(a, K-N)
	y := make([]int, N+1)
	y[0] = 0
	for i := 1; i < N+1; i++ {
		y[i] = y[i-1] + x[i-1]
	}
	// out(y)
	s := 0
	for i := 0; i < N+1; i++ {
		s += y[N-i] * a[0][i]
		s %= mod
		if s < 0 {
			s += mod
		}
	}

	out(f, s)
}

func main() {
	sc.Split(bufio.ScanWords)
	N, K := getInt(), getInt()
	a := getInts(N)
	if K > 1000000 {
		solve2(N, K, a)
	} else {
		solve1(N, K, a)
	}
}

yukicoder

結果

テストケース

ソースコード