题解 [ABC020D] LCM Rush

该博客是 [ABC020D] LCM Rush 的题解。题目要求给定 n、k,求 i 从 1 到 n 的 lcm(i,k) 之和。通过公式转换,将问题转化为给定 m、p 求 i 从 1 且与 p 互质到 m 的和,可使用容斥原理求解,还给出了具体的容斥计算方式。

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

题解 [ABC020D] LCM Rush

题意:给 n,kn,kn,k,求 ∑i=1nlcm⁡(i,k)\sum\limits_{i=1}^n \operatorname{lcm}(i,k)i=1nlcm(i,k)1≤n,k≤1091\leq n,k\leq 10^91n,k109

∑i=1nlcm⁡(i,k)=∑i=1nikgcd⁡(i,k)=k∑i=1nigcd⁡(i,k) \begin{aligned} \sum_{i=1}^n \operatorname{lcm}(i,k)&=\sum_{i=1}^n \dfrac{ik}{\gcd(i,k)}\\ &=k\sum_{i=1}^n \dfrac{i}{\gcd(i,k)} \end{aligned} i=1nlcm(i,k)=i=1ngcd(i,k)ik=ki=1ngcd(i,k)i

枚举 g=gcd⁡(i,k)g=\gcd(i,k)g=gcd(i,k)

∑i=1nlcm⁡(i,k)=k∑g∣k∑i=1,(i,kg)=1⌊ng⌋i \begin{aligned} \sum_{i=1}^n \operatorname{lcm}(i,k)&=k\sum_{g|k}\sum_{i=1,(i,\frac kg)=1}^{\lfloor\frac ng\rfloor}i \end{aligned} i=1nlcm(i,k)=kgki=1,(i,gk)=1gni

现在问题变成了,给定 m,pm, pm,p,求 ∑i=1,(i,p)=1mi\sum\limits_{i=1,(i,p)=1}^mii=1,(i,p)=1mi。可以容斥求解。

具体地,枚举所有 q∣pq|pqpqqq 每个质因子次数不大于 111,则

∑i=1,(i,p)=1mi=∑qq×(∑i=1⌊mq⌋i)×valq \sum\limits_{i=1,(i,p)=1}^mi=\sum_q q\times \left(\sum_{i=1}^{\lfloor\frac mq\rfloor}i\right)\times val_q i=1,(i,p)=1mi=qq×i=1qmi×valq

其中当 qqq 的质因子个数为偶数时 valq=1val_q=1valq=1,否则 valq=−1val_q=-1valq=1。(111 不算质因子)

/*
    name: LCM Rush
    id:   AT_abc020_d
    date: 2023/01/25
*/

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const ll P = 1e9 + 7;
int n, k, top;
ll ans = 0;

int calc(int p){
	int sum = 1;
	for(int i = 2; i * i <= p; ++ i){
		if(p % i == 0){
			++ sum;
			p /= i;
		}
		if(p % i == 0){
			return 0;	
		}
	}
	if(p != 1){
		++ sum;
	}
	return (sum&1) ? 1 : -1;
}

#define adall(k) ((ll)(k) * ((k)+1) / 2)

ll solve(int g){
	int p = k / g, m = n / g;
	ll cur = 0;
	for(int q = 1; q * q <= p; ++ q){
		if(p % q == 0){
			cur += calc(q) * q * adall((int)m/q);
			int r = p/q;
			if(r != q){
				cur += calc(r) * r * adall((int)m/r);
			}
		}
	}
	return cur % P;
}

int main(){
	scanf("%d%d", &n, &k);
	for(int g = 1; g * g <= k; ++ g){
		if(k % g != 0){
			continue;
		}
		ans += solve(g);
		if(g * g != k){
			ans += solve(k / g);
		}
		ans %= P;
	}
	printf("%lld\n", ans * k % P);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值