#-----------
# 素数のリストを生成、テキスト形式で保存
#-----------
# 環境変数、初期変数の設定
sfold="D:\\R\\makeprime" #データ保存フォルダ
maxval=10000000 #探索範囲最大値
viewval=as.integer(maxval/10) #画面表示判定用
sfile=paste("primelist",maxval,".txt",sep="") #保存ファイル名
#----------
setwd(sfold) #カレントディレクトリの変更
arr_p  <- numeric(1) #素数を格納する配列（0を並べた配列を用意）
pn=1  #n番目の素数
arr_p[pn]<-2  #最初の素数は2
i<-3 #素数判定の探索初期値
while(i<maxval){ #探索範囲
  judmax=as.integer(sqrt(i)) #割れるかの判定範囲
  #cat(i , judmax,"\n")
  j=1 #素数で回すための変数初期化
  flag=0  #判定フラグ（0:割り切れない、1:割り切れた）
  while(arr_p[j]<=judmax  &&  flag==0){ #チェックすべき素数はjudmax以下で、まだ割切れていない
    if(i%%arr_p[j]==0){ #割り切れたら
      flag=1 #フラグを変更→ループから抜ける
    }
    j=j+1 #次の素数
  }
  if(flag==0){ #すべて割り切れなかったら新しい素数
    pn=pn+1
    arr_p[pn]=i
  }
  if((i+1)%%viewval==0){ #画面表示の条件
    cat("p(",pn,")=",arr_p[pn],"\n")
  }
  i=i+2 #次のチェック数字（奇数だけチェック）
}
write.table(arr_p,sfile,quote=FALSE,row.names=FALSE,col.names=FALSE) #データファイル保存
rm(list=ls(all=TRUE)) #全体の初期化