我来帮你一下吧，能跑出结果，但不一定是写得最简洁漂亮的：

for (i in 1:10000){
  a <- seq(1:10000)
  a[i+1] <- (4*a[i]-1)/5
  a[i+2] <- (4*a[i+1]-1)/5
  a[i+3] <- (4*a[i+2]-1)/5
  a[i+4] <- (4*a[i+3]-1)/5
  if(a[i+1]%%1==0 & a[i+2]%%1==0 & a[i+3]%%1==0 & a[i+4]%%1==0) {
    print(i)
  }
}


搜索范围1:10000（其实1:1000就够了），你还可以自己改。结果的话，你运行之后就知道第一个猴子拿的桃子数，剩下的就简单了。总桃子数是3121，五猴分别拿了：624，499，399，319，255。最后还剩1020。

改良一下，一次给出全部你要的结果：
store <- rep(0,7)
total <- 0
remainder <- 0
for (i in 1:1000){
  a <- seq(1:1000)
  a[i+1] <- (4*a[i]-1)/5
  a[i+2] <- (4*a[i+1]-1)/5
  a[i+3] <- (4*a[i+2]-1)/5
  a[i+4] <- (4*a[i+3]-1)/5
  if(a[i+1]%%1==0 & a[i+2]%%1==0 & a[i+3]%%1==0 & a[i+4]%%1==0) {
    print(i)
    total <- a[i]*5 + 1
    remainder <- total-a[i]-a[i+1]-a[i+2]-a[i+3]-a[i+4]-5
    store <- cbind(total,a[i],a[i+1],a[i+2],a[i+3],a[i+4],remainder)
  }

}
store

真厉害，非常感谢，认真学习。外加膜拜！

不客气啊。很高兴能帮上忙。欢迎常来论坛逛逛。

请问If条件语句里的命令具体是什么意思【 if(a[i+1]%%1==0 & a[i+2]%%1==0 & a[i+3]%%1==0 & a[i+4]%%1==0)】，是满足a[i+n]项是整数就执行if条件下的语句？？？？那我改成 if（is.integer(a[i+1]) & is.integer(a[i+2]) &is.integer(a[i+3]) & is.integer(a[i+4]))后为何无效？新手入门，许多东西不了解，希望赐教。