【数学検定】数学検定（数検）総合スレッド Part.12

>>193
> S <- function(n){
+ s=numeric()
+ s[1]=3
+ if(n > 1){
+ for(i in 1:(n-1)) s[i+1]=3*s[i]+3^(i+1)
+ }
+ return(s[n])
+ }
> a <- function(n){
+ if(n==1) return(3)
+ else return(S(n)-S(n-1))
+ }
> sapply(1:20,a)
[1] 3 15 63 243 891
[6] 3159 10935 37179 124659 413343
[11] 1358127 4428675 14348907 46235367 148272039
[16] 473513931 1506635235 4778186031 15109399071 47652720147
>