8.12.24 ACM-ICPC数学数论 Powerful Number 筛

本文链接：https://blog.csdn.net/tang7mj/article/details/139204910

8.12.23 ACM-ICPC数学数论杜教筛

杜教筛是一种用于求解数论函数前缀和的高效算法，尤其适用于积性函数。通过构造合适的数论函数 𝑔g 和递推关系，可以在低于线性时间复杂度内计算出目标函数的前缀和。

算法思想

杜教筛的核心在于构造一个数论函数 𝑔g 满足特定性质，并利用狄利克雷卷积性质来实现递推求解。对于任意数论函数 𝑔g，必满足：

其中 𝑓∗𝑔f∗g 为数论函数 𝑓f 和 𝑔g 的狄利克雷卷积。

通过将双重求和变换为递推关系，可以得到：

进一步化简得到递推式：

𝑆(𝑛)=1𝑔(1)(∑𝑖=1𝑛(𝑓∗𝑔)(𝑖)−∑𝑖=2𝑛𝑔(𝑖)𝑆(⌊𝑛𝑖⌋))S(n)=g(1)1(i=1∑n(f∗g)(i)−i=2∑ng(i)S(⌊in⌋))

复杂度分析

例题分析

问题一：求莫比乌斯函数和欧拉函数前缀和

代码实现：

#include <algorithm>
#include <cstdio>
#include <cstring>
#include <map>
using namespace std;
const int maxn = 2000010;
long long T, n, pri[maxn], cur, mu[maxn], sum_mu[maxn];
bool vis[maxn];
map<long long, long long> mp_mu;

long long S_mu(long long x) {  // 求mu的前缀和
  if (x < maxn) return sum_mu[x];
  if (mp_mu[x]) return mp_mu[x];  // 如果map中已有该大小的mu值，则可直接返回
  long long ret = 1;
  for (long long i = 2, j; i <= x; i = j + 1) {
    j = x / (x / i);
    ret -= S_mu(x / i) * (j - i + 1);
  }
  return mp_mu[x] = ret;  // 路径压缩，方便下次计算
}

long long S_phi(long long x) {  // 求phi的前缀和
  long long ret = 0;
  long long j;
  for (long long i = 1; i <= x; i = j + 1) {
    j = x / (x / i);
    ret += (S_mu(j) - S_mu(i - 1)) * (x / i) * (x / i);
  }
  return (ret - 1) / 2 + 1;
}

int main() {
  scanf("%lld", &T);
  mu[1] = 1;
  for (int i = 2; i < maxn; i++) {  // 线性筛预处理mu数组
    if (!vis[i]) {
      pri[++cur] = i;
      mu[i] = -1;
    }
    for (int j = 1; j <= cur && i * pri[j] < maxn; j++) {
      vis[i * pri[j]] = true;
      if (i % pri[j])
        mu[i * pri[j]] = -mu[i];
      else {
        mu[i * pri[j]] = 0;
        break;
      }
    }
  }
  for (int i = 1; i < maxn; i++)
    sum_mu[i] = sum_mu[i - 1] + mu[i];  // 求mu数组前缀和
  while (T--) {
    scanf("%lld", &n);
    printf("%lld %lld\n", S_phi(n), S_mu(n));
  }
  return 0;
}