欧拉phi函数

1.小于n且与n互素的整数个数。

给出正整数n的唯一分解式n=p1^a1*p2^a2*p3^a3...pk^ak，求1,2,3...n中与n互素的数的个数。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

int euler_phi(int n)
{
    int m=(int)sqrt(n+0.5);
    int ans=n;
    for(int i=2;i<=m;i++)
    {
        if(n%i==0)
        {
            ans=(ans/i)*(i-1);
            while(n%i==0) n/=i;
        }
    }
    if(n>1) ans=ans/n*(n-1);
    return ans;
}

(算法竞赛入门经典第二版P321)

2. 1-n中所有数的欧拉phi函数值

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

void phi_table(int n,int *phi)
{
    memset(phi,0,sizeof(phi));
    phi[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!phi[i])
        {
            for(int j=i;j<=n;j+=i)
            {
                if(!phi[j]) phi[j]=j;
                phi[j]=phi[j]/i*(i-1);
            }
        }
    }
}

uva10820

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int INF=0x3f3f3f3f;
const int MAXN=5e4+5;
int phi[MAXN];
int n;
ll sum[MAXN];

void phi_table()
{
	phi[1]=1;
	for(int i=2;i<MAXN;i++)
	{
		if(!phi[i])
		{
			for(int j=i;j<MAXN;j+=i)
			{
				if(!phi[j]) phi[j]=j;
				phi[j]=phi[j]/i*(i-1);
			}
		}
	}
}

int main()
{
	phi_table();
	for(int i=2;i<MAXN;i++)
	{
		sum[i]=sum[i-1]+phi[i];
	}
	while(~scanf("%d",&n) && n)
	{
		if(n==1)
		{
			printf("1\n");
			continue;
		 } 
		printf("%lld\n",sum[n]*2+1);
	}
	return 0;
}