- 461 名前：１３２人目の素数さん mailto:sage [2020/11/26(木) 21:03:20.05 ID:nTBSTLPq.net]
- 暇つぶしに反例探しのプログラムを組んで処理が終わらないことを体感してみた。
f <- function(x) length(unique(x))==5 # 異なる実数
g <- function(x){ # (1)どの実数も残りの4つの和より小さい
flg=FALSE
for(i in 1:5){
if(x[i] < sum(x[-i])){
flg=TRUE
}else{
flg=FALSE
break
}
}
return(flg)
}
# (2)任意に２数a,bを選んだときa<=bなら2a<=bである
'%=>%' = function(P,Q) !(P & !Q) # PならばQ
pm=gtools::permutations(5,2)
h <- function(x){
sub <- function(i){
a=x[i[1]]
b=x[i[2]]
(a <= b) %=>% (2*a <= b)
}
all(apply(pm,1,sub))
}
sim <- function(x) f(x) & g(x) & h(x)
flg=FALSE
while(!flg){
x=runif(5)
flg=sim(x)
}
|
 |
