CodeForces 1717E【线性筛】

最新推荐文章于 2024-11-17 15:19:52 发布

willem248

最新推荐文章于 2024-11-17 15:19:52 发布

阅读量203

点赞数

分类专栏： CodeForces 文章标签： c++ 算法图论

转载请注明出处，谢谢。

本文链接：https://blog.csdn.net/lxylluvio/article/details/126677171

版权

CodeForces 专栏收录该内容

17 篇文章 1 订阅

订阅专栏

本文介绍了如何使用埃拉托斯特尼筛法求解给定整数n下，满足a+b+c=n且lcm(gcd(a,b),c)的三元组(a,b,c)的和。通过枚举gcd(a,b)和a+b的组合，并利用欧拉函数计算可能的组合数量，实现了高效的算法实现。

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

题意：给定 $n$ ，求所有三元组 $(a, b, c)$ 满足 $a + b + c = n$ 的 $\operatorname{lcm}(\gcd(a,b),c)$ 的和。

思想：

第一层循环枚举 $\gcd(a,b)$ 的值。

然后枚举 $a + b$ 的值。

容易发现， $a + b$ 的值一定是 $k\times \gcd (a,b)$ ， $k\in N^*$ 。

也就是说， $(a, b)$ 的取值方法一共有 $\varphi(\lfloor\frac{n-c}{\gcd(a,b)}\rfloor)$ 种，所以考虑进行类似于埃式筛法的算法，每一次更新的时候将答案加上 $\operatorname{lcm}(\gcd(a,b),c)\times \varphi(\lfloor\frac{n-c} {\gcd(a,b)}\rfloor)$ 即可。

// 这回只花了45min就打完了。
// 真好。记得多手造几组。最好有暴力对拍。

#include <bits/stdc++.h>

#define int long long

using namespace std;

const int N = 2e5 + 10, mod = 1e9 + 7;
int phi[N]; // debug 这个地方要开 long long
bool not_prime[N];
int cnt, primes[N], n, ans = 0;

void sieve()
{
  phi[1] = 1;
  for (int i = 2; i < N; i ++)
  {
    if (!not_prime[i])
    {
      primes[++ cnt] = i;
      phi[i] = i - 1;
    }
    for (int j = 1; i * primes[j] < N; j ++)
    {
      not_prime[i * primes[j]] = true; // debug i*primes[j] --> j*primes[j]
      if (i % primes[j])
        phi[i * primes[j]] = phi[i] * (primes[j] - 1);
      else
      {
        phi[i * primes[j]] = phi[i] * primes[j];
        break;
      }
    }
  }
}

signed main()
{
  cin >> n;
  sieve();
  for (int i = 1; i < n; i ++)
    for (int j = 2 * i; j < n; j += i) // debug =n 的话没有法子选了
      (ans += phi[j / i] * lcm(i, n - j) % mod) %= mod;
  cout << ans << '\n';
  return 0;
}

// 足下千里，移步换景，寰宇纷呈万花筒