#----------- # ユークリッドの互助法により最大公約数を求める #----------- # 環境変数、初期変数の設定 options(scipen=100) #指数表記回避 #----------- gcd=function(x,y){ if(x>y){ #数が大きいほうをaにする a=x b=y }else{ a=y b=x } flag=0 #フラグ用(念のため) while(flag==0){ r=a%%b #余りを求める if(r==0){ #余りが0のときbがgcd flag=1 return(b) }else{ #余りが0でない場合は数字を置き換え a=b b=r } } }