//============================================================================
// Name : POJ_2407_number_theroy.cpp
// Author : tiger
// Version :
// Copyright : 欧拉函数
// Description : 欧拉函数是积性函数——若m,n互质,φ(mn)=φ(m)φ(n),
// 若n是质数p的k次幂,φ(n)=p^k-p^(k-1)=(p-1)p^(k-1),因为除了p的倍数外,其他数都跟n互质。
//============================================================================
#include <iostream>
#include <cmath>
using namespace std;
#define size 5000000
bool isprime[size];
int prime[size/2];
int counts;
int ans[size/2];
//筛选法素数表
void initprime()
{
int i,j;
counts = 0;
prime[counts++] = 2;
memset(isprime,true,sizeof(isprime));
for(i = 3; i*i<size;i+=2)
{
if(isprime[i])
{
prime[counts++] = i;
for(j = i*i; j<size; j+=i)
isprime[j] = false;
}
}
}
int euler(int n)
{
memset(ans,0,sizeof(ans));
int i;
int ret = 1;
//素因子分解
for(i = 0; i < counts && n>1 ;i++)
{
while(n % prime[i] == 0)
{
ans[i]++;
n /= prime[i];
}
if(ans[i]) ret *= (int)ceil(pow((double)(prime[i]),(double)(ans[i])-1) )*(prime[i]-1);
}
if(n>1)
ret *= n-1;
return ret;
}
int main() {
int n;
freopen("in","r",stdin);
initprime();
while(scanf("%d",&n) && n)
{
printf("%d/n",euler(n));
}
return 0;
}
POJ_2407_欧拉函数
最新推荐文章于 2019-11-16 11:02:04 发布