(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;
}
时间: 2024-10-28 22:51:36

(hdu step 2.1.6)找新朋友(欧拉函数的简单使用:求与n互质的元素的个数)的相关文章

hdu 1286 找新朋友 (欧拉函数)

Description 新年快到了,“猪头帮协会”准备搞一个聚会,已经知道现有会员N人,把会员从1到N编号,其中会长的号码是N号,凡是和会长是老朋友的,那么该会员的号码肯定和N有大于1的公约数,否则都是新朋友,现在会长想知道究竟有几个新朋友?请你编程序帮会长计算出来. Input 第一行是测试数据的组数CN(Case number,1<CN<10000),接着有CN行正整数N(1<n<32768),表示会员人数. Output 对于每一个N,输出一行新朋友的人数,这样共有CN行输出

hdu-1286 找新朋友(欧拉函数,水题)

题目链接: 找新朋友 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 10120    Accepted Submission(s): 5344 Problem Description 新年快到了,“猪头帮协会”准备搞一个聚会,已经知道现有会员N人,把会员从1到N编号,其中会长的号码是N号,凡是和会长是老朋友的,那么该会员的号码肯定和N有

hdu1286找新朋友 欧拉函数模板题

#include<cstdio> #include<cstring> #include<iostream> using namespace std ; int Euler(int n) { int rea = n ; for(int i = 2;i*i <= n;i++) { if(n%i == 0) rea -= rea/i ; while(n%i == 0) n/=i ; } if(n>1) rea-=rea/n ; return rea ; } int

杭电ACM1286——找新朋友~欧拉函数的应用

欧拉函数 欧拉函数,对于正整数n,欧拉函数就是小于或者等于n的数中与n互质的数的数目. 通式为: F(x) = x * (1 - 1 / p1) * (1 - 1 / p2) * ....... * (1 - 1 / pn):p1,p2,p3.....pn为x的质因子.每一质因子只出现一次.即p1 ≠ p2 ≠ ....pn: 例如12 = 2 * 2 * 3:2只能算一次. 有了这个欧拉通式,便可以快速的求解这一题,该开始我是用普通的方法做的,很明显,超时了,一次输入的数很多个的时候,便容易超

HDU1286新朋友欧拉函数版

找新朋友 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 10401    Accepted Submission(s): 5493 Problem Description 新年快到了,“猪头帮协会”准备搞一个聚会,已经知道现有会员N人,把会员从1到N编号,其中会长的号码是N号,凡是和会长是老朋友的,那么该会员的号码肯定和N有大于1的公约

HDU 1286 找新朋友 (欧拉phi函数打表)

题意:你懂得. 析:一看这个题应该是欧拉phi函数,也就说欧拉phi函数是指求从 1 到 n 中与 n 互素的数的个数,这个题很明显是这个意思嘛,不多说了. 代码如下: #include <iostream> #include <cstdio> #include <cstring> using namespace std; const int maxn = 32768 + 5; int phi[maxn+5]; void phi_table(int n){ memset

(hdu step 7.2.2)GCD Again(欧拉函数的简单应用——求[1,n)中与n不互质的元素的个数)

题目: GCD Again Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 125 Accepted Submission(s): 84   Problem Description Do you have spent some time to think and try to solve those unsolved problem afte

(hdu step 7.2.1)The Euler function(欧拉函数模板题——求phi[a]到phi[b]的和)

题目: The Euler function Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 166 Accepted Submission(s): 96   Problem Description The Euler function phi is an important kind of function in number theory

HDU ACM 1286-找新朋友-欧拉phi函数

φ函数的值,通式:φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)-..(1-1/pn),其中p1, p2--pn为x的所有质因数,x是不为0的整数.φ(1)=1(唯一和1互质的数(小于等于1)就是1本身). #include<iostream> #include<cmath> using namespace std; int Get(int n) { int res,i; if(n==0) return 0; else { res=n; for(i=2