SPOJ - GCDEX GCD Extreme 题解

最新推荐文章于 2023-08-25 11:46:11 发布

Gary_2005

最新推荐文章于 2023-08-25 11:46:11 发布

阅读量129

点赞数

分类专栏：题解

本文链接：https://blog.csdn.net/qq_42925924/article/details/108176180

版权

题解专栏收录该内容

119 篇文章 1 订阅

订阅专栏

题意：
有q此询问每次询问一个N，求G的值。

G = 0;
for (i = 1; i < N; i++)
  for (j = i+1; j <= N; j++) 
    G += gcd(i, j);

首先N的范围很小一看就知道是预处理询问。
$answer_i$ 表示N=i时的答案。 $answer_i=answer_{i-1}+\sum_{j=1}^{i-1} gcd(j,i)$ 。
这时就是需要求对于一个i，求出 $\sum_{j=1}^{i-1} gcd(j,i)$ 可以发现gcd(j,i)时i的一个因子，所以只需要枚举i的因子g。然后就是在 $[1, i - 1]$ 中，有多少j满足 $g|j,gcd(\frac{j}{g},\frac{j}{g})=1$ 也就是和 $\frac{j}{g}$ 互质的数的个数，显然这就是 $\phi$ 的定义了。

/*
{

AuThOr Gwj
*/
#pragma GCC optimize(2)
#include<bits/stdc++.h>
#define rb(a,b,c) for(int a=b;a<=c;++a)
#define rl(a,b,c) for(int a=b;a>=c;--a)
#define LL long long
#define IT iterator
#define PB push_back
#define II(a,b) make_pair(a,b)
#define FIR first
#define SEC second
#define FREO freopen("check.out","w",stdout)
#define rep(a,b) for(int a=0;a<b;++a)
#define SRAND mt19937 rng(chrono::steady_clock::now().time_since_epoch().count())
#define random(a) rng()%a
#define ALL(a) a.begin(),a.end()
#define POB pop_back
#define ff fflush(stdout)
#define fastio ios::sync_with_stdio(false)
#define R(a) cin>>a
#define R2(a,b) cin>>a>>b
using namespace std;
const int INF=0x3f3f3f3f;
typedef pair<int,int> mp;
/*}
*/
const int MAXN=1e6;
int phi[MAXN+2];
LL rest[MAXN+2];
void phi_table(int n) {
  for (int i = 2; i <= n; i++) phi[i] = 0;
  phi[1] = 1;
  for (int i = 2; i <= n; i++)
    if (!phi[i])
      for (int j = i; j <= n; j += i) {
        if (!phi[j]) phi[j] = j;
        phi[j] = phi[j] / i * (i - 1);
      }
}

int main(){
	fastio;
	phi_table(MAXN);
	rb(i,1,MAXN){
		for(int j=i+i;j<=MAXN;j+=i){
			rest[j]+=phi[j/i]*i;
		}
	} 
	rb(i,1,MAXN)
		rest[i]+=rest[i-1];
	while(1){
		int ai;
		R(ai);
		if(!ai) break;
		cout<<rest[ai]<<endl; 
	}
	return 0;
}