hdu 找新朋友

来源:互联网 发布:网络基础技术视频教程 编辑:程序博客网 时间:2024/05/16 00:56
题目:
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
22560824027
 
Sample Output
768016016
 
Author
SmallBeer(CML)

算法:
1.预处理求出所有数的质数
2.根据欧拉函数E(N) = N * (1- 1/P1)(1-1/P2)...(1-1/Pm) (P1,P2,..Pm都是N的质数因子),求出E(N)

#include <iostream>#include <algorithm>using namespace std; #define maxN 32768int nums[ maxN ] = { 0 };int primes[ maxN ][ 6 ] = { 0 };void pre( ){for( int i = 2; i < maxN; ++i ){if( !nums[ i ] ){for( int j = i; j < maxN; j += i ){primes[ j ][ nums[ j ]++ ] = i;}}}}int main( void ){pre( );int kase;scanf( "%d", &kase );while( kase-- ){int n;scanf( "%d", &n );int sum = n;int num = nums[ n ];//cout << "num: " << num << endl;for( int i = 0; i < nums[ n ]; ++i ){//cout << primes[ n ][ i ] << " factor " << endl;sum = sum * ( primes[ n ][ i ] - 1 ) / primes[ n ][ i ];//cout << "sum: " << sum << endl;}printf( "%d\n", sum );}}