51nod 1220 约数之和【莫比乌斯反演+杜教筛】
首先由这样一个式子:\( d(ij)=\sum_{p|i}\sum_{q|j}[gcd(p,q)==1]\frac{pj}{q} \)大概感性证明一下吧我不会证
然后开始推:
\[ \sum_{i=1}^{n}\sum_{j=1}^{n}\sum_{p|i}\sum_{q|j}[gcd(p,q)==1]\frac{pj}{q} \]
\[ \sum_{p=1}^{n}\sum_{q=1}^{n}[gcd(p,q)==1]\sum_{p|i}\sum_{q|j}\frac{pj}{q} \]
\[ \sum_{p=1}^{n}p\sum_{q=1}^{n}[gcd(p,q)==1]\left \lfloor \frac{n}{p} \right \rfloor\sum_{j=1}^{\left \lfloor \frac{n}{q} \right \rfloor}j \]
方便起见设\( f(n)=\sum_{i=1}^{n}i \)
\[ \sum_{p=1}^{n}p\sum_{q=1}^{n}\sum_{k|p,k|q}\mu(k)\left \lfloor \frac{n}{p} \right \rfloor f(\left \lfloor \frac{n}{q} \right \rfloor) \]
\[ \sum_{k=1}^{n}\mu(k)\sum_{k|p}p\left \lfloor \frac{n}{p} \right \rfloor\sum_{k|q}f(\left \lfloor \frac{n}{q} \right \rfloor) \]
\[ \sum_{k=1}^{n}\mu(k)\sum_{i=1}^{\left \lfloor \frac{n}{k} \right \rfloor}ik\left \lfloor \frac{n}{ik} \right \rfloor\sum_{j=1}^{\left \lfloor \frac{n}{k} \right \rfloor}f(\left \lfloor \frac{n}{jk} \right \rfloor) \]
\[ \sum_{k=1}^{n}\mu(k)k\sum_{i=1}^{\left \lfloor \frac{n}{k} \right \rfloor}i\left \lfloor \frac{n}{ik} \right \rfloor\sum_{j=1}^{\left \lfloor \frac{n}{k} \right \rfloor}f(\left \lfloor \frac{n}{jk} \right \rfloor) \]
这个样子显然可以用杜教筛了,但是注意到后面有两个求和式,可能会增大常数(但是也不会T啦),所以考虑这两个求和式的关系:
\[ \sum_{i=1}^{n}f(\left \lfloor \frac{n}{i} \right \rfloor) \]
\[ =\sum_{i=1}^{n}\sum_{j=1}^{\left \lfloor \frac{n}{i} \right \rfloor}j \]
\[ =\sum_{i=1}^{n}\sum_{j=1}^{\left \lfloor \frac{n}{i} \right \rfloor}j \]
\[ =\sum_{j=1}^{n}j\left \lfloor \frac{n}{j} \right \rfloor \]
所以这两个式子是一样的!于是就变成了:
\[ \sum_{k=1}^{n}\mu(k)k(\sum_{j=1}^{\left \lfloor \frac{n}{k} \right \rfloor}f(\left \lfloor \frac{n}{jk} \right \rfloor))^2 \]
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int N=1000005,inv2=500000004,mod=1e9+7;
int m,mb[N],q[N],tot;
long long n,s[N],ans,ha[N];
bool v[N];
long long slv(long long n)
{return n*(n+1)%mod*inv2%mod;
}
long long wk(long long x)
{if(x<=m)return s[x];//cout<<x<<endl;if(ha[n/x])return ha[n/x];long long re=1ll;for(int i=2,la;i<=x;i=la+1){la=x/(x/i);re=(re-(slv(la)-slv(i-1))*wk(x/i)%mod)%mod;}return ha[n/x]=re;
}
long long clc(long long n)
{long long re=0ll;for(int i=1,la;i<=n;i=la+1){la=n/(n/i);re=(re+(la-i+1)*slv(n/i)%mod)%mod;}return re;
}
int main()
{scanf("%lld",&n);m=(int)ceil(pow((int)n,2.0/3));mb[1]=1;for(int i=2;i<=m;i++){if(!v[i]){q[++tot]=i;mb[i]=-1;}for(int j=1;j<=tot&&q[j]*i<=m;j++){int k=i*q[j];v[k]=1;if(i%q[j]==0){mb[k]=0;break;}mb[k]=-mb[i];}}for(int i=1;i<=m;i++)s[i]=(s[i-1]+i*mb[i])%mod;//cout<<wk(102)<<" "<<wk(101)<<endl;for(int i=1,la;i<=n;i=la+1){la=n/(n/i);long long ml=clc(n/i);//if(i!=1)cout<<i-1<<" "<<n/(i-1)<<endl<<la<<" "<<n/la<<endl;ans=(ans+(wk(la)-wk(i-1))*ml%mod*ml%mod)%mod;}printf("%lld",(ans%mod+mod)%mod);return 0;
}
转载于:https://www.cnblogs.com/lokiii/p/8334216.html
51nod 1220 约数之和【莫比乌斯反演+杜教筛】相关推荐
- [复习]莫比乌斯反演,杜教筛,min_25筛
[复习]莫比乌斯反演,杜教筛,min_25筛 莫比乌斯反演 做题的时候的常用形式: \[\begin{aligned}g(n)&=\sum_{n|d}f(d)\\f(n)&=\sum_ ...
- 【bzoj4176】Lucas的数论 莫比乌斯反演+杜教筛
题目描述 去年的Lucas非常喜欢数论题,但是一年以后的Lucas却不那么喜欢了. 在整理以前的试题时,发现了这样一道题目"求Sigma(f(i)),其中1<=i<=N" ...
- 51nod1220-约数之和【莫比乌斯反演,杜教筛】
正题 题目链接:http://www.51nod.com/Challenge/Problem.html#problemId=1220 题目大意 给出nnn,求∑i=1n∑j=1nσ(i∗j)\sum_ ...
- matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...
数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...
- 牛客练习赛84F-牛客推荐系统开发之下班【莫比乌斯反演,杜教筛】
正题 题目链接:https://ac.nowcoder.com/acm/contest/11174/F 题目大意 给出n,kn,kn,k求 ∑i1=1n∑i2=1n...∑ik=1ngcd(fi1,f ...
- 【LOJ#572】Misaka Network 与求和(莫比乌斯反演/杜教筛/min_25筛)
[LOJ#572]Misaka Network 与求和 https://www.cnblogs.com/cjyyb/p/10170630.html 看到次大质因子就可以想到是min_25筛了,然后只需 ...
- 牛客P19836 裴蜀定理+莫比乌斯反演+杜教筛
题意: 一开始一个人在原点,它拥有 n n n个整数 x i x_{i} xi,并且 x i ∈ [ 1 , m ] x_{i}\in[1,m] xi∈[1,m],他每次可以选择一个 x i x_ ...
- CCPC-Wannafly Winter Camp Day3 (Div2, onsite) F 小清新数论 欧拉函数的利用 莫比乌斯反演 杜教筛
F - 小清新数论 做法一:欧拉函数 #include<stdio.h> #include<bits/stdc++.h> using namespace std; #defin ...
- Loj #572. 「LibreOJ Round #11」Misaka Network 与求和(莫比乌斯反演 + 杜教筛 + min_25筛(递推版))
直接反演一下:∑i=1n∑i=1nf(gcd(i,j))k\sum_{i = 1}^n\sum_{i = 1}^nf(gcd(i,j))^ki=1∑ni=1∑nf(gcd(i,j))k=∑d=1n ...
最新文章
- [20180806]tune2fs调整保留块百分比.txt
- 原生js设置div隐藏或者显示_10种JS控制DIV的显示隐藏代码
- 【Python】Paramiko模块在Windows10中import ssl报错的处理方法
- Linux mcheck机制检测内存溢出、内存越界
- .net与.net core学习目录
- 20190908:(leetcode习题)最大子序和
- 八皇后(N皇后)问题算法程序(回溯法)
- iOS开发之旧版本Xcode下载
- PostgreSQL on Docker部署是的shm配置注意
- 数字图像处理实验一 直方图均衡化与直方图规定化
- 【总结】计算机网络常用协议总结------结合cpri
- 高端游戏计算机电源外形规格,游戏新高度!实述高端电源选购技巧
- Java注解中produces啥意思_SpringBoot常用注解
- 森林中的兔子java
- 学计算机土味情话,很套路的土味情话
- [英语阅读]你的英文名特别吗
- 味知香VS千味央厨,谁是年夜饭C位?
- 计算机内存tm,Intel(R) Core(TM) i3-3110M CPU @ 2.40GHz支持多大内存
- leetcode 笨阶乘
- 三星4k3d电视测试软件,技术篇 | 央视4K超高清频道怎么看?三星电视让您大饱眼福!...
热门文章
- Redis 笔记(05)— hash 类型(设置获取单个hash值、判断hash键是否存在、设置获取多个hash值、获取所有hash的键、获取所有hash的值、删除hash指定键)
- centos 安装 node12以及yarn
- mybatis配置文件解析
- Ubuntu使用wget下载zabbix的扩展源
- 使用哈工大LTP进行文本命名实体识别并保存到txt
- Python学习--not语句
- 接触追踪解决方案建立在UWB而不是蓝牙上
- 分层条件关系网络在视频问答VideoQA中的应用:CVPR2020论文解析
- 2021年大数据HBase(八):Apache Phoenix的基本介绍
- android数据库isnull,Android中SQLite数据库知识点总结