bzoj2818

线性筛+线性欧拉函数。

参考：http://blog.csdn.net/popoqqq/article/details/38760601

（废柴也是木头!只要能烧!  ——江南 《龙族·火之晨曦》）

2015.10.29：

仅以此书献给所有有梦想的衰小孩，如果你知道去哪，全世界都会为你让路。 ——江南 《龙族》

（矫情，删了 ）

//好搜百科：互质（relatively prime）又叫互素。若N个整数的最大公因数是1，则称这N个整数互质。
#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define N 10000010

bool isprime[N];
int primenum[N];
int cou;
int ef[N];
long long int efsum[N];

void linear_sieve(){//1余1互素没有计数计进去
    cou=0;
    memset(isprime,true,sizeof(isprime));
    isprime[0]=false;
    isprime[1]=false;
    for(int i=2;i<N;i++){
        if(isprime[i]){
            primenum[cou++]=i;
            ef[i]=i-1;
        }

        for(int j=0;j<cou&&primenum[j]*i<N;j++){
            isprime[i*primenum[j]]=false;
            if(i%primenum[j]==0){
                ef[i*primenum[j]]=primenum[j]*ef[i];
                break;
            }
            else{
                ef[i*primenum[j]]=(primenum[j]-1)*ef[i];//所以求欧拉函数时不一定要按照因数分解一样分解的那么彻底，只要分解成两两互素就可以了，有时这样用着更方便。
            }
        }
    }

    ef[0]=0;
    ef[1]=0;
    efsum[0]=0;
    for(int i=1;i<N;i++){
        efsum[i]=efsum[i-1]+ef[i];
    }

    return;
}

int main(){
    linear_sieve();

    int n;

    while(scanf("%d",&n)!=EOF){
        long long int ans=0;
        for(int i=0;i<cou&&primenum[i]<n;i++){
            ans=ans+efsum[n/primenum[i]]*2+1;
        }
        printf("%lld\n",ans);
    }

    return 0;
}