3semi20250115

コード

\[ p=c\sqrt{\frac{\log{n}}{n}} \]

library(igraph)

n=100 #頂点数
N=400 #サンプル数
cs=seq(1,1.8,0.1) 

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)<=2){
      count=count+1
    }
  }
  ratio=c(ratio,count/N)  
}

グラフ

barplot(ratio,names.arg = cs)