【HDU】-1286-找新朋友（欧拉函数）

找新朋友

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

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
  
  
  
  


  
  

题解：求一个数与他互质数的个数，正好用欧拉函数。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define CLR(a,b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define M 32768+10
int ant[M];
void Eurla()	//打表，求出所有的欧拉函数 
{
	for(int i=2;i<M;i++)
	{
		if(!ant[i])
		{
			for(int j=i;j<M;j+=i)
			{
				if(!ant[j])
					ant[j]=j;
				ant[j]=ant[j]/i*(i-1);
			}
		}
	}
}
int main()
{
	Eurla();
	int u;
	scanf("%d",&u);
	while(u--)
	{	
		int t;
		scanf("%d",&t);
		printf("%d\n",ant[t]);
	}
	return 0;
}

2017.3.10 更新 

求一个数N的欧拉函数值公式为：N*(M1-1)/M1*(M2-1)/M2......*(Mn-1)/M

其中M1...Mn为N的质因数

#include <cstdio>
#include <stack>
#include <queue>
#include <cmath>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
#define CLR(a,b) memset(a,b,sizeof(a))
#define INF 0x3f3f3f3f
#define LL long long
int Eular(int n)
{
	int ans=n;
	int endd=sqrt(n);
	for(int i=2;i<=endd;i++)
	{
		if(n%i==0)
		{
			ans=ans*(i-1)/i;
			while(n%i==0)
				n/=i;
		}
	}
	if(n!=1) 	
		ans=ans*(n-1)/n;
	return ans;
}
int main()
{
	int u;
	scanf("%d",&u);
	while(u--)
	{
		int n;
		scanf("%d",&n);
		printf("%d\n",Eular(n));
	}
	return 0;
}