HDU 6428 Problem C. Calculate（积性函数）

Problem C. Calculate

ϕ=ϕ∗ϵ=ϕ∗μ∗Iϕ(n)=∑d∣n(ϕ∗μ)(d)设g(n)=∑d∣n(ϕ∗μ)(d)∑i=1A∑j=1B∑k=1Cϕ(gcd(i,j2,k3))∑i=1A∑j=1B∑k=1C∑d∣i,d∣j2,d∣k3(ϕ∗μ)(d)∑d=1A(ϕ∗μ)(d)∑i=1A∑j=1B∑k=1C[d∣i,d∣j2,d∣k3]\phi = \phi * \epsilon = \phi * \mu * I\\ \phi(n) = \sum_{d \mid n} (\phi*\mu)(d)\\ 设g(n) = \sum_{d \mid n}(\phi * \mu)(d)\\ \sum_{i = 1} ^{A} \sum_{j = 1} ^{B} \sum_{k = 1} ^{C} \phi(gcd(i, j ^ 2, k ^ 3))\\ \sum_{i = 1} ^{A} \sum_{j = 1} ^{B} \sum_{k = 1} ^{C} \sum_{d \mid i, d \mid j ^ 2, d \mid k ^ 3} (\phi *\mu)(d)\\ \sum_{d = 1} ^{A} (\phi *\mu)(d) \sum_{i = 1} ^{A} \sum_{j = 1} ^{B} \sum_{k = 1} ^{C}[d \mid i, d \mid j ^ 2, d \mid k ^ 3]\\ ϕ=ϕ∗ϵ=ϕ∗μ∗Iϕ(n)=d∣n∑(ϕ∗μ)(d)设g(n)=d∣n∑(ϕ∗μ)(d)i=1∑Aj=1∑Bk=1∑Cϕ(gcd(i,j2,k3))i=1∑Aj=1∑Bk=1∑Cd∣i,d∣j2,d∣k3∑(ϕ∗μ)(d)d=1∑A(ϕ∗μ)(d)i=1∑Aj=1∑Bk=1∑C[d∣i,d∣j2,d∣k3]
设x=∏pikix = \prod p_i ^{k_i}x=∏piki，唯一分解，如果x∣ytx \mid y ^ tx∣yt，则有∏pi⌈kit⌉∣y\prod p_i ^{\lceil \frac{k_i}{t} \rceil} \mid y∏pi⌈tki⌉∣y，设ft(n)=∏pjkjtf_t(n) = \prod p_j ^ \frac{k_j}{t}ft(n)=∏pjtkj。

上式可得：
∑d=1A(ϕ∗μ)(d)⌊df1(d)⌋⌊df2(d)⌋⌊df3(d)⌋\sum_{d = 1} ^{A}(\phi * \mu)(d) \lfloor \frac{d}{f_1(d)} \rfloor \lfloor \frac{d}{f_2(d)} \rfloor \lfloor \frac{d}{f_3(d)} \rfloor d=1∑A(ϕ∗μ)(d)⌊f1(d)d⌋⌊f2(d)d⌋⌊f3(d)d⌋

对于h=(ϕ∗μ)(d),f1,f2,f3h = (\phi * \mu)(d), f_1, f_2, f_3h=(ϕ∗μ)(d),f1,f2,f3都具有积性，所以可以线性筛预处理，最后O(A)O(A)O(A)统计答案。

