library(igraph)Warning: package 'igraph' was built under R version 4.3.3
n=100 #頂点数
N=200 #サンプル数
cs=seq(0.2,0.7,0.01)
ratio=c()
for (c in cs){
count=0
for (i in 1:N){
p=c*sqrt(log(n)/n)
g=sample_gnp(n,p)
if(diameter(g)<=4){
count=count+1
}
}
ratio=c(ratio,count/N)
}