*HDU 1286，2824欧拉函数

//
#include<iostream>
#include<string>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<vector>
#include<iomanip>
using namespace std;
int main()
{
    int t,n,sum;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        sum=n;
        for(int i=;i<=sqrt(n);i++)
        {
            if(n%i==)
            {
                sum=sum/i*(i-);
                while(n%i==)
                n/=i;
            }
        }
        if(n>)
        sum=sum/n*(n-);
     printf("%d\n",sum);
    }
    return ;
}
//
#include<cstdio>
#include<cstring>
#define size 3000001
int euler[size];
void Init()
{
    memset(euler,,sizeof(euler));
    euler[]=;
    for(int i=; i<size; i++)
        if(!euler[i])
            for(int j=i; j<size; j+=i)
            {
                if(!euler[j])
                    euler[j]=j;
                euler[j]=euler[j]/i*(i-);//先进行除法是为了防止中间数据的溢出
            }
}
int main()
{
    int a,b;
    Init();
    while(scanf("%d%d",&a,&b)!=EOF)
    {
        long long ans=;
        for(int i=a;i<=b;i++)
        ans+=euler[i];
        printf("%lld\n",ans);
    }
    return ;
}
/*******************************************************/
模板：
()直接求小于或等于n,且与n互质的个数:
int Euler(int n)
{
    int ret=n;
    for(int i=;i<=sqrt(n);i++)
    if(n%i==)
    {
        ret=ret/i*(i-);//先进行除法防止溢出(ret=ret*(1-1/p(i)))
        while(n%i==)
        n/=i;
    }
    if(n>)
    ret=ret/n*(n-);
    return ret;
}
/********************************************************/
筛选模板:求[,n]之间每个数的质因数的个数
#define size 1000001
int euler[size];
void Init()
{
    memset(euler,,sizeof(euler));
    euler[]=;
    for(int i=;i<size;i++)
    if(!euler[i])
    for(int j=i;j<size;j+=i)
    {
        if(!euler[j])
        euler[j]=j;
        euler[j]=euler[j]/i*(i-);//先进行除法是为了防止中间数据的溢出
    }
}
/*****************************************************/
//比上面更快的方法
#include<cstdio>
using namespace std;
const int N = 1e6+ ;
int phi[N], prime[N];
int tot;//tot计数，表示prime[N]中有多少质数
void Euler(){
    phi[] = ;
    for(int i = ; i < N; i ++){
        if(!phi[i]){
            phi[i] = i-;
            prime[tot ++] = i;
        }
        for(int j = ; j < tot && 1ll*i*prime[j] < N; j ++){
            if(i % prime[j]) phi[i * prime[j]] = phi[i] * (prime[j]-);
            else{
                phi[i * prime[j] ] = phi[i] * prime[j];
                break;
            }
        }
    }
}int main(){
    Euler();
}