欧拉函数的计算

#欧拉函数是小于n的正整数中与n的数的数目

#判里两个数是不是互质

is.prime <- function(n,m)

{
  flag <- 1
  temp1 <- min(m,n)
  if(temp1 ==1)
  {
    flag <- 1
  }
  else if(temp1 == 2)
  {
    if(max(n,m)%%2 == 0)
    {
      flag <- 0
    }
  }
  else
  {
    for(i in 2:temp1)
    {
      if(m%%i==0 & n%%i==0)
      {
        flag <- 0
      }
    }
  }
  return(flag)
}

#计算具体的数值

Euler <- function(n)
{
  re <- 0
  if(n ==2|n==1)
  {
    re <- 1
  }
  else
  {
    for(i in 1:(n-1))
    {
      temp1 <- is.prime(i,n)
      re <- re + temp1
    }
  }
  return(re)
}

n <- 100
a <- numeric(n)
for(i in 1:n)
{
  a[i] <- Euler(i)
}
plot(a)