找新朋友
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