>>93
コードにコメントいれて総当たりした方法を解説

f <- function(nmax,N=13,showALL=FALSE){
library(numbers) # 約数や最小公倍数を計算するライブラリ
y=divisors(nmax) # nmaxまでの約数の数列y
y=y[y>(N-1)] # N(=13)以上の約数のみ
re=NULL # 答の格納場所
for(i in 1:length(y)){ # N以上の約数y[i]について
if((2*y[i]) %in% y){ # y[i]の2倍の数がyに含まれれば
i2=which(y==2*y[i]) # 何番目かをi2に収納
if((i2-i)==(N-1)){ # i2とiの差がN-1(=12)であれば
re=y[i:i2] # 答として格納i番目からi2番目を
cat(re,':','LCM =',mLCM(re),'\n') # その最小公倍数を返す
if(!showALL) break # showALLでなければ1つ表示してループからでる
}
}
}
invisible(re) # 答を返す
}