F - GCD - Extreme (II)欧拉函数

https://cn.vjudge.net/problem/18553/origin

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

ll a[maxn];
void make()
{
      for(ll i=0;i<maxn-3;i++)
            a[i]=i;
      for(ll i=2;i<maxn-4;i++)
      {
            if(a[i]==i)
            {
                  for(ll j=i;j<maxn-4;j+=i)
                  {
                        
                        a[j]=a[j]*(i-1)/i;
                  }
            }
      }
      // for(int i=1;i<100;i++)
      // printf("%d    %d\n",a[i],i );
      // cout<<endl;
}
long long  b[maxn]={0};
void houxu()
{
      for(int i=1;i<maxn-4;i++)///当gcd为i的时候
      {
            for(int j=i*2;j<maxn-4;j+=i)
            {
                  b[j]=b[j]+i*a[j/i];
            }
      }
      for(int i=1;i<maxn-4;i++)
      {
            b[i]+=b[i-1];
      }
}
int main()
{
      ll t;
      make();
      houxu();
      while (scanf("%lld",&t),t) {
            printf("%lld\n",b[t] );
      }

      return 0;
}