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