>>10
プログラムを組んで体感してみた。

> # 5^n の mod 2^m の周期が 2^(m-2) であることを示す
> library(gmp)
> f <- function(m=4,ps=5){ # ps周期表示
+ M=numeric()
+ for(i in 1:(2^(m-2)*ps)){
+ n=as.bigz(i)
+ M[i] = asNumeric(mod.bigz(as.bigz(5^n),2^m))
+ }
+ matrix(M,nrow=ps,byrow=TRUE)
+ }
> f(4)
[,1] [,2] [,3] [,4]
[1,] 5 9 13 1
[2,] 5 9 13 1
[3,] 5 9 13 1
[4,] 5 9 13 1
[5,] 5 9 13 1
> f(5)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,] 5 25 29 17 21 9 13 1
[2,] 5 25 29 17 21 9 13 1
[3,] 5 25 29 17 21 9 13 1
[4,] 5 25 29 17 21 9 13 1
[5,] 5 25 29 17 21 9 13 1
> f(6)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16]
[1,] 5 25 61 49 53 9 45 33 37 57 29 17 21 41 13 1
[2,] 5 25 61 49 53 9 45 33 37 57 29 17 21 41 13 1
[3,] 5 25 61 49 53 9 45 33 37 57 29 17 21 41 13 1
[4,] 5 25 61 49 53 9 45 33 37 57 29 17 21 41 13 1
[5,] 5 25 61 49 53 9 45 33 37 57 29 17 21 41 13 1