E - Madoka and The Best University（数论/欧拉函数）

对方正在debug

已于 2022-09-04 20:22:20 修改

阅读量295

点赞数

分类专栏： Codeforces 文章标签：算法

于 2022-09-04 09:54:33 首次发布

本文链接：https://blog.csdn.net/weixin_43918473/article/details/126685672

版权

Codeforces 专栏收录该内容

228 篇文章 7 订阅

订阅专栏

这篇博客介绍了如何利用欧拉函数解决一道数学问题，即给定正数n，计算∑lcm(c,gcd(a,b))的值，其中a+b+c=n。通过暴力枚举和欧拉函数的性质，可以计算出满足条件的组合数。博客还提供了C++代码实现，并邀请读者一起刷题讨论。

摘要由CSDN通过智能技术生成

题目
 参考

题意

给定正数n，求 $\sum lcm(c,gcd(a,b)),a+b+c=n,1<=a,b,c$

思路

暴力枚举c，对于固定的c，a+b=n-c
因为gcd(a,b)=gcd(a,a+b)=gcd(a,n-c)
枚举n-c的所有因子d，对于固定的d
有d=gcd(a,b)=gcd(a,n-c),d|n，要求满足能被d整除且与n-c互素的a的个数，
实际上，就是欧拉值 $\phi(n-c/d)$ 。

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int maxn = 100010;
const int mod = 1e9 + 7;

int n;

int phi[maxn]; // 欧拉函数 <= x中与x互素的个数 
int vis[maxn];
int pri[maxn];
int tot = 0;
// 素数筛 求欧拉函数 
/*
素数p的欧拉值 
phi(p) = p - 1;

合数x的欧拉值 
x = p1^k1 * p2^k2 * ... pr^kr. p1,p2,...,pr为素数且p1<p2<...<pr
phi(x) = p1^(k1-1)*(p1-1) *  p2^(k2-1)*(p2-1) * ... * pr^(kr-1)*(pr-1).
*/
void prime_init(int n) {
    phi[1] = 0;
    for (int i = 2; i <= n; ++i) {
        if (!vis[i]) {
        	pri[++tot] = i;
			phi[i] = i - 1;// 素数x的欧拉函数值为x-1 
		}
        for (int j = 1; j <= tot && i * pri[j] <= n; ++j) {
            vis[i * pri[j]] = 1;
            if (i % pri[j] == 0) {
            	// 根据欧拉函数的定义 
				// phi(x*p) = phi(x) * p, x % p == 0,p为素数 
                phi[i * pri[j]] = phi[i] * pri[j];
                break;
            }
            // 根据欧拉函数的定义 
			// phi(x*p) = phi(x) * (p-1), gcd(x,p)==1,p为素数 
            phi[i * pri[j]] = phi[i] * (pri[j] - 1);
        }
    }
}
int lcm(int x, int y) {
	return 1LL * x * y / __gcd(x, y) % mod;
}
void solve() {
	scanf("%d", &n);
	ll res = 0;
	for (int c = 1; c <= n - 2; ++c) {
		int tmp = n - c;
		for (int d = 1; d <= tmp / d; ++d) {
			if (tmp % d != 0) {
				continue;
			}
			res += 1LL * lcm(c, d) * phi[tmp / d] % mod, res %= mod;
			if (d != tmp / d)
				res += 1LL * lcm(c, tmp / d) * phi[d] % mod, res %= mod;
		}
	}
	printf("%lld\n", res);
}

int main() {
	int t;
	prime_init(maxn - 5);
//	scanf("%d", &t);
	t = 1; 
	while (t--) {
		solve();
	}
}