G. GCD Festival（莫比乌斯、欧拉函数）

G. GCD Festival

∑i=1n∑j=1ngcd⁡(ai,aj)gcd⁡(i,j)∑d=1nd∑i=1nd∑j=1ndgcd⁡(aid,ajd)[gcd⁡(i,j)=1]∑d=1nd∑k=1ndμ(k)∑i=1nkd∑j=1nkdgcd⁡(aikd,ajkd)T=kd∑T=1n∑i=1nT∑j=1nTgcd⁡(aiT,ajT)∑d∣Tdμ(Td)∑T=1nϕ(T)∑i=1nT∑j=1nTgcd⁡(aiT,ajT)\sum_{i = 1} ^{n} \sum_{j = 1} ^{n} \gcd(a_i, a_j) \gcd(i, j)\\ \sum_{d = 1} ^{n} d \sum_{i = 1} ^{\frac{n}{d}} \sum_{j = 1} ^{\frac{n}{d}} \gcd(a_{id}, a_{jd})[\gcd(i, j) = 1]\\ \sum_{d = 1} ^{n} d \sum_{k = 1} ^{\frac{n}{d}} \mu(k) \sum_{i = 1} ^{\frac{n}{kd}} \sum_{j = 1} ^{\frac{n}{kd}} \gcd(a_{i kd}, a_{jkd})\\ T = kd\\ \sum_{T = 1} ^{n} \sum_{i = 1} ^{\frac{n}{T}} \sum_{j = 1} ^{\frac{n}{T}} \gcd(a_{iT}, a_{jT}) \sum_{d \mid T} d \mu(\frac{T}{d})\\ \sum_{T = 1} ^{n} \phi(T) \sum_{i = 1} ^{\frac{n}{T}} \sum_{j = 1} ^{\frac{n}{T}} \gcd(a_{iT}, a_{jT})\\ i=1∑nj=1∑ngcd(ai,aj)gcd(i,j)d=1∑ndi=1∑dnj=1∑dngcd(aid,ajd)[gcd(i,j)=1]d=1∑ndk=1∑dnμ(k)i=1∑kdnj=1∑kdngcd(aikd,ajkd)T=kdT=1∑ni=1∑Tnj=1∑Tngcd(aiT,ajT)d∣T∑dμ(dT)T=1∑nϕ(T)i=1∑Tnj=1∑Tngcd(aiT,ajT)
我们考虑设f(n,T)=∑i=1nT∑j=1nTgcd⁡(aiT,ajT)f(n, T) = \sum\limits_{i = 1} ^{\frac{n}{T}} \sum\limits_{j = 1} ^{\frac{n}{T}} \gcd(a_{iT}, a_{jT})f(n,T)=i=1∑Tnj=1∑Tngcd(aiT,ajT)，g(x)g(x)g(x)为i∈[T,2T,…,nTT]i \in [T, 2T, \dots, \frac{n}{T} T]i∈[T,2T,…,TnT]时xxx的出现次数。
f(n,T)=∑i=1m∑j=1mg(i)g(j)gcd⁡(i,j),(m=105)∑d=1md∑i=1md∑j=1mdg(id)g(jd)[gcd⁡(i,j)=1]∑d=1md∑k=1ndμ(k)(∑i=1mkdg(ikd))2T=kd∑T=1mϕ(T)(∑i=1mTg(iT))2f(n, T) = \sum_{i = 1} ^{m} \sum_{j = 1} ^{m} g(i) g(j) \gcd(i, j), (m = 10 ^ 5)\\ \sum_{d = 1} ^{m} d \sum_{i = 1} ^{\frac{m}{d}} \sum_{j = 1} ^{\frac{m}{d}} g(id) g(jd) [\gcd(i, j) = 1]\\ \sum_{d = 1} ^{m} d \sum_{k = 1} ^{\frac{n}{d}} \mu(k) \left( \sum_{i = 1} ^{\frac{m}{kd}} g(ikd) \right) ^ 2\\ T = kd\\ \sum_{T = 1} ^{m} \phi(T) \left( \sum_{i = 1} ^{\frac{m}{T}} g(iT) \right) ^ 2\\ f(n,T)=i=1∑mj=1∑mg(i)g(j)gcd(i,j),(m=105)d=1∑mdi=1∑dmj=1∑dmg(id)g(jd)[gcd(i,j)=1]d=1∑mdk=1∑dnμ(k)⎝⎛i=1∑kdmg(ikd)⎠⎞2T=kdT=1∑mϕ(T)⎝⎛i=1∑Tmg(iT)⎠⎞2
考虑重新定义g(n)g(n)g(n)表示为是nnn的倍数的数字有多少个，则上式可以直接写成：
∑T=1mϕ(T)g(T)2\sum_{T = 1} ^{m} \phi(T) g(T) ^ 2\\ T=1∑mϕ(T)g(T)2
由此我们可以在O(nlog⁡2n)O(n \log ^ 2n)O(nlog2n)的时间内完成这题。