#include <bits/stdc++.h>using namespace std;typedef unsigned uint;const int N = 1e7 + 10;uint prime[N], h[N], f1[N], f2[N], f3[N], num[N], primes[N], cnt;bool st[N];void init() {h[1] = f1[1] = f2[1] = f3[1] = primes[1] = 1;for (int i = 2; i < N; i++) {if (!st[i]) {prime[++cnt] = i;h[i] = i - 2;f1[i] = i;f2[i] = i;f3[i] = i;primes[i] = i;num[i] = 1;}for (int j = 1; j <= cnt && 1ll * i * prime[j] < N; j++) {st[i * prime[j]] = 1;if (i % prime[j] == 0) {primes[i * prime[j]] = primes[i] * prime[j];num[i * prime[j]] = num[i] + 1;f1[i * prime[j]] = f1[i] * prime[j];f2[i * prime[j]] = f2[i / prime[j]] * prime[j];if (i / prime[j] % prime[j] == 0) {f3[i * prime[j]] = f3[i / prime[j] / prime[j]] * prime[j];}else {f3[i * prime[j]] = f3[i];}if(i * prime[j] == primes[i * prime[j]]) {h[i * prime[j]] = i * prime[j] - 2 * i + i / prime[j];}else {h[i * prime[j]] = h[i / primes[i]] * h[primes[i * prime[j]]];}break;}h[i * prime[j]] = h[i] * h[prime[j]];f1[i * prime[j]] = f1[i] * f1[prime[j]];f2[i * prime[j]] = f2[i] * f2[prime[j]];f3[i * prime[j]] = f3[i] * f3[prime[j]];primes[i * prime[j]] = prime[j];num[i * prime[j]] = 1;}}
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);init();int T;scanf("%d", &T);while (T--) {uint A, B, C, ans = 0;scanf("%u %u %u", &A, &B, &C);for (int d = 1; d <= (int)A; d++) {ans += h[d] * (A / f1[d]) * (B / f2[d]) * (C / f3[d]);}printf("%u\n", ans & ((1 << 30) - 1));}return 0;
}

HDU 6428 Problem C. Calculate（积性函数）相关推荐

CCPC-2017 杭州站B丨HDU - 6265丨数论丨积性函数丨欧拉函数丨狄利克雷卷积丨思维变换
[参考博客]@WJHKDGHP ccpc2017杭州站 B [参考博客]@灬从此以后灬 2017 CCPC 杭州 HDU6265B 积性函数特别感谢以上两位博主,让我看懂了许多细节. 交题网址(HD ...
POJ-2480 Longge's problem 积性函数
题目链接:http://poj.org/problem?id=2480 题意:多次求sigma(gcd(i,n), 1<=i<=n<2^32) 这题不能直接搜了,需要考虑函数的性质. ...
POJ 2480 Longge#39;s problem 积性函数
题目来源:POJ 2480 Longge's problem 题意:求i从1到n的gcd(n, i)的和思路:首先假设m, n 互质 gcd(i, n*m) = gcd(i, n)*gcd(i, m ...
HDU 6134 2017 多校训练：Battlestation Operational（莫比乌斯反演+积性函数）
实在太长了直接放题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=6134 这题就是求考虑当Gcd(i, j)==1时,除了j为1的情况,其它时候i/j一 ...
《算法竞赛中的初等数论》（三）正文 0x30 积性函数（ACM / OI / MO）（十五万字符数论书）
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划写在最前面:本文部分内容来自网上各大博客或是各类图书,由我个人整理,增加些许见解,仅做学习交流使用,无 ...
【数学专题】莫比乌斯反演与积性函数
整理的算法模板合集: ACM模板点我看算法全家桶系列!!! 实际上是一个全新的模板整合计划目录莫比乌斯反演 AcWing 2702. problem b AcWing 1358. 约数个数和(莫 ...
浅谈积性函数求前缀和
转载至https://blog.csdn.net/skywalkert/article/details/50500009 前置技能积性函数的定义若f(n)f(n)的定义域为正整数域,值域为复数,即 ...
2018ACM-ICPC南京赛区网络赛: J. Sum（积性函数前缀和）
J. Sum A square-free integer is an integer which is indivisible by any square number except 11. For ...
积性函数性质-POJ2480
https://vj.xtuacm.cf/contest/view.action?cid=146#problem/D 参考:http://www.cnblogs.com/xiaowuga/p/7161 ...

HDU 6428 Problem C. Calculate（积性函数）

Problem C. Calculate

HDU 6428 Problem C. Calculate（积性函数）相关推荐

最新文章

热门文章