问题都是求 ∑ i = 1 n ∑ j = i + 1 n gcd ( i , j ) \sum\limits_{i=1}^n\sum\limits_{j=i+1}^n \gcd(i,j) i=1∑nj=i+1∑ngcd(i,j)
总结一下套路
首先第一题直接 O ( n 2 ) O(n^2) O(n2) 暴力就能解决
int n, ans ;
signed main(){
while (scanf("%d", &n) != EOF && n) {
ans = 0 ;
rep(i, 1, n) rep(j, i + 1, n) ans += __gcd(i, j) ;
printf("%d\n", ans) ;
}
return 0 ;
}
对于第二题, n < = 200000 n <= 200000 n<=200000
我们枚举 gcd \gcd gcd 为 k k k,现在的问题就是有多少对, 我们记做 f ( k ) f(k) f(k), 那么答案就是 ∑ i = 1 n ( f ( i ) − 1 ) ∗ i \sum\limits_{i=1}^n(f(i)-1)*i i=1∑n(f(i)−1)∗i,注意要 − 1 -1 −1 因为不存在 gcd ( i , i ) = i \gcd(i,i)=i gcd(i,i)=i 的情况
f ( k ) f(k) f(k) 怎么求?
若 a,b 互质, 则 gcd ( a k , b k ) = k \gcd(ak,bk)=k gcd(ak,bk)=k
枚举 a , b a,b a,b 中较大的一个,记做 i i i, 那么另一个数有 φ ( i ) \varphi(i) φ(i)种可能,所以 f ( k ) = ∑ i = 1 n / k φ ( i ) f(k)=\sum\limits_{i=1}^{n/k}\varphi(i) f(k)=i=1∑n/kφ(i),用前缀和处理
复杂度 O ( n ) O(n) O(n)
int n, cnt ;
ll phi[N], f[N], s[N], prime[N / 10] ;
ll ans ;
void sieve(int n) {
phi[1] = 1 ;
rep(i, 2, n) {
if (!f[i]) {
prime[++cnt] = i ;
phi[i] = i - 1 ;
}
for (int j = 1; j <= cnt && prime[j] * i <= n; j++) {
f[i * prime[j]] = 1 ;
if (i % prime[j]) phi[i * prime[j]] = phi[i] * (prime[j] - 1) ;
else {
phi[i * prime[j]] = phi[i] * prime[j] ;
break ;
}
}
}
}
void init(int n) {
clr(f) ;
rep(i, 1, n)
for (int j = 2 * i; j <= n; j += i)
f[j] += i * phi[j / i] ;
rep(i, 1, n) s[i] = s[i - 1] + f[i] ;
}
signed main() {
sieve(N - 10) ;
init(N - 10) ;
while (scanf("%d", &n) != EOF && n) printf("%lld\n", s[n]) ;
}
对于第三题, n < = 4000000 n<=4000000 n<=4000000
其实本质就是 Z A P ZAP ZAP 那题一样的,直接莫比乌斯反演+整除分块就行了
时间复杂度 O ( n ) O(\sqrt n) O(n)