bzoj 2154 Crash的数字表格(莫比乌斯反演及优化)
Description
今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple)。对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数。例如,LCM(6, 8) = 24。回到家后,Crash还在想着课上学的东西,为了研究最小公倍数,他画了一张N*M的表格。每个格子里写了一个数字,其中第i行第j列的那个格子里写着数为LCM(i, j)。一个4*5的表格如下: 1 2 3 4 5 2 2 6 4 10 3 6 3 12 15 4 4 12 4 20 看着这个表格,Crash想到了很多可以思考的问题。不过他最想解决的问题却是一个十分简单的问题:这个表格中所有数的和是多少。当N和M很大时,Crash就束手无策了,因此他找到了聪明的你用程序帮他解决这个问题。由于最终结果可能会很大,Crash只想知道表格里所有数的和mod 20101009的值。
Input
输入的第一行包含两个正整数,分别表示N和M。
Output
输出一个正整数,表示表格中所有数的和mod 20101009的值。
Sample Input
Sample Output
【数据规模和约定】
100%的数据满足N, M ≤ 107。
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std; typedef long long ll;
const int N = 1e7+;
const int MOD = ; ll mu[N],su[N],sz,np[N];
int n,m,mx; void get_mu()
{
int i,j;
mu[]=;
for(i=;i<=mx;i++) {
if(!np[i])
su[++sz]=i,mu[i]=-;
for(j=;j<=sz&&i*su[j]<=mx;j++) {
np[i*su[j]]=;
if(i%su[j]==)
mu[i*su[j]]=;
else
mu[i*su[j]]=-mu[i];
}
}
for(i=;i<=mx;i++)
mu[i]=(mu[i-]+(ll)(mu[i]*i*i)%MOD)%MOD;
}
ll t(ll x,ll y)
{
return ((ll)(x*(x+)/%MOD)*(ll)(y*(y+)/%MOD)%MOD);
}
ll F(int n,int m)
{
int i,last; ll ans=;
for(i=;i<=n;i=last+) {
last=min(n/(n/i),m/(m/i));
ans=(ans+(ll)(mu[last]-mu[i-])*t(n/i,m/i)%MOD)%MOD;
}
return ans;
}
int main()
{
//freopen("in.in","r",stdin);
//freopen("out.out","w",stdout);
scanf("%d%d",&n,&m);
if(n>m) swap(n,m);
mx=n;
get_mu();
int last; ll ans=;
for(int i=;i<=n;i=last+) {
last=min(n/(n/i),m/(m/i));
ans=(ans+((ll)(i+last)*(last-i+)/*F(n/i,m/i)%MOD))%MOD;
}
printf("%lld",(ans+MOD)%MOD);
return ;
}
【优化】
我是图片的搬运工
(H打错为F 233)
积性函数的约数和也是积性函数,H(D) 可以用线性筛求,然后求下前缀和就好了。
至此为止,可以解决2693的多查询问题了。
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std; typedef long long ll;
const int N = 1e7+;
const int MOD = ; ll sum[N],su[N],sz,np[N];
int n,m,mx; void get_mu()
{
int i,j;
sum[]=;
for(i=;i<=mx;i++) {
if(!np[i]) {
su[++sz]=i,
sum[i]=(i-(ll)i*i)%MOD;
}
for(j=;j<=sz&&i*su[j]<=mx;j++) {
np[i*su[j]]=;
if(i%su[j]==)
sum[i*su[j]]=(su[j]*sum[i])%MOD;
else
sum[su[j]*i]=(sum[su[j]]*sum[i])%MOD;
}
}
for(i=;i<=mx;i++)
sum[i]=(sum[i]+sum[i-])%MOD;
}
ll S(ll x,ll y)
{
return ((ll)(x*(x+)/%MOD)*(ll)(y*(y+)/%MOD)%MOD);
} int main()
{
//freopen("in.in","r",stdin);
//freopen("out.out","w",stdout);
mx=1e7;
get_mu();
int T; scanf("%d",&T);
while(T--) {
scanf("%d%d",&n,&m);
if(n>m) swap(n,m);
mx=n; get_mu();
int last; ll ans=;
for(int i=;i<=n;i=last+) {
last=min(n/(n/i),m/(m/i));
ans=(ans+(ll)S(n/i,m/i)*(sum[last]-sum[i-])%MOD)%MOD;
}
printf("%lld\n",(ans+MOD)%MOD);
}
return ;
}
bzoj 2154 Crash的数字表格(莫比乌斯反演及优化)的更多相关文章
- [BZOJ 2154]Crash的数字表格(莫比乌斯反演+数论分块)
[BZOJ 2154]Crash的数字表格(莫比乌斯反演+数论分块) 题面 求 \[\sum_{i=1}^{n} \sum_{j=1}^{m} \mathrm{lcm}(i,j)\] 分析 \[\su ...
- [bzoj 2693] jzptab & [bzoj 2154] Crash的数字表格 (莫比乌斯反演)
题目描述 TTT组数据,给出NNN,MMM,求∑x=1N∑y=1Mlim(x,y)\sum_{x=1}^N\sum_{y=1}^M lim(x,y)\newlinex=1∑Ny=1∑Mlim(x, ...
- BZOJ 2154: Crash的数字表格 [莫比乌斯反演]
2154: Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 2924 Solved: 1091[Submit][Status][ ...
- BZOJ 2154 Crash的数字表格 ——莫比乌斯反演
求$\sum_{i=1}^n\sum_{j=1}^n lcm(i,j)$ 枚举因数 $ans=\sum_{d<=n} F(d) * d$ $F(d)$表示给定范围内两两$\sum_{gcd(i, ...
- 【BZOJ】2154: Crash的数字表格 莫比乌斯反演
[题意]给定n,m,求Σlcm(i,j),1<=i<=n,1<=j<=m,n,m<=10^7. [算法]数论(莫比乌斯反演) [题解] $$ans=\sum_{i\leq ...
- 【莫比乌斯反演】关于Mobius反演与lcm的一些关系与问题简化(BZOJ 2154 crash的数字表格&&BZOJ 2693 jzptab)
BZOJ 2154 crash的数字表格 Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b ...
- Bzoj 2154: Crash的数字表格(积性函数)
2154: Crash的数字表格 Time Limit: 20 Sec Memory Limit: 259 MB Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least ...
- 【刷题】BZOJ 2154 Crash的数字表格
Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数.例如 ...
- 【bzoj2154】Crash的数字表格 莫比乌斯反演
题目描述 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b)表示能同时被a和b整除的最小正整数.例如,LCM(6, ...
- ●BZOJ 2154 Crash的数字表格
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=2154 题解: 莫比乌斯反演. 题意还是很清楚的,就不赘述了. 显然有 $ANS=\sum_{ ...
随机推荐
- 2336: [HNOI2011]任务调度 - BZOJ
一道随机算法的题目 随便用什么随机算法 首先我们可以想到枚举类型3的最终类型,然后再做 先贪心出一个较优的序列,首先我们知道肯定是在A机器上先做完类型1的事件再做类型2的事件,机器B也类似,因为这些没 ...
- Hibernate应用SQL查询返回实体类型
Hibernate应用SQL查询返回实体类型 Hibernate使用SQL查询返回实体类型 以前,使用SQL查询,结果放在 RS 结果集中,还要去转换影射到Java类中.Hibernate中,可以自动 ...
- java模拟OSUnMapTbl[]
问题描述: 任务就绪表,记录当前就绪的任务,就绪表中把64个优先级的任务分成8组,优先级的1-3bit表示OSRdyTbl[]中组别OSRedyGrp,优先级的4-6bit表示每组中就绪任务的位置,当 ...
- Javascript和ECMAScript二三事
来自<javascript高级程序设计 第三版:作者Nicholas C. Zakas>的学习笔记(一) Javascript是一种专为与网页交互而设计的脚本语言,由下列三个不同部分组成: ...
- 可编辑的select框的实现(实用版)
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> <HTML> <HEA ...
- PHP MSSQL数据操作PDO API
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 3 ...
- 数据聚合 & 分组:新一代系统监控的核心功能
遥想 2015 年 8 月 17 日,Cloud Insight 还在梳理功能原型,畅想 Cloud Insight 存在的意义:为什么阿里云用户需要使用 Cloud Insight 来加强管理. 而 ...
- c++调用matlab生成的Dll动态连接库
点击打开链接http://download.csdn.net/detail/nuptboyzhb/4228429 c++调用matlab生成的Dll动态连接库 实验平台: matlab 7.0(R ...
- [itint5]支持删除的后继查询
http://www.itint5.com/oj/#49 这一题一开始想到是用HashSet+链表来做,链表记录prev和next.这样也可以,后来看到都是连续的整数,而且交流了一下觉得可以用类似并查 ...
- 数据库 MySQL Jdbc JDBC的六个固定步骤
*0 案例: a)在JavaScript中使用正则表达式,在JS中正则表达式的定界符是:// var regexp = /^[0-9]+$/; if(regexp.test(nu ...