#include <bits/stdc++.h>using namespace std;const int N = 1e5 + 10, mod = 1e9 + 7;int prime[N], phi[N], a[N], n, cnt;int sum[N], m;bool st[N];vector<int> fac[N];inline int add(int x, int y) {return x + y < mod ? x + y : x + y - mod;
}inline int sub(int x, int y) {return x >= y ? x - y : x - y + mod;
}void init() {phi[1] = 1;for (int i = 2; i < N; i++) {if (!st[i]) {prime[++cnt] = i;phi[i] = i - 1;}for (int j = 1; j <= cnt && 1ll * i * prime[j] < N; j++) {st[i * prime[j]] = 1;if (i % prime[j] == 0) {phi[i * prime[j]] = phi[i] * prime[j];break;}phi[i * prime[j]] = phi[i] * (prime[j] - 1);}}for (int i = 1; i < N; i++) {for (int j = i; j < N; j += i) {fac[j].push_back(i);}}
}int f(int n, int T) {int ans = 0;for (int i = T; i <= n; i += T) {for (auto it : fac[a[i]]) {ans = sub(ans, 1ll * phi[it] * sum[it] % mod * sum[it] % mod);sum[it]++;ans = add(ans, 1ll * phi[it] * sum[it] % mod * sum[it] % mod);}}for (int i = T; i <= n; i += T) {for (auto it : fac[a[i]]) {sum[it]--;}}return ans;
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);init();scanf("%d", &n);for (int i = 1; i <= n; i++) {scanf("%d", &a[i]);}int ans = 0;for (int T = 1; T <= n; T++) {ans = add(ans, 1ll * phi[T] * f(n, T) % mod);}printf("%d\n", ans);return 0;
}

G. GCD Festival（莫比乌斯、欧拉函数）相关推荐

数论 GCD 最大公约数欧拉函数经典题洛谷 CF1295D Same GCDs Codeforces1295D
前言两个月了,我终于更了-- 这两个月忙(chen)于(mi)内(xiang)卷(le),现在终于出新文章啦,(也算兑现了当初的出数论题文章的承诺)~ 不说废话了,今天给大家介绍一道CF/洛谷上的 ...
GCD - Extreme（欧拉函数变形）
题目链接:https://vjudge.net/problem/UVA-11426 题目大意: 给出整数n∈[2,4000000],求解∑gcd(i,j),其中(i,j)满足1≤i<j≤n. 的 ...
poj2154Color polya定理+欧拉函数优化
没想到贱贱的数据居然是错的..搞得我调了一中午+晚上一小时(哦不d飞LJH掉RP毕竟他是BUFF)结果重判就对了五次.. 回归正题,这题傻子都看得出是polya定理(如果你不是傻子就看这里),还没有翻 ...
POJ 2480 (约数+欧拉函数)
题目链接: http://poj.org/problem?id=2480 题目大意:求Σgcd(i,n). 解题思路: 如果i与n互质,gcd(i,n)=1,且总和=欧拉函数phi(n). 如果i与n ...
matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...
数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...
欧拉函数+狄利克雷卷积+莫比乌斯函数+莫比乌斯反演+整除分块+杜教筛
Powered by:NEFU AB-IN 文章目录欧拉函数狄利克雷卷积莫比乌斯函数莫比乌斯反演 P3455 [POI2007]ZAP-Queries 整除分块 P2522 [HAOI2011 ...
专题·莫比乌斯函数与欧拉函数【including 整除分块，积性函数，狄利克雷卷积，欧拉函数，莫比乌斯函数，莫比乌斯反演
初见安~又是好久没写博客了--加上CSP才炸了一波. 目录一.整除分块题解二.积性函数三.狄利克雷卷积四.欧拉函数五.莫比乌斯函数(mu) 六.莫比乌斯反演一.整除分块看个例题:洛谷P ...
数学/数论专题：莫比乌斯函数与欧拉函数
数学/数论专题:莫比乌斯函数与欧拉函数(进阶) 0. 前言 1. 前置知识 2. 正文 3. 总结 4. 参考资料 0. 前言本篇文章会从狄利克雷卷积的角度,讨论莫比乌斯函数与欧拉函数的相关性质. ...
Educational Codeforces Round 81 (Rated for Div. 2) D. Same GCDs 欧拉函数\莫比乌斯
传送门文章目录题意: 思路: 题意: 给定a,ma,ma,m,求满足gcd(a,m)=gcd(a+x,m)gcd(a,m)=gcd(a+x,m)gcd(a,m)=gcd(a+x,m)的xxx的个数 ...

G. GCD Festival（莫比乌斯、欧拉函数）

G. GCD Festival

G. GCD Festival（莫比乌斯、欧拉函数）相关推荐

最新文章

热门文章