- 651 名前：１３２人目の素数さん mailto:sage [2021/01/28(木) 07:45:46.75 ID:7P00nNRm.net]
- >>624
数値を変えても数えられるように関数化
f <- function(a,b,n){
sub <- function(x,y) a*x + b*y <= n
x=0:ceiling(n/a)
y=0:ceiling(n/b)
xy=expand.grid(x,y)
sum(mapply(sub,xy[,1],xy[,2]))
}
f(a=3,b=2,n=2008)
f(a=3,b=7,n=2021)
結果
> f(a=3,b=2,n=2008)
[1] 337010
> f(a=3,b=7,n=2021)
[1] 97778
