(hdu step 2.1.6)觅新朋友(欧拉函数的简单使用:求与n互质的元素的个数)
(hdu step 2.1.6)找新朋友(欧拉函数的简单使用:求与n互质的元素的个数)
题目:
找新朋友 |
| Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) |
| Total Submission(s): 2788 Accepted Submission(s): 1307 |
|
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 |
|
Author
SmallBeer(CML)
|
|
Source
杭电ACM集训队训练赛(VII)
|
|
Recommend
lcy
|
题目分析:
本题如果是单纯的暴力(就是用两个for循环求出数据规模里面从1~i中与i互质的元素的个数),超时超的很厉害(n本身就大,而且还用到了>O(n^2)的暴力算法)。这道题实际上是“求出从1~n中与n互质的元素的个数”的问题。这宗问题在数论中有专门的知识来解决——欧拉函数。
以下介绍一下这道题用到的基本知识:
1)积性函数:满足f(ab) = f(a)*f(b)的函数就积性函数,其中a与b是互质的任意两个数。若a,b为任意两个整数,则此时的f成为完全积性函数。对于积性函数,f(1)=1 。
若n为素数(质数),f(n) = n-1;
若n为合数,f(n) < n-1;
代码如下:
/*
* f1.cpp
*
* Created on: 2015年1月31日
* Author: Administrator
*/
#include <iostream>
#include <cstdio>
#include <cstring>
const int maxn = 32770;
bool u[maxn];
int su[maxn];
int f[maxn];
/**
* 欧拉函数:
* 用于求解:
* 求出从1~n中与n互质的元素的个数
*/
void oula(){
memset(u,true,sizeof(u));
f[1] = 1;
int i;
int j;
int num = 1;
for(i = 2 ; i < maxn ; ++i){
if(u[i] == true){
su[num++] = i;
f[i] = i-1;
}
for(j = 1 ; j < num ; ++j){
if(i*su[j] > maxn){
break;
}
u[i*su[j]] = false;
if(i%su[j]){
f[i*su[j]] = f[i] * (su[j]-1);
}else{
f[i*su[j]] = f[i] * su[j];
break;
}
}
}
}
int main(){
oula();
int t;
scanf("%d",&t);
while(t--){
int n;
scanf("%d",&n);
printf("%d\n",f[n]);
}
return 0;
}