[lydsy2005]能量采集

最新推荐文章于 2019-01-02 20:45:35 发布

DraZxlNDdt

最新推荐文章于 2019-01-02 20:45:35 发布

阅读量590

点赞数

分类专栏：数论筛法欧拉函数莫比乌斯反演

本文链接：https://blog.csdn.net/drazxlnddt/article/details/51957856

版权

数论同时被 3 个专栏收录

6 篇文章 0 订阅

订阅专栏

莫比乌斯反演

3 篇文章 0 订阅

订阅专栏

欧拉函数

2 篇文章 0 订阅

订阅专栏

原题地址

前文：很久没有写博客了呢，一直没有时间（在浪）

题意不提。

最后要求的式子

我们单独看这坨cegma

我们可以在o（n）的时间筛出phi（d）。

总复杂度o（n）

代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cctype>
#include<string>
#include<cstring>

const int N = 1e5 + 7;

int n, m, prime[N], p, phi[N];
long long ans;
bool c[N];

int main () {
	scanf ("%d%d", &n, &m);
	if (n > m) n ^= m ^= n ^= m;
	phi[1] = 1;
	for (int i = 2; i <= n; ++i) {
		if (!c[i]) prime[++p] = i, phi[i] = i - 1;
		for (int j = 1; j <= p && (long long) i * prime[j] <= n; ++j) {
			c[i * prime[j]] = true;
			if (i % prime[j] == 0) {
				phi[i * prime[j]] = phi[i] * prime[j];
				break;
			} else {
				phi[i * prime[j]] = phi[i] * (prime[j] - 1);
			}
		}
	}
	for (int i = 1; i <= n; ++i) {
		ans += (long long) (n / i) * (m / i) * phi[i];
	}
	printf ("%lld\n", ans * 2 - (long long) n * m);
	return 0;
}