library(gmp)func = function(n) { a = factorize(n) # 素因数分解 u = unique(a) # ユニークな素因数 b = numeric(length(u)) # 同じ素因数の個数 for (i in 1:length(a)) { suf = which(as.character(a[i]) == u) b[suf] = b[suf]+1 } return(prod(b+1)) #…