poj2478求（2~n）之间的欧拉函数之和

最新推荐文章于 2020-07-17 07:59:37 发布

&小鹏鹏

最新推荐文章于 2020-07-17 07:59:37 发布

阅读量405

点赞数

本文链接：https://blog.csdn.net/u012785169/article/details/44839031

版权

/*
求Fi集合中元素的个数，其中Fi集合的元素满足下列条件
形如a/b，且0<a<b<=i, gcd(a,b)=1
欧拉函数表示{0……n}中与n互质的整数的个数，1,2的欧拉函数是1；
当n为素数是欧拉函数是n-1；
当n是合数是欧拉函数是<n-1;
若m，n互质时则，m*n的欧拉函数的值是m的欧拉函数值乘以n的欧拉函数值；
质数p的正整数k次幂的欧拉函数是P^k-P^(k-1)；=（p-1)*p^(k-1);

*/
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int phi[1000005];
int main()
{
    int i , j ;
    memset ( phi , 0 ,sizeof ( phi ) ) ;
    phi[1]=1;
    for ( i = 2 ; i <= 1000000 ; i ++ )
    {//筛选求phi
        if ( ! phi [i] )
        {
            for ( j = i ; j <= 1000000 ; j += i )
            {
                if ( ! phi [j] )
                    phi [j ] = j ;
                phi [j] = phi [j] / i * ( i - 1 ) ;
            }
        }
    }


    int n ;
    while ( cin >> n  , n )
    {
        __int64 sum = 0 ;
        for ( i = 2 ; i <= n ; i ++ )
            sum += phi [i] ;
        cout<<sum<<endl;
    }
    return 0;
}
/*
/*******************************************************************************
 # Author : Neo Fung
 # Email : neosfung@gmail.com
 # Last modified: 2012-07-17 18:57
 # Filename: POJ2478 Farey Sequence.cpp
 # Description : 
 ******************************************************************************/
#ifdef _MSC_VER
#define DEBUG
#define _CRT_SECURE_NO_DEPRECATE
#endif

#include <fstream>
#include <stdio.h>
#include <iostream>
#include <string.h>
#include <string>
#include <limits.h>
#include <algorithm>
#include <math.h>
#include <numeric>
#include <functional>
#include <ctype.h>
using namespace std;

const int kMAX=1000010;
const double kEPS=10E-6;
int prime[kMAX]={0},num_prime=0;//num_pirme记录素数个数
bool is_prime[kMAX];

void GetPrime(const int m)
{
  memset(is_prime,true,sizeof(is_prime));
  for(int i=2;i<m;i++)
  {
    if(is_prime[i])
      prime[num_prime++]=i;
    for(int j=0;j<num_prime && i*prime[j]<m;j++)
    {
      is_prime[i*prime[j]]=false;//合数标为1，同时，prime[j]是合数i*prime[j]的最小素因子
      if(!(i%prime[j]))//即比一个合数大的质数和该合数的乘积可用一个更大的合数和比其小的质数相乘得到
        break;
    }
  }	
}

long long dp[kMAX]={0ll};
void Solve(void)
{
	GetPrime(kMAX);
	for(int i=2;i<kMAX;++i)
	{
		if(is_prime[i])
			dp[i]=i-1;
		else
		{
			for(int j=0;j<num_prime;++j)
				if(i%prime[j]==0)
				{
					if((i/prime[j])%prime[j]==0)
						dp[i]=dp[i/prime[j]]*prime[j];
					else
						dp[i]=dp[i/prime[j]]*(prime[j]-1);
					break;
				}
		}
	}
	for(int i=3;i<kMAX;++i)
		dp[i]+=dp[i-1];
}

int main(void)
{
#ifdef DEBUG  
  freopen("../stdin.txt","r",stdin);
  freopen("../stdout.txt","w",stdout); 
#endif  

  int n,ncase=1;

  Solve();
  while(~scanf("%d",&n) && n)
  {
		printf("%lld\n",dp[n]);
  }

  return 0;
}
*/