欧拉函数FZU1969

（转）

/********************
题意：G=0;
            for(i=1;i<N;i++)
                  for(j=i+1;j<=N;j++)
                          G+=gcd(i,j);
求G

UESTC1723 吴神的烦劳
http://acm.uestc.edu.cn/problem.php?pid=1723
题目开始给出的公式可以转换为是给出一个n，求
for(m=1;m<=n;m++)
每一个m值与它（不包括自身）之前所有数的公约数之和。
显然直接对每一个m从1-(m-1)循环每次求出最大公约数的时间复杂度是O(n^2)肯定TLE.。
由于m从1循环到n，而求m与它之前的数最大公约数又是从1循环到(m-1)，其中有很多公约
数的重复的计算，那么可以换个角度考虑，考虑一个数d它可以做哪些数的最大公约数，则
枚举公约数d （1~n），d一共出现了sigma(phi(i))-1  (1<=i<=[N/d])次 ，
而sigma(phi(i))通过预处理可以O(1)得到。所以时间复杂度从O（n^2）降到了O（n），
UESTC  OJ上面此题说明了case最多只有100组，那么O(n)的复杂度便可以过了，
500+ms过了
时间复杂度O（n）

FZU1969
http://acm.fzu.edu.cn/problem.php?pid=1969
在这道题目中数据较多，在O（case*n）的情况下容易超时，因此不能用UESTC的代码

模拟素数筛选的过程，d可以作为后面某些数的最大公约数，则将d*phi[n/d]与n/d*phi[d]全都更新保存到
相应的数组里面，用O(sqrt(n))的效率预处理出1000000以内所有数的G值，最后O（1）
查询便可
时间复杂度：O(sqrt(n))

********************/
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define maxn 1000001
using namespace std;

int phi[maxn];
long long ans[maxn];
void euler()
{
    for(int i=1;i<=maxn;i++)phi[i]=i;
    for(int i=2;i<=maxn;i+=2)phi[i]/=2;
    for(int i=3;i<=maxn;i+=2)
    {
        if(phi[i]==i)
        {
            for(int j=i;j<=maxn;j+=i)
            {
                phi[j]=phi[j]/i*(i-1);
            }
        }
    }
    for(int i=2;i<=maxn;i++)
          ans[i]=phi[i];
    for(int i=2;i<=1000;i++)
    {
        ans[i*i]+=i*phi[i];
        for(int j=i*i+i,k=i+1;j<=maxn;j+=i,k++)
        {
            ans[j]+=phi[i]*k+phi[k]*i;
        }
    }
    for(int i=1;i<=maxn;i++)
          ans[i]+=ans[i-1];
}

int main()
{
    euler();
    int n;
    while(cin>>n,n)
    {
        cout<<ans[n]<<endl;
    }
    return 0;
}
/*
int main()//UESTC1723的主函数
{
    euler();
    int n,t;
    cin>>t;
    while(t--)
    {
        cin>>n;
        cout<<ans[n]<<endl;
    }
    return 0;
}
*/
/**
附O（n）的code
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define maxn 1000001
using namespace std;

int phi[maxn];
long long sum[maxn];
void euler()
{
    memset(sum,0,sizeof(sum));
    for(int i=1; i<=maxn; i++)phi[i]=i;
    for(int i=2;i<=maxn;i+=2)phi[i]/=2;
    for(int i=3; i<=maxn; i+=2)
    {
        if(phi[i]==i)
        {
            for(int j=i; j<=maxn; j+=i)
            {
                phi[j]=phi[j]/i*(i-1);
            }
        }
    }
    for(int i=1;i<=maxn;i++)
            sum[i]=sum[i-1]+phi[i];
}

long long ans;
int main()
{
    int t,n;
    cin>>t;
    euler();
    while(t--)
    {
        ans=0;
        cin>>n;
        for(int i=1;i<=n;i++)
        {
            ans+=i*(sum[n/i]-1);
        }
        cout<<ans<<endl;
    }
    return 0;
}
********************/