(hdu step 2.1.6)找新朋友(欧拉函数的简单使用:求与n互质的元素的个数)
2015-01-31 11:56
423 查看
题目:
题目分析:
本题如果是单纯的暴力(就是用两个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
);
}
return 0;
}
找新朋友 |
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) |
Total Submission(s): 2788 Accepted Submission(s): 1307 |
[align=left]Problem Description[/align]新年快到了,“猪头帮协会”准备搞一个聚会,已经知道现有会员N人,把会员从1到N编号,其中会长的号码是N号,凡是和会长是老朋友的,那么该会员的号码肯定和N有大于1的公约数,否则都是新朋友,现在会长想知道究竟有几个新朋友?请你编程序帮会长计算出来。 |
[align=left]Input[/align]第一行是测试数据的组数CN(Case number,1<CN<10000),接着有CN行正整数N(1<n<32768),表示会员人数。 |
[align=left]Output[/align]对于每一个N,输出一行新朋友的人数,这样共有CN行输出。 |
[align=left]Sample Input[/align]2 25608 24027 |
[align=left]Sample Output[/align]7680 16016 |
[align=left]Author[/align]SmallBeer(CML) |
[align=left]Source[/align]杭电ACM集训队训练赛(VII) |
[align=left]Recommend[/align]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
);
}
return 0;
}
相关文章推荐
- (hdu step 7.2.2)GCD Again(欧拉函数的简单应用——求[1,n)中与n不互质的元素的个数)
- (hdu step 7.2.2)GCD Again(欧拉函数的简单应用——求[1,n)中与n不互质的元素的个数)
- HDU ACM Step 2.1.5 找新朋友(欧拉函数)
- 找新朋友 hdu 1286 欧拉函数的简单运用
- (step8.2.2)hdu 2509(Be the Winner——简单博弈)
- [Step By Step]使用SAP Business Objects Data Services将Mysql中的数据导入到SAP HANA中,并使用简单的Transformation
- 母函数简介及简单使用(附HDU 1284 钱币兑换问题 解题报告)
- hdu 1286 找新朋友(欧拉函数)
- 【Azure Services Platform Step by Step-第10篇】使用Blob Storage搭建简单网络硬盘
- (step7.2.1)hdu 1395(2^x mod n = 1——简单数论)
- HDU 1286.找朋友(简单的筛选法)
- SQL Server BI Step by Step 2--- 使用SSIS进行简单的数据导入导出
- (step3.2.1)hdu 1003(Max Sum——简单DP)
- (step3.2.2)hdu 1087(Super Jumping! Jumping! Jumping!——简单DP)
- (step7.2.1)hdu 1395(2^x mod n = 1——简单数论)
- hdu 1286 找新朋友 (容斥原理 || 欧拉函数)
- hdu 1286 找新朋友(欧拉函数)
- SQL Server BI Step by Step 2--- 使用SSIS进行简单的数据导入导出
- hdu 1286找新朋友 (简单数学题)
- hdu 1286 找新朋友(欧拉函数)