hdoj 1286 找新朋友 【欧拉函数】

找新朋友

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 9370    Accepted Submission(s): 4961

Problem Description

新年快到了，“猪头帮协会”准备搞一个聚会，已经知道现有会员N人，把会员从1到N编号，其中会长的号码是N号，凡是和会长是老朋友的，那么该会员的号码肯定和N有大于1的公约数，否则都是新朋友，现在会长想知道究竟有几个新朋友？请你编程序帮会长计算出来。

 

Input

第一行是测试数据的组数CN（Case number，1<CN<10000），接着有CN行正整数N（1<n<32768），表示会员人数。

 

Output

对于每一个N，输出一行新朋友的人数，这样共有CN行输出。

 

Sample Input

  
  

   
   2
25608
24027
  
  

 

Sample Output

  
  

   
   7680
16016
  
  
  
  

   
    
  
  
  
  

   
    
  
  
  
  

   
    
  
  
  
  

   
   分析:
  
  
  
  

   
   
(1)   p^k型的欧拉函数的证明：
   
   
对于给定的一个素数p: φ(p)=p-1 那么容易证明φ(n)=p^k-p^(k-1)
   
   
已知少于或等于p^k的正整数的个数为p^k-1,其中和p^k不互质的正整数有{ p×1,p×2,...,p×(p^(k-1)-1)}，共计p^(k-1)-1个
   
   
故: φ(n) = p^k-1-(p^(k-1)-1)=p^k-p^(k-1)。
代码：
  
  
  
  

   
   
/*
打表筛选法。 

*/
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#define maxn 32780
using namespace std;
int oula[maxn];
int n;

void eular()
{
  memset(oula,0,sizeof(oula));
    oula[1]=1;
    for(int i=2;i<maxn;i++)
    if(!oula[i])
    for(int j=i;j<maxn;j+=i)
    {
     if(!oula[j])
     oula[j]=j;
     oula[j]=oula[j]/i*(i-1);
    }
}

int main()
{
    eular();
    int t;
    scanf("%d",&t);
    while(t--)
    {
      scanf("%d",&n);
      printf("%d\n",oula[n]);
    }
    system("pause");
     return 0;
}
/* 
直接求解法 数据范围在（0~32768之间） 
*/
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;

int eular(int n)
{
  int ret=n;
  for(int i=2;i<=sqrt(n);i++)
  {
    if(n%i==0)
    {
     ret=ret/i*(i-1);
     while(n%i==0)
     n/=i;
    }
  }
  if(n>1)
  ret=ret/n*(n-1);
  return ret;
}

int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
       int n;
      scanf("%d",&n);
      printf("%d\n",eular(n));
    }
    system("pause");
  return 0;
}