21/02/19 08:00:47.01 OB1eKa3C.net
f <- function(n){
a=numeric()
i=1
a[i]=n
while(a[i]!=1){
a[i+1]=ifelse(a[i]%%3,a[i]+2,a[i]%/%3)
i=i+1
}
c(i,a[i])
}
f=Vectorize(f)
x=0:500000*2+1
y=t(f(x))
y[y[,1]==max(y[,1]),]
x[y[,1]==max(y[,1])]
21/02/19 08:00:47.01 OB1eKa3C.net
f <- function(n){
a=numeric()
i=1
a[i]=n
while(a[i]!=1){
a[i+1]=ifelse(a[i]%%3,a[i]+2,a[i]%/%3)
i=i+1
}
c(i,a[i])
}
f=Vectorize(f)
x=0:500000*2+1
y=t(f(x))
y[y[,1]==max(y[,1]),]
x[y[,1]==max(y[,1])]