HDU6715 算术(莫比乌斯反演)
莫比乌斯反演的变形。
对 \(\mu(lcm(i,j))\) 变换,易得 \(\mu(lcm(i,j)) = \mu(i)\cdot\mu(j)\cdot \mu(gcd(i,j))\) 。那么有:
\sum_{i=1}^{n} \sum_{j=1}^{m} \mu(lcm(i,j)) &= \sum_{i=1}^{n}\mu(i) \sum_{j=1}^{m}\mu(j)\cdot \mu(gcd(i,j)) \\
&= \sum_{d=1}^{\min(n,m)}\sum_{i=1}^{n/d}\sum_{j=1}^{m/d}\mu(id)\mu(jd)\mu(d)[gcd(i,j)=1]
\end{split}\]
由于莫比乌斯函数的性质 \(\sum_{d\ |\ n}\mu(d)=[n=1]\) ,我们有:
\text{上式} = \sum_{d=1}^{\min(n,m)}\sum_{d_1 = 1}^{\min(n,m)/d}\sum_{i=1}^{n/dd_1}\sum_{j=1}^{m/dd_1}\mu(idd_1)\mu(jdd_1)\mu(d)\mu(d_1)
\end{split}\]
我们令 \(T = dd_1\) ,有:
\]
令 \(f(T) = \sum_{d|T} \mu(d)\mu(T/d)\) 。
令 \(g(T,N,M)=\sum_{i=1}^{N}\sum_{j=1}^{M}\mu(iT)\mu(jT)=(\sum_{i=1}^{N}\mu(iT))\cdot(\sum_{j=1}^{M}\mu(jT))\) 。
那么我们有:
\]
\(g(T,n/T,m/T)\) 可以在 \(O(n/T+m/T)\) 的时间内计算。
时间复杂度为 \(O(n\log n)\) 。
#include<stdio.h>
#include<algorithm>
using namespace std;
const int maxn = 1000005;
int t, n, m, tot;
int mu[maxn], check[maxn], prime[maxn];
long long ans, f[maxn];
void init()
{
mu[1] = 1; tot = 0;
for(int i = 2; i <= 1000000; i++){
if(check[i] == 0){
mu[i] = -1;
prime[++tot] = i;
}
for(int j = 1; j <= tot && prime[j] * i <= 1000000; j++){
check[i * prime[j]] = 1;
if(i % prime[j] == 0) break;
mu[i * prime[j]] = -1 * mu[i];
}
}
for(int i = 1; i <= 1000000; i++){
for(int j = 1; j * i <= 1000000; j++){
f[i * j] += mu[i] * mu[j];
}
}
}
int main()
{
init();
for(scanf("%d", &t); t--;){
scanf("%d%d", &n, &m);
ans = 0;
for(int T = 1; T <= min(n, m); T++){
if(f[T]){
long long g1 = 0, g2 = 0;
for(int i = 1; i * T <= n; i++) g1 += mu[i * T];
for(int i = 1; i * T <= m; i++) g2 += mu[i * T];
ans += f[T] * g1 * g2;
}
}
printf("%lld\n", ans);
}
return 0;
}
HDU6715 算术(莫比乌斯反演)的更多相关文章
- hdu1695 GCD(莫比乌斯反演)
题意:求(1,b)区间和(1,d)区间里面gcd(x, y) = k的数的对数(1<=x<=b , 1<= y <= d). 知识点: 莫比乌斯反演/*12*/ 线性筛求莫比乌 ...
- BZOJ 2154: Crash的数字表格 [莫比乌斯反演]
2154: Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 2924 Solved: 1091[Submit][Status][ ...
- BZOJ2301: [HAOI2011]Problem b[莫比乌斯反演 容斥原理]【学习笔记】
2301: [HAOI2011]Problem b Time Limit: 50 Sec Memory Limit: 256 MBSubmit: 4032 Solved: 1817[Submit] ...
- Bzoj2154 Crash的数字表格 乘法逆元+莫比乌斯反演(TLE)
题意:求sigma{lcm(i,j)},1<=i<=n,1<=j<=m 不妨令n<=m 首先把lcm(i,j)转成i*j/gcd(i,j) 正解不会...总之最后化出来的 ...
- 莫比乌斯函数筛法 & 莫比乌斯反演
模板: int p[MAXN],pcnt=0,mu[MAXN]; bool notp[MAXN]; void shai(int n){ mu[1]=1; for(int i=2;i<=n;++i ...
- 【BZOJ-2440】完全平方数 容斥原理 + 线性筛莫比乌斯反演函数 + 二分判定
2440: [中山市选2011]完全平方数 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2371 Solved: 1143[Submit][Sta ...
- POI2007_zap 莫比乌斯反演
题意:http://hzwer.com/4205.html 同hdu1695 #include <iostream> #include <cstring> #include & ...
- hdu.5212.Code(莫比乌斯反演 && 埃氏筛)
Code Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others) Total Submi ...
- CSU 1325 莫比乌斯反演
题目大意: 一.有多少个有序数对(x,y)满足1<=x<=A,1<=y<=B,并且gcd(x,y)为p的一个约数: 二.有多少个有序数对(x,y)满足1<=x<=A ...
随机推荐
- linux-memcache安装及memcached memcache扩展
linux memcached安装yum -y install libevent libevent-deve yum list memcached yum -y install memcached m ...
- 管道(Pipe)----计算机进程间通信
参至他人博客:https://blog.csdn.net/u011583316/article/details/83419805
- 移动端和pc端公用样式及rem布局
一:移动端准备工作<meta name="viewport" content="width=device-width, initial-scale=1.0, max ...
- [七月挑选]优化hexo目录,使本地图片能显示出来
title: 优化hexo目录,使本地图片能显示出来 查看了一下从此蜕变作者的Hexo中添加本地图片,提炼了一些能优化本地图片存放及编写是图片查看的问题. 1.修改配置文件_config.yml 里的 ...
- mpg123 - 播放 MPEG 1.0/2.0 Layer-1, -2, -3 音频文件
语法 mpg123 [ -tscvqy01m24 ][ -b size ][ -k num ][ -n num ][ -f factor ][ -r rate ][ -g gain ][ -a dev ...
- 长沙理工大学第十二届ACM大赛-重现赛 D 小M和天平 (01背包)
链接:https://ac.nowcoder.com/acm/contest/1/D 来源:牛客网 小M和天平 时间限制:C/C++ 2秒,其他语言4秒 空间限制:C/C++ 131072K,其他语言 ...
- linux c++下遍历文件
https://blog.csdn.net/u013617144/article/details/44807333
- Mybatis SQL 使用JAVA 静态资源
常量:${@com.htsc.backtest.component.Global@PAGE_SIZE} 静态方法:${@com.htsc.backtest.component.Global@doMet ...
- xtrabackup备份mysql以及创建自动定时任务
xtrabackup的安装 安装依赖关系 wget -O /etc/yum.repos.d/epel.repo http://mirrors.aliyun.com/repo/epel-7.repo y ...
- DevExpress ASP.NET Core Controls 2019发展蓝图(No.5)
本文主要为大家介绍DevExpress ASP.NET Core Controls 2019年的官方发展蓝图,更多精彩内容欢迎持续收藏关注哦~ [DevExpress ASP.NET Controls ...