1190 最小公倍数之和 V2

1190 最小公倍数之和 V2

∑i=ablcm(i,b)∑i=abi×bgcd⁡(i,b)b∑d∣b∑i=⌈ad⌉bdi[gcd(i,bd)=1]b∑d∣b∑k∣bdμ(k)k∑i=⌈⌈ad⌉k⌉abkib∑T∣n∑i=⌈aT⌉bTi∑k∣Tμ(k)kb∑T∣n(bT+⌈aT⌉)×(bT−⌈aT⌉+1)2∑k∣Tμ(k)k设f(n)=∑d∣nμ(d)d,f(1)=1,f(p)=1−p,f(pk)=1−p,且为积性函数\sum_{i = a} ^{b} lcm(i, b)\\ \sum_{i = a} ^{b} \frac{i \times b}{\gcd(i, b)}\\ b \sum_{d \mid b} \sum_{i = \lceil \frac{a}{d} \rceil} ^{\frac{b}{d}} i [gcd(i, \frac{b}{d}) = 1]\\ b \sum_{d \mid b} \sum_{k \mid \frac{b}{d}} \mu(k) k \sum_{i = \lceil \frac{ \lceil \frac{a}{d} \rceil}{k} \rceil} ^{\frac{a}{bk}}i\\ b \sum_{T \mid n} \sum_{i = \lceil \frac{a}{T} \rceil} ^{\frac{b}{T}} i \sum_{k \mid T} \mu(k) k\\ b \sum_{T \mid n} \frac{ (\frac{b}{T} + \lceil \frac{a}{T} \rceil) \times ( \frac{b}{T} - \lceil \frac{a}{T} \rceil + 1)}{2} \sum_{k \mid T} \mu(k) k\\ 设f(n) = \sum_{d \mid n} \mu(d) d,f(1) = 1, f(p) = 1 - p, f(p ^ k) = 1- p,且为积性函数\\ i=ablcm(i,b)i=abgcd(i,b)i×bbdbi=dadbi[gcd(i,db)=1]bdbkdbμ(k)ki=kdabkaibTni=TaTbikTμ(k)kbTn2(Tb+Ta)×(TbTa+1)kTμ(k)kf(n)=dnμ(d)d,f(1)=1,f(p)=1p,f(pk)=1p,

所以最后我们只要2prime_facofb2 ^{prime\_fac\ of\ b}2prime_fac of b的复杂度去枚举ddd即可。

