【CF 235E】Number Challenge

最新推荐文章于 2024-05-02 13:48:13 发布

原创最新推荐文章于 2024-05-02 13:48:13 发布 · 824 阅读

1 ·

CC 4.0 BY-SA版权

数学同时被 2 个专栏收录

20 篇文章

订阅专栏

codeforces

6 篇文章

订阅专栏

本文详细介绍了一个基于数论的算法实现，该算法通过预处理获得特定数值范围内所有整数的μ函数值，并结合输入参数计算最终结果。算法涉及了数学中的筛法、倍数关系和组合计数等内容。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

求 $\text{[math]}$ 1 ≤ a, b, c ≤ 2000

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#define Rep(i, x, y) for (int i = x; i <= y; i ++)
#define Dwn(i, x, y) for (int i = x; i >= y; i --)
#define RepE(i, x) for(int i = pos[x]; i; i = g[i].nex)
using namespace std;
typedef long long LL;
const int N = 2005, M = N * N, mod = 1073741824;
int A, B, C, mu[M], pri[M], pz, num[M], ab;
LL f[M], ans; bool c[M];
void Pre() {
	mu[1] = 1;
	Rep(i, 2, C) {
		if (!c[i]) pri[++ pz] = i, mu[i] = -1;
		Rep(j, 1, pz) {
			int k = pri[j] * i;
			if (k > ab) break ;
			c[k] = 1, mu[k] = -mu[i];
			if (i % pri[j] == 0) { mu[k] = 0; break ; }
		}
	}
	Rep(e, 1, C) {
		LL k = 0;
		Rep(d, 1, C) k += C / (e * d);
		for (int g = e; g <= ab; g += e) f[g] += mu[e] * k;
	}
	Rep(g, 1, ab) f[g] %= mod;
}
int main()
{
	scanf ("%d%d%d", &A, &B, &C);
	Rep(i, 1, A) Rep(j, 1, B) num[i * j] ++;
	ab = A * B;
	Pre(); // puts("fin");
	Rep(g, 1, ab) if (num[g]) {
		LL k1 = 0;
		Rep(j, 1, ab / g) (k1 += num[j * g]) %= mod;
		(ans += f[g] * k1) %= mod;
	}
	cout << ((ans + mod) % mod) << endl;

	return 0;
}