●BZOJ 3529 [Sdoi2014]数表
题链:
http://www.lydsy.com/JudgeOnline/problem.php?id=3529
题解:
莫比乌斯反演。
按题目的意思,令$f(i)$表示i的所有约数的和,就是要求:
$ANS=\sum f(gcd(i,j)),满足1 \leq i \leq n,1 \leq j \leq m,且 f(gcd(i,j))\leq a$
首先 $f(i)$ 应该还是比较好推的,利用其为积性函数的特点,可以在线性筛时完成计算。
令$g[k]$表示$gcd(i,j)=k$的$(i,j)$的对数
$G[k]$表示$gcd(i,j)=\lambda k$的$(i,j)$的对数,其值$G[k]=\lfloor \frac{n}{k} \rfloor \lfloor \frac{m}{k} \rfloor$
那么显然,$G[k]$为$g[k]$的倍数关系和函数,
即满足$G[k]=\sum_{k|d} g[d]$
则由莫比乌斯反演得:
$g[k]=\sum_{k|d}\mu(\frac{d}{k})G[d]$
$\quad\quad=\sum_{k|d}\mu(\frac{d}{k})\lfloor \frac{n}{d} \rfloor \lfloor \frac{m}{d} \rfloor$
那么现在,直接从gcd的值的角度出发,ANS可以写成如下形式:
$ANS=\sum_{i=1}^{min(n,m)}f(i)g(i)$
$\quad\quad=\sum_{i=1}^{min(n,m)}f(i)\sum_{i|d}\mu(\frac{d}{i})\lfloor \frac{n}{d} \rfloor \lfloor \frac{m}{d} \rfloor$
然后再化一下:
$\quad\quad=\sum_{d=1}^{min(n,m)}\lfloor \frac{n}{d} \rfloor \lfloor \frac{m}{d} \rfloor\sum_{i|d}f(i)\mu(\frac{d}{i})$
令 $w(d)=\sum_{i|d}f(i)\mu(\frac{d}{i})$
那么$ANS=\sum_{d=1}^{min(n,m)}\lfloor \frac{n}{d} \rfloor \lfloor \frac{m}{d} \rfloor w(d)$
如果不考虑题目中$f(gcd(i,j))\leq a$的限制
我们就可以枚举每个i,然后把其倍数$d=\lambda i$的$w(d)+=f(i)\mu(\frac{d}{i})$
以此计算出所有的w(d),复杂度为O(Nlog_2N)的。
然后那个求ANS的式子就可以运用向下取整的特性进行分块计算,就可以达到每个询问$O(\sqrt N)$的复杂度。
再来考虑有a的限制时的做法,(其实也不麻烦)
离线询问,按a从小到大排序,
同时把f(i)按从小到大排序,
一次遍历每个询问,并把$f(i)$小于当前询问的$a$的$i$按之前的做法:枚举倍数,加入对应的$w(d)$。
但是为了维护好前缀和,以便使用分块计算,
所以用树状数组维护,即把值$f(i)\mu(\frac{d}{i})$加入到树状数组里面。
然后同样的用树状数组查询前缀和就可以继续对当前询问进行分块计算了。
复杂度:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define MAXN 100050
using namespace std;
struct BIT{
int val[MAXN],N;
void Reset(int n){memset(val,0,sizeof(val));N=n;}
int Lowbit(int p){return p&-p;}
void Modify(int p,int v){
while(p<=N) val[p]+=v,p+=Lowbit(p);
}
int Query(int p,int ret=0){
while(p) ret+=val[p],p-=Lowbit(p);
return ret;
}
}DT;
struct Question{
int n,m,a,id;
friend bool operator < (Question A,Question B){
return A.a<B.a;
}
}Q[MAXN];
int f[MAXN],fi[MAXN],mu[MAXN],ANS[MAXN];
bool cmp(int i,int j){
return f[i]<f[j];
}
void Sieve(){
static bool np[MAXN];
static int prime[MAXN],pnt,tmp;
mu[1]=f[1]=fi[1]=1;
for(int i=2;i<=100000;i++){
fi[i]=i;
if(!np[i]) prime[++pnt]=i,mu[i]=-1,f[i]=1+i;
for(int j=1;j<=pnt&&i<=100000/prime[j];j++){
np[i*prime[j]]=1;
if(i%prime[j]){
mu[i*prime[j]]=-mu[i];
f[i*prime[j]]=f[i]*f[prime[j]];
}
else{
mu[i*prime[j]]=0;
tmp=i; while(!(tmp%prime[j])) tmp=tmp/prime[j];
f[i*prime[j]]=f[i]*prime[j]+f[tmp];
}
if(i%prime[j]==0) break;
}
}
sort(fi+1,fi+100000+1,cmp);
}
int main(){
Sieve(); int Case,ret,mini,n,m;
scanf("%d",&Case);
DT.Reset(100000);
for(int i=1;i<=Case;i++)
scanf("%d%d%d",&Q[i].n,&Q[i].m,&Q[i].a),Q[i].id=i;
sort(Q+1,Q+Case+1);
for(int q=1,p=1;q<=Case;q++){
while(p<=100000&&f[fi[p]]<=Q[q].a){
for(int d=fi[p];d<=100000;d+=fi[p])
DT.Modify(d,f[fi[p]]*mu[d/fi[p]]);
p++;
}
ret=0; n=Q[q].n; m=Q[q].m; mini=min(n,m);
for(int i=1,last;i<=mini;i=last+1){
last=min(n/(n/i),m/(m/i));
ret+=(DT.Query(last)-DT.Query(i-1))*(n/i)*(m/i);
}
ANS[Q[q].id]=ret;
}
for(int i=1;i<=Case;i++)
printf("%d\n",ANS[i]&(int)((1ll<<31)-1));
return 0;
}
●BZOJ 3529 [Sdoi2014]数表的更多相关文章
- BZOJ 3529: [Sdoi2014]数表 [莫比乌斯反演 树状数组]
3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1399 Solved: 694[Submit][Status] ...
- bzoj 3529 [Sdoi2014]数表(莫比乌斯反演+BIT)
Description 有一张N×m的数表,其第i行第j列(1 < =i < =礼,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a ...
- 【刷题】BZOJ 3529 [Sdoi2014]数表
Description 有一张n×m的数表,其第i行第j列(1<=i<=n,1<=j<=m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a的数之和. In ...
- BZOJ 3529 [Sdoi2014]数表 (莫比乌斯反演+树状数组+离线)
题目大意:有一张$n*m$的数表,第$i$行第$j$列的数是同时能整除$i,j$的所有数之和,求数表内所有不大于A的数之和 先是看错题了...接着看对题了发现不会做了...刚了大半个下午无果 看了Po ...
- bzoj 3529: [Sdoi2014]数表
#include<cstdio> #include<iostream> #include<algorithm> #define M 200009 //#define ...
- BZOJ 3529 [Sdoi2014]数表 ——莫比乌斯反演 树状数组
$ans=\sum_{i=1}^n\sum_{j=1}^n\sigma(gcd(i,j))$ 枚举gcd为d的所有数得到 $ans=\sum_{d<=n}\sigma(d)*g(d)$ $g(d ...
- BZOJ 3259 [Sdoi2014]数表 (莫比乌斯反演 + 树状数组)
3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 2321 Solved: 1187[Submit][Status ...
- 3529: [Sdoi2014]数表 - BZOJ
Description 有一张N×m的数表,其第i行第j列(1 < =i < =n,1 < =j < =m)的数值为能同时整除i和j的所有自然数之和.给定a,计算数表中不大于a ...
- 【BZOJ】3529: [Sdoi2014]数表
题意:求 $$\sum_{i=1}^{n} \sum_{j=1}^{m} \sum_{d|(i, j)} d 且 (\sum_{d|(i, j)} d)<=a$$ n, m<=1e5,q次 ...
随机推荐
- C语言第七次作业
一.PTA实验作业 题目1:求整数序列中出现次数最多的数 1.本题PTA提交列表 2.设计思路 定义一个整型数组a[1001],i,j 为循环变量,N,定义数组b[1001]={0} 输入N for( ...
- 2017-2018-1 20155215 第九周 加分项 PWD命令的实现
1 学习pwd命令 Linux中用 pwd 命令来查看"当前工作目录"的完整路径. 简单得说,每当你在终端进行操作时,你都会有一个当前工作目录. 在不太确定当前位置时,就会使用pw ...
- 2018上c语言第0次作业
随笔: 1.翻阅邹欣老师博客关于师生关系博客,并回答下列问题,每个问题的答案不少于500字: (1)最理想的师生关系是健身教练和学员的关系,在这种师生关系中你期望获得来自老师的哪些帮助? 答:对此问题 ...
- Linux下进程间通信--消息队列
消息队列的定义遍地都是,不想移驾,请看下文: 一.定义: 消息队列提供了一种从一个进程向另一个进程发送一个数据块的方法. 每个数据块都被认 为是有一个类型,接收者进程接收的数据块可以有不同的类型值.我 ...
- hashlib 加密
import hashlib def md5(args): hash = hashlib.md5(bytes('aaadf',encoding='utf-8')) hash.update(bytes( ...
- HTML 字符集
在 HTML 中,正确的字符编码是什么? HTML5 中默认的字符编码是 UTF-8. 这并非总是如此.早期网络的字符编码是 ASCII 码.后来,从 HTML 2.0 到 HTML 4.01,I ...
- New UWP Community Toolkit - DropShadowPanel
概述 UWP Community Toolkit 中有一个为 Frmework Element 提供投影效果的控件 - DropShadowPanel,本篇我们结合代码详细讲解 DropShado ...
- hadoop2.6.0实践:A03 例子验证
[hadoop@LexiaofeiN1 ~]$ hdfs dfs -ls /output/grep[hadoop@LexiaofeiN1 ~]$ hdfs dfs -rm -R /output/gre ...
- springmvc4开发rest
Spring MVC 4 RESTFul Web Services CRUD Example+RestTemplate Created on: August 11, 2015 | Last upd ...
- UnicodeEncodeError: 'gbk' codec can't encode character '\xa0' in position 1987: illegal multibyte sequence
在爬取 url = "http://stats.meizhou.gov.cn/show/index/1543/1689" 时出现了问题: UnicodeEncodeError: ' ...