HDU 1286 找新朋友

 
找新朋友

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)

Total Submission(s): 1487 Accepted Submission(s): 744

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的公因子，所以假设会长的号为n，新会员的号为m，那么m与n互质（即没有公因子）

•1.暴力的方法不行，计算一下复杂度，是O(nlogn),考虑到case数目，TLE一定的。

•2.容斥原理。基本想法是将n做一个素因子分解，考虑每个素因子在n之内可以整除的数的个数，采用容斥原理，计算出和n互素的个数。时间复杂度应该可以，程序复杂。

•3.euler函数。这个题就是裸的euler函数的应用。

         欧拉公式（数论）

欧拉φ函数：φ(n)是所有小于n的正整数里，和n互素的整数的个数。n是一个正整数。

•Phi（n）=n(1-1/p1) (1-1/p2)….. (1-1/pk)

•其中p1， p2 ，pk是n的所有素数因子

•Phi（n）：所有小于等于n的且与n互素的数的个数

第一种：

#include<iostream>
#include<cstring>
#include<cmath>
#define maxn 1000005
using namespace std;
int s[maxn];//用于储存素数因子
int prime(int n)//求n的素数因子 
{
    int cnt=0;
    if(n%2==0) s[cnt++]=2;
    while(n%2==0) n/=2;
    for(int i=3;i*i<=n;i+=2)
    {
        if(n%i==0)
        {
            s[cnt++]=i;
            n/=i;
        }
        while(n%i==0)
        n/=i;
    }
    if(n!=1)
    {
        s[cnt++]=n;
    }
    return cnt;
} 
//计算欧拉函数 
int euler(int n)
{
    int cnt=prime(n);
    double f=n;
    for(int i=0;i<cnt;i++)
    {
        f*=1.0*(s[i]-1)/s[i];
    }
    return (int)(f+0.5);
}
int main()
{
    int cn,n;
    cin>>cn;
    while(cn--)
    {
        memset(s,0,sizeof(s));
        cin>>n;
        cout<<euler(n)<<endl;
    }
    return 0;
}

第二种：先将素数打表，再挑出因子求欧拉函数 

#include<bits/stdc++.h>
using namespace std;
int p[16384]={0};//已经将偶数排除 
int prm[3600],pn;
//先素数打表
int sieve()
{
	int i,j,k;
	for(i=0;i<90;i++)
	{
		if(!p[i])
		{
			for(k=(i<<1)+3,j=k*i+i+k;j<16384;j+=k)
			{
				p[j]=1;
			}
		}
	}
	for(prm[i=j=0]=2;i<16384;++i)
	{
		if(!p[i])
		prm[++j]=(i<<2)+3;
	}
		return j+1;
} 
int phi(int n)
{
	for(int i=0;a=t=n;i<pn&&prm[i]*prm[i]<=t;i++)
	{
		if(t%prm[i])
		continue;
		a-=a/prm[i];
		while(t%prm[i]==0)
		t/=prm[i];//算法加速 
	}
	if(t!=1) a-=a/t;
	return a; 
}
int main()
{
	int t,n;
	pn=sieve();
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d",&n);
		printf("%d\n",phi(n));
	}
	return 0;
 } 

第三种：主要还是理解一下欧拉公式

#include<iostream>
#include<cmath>
using namespace std;
 
//初等数论的欧拉公式
//返回小于n的所有互质的元素的个数
int oula(int n)  
{   
	int t=(int) sqrt((float) n);
	int ret=1,i;
	for(i=2;i<=t;i++)    
		if(n%i==0)    
		{              
			n/=i;
			ret*=i-1;              
			while(n%i==0)  
			{
				n/=i;
				ret*=i;
			}
		}    
	if(n>1)           
		ret*=n-1;    
	return ret;
}
int main()
{    
	int n,t;    
	scanf("%d",&n);    
	while(n--)    
	{           
		scanf("%d",&t);           
		printf("%d\n",oula(t));    
	}    
	return 0;
}

 2019.8.6