Akiyah

# http://r-statistics-fan.hatenablog.com/entry/2014/09/22/210455
is.sosu <- function(x){   #素数判定の自作関数作成。（2は判定できない）
  sum(x %% 2:sqrt(x) == 0) == 0
}

NUMS <- c(1,3,7,9)

result <- list()

Akiyah

●をクリックすると消えて右と下に分裂する。
左上6マスの白いブロックから●を追い出せ。

参考
バナッハ‐タルスキの密室: 『数学者シャーロック・ホームズ』増補・改題
http://www.amazon.co.jp/dp/4535787409

Akiyah

u <- c()
for (i in 1:10000) {
  j <- sample(data, 10000)
  u <- c(u, mean(j))
}
 
sqrt(var(u)) #=> 10.06759
 
quantile(u, c(0.05, 0.95))
#=>    5%   95% 

Akiyah

set.seed(0) # 同じ値を再現できるようにしておく

data <- c(
  rep(    100000,     10),
  rep(     10000,    100),
  rep(      1000,   1000),
  rep(       100,  10000),
  rep(         0, 100000 - 11110)
)

Akiyah

* {
  margin: 0;
  padding: 0;
  border: 0;
}

body {
  background: #ddd;
  font: 30px sans-serif;
}

Akiyah

* {
  margin: 0;
  padding: 0;
  border: 0;
}

body {
  background: #ddd;
  font: 30px sans-serif;
}

Akiyah

# plus
"%+%" <- function(f, g) {
  return(paste0("(", f, ") * (", g,")"))
}

# minus
m <- function(f) {
  paste0( "-(", f, ")" )
}

Akiyah

s <- 100
index_h <- 500^(((1:s)-1)/(s-1))
index <- as.integer(c(index_h, rev(index_h)))
for (i in 1:length(index)) {
  png(sprintf("plot_%04d.png", i), width=500, height=500)
  np <- index[i]
  n <- 1000
  p <- np/n
  q <- 1 - p
  x <- 0:n

Akiyah

jibanyan_equation <- function(x,y) {
  min(max(min(1-(x/108)^2-(y/94)^2,y),min(1-((abs(x)-119)/103)^2-((y-56)/86)^2,1-((abs(x)-15)/77)^2-((y-119)/100)^2),1-((abs(x)-42)/66)^2-(y/55)^2,min(55+y,51-abs(x),-y)),3*abs(y-100)-2*(x-75)) * 
  min(min(max(min(1-(x/106)^2-(y/92)^2,y),min(1-((abs(x)-119)/101)^2-((y-56)/84)^2,((abs(x)-99)/40)^2+((y-54)/86)^2-1,92-abs(x)),1-((abs(x)-42)/64)^2-(y/53)^2),min(((abs(x)-52)/26)^2+((y+28)/26)^2-1,((abs(x)-51)/13)^2+(y/13)^2-1,max(abs(x)-51,y))),abs(x/51+10/51*sin(abs(y/61.2)^(1.2)*pi*(7/2)))^(2/3)+abs(y/61.2)^(2/3)-1) *
  min(1-(x/32)^2-((y+30)/32)^2,1-((abs(x)+5)/22)^2-((y-18)/22)^2) *
  min(1-((abs(x)-18)/20)^2-((y+10)/20)^2,((abs(x)-20)/22)^2+((y+7)/20)^2-1) *
  (1-((abs(x)-51)/11)^2-(y/11)^2)
}

x <- seq(-150, 150, length=400)
z <- outer(x,x, Vectorize(jibanyan_equation))

Akiyah

"+" <- function(e1, e2=NA) {
    if (is.character(e1) & is.character(e2)) {
        return(paste0(e1, e2))
    }
    ifelse(is.na(e2), base::"+"(e1), base::"+"(e1, e2))
}

1+2+3
#=> [1] 6
"f"+"g"+"h"