/*Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>#define mp make_pair
#define pb push_back
#define endl '\n'
#define mid (l + r >> 1)
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
#define ls rt << 1
#define rs rt << 1 | 1using namespace std;typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;const double pi = acos(-1.0);
const double eps = 1e-7;
const int inf = 0x3f3f3f3f;inline ll read() {ll f = 1, x = 0;char c = getchar();while(c < '0' || c > '9') {if(c == '-')    f = -1;c = getchar();}while(c >= '0' && c <= '9') {x = (x << 1) + (x << 3) + (c ^ 48);c = getchar();}return f * x;
}const int N = 1e5 + 10, mod = 1e9 + 7;int prime[N], cnt;bool st[N];void init() {for(int i = 2; i < N; i++) {if(!st[i]) prime[cnt++] = i;for(int j = 0; j < cnt && i * prime[j] < N; j++) {st[i * prime[j]] = 1;if(i % prime[j] == 0) break;}}
}ll a, b, ans;int fac[50], num[50], tot;void dfs(int step, ll T, ll g){if(step == tot + 1){ans = (ans + ((1 + b / T ) * (b / T) / 2 - (1 + a / T) * (a / T) / 2) % mod * g % mod) % mod;return ;}ll sum = 1;dfs(step + 1, T * sum, g);for(int i = 1; i <= num[step]; i++){sum = sum * fac[step];dfs(step + 1, T * sum, g * (1 - fac[step]));}
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);init();int T = read();while(T--) {a = read() - 1, b = read(); ll n = b;tot = 0;for(int i = 0; prime[i] * prime[i] <= n; i++) {if(n % prime[i] == 0) {tot++;fac[tot] = prime[i], num[tot] = 0;while(n % prime[i] == 0) {n /= prime[i];num[tot]++;}}}if(n != 1) {tot++;fac[tot] = n, num[tot] = 1;}ans = 0;dfs(1, 1, 1);printf("%lld\n",(ans * b % mod + mod) % mod);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/313833.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

有关 VS Code 的五大谣言,背后的真相到底是如何的?

2015 年 4 月 29 日&#xff0c;在微软 Build 2015 大会上&#xff0c;微软发布了 Visual Studio Code 第一个预览版本。随着 Visual Studio Code 有了越来越多的使用者&#xff0c;随之而来的各类谣言也层出不穷。让我们就来看看有哪些与 VS Code 相关的谣言&#xff0c;背后的…

Mysql数据库锁机制

一&#xff1a;概念介绍 MySQL数据库锁管理机制&#xff1a; SQL层实现的锁机制    Meta-data元数据锁&#xff1a;在table cache缓存里实现的&#xff0c;为DDL&#xff08;Data Definition Language&#xff09;提供隔离操作。一种特别的meta-data元数据类型&#xff0c;…

HDU 6340 Problem I. Delightful Formulas(伯努利数 + 积性函数反演)

Problem I. Delightful Formulas 大概就是照着题解抄了一遍吧&#xff0c;这道题太神仙了…… aiik,si∑j1iajcalc∑i1nsi[gcd⁡(i,n)1]∑d∣nμ(d)∑i1ndsida_i i ^ k, s_i \sum_{j 1} ^{i} a_j\\ calc\ \sum_{i 1} ^{n} s_i[\gcd(i, n) 1]\\ \sum_{d \mid n} \mu(d) \s…

干货|亲测有效的N倍学习效果笔记法

这里是Z哥的个人公众号每周五11&#xff1a;45 按时送达当然了&#xff0c;也会时不时加个餐&#xff5e;我的第「108」篇原创敬上大家好&#xff0c;我是Z哥。先祝大家中秋快乐。我猜你现在心情不错&#xff0c;毕竟小长假的第一天才开始&#xff0c;后面还有60个小时的假期&a…

Java偏向锁、轻量级锁、重量级锁

先Mark&#xff0c;后补充 参照&#xff1a; https://www.infoq.cn/article/java-se-16-synchronized http://www.cnblogs.com/paddix/p/5405678.html http://www.cnblogs.com/lzh-blogs/p/7477157.html

.NET Core 3.0 可卸载程序集原理简析

文章转载授权级别&#xff1a;A 预计阅读时间&#xff1a;8分钟 损失发量&#xff1a;不好统计因为最近在群里被问到如何理解 .NET Core 3.0 可卸载程序集&#xff0c;所以就写了这篇简单的分析。因为时间实在很少&#xff0c;这篇文章只简单的罗列了相关的代码&…

P4331 [BalticOI 2004]Sequence 数字序列(左偏树)

P4331 [BalticOI 2004]Sequence 数字序列 给定一个序列整数a1,a2,a3,…,an−1,ana_1, a_2, a_3, \dots, a_{n - 1}, a_na1​,a2​,a3​,…,an−1​,an​&#xff0c;要找一个整数序列bbb&#xff0c;满足b1<b2<b3<⋯<bn−1<bnb_1 < b_2 < b_3< \dots&…

.NetCore技术研究-ConfigurationManager在单元测试下的坑

最近在将原有代码迁移.NET Core, 代码的迁移基本很快&#xff0c;当然也遇到了不少坑&#xff0c;重构了不少&#xff0c;后续逐步总结分享给大家。今天总结分享一下ConfigurationManager遇到的一个问题。先说一下场景&#xff1a;迁移.NET Core后&#xff0c;已有的配置文件&a…

E 速度即转发(牛客挑战赛48)(树套树)

速度即转发 给定一个长度为nnn的数组aaa&#xff0c;里面元素为a1,a2,a3,…,an−1,ana_1, a_2, a_3, \dots, a_{n - 1}, a_na1​,a2​,a3​,…,an−1​,an​。 有两种操作&#xff1a; 修改apka_p kap​k。给定l,r,kl, r, kl,r,k&#xff0c;设S(x)∑ilrmax(ai−x,0)S(x) …

volatile实现原理

先Mark&#xff0c;后续完成 https://segmentfault.com/a/1190000017255405 http://ifeve.com/volatile/ http://cmsblogs.com/?hmsrtoutiao.io&p2092&utm_mediumtoutiao.io&utm_sourcetoutiao.io https://my.oschina.net/u/2288283/blog/656572 https://blog.…

分析一次double强转float的翻车原因

人逢喜事精神爽,总算熬到下班撩~~正准备和同事打个招呼回家,被同事拖住问了.?‍♂️: 你们组做的那块代码,把double类型数据成float有问题啊?.?‍♀️: 嗯?不对是正常啊,float精度是没有double高,但float能保存到小数点后好多位,对我们来说完全够用了!?‍♂️: 不是啊,这不…

L 苍天阻我寻你,此情坚贞如一(西南科技大学2021届新生赛)(线段树)

苍天阻我寻你&#xff0c;此情坚贞如一 给定两个长度为nnn的数组a,ba, ba,b&#xff0c;满足−103≤ai,bi≤103-10 ^ 3 \leq a_i, b_i \leq 10 ^ 3−103≤ai​,bi​≤103&#xff0c;每个数字xxx表示向前走xxx步&#xff0c;如果是负数则后退嘛&#xff0c;设小AAA执行aaa数组…

Java并发常用方法 sleep 和 wait

一&#xff1a;sleep 和 wait sleep()方法&#xff1a; 功能&#xff1a;让当前线程休眠指定时间&#xff0c;休眠时间的准确性依赖于系统时钟和CPU调度机制是Thread类的静态方法可在任何地方调用&#xff0c;需要处理InterruptedException当前线程调用sleep()方法&#xff0…

.NET Core 3.0之深入源码理解Host(二)

写在前面停了近一个月的技术博客&#xff0c;随着正式脱离996的魔窟&#xff0c;接下来也正式恢复了。本文从源码角度进一步讨论.NET Core 3.0 中关于Host扩展的一些技术点&#xff0c;主要内容是关于创建Long Run Program的创建与守护。关于Host&#xff0c;我们最容易想到的就…

Finding Hotels(牛客国庆集训派对Day7 )(2016ICPC青岛K)(K-D Tree)

Finding Hotels 给定二维平面上nnn个点&#xff0c;每个点描述为x,y,cx, y, cx,y,c&#xff0c;x,yx, yx,y为坐标&#xff0c;ccc为该点的价值&#xff0c; 有mmm个询问&#xff0c;每次询问给x,y,cx, y, cx,y,c&#xff0c;要求&#xff0c;点的价值小于等于ccc的条件下&…

Java多线程常用方法 wait 和 notify

一&#xff1a;从一道面试题说起 启动两个线程, 一个输出 1,3,5,7…99, 另一个输出 2,4,6,8…100 最后 STDOUT 中按序输出 1,2,3,4,5…100 要求用 Java 的 wait notify 机制来实现 解法&#xff1a; public class Test {private static Object lock new Object();private st…

dotNET Core实现分布式环境下的流水号唯一

业务背景在管理系统中&#xff0c;很多功能模块都会涉及到各种类型的编号&#xff0c;例如&#xff1a;流程编号、订单号、合同编号等等。编号各有各自的规则&#xff0c;但通常有一个流水号来确定编号的唯一性&#xff0c;保证流水号的唯一&#xff0c;在不同的环境中实现方式…

单位根反演小记

单位根反演 一个等式&#xff1a;[n∣a]1n∑k0n−1wnak[n \mid a] \frac{1}{n} \sum\limits_{k 0} ^{n - 1}w_n ^{ak}[n∣a]n1​k0∑n−1​wnak​ 证明&#xff1a; wnaw_n ^ awna​是nnn次单位根的aaa次方&#xff0c;所以这里是一个公比为wnaw_n ^ awna​的等比数列&…

.NET Core3发布Json API

我们给DNC3&#xff08;.NET Core 3&#xff09;上了一个新包&#xff0c;叫做System.Text.Json(点我下载)&#xff0c;支持读写器&#xff0c;DOM&#xff08;文档对象模型&#xff09;&#xff0c;和序列化&#xff0c;在这篇博文里&#xff0c;我会告诉大家为什么要做这个&a…

Java ThreadLocal

** 一&#xff1a;ThreadLocal的简要介绍及使用 ** Java中的ThreadLocal类允许我们创建只能被同一个线程读写的变量。因此&#xff0c;如果一段代码含有一个ThreadLocal变量的引用&#xff0c;即使两个线程同时执行这段代码&#xff0c;它们也无法访问到对方的ThreadLocal变…