51nod 1188 最大公约数之和 V2

1188 最大公约数之和 V2

1. 2 秒
2.  
3. 262,144 KB
4.  
5. 160 分
6.  
7. 6 级题

给出一个数N，输出小于等于N的所有数，两两之间的最大公约数之和。

 

 

相当于计算这段程序（程序中的gcd(i,j)表示i与j的最大公约数）：

 

G=0;

for(i=1;i<N;i++)

for(j=i+1;j<=N;j++)

{

    G+=gcd(i,j);

}

 收起

输入

第1行：1个数T，表示后面用作输入测试的数的数量。(1 <= T <= 50000)
第2 - T + 1行：每行一个数N。(2 <= N <= 5000000)

输出

共T行，输出最大公约数之和。

输入样例

3
10
100
200000

输出样例

67
13015
143295493160

 

最这个题之前，首先要做一下上面的那个题，https://blog.csdn.net/qq_38185591/article/details/88297751

即51nod 1040 最大公约数之和  https://www.51nod.com/Challenge/Problem.html#!#problemId=1040

可以发现，斜率为1的那条线上都是从1到n之间组合gcd

1,2

1,3 2,3

1,4 2,4 3,4

 

那么ans[i] = ans[i - 1] + sum(i) （sum（i）代表的是1到i - 1与i的gcd的和）

在求sum的时候会发现，需要sqrt(n)的时间，所以现在n*sqrt(n)的复杂度我们是接收不了的。

开始考虑其他的

 

代码先贴，打完比赛再详细解释：

#include <bits/stdc++.h>
#include <time.h>
#define first fi
#define second se


using namespace std;

typedef long long ll;
typedef double db;
int xx[4] = {1,-1,0,0};
int yy[4] = {0,0,1,-1};
const double eps = 1e-9;
typedef pair<int,int>  P;
const int maxn = 5000010 ;
const ll mod = 1e9 + 7;
inline int sign(db a) {
    return a < -eps ? -1 : a > eps;
}
inline int cmp(db a,db b) {
    return sign(a - b);
}
ll mul(ll a,ll b,ll c) { ll res = 1; while(b) {  if(b & 1) res *= a,res %= c;  a *= a,a %= c,b >>= 1;  }  return res;}
//ll phi(ll x) {  ll res = x;  for(ll i = 2; i * i <= x; i++) { if(x % i == 0) res = res / i * (i - 1);   while(x % i == 0) x /= i;   }  if(x > 1) res = res / x  * (x - 1);    return res;}
ll c,n,k;
ll ans[maxn + 20];




int main() {
    long long a[maxn + 20];
memset(a,0,sizeof(a));
a[1]=1;
for(int i=2; i<maxn; i++)
    if(!a[i])
        for(int j=i; j<maxn; j+=i) {
            if(!a[j])
                a[j]=j;
            a[j]=a[j]/i*(i-1);
        }

    for(ll i = 1;i <= maxn;i++){
        ll anss = 0;
//        for(ll j = i + i;j < maxn;j += i){
//                    ans[j] +=  phi[j / i] * i;
//        }
        for(ll j = 2;i * j < maxn ;j++){
                    ans[i * j] += a[j] * i;
        }
    }
    for(int i = 1;i <= maxn;i++){
        ans[i] += ans[i - 1];
    }
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        printf("%lld\n",ans[n]);
//        cout << ans[n] << endl;
    }
//    cout << "time: " << (long long)clock() * 1000 / CLOCKS_PER_SEC << " ms" << endl;
    return 0;
}
//
//