[51Nod 1238] 最小公倍数之和 (恶心杜教筛)
题目描述
求∑i=1N∑j=1Nlcm(i,j)\sum_{i=1}^N\sum_{j=1}^Nlcm(i,j)i=1∑Nj=1∑Nlcm(i,j)
2<=N<=10102<=N<=10^{10}2<=N<=1010
题目分析
这道题题面跟[bzoj 2693] jzptab & [bzoj 2154] Crash的数字表格一样,然而数据范围加强到了101010^{10}1010,莫比乌斯反演不行了了,所以我们看看怎样玄学杜教筛
Ans=∑i=1n∑j=1nlcm(i,j)=2∑i=1n∑j=1ilcm(i,j)−n(n+1)2Let s(n)=∑i=1n∑j=1ilcm(i,j),f(n)=∑i=1nlcm(i,n)∴f(n)=∑i=1nin(i,n)=n∑i=1ni(i,n)=n∑d∣n∑i=1n[(i,n)==d]id=n∑d∣n∑i=1nd[(i,nd)==1]i=n∑d∣n∑i=1d[(i,d)==1]i=n∑d∣nφ(d)d+[d==1]2\large Ans=\sum_{i=1}^n\sum_{j=1}^nlcm(i,j)=2\sum_{i=1}^n\sum_{j=1}^ilcm(i,j)-\frac{n(n+1)}{2}\\Let~s(n)=\sum_{i=1}^n\sum_{j=1}^ilcm(i,j),f(n)=\sum_{i=1}^nlcm(i,n)\\ \therefore f(n)=\sum_{i=1}^n\frac{in}{(i,n)}=n\sum_{i=1}^n\frac i{(i,n)}\\=n\sum_{d|n}\sum_{i=1}^n[(i,n)==d]\frac id\\=n\sum_{d|n}\sum_{i=1}^{\frac nd}[(i,\frac nd)==1]i\\=n\sum_{d|n}\sum_{i=1}^d[(i,d)==1]i\\=n\sum_{d|n}\frac{\varphi(d)d+[d==1]}2Ans=i=1∑nj=1∑nlcm(i,j)=2i=1∑nj=1∑ilcm(i,j)−2n(n+1)Let s(n)=i=1∑nj=1∑ilcm(i,j),f(n)=i=1∑nlcm(i,n)∴f(n)=i=1∑n(i,n)in=ni=1∑n(i,n)i=nd∣n∑i=1∑n[(i,n)==d]di=nd∣n∑i=1∑dn[(i,dn)==1]i=nd∣n∑i=1∑d[(i,d)==1]i=nd∣n∑2φ(d)d+[d==1]
此处有一个常识
∑i=1ni[(i,n)==1]=φ(n)n+[n==1]2\sum_{i=1}^ni[(i,n)==1]=\frac {\varphi(n)n+[n==1]}2i=1∑ni[(i,n)==1]=2φ(n)n+[n==1]
- 证明如下
- 当n>1n>1n>1时,若(i,n)=1  ⟺  (n−i,n)=1(i,n)=1\iff(n-i,n)=1(i,n)=1⟺(n−i,n)=1,所以与nnn互质的数是成对出现,且他们的和为nnn
- 再加之n=1n=1n=1的特殊情况,可得
∑i=1ni[(i,n)==1]=φ(n)n+[n==1]2\large \sum_{i=1}^ni[(i,n)==1]=\frac {\varphi(n)n+[n==1]}2∑i=1ni[(i,n)==1]=2φ(n)n+[n==1]
继续
∴f(n)=n⋅1+∑d∣nφ(d)d2s(n)=∑i=1nf(i)=∑i=1ni(1+∑d∣iφ(d)d)2=n(n+1)2+∑i=1ni∑d∣iφ(d)d2=n(n+1)2+∑d=1nφ(d)d∑d∣ii2=n(n+1)2+∑d=1nφ(d)d2∑i=1⌊nd⌋i2=n(n+1)2+∑i=1ni∑d=1⌊ni⌋φ(d)d22Ans=2s(n)−n(n+1)2=∑i=1ni∑d=1⌊ni⌋φ(d)d2Let h(d)=φ(d)d2,g(n)=∑d=1nh(d)n=∑d∣nφ(d)n3=∑d∣nφ(d)n2=∑d∣nφ(d)d2(nd)2=∑d∣nh(d)(nd)2∑i=1ni3=∑i=1n∑d∣nh(d)(id)2=∑d=1nh(d)∑d∣i(id)2=∑d=1nh(d)∑i=1⌊nd⌋i2=∑i=1ni2∑d=1⌊ni⌋h(d)=∑i=1ni2g(⌊ni⌋)g(n)=∑i=1ni3−∑i=2ni2g(⌊ni⌋)=(n(n+1)2)2−∑i=2ni2g(⌊ni⌋)
\large \therefore f(n)=n\cdot\frac {1+\sum_{d|n}\varphi(d)d}2\\s(n)=\sum_{i=1}^nf(i)=\frac{\sum_{i=1}^ni(1+\sum_{d|i}\varphi(d)d)}2\\=\frac{\frac{n(n+1)}2+\sum_{i=1}^ni\sum_{d|i}\varphi(d)d}2\\=\frac{\frac{n(n+1)}2+\sum_{d=1}^n\varphi(d)d\sum_{d|i}i}2\\=\frac{\frac{n(n+1)}2+\sum_{d=1}^n\varphi(d)d^2\sum_{i=1}^{\lfloor\frac nd\rfloor}i}2\\=\frac{\frac{n(n+1)}2+\sum_{i=1}^ni\sum_{d=1}^{\lfloor\frac ni\rfloor}\varphi(d)d^2}2\\Ans=2s(n)-\frac{n(n+1)}2=\sum_{i=1}^ni\sum_{d=1}^{\lfloor\frac ni\rfloor}\varphi(d)d^2\\Let~h(d)=\varphi(d)d^2,g(n)=\sum_{d=1}^nh(d)\\n=\sum_{d|n}\varphi(d)\\n^3=\sum_{d|n}\varphi(d)n^2\\=\sum_{d|n}\varphi(d)d^2(\frac nd)^2\\=\sum_{d|n}h(d)(\frac nd)^2\\\sum_{i=1}^ni^3=\sum_{i=1}^n\sum_{d|n}h(d)(\frac id)^2\\=\sum_{d=1}^nh(d)\sum_{d|i}(\frac id)^2\\=\sum_{d=1}^nh(d)\sum_{i=1}^{\lfloor\frac nd\rfloor}i^2\\=\sum_{i=1}^ni^2\sum_{d=1}^{\lfloor\frac ni\rfloor}h(d)\\=\sum_{i=1}^ni^2g(\lfloor\frac ni\rfloor)\\g(n)=\sum_{i=1}^ni^3-\sum_{i=2}^ni^2g(\lfloor\frac ni\rfloor)\\=(\frac{n(n+1)}2)^2-\sum_{i=2}^ni^2g(\lfloor\frac ni\rfloor)∴f(n)=n⋅21+∑d∣nφ(d)ds(n)=i=1∑nf(i)=2∑i=1ni(1+∑d∣iφ(d)d)=22n(n+1)+∑i=1ni∑d∣iφ(d)d=22n(n+1)+∑d=1nφ(d)d∑d∣ii=22n(n+1)+∑d=1nφ(d)d2∑i=1⌊dn⌋i=22n(n+1)+∑i=1ni∑d=1⌊in⌋φ(d)d2Ans=2s(n)−2n(n+1)=i=1∑nid=1∑⌊in⌋φ(d)d2Let h(d)=φ(d)d2,g(n)=d=1∑nh(d)n=d∣n∑φ(d)n3=d∣n∑φ(d)n2=d∣n∑φ(d)d2(dn)2=d∣n∑h(d)(dn)2i=1∑ni3=i=1∑nd∣n∑h(d)(di)2=d=1∑nh(d)d∣i∑(di)2=d=1∑nh(d)i=1∑⌊dn⌋i2=i=1∑ni2d=1∑⌊in⌋h(d)=i=1∑ni2g(⌊in⌋)g(n)=i=1∑ni3−i=2∑ni2g(⌊in⌋)=(2n(n+1))2−i=2∑ni2g(⌊in⌋)然后就是杜教筛的形式了,上杜教筛即可.先预处理出小范围的ggg然后较大的就用杜教筛计算
又Ans=∑i=1ni⋅g(⌊ni⌋)\large Ans=\sum_{i=1}^ni\cdot g(\lfloor\frac ni\rfloor)Ans=∑i=1ni⋅g(⌊in⌋)
因为ggg函数求解时是用的记忆化,所以在外面套上一层分块优化不会影响ggg函数的时间复杂度,所以复杂度为Θ(n23)\Theta(n^{\frac 23})Θ(n32)
AC code
#include <cstdio>
#include <map>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long LL;
const int mod = 1e9+7;
const int MAXN = 5e6+1;
const int inv2 = 500000004;
const int inv3 = 333333336;
map<LL, LL> G; LL g[MAXN];
int Prime[MAXN], Cnt, phi[MAXN];
bool IsnotPrime[MAXN];
void init()
{
phi[1] = 1;
for(int i = 2; i < MAXN; ++i)
{
if(!IsnotPrime[i]) Prime[++Cnt] = i, phi[i] = i-1;
for(int j = 1; j <= Cnt && i * Prime[j] < MAXN; ++j)
{
IsnotPrime[i * Prime[j]] = 1;
if(i % Prime[j] == 0)
{
phi[i * Prime[j]] = phi[i] * Prime[j];
break;
}
phi[i * Prime[j]] = phi[i] * phi[Prime[j]];
}
}
for(int i = 1; i < MAXN; ++i) g[i] = (g[i-1] + 1ll * phi[i] * i % mod * i % mod) % mod;
}
inline LL sum2(LL i) { return (i%mod) * ((i+1)%mod) % mod * ((2*i+1)%mod) % mod * inv2 % mod * inv3 % mod; }
inline LL calc(LL n)
{
if(n < MAXN) return g[n];
if(G.count(n)) return G[n];
LL ret = (n%mod) * ((n+1)%mod) % mod * inv2 % mod;
ret = ret * ret % mod;
for(LL i = 2, j; i <= n; i=j+1)
{
j = n/(n/i);
ret = (ret - (sum2(j)-sum2(i-1)) % mod * calc(n/i) % mod) % mod;
}
return G[n] = ret;
}
inline LL sum(LL i, LL j) { return ((i+j)%mod) * ((j-i+1)%mod) % mod * inv2 % mod; }
inline LL solve(LL n)
{
LL ret = 0;
for(LL i = 1, j; i <= n; i=j+1)
{
j = n/(n/i);
ret = (ret + sum(i, j) * calc(n/i) % mod) % mod;
}
return ret;
}
int main ()
{
init(); LL n;
scanf("%lld", &n);
printf("%lld\n", (solve(n)+mod)%mod);
}
[51Nod 1238] 最小公倍数之和 (恶心杜教筛)的更多相关文章
- 51NOD 1238 最小公倍数之和 V3 [杜教筛]
1238 最小公倍数之和 V3 三种做法!!! 见学习笔记,这里只贴代码 #include <iostream> #include <cstdio> #include < ...
- 【51nod】1238 最小公倍数之和 V3 杜教筛
[题意]给定n,求Σi=1~nΣj=1~n lcm(i,j),n<=10^10. [算法]杜教筛 [题解]就因为写了这个非常规写法,我折腾了3天…… $$ans=\sum_{i=1}^{n}\s ...
- 51 Nod 1238 最小公倍数之和 V3 杜教筛
题目链接:http://www.51nod.com/Challenge/Problem.html#!#problemId=1238 题意:求$\sum_{i=1}^{n}\sum_{j=1}^{n}l ...
- [51Nod1238]最小公倍数之和 V3[杜教筛]
题意 给定 \(n\) ,求 \(\sum_{i=1}^n \sum_{j=1}^n lcm(i,j)\). \(n\leq 10^{10}\) 分析 推式子 \[\begin{aligned} an ...
- 51nod 1238 最小公倍数之和 V3
51nod 1238 最小公倍数之和 V3 求 \[ \sum_{i=1}^N\sum_{j=1}^N lcm(i,j) \] \(N\leq 10^{10}\) 先按照套路推一波反演的式子: \[ ...
- 51nod 1238 最小公倍数之和 V3 【欧拉函数+杜教筛】
首先题目中给出的代码打错了,少了个等于号,应该是 G=0; for(i=1;i<=N;i++) for(j=1;j<=N;j++) { G = (G + lcm(i,j)) % 10000 ...
- 51Nod 1238 - 最小公倍数之和 V3(毒瘤数学+杜教筛)
题目 戳这里 推导 ∑i=1n∑j=1nlcm(i,j)~~~\sum_{i=1}^{n}\sum_{j=1}^{n}lcm(i,j) ∑i=1n∑j=1nlcm(i,j) =∑i=1n∑j= ...
- 51nod 1244 莫比乌斯函数之和 【杜教筛】
51nod 1244 莫比乌斯函数之和 莫比乌斯函数,由德国数学家和天文学家莫比乌斯提出.梅滕斯(Mertens)首先使用μ(n)(miu(n))作为莫比乌斯函数的记号.具体定义如下: 如果一个数包含 ...
- 51nod 1244 莫比乌斯函数之和(杜教筛)
[题目链接] http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1244 [题目大意] 计算莫比乌斯函数的区段和 [题解] 利 ...
随机推荐
- LINUX驱动笔记 目录
笔记参考了宋宝华老师的<Linux设备驱动开发详解:基于最新的Linux 4.0内核>以及韦东山老师的嵌入式驱动教程 笔记开发环境: 单板:第一章到第八章使用TINY4412-1611:第 ...
- python3 内置方法 字符串转换为字典
内置方法:eval()将字符串转换为字典代码: str = '''{'backend':'www.oldboy.org', 'record':{ 'server':'122.111.2.23', 'w ...
- UnknownError: session deleted because of page crash from tab crashed
一.问题 在docker上跑Selenium+ChromeDriver+Chrome无头模式报错: UnknownError: unknown error: session deleted becau ...
- The Preliminary Contest for ICPC Asia Xuzhou 2019 E XKC's basketball team [单调栈上二分]
也许更好的阅读体验 \(\mathcal{Description}\) 给n个数,与一个数m,求\(a_i\)右边最后一个至少比\(a_i\)大\(m\)的数与这个数之间有多少个数 \(2\leq n ...
- javascript:void(0);的含义以及使用场景
一.含义: javascript:是伪协议,表示内容通过javascript执行. void(0)表示不作任何操作. 二.使用场景 1.href=”javascript:void(0);” 作用:为了 ...
- ESlint 格式化代码 备忘
vscode 代码格式化配置 vscode 菜单 文件->首选项->设置 --->进入扩展查找到ESlint,点击任一选项中的[在setting.json中配置],复制以下代码 { ...
- Go net/http,web server
net/http 包实现 HTTP Server Go 中,实现一个最简单的 http server 非常容易,代码如下: package main import ( "fmt" ...
- Python爬虫之BeautifulSoap的用法
1. Beautiful Soup的简介 简单来说,Beautiful Soup是python的一个库,最主要的功能是从网页抓取数据.官方解释如下: Beautiful Soup提供一些简单的.pyt ...
- 解决在web.xml中配置server服务器启动失败问题
一.问题"Server Tomacat v8.5 Server at locallhost failed to start" 二.解决方法:删除注释@webServlet 三.分析 ...
- [iOS]UIButton内、外边距设置
- (void)viewDidLoad { [super viewDidLoad]; /* UIButton设置对应的边距image跟title的边距属性 ...