BZOJ 3529 【SDOI2014】 数表
题目链接:数表
我们一起来膜PoPoQQQ大爷的题解吧Orz
首先我们来考虑没有\(a\)的限制该怎么做。显然交换\(n\),\(m\)答案不变,所以后面默认\(n \le m\)。
我们定义两个函数:
\[f(x)=\sum_{d|x}d\]
\[g(x)=\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=x]\]
那么显然有:\[ans=\sum_{i=1}^nf(i)g(i)\]
\(g\)函数我们可以考虑化简一下:
\begin{aligned}
g(x)&=\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=x] \\
&=\sum_{i=1}^{\lfloor \frac{n}{x} \rfloor}\sum_{j=1}^{\lfloor \frac{m}{x} \rfloor}\sum_{d|i,d|j}\mu(d) \\
&=\sum_{d=1}^{\lfloor \frac{n}{x} \rfloor}\mu(d)\lfloor \frac{n}{dx} \rfloor\lfloor \frac{m}{dx} \rfloor \\
&=\sum_{x|d}\mu(\frac{d}{x})\lfloor \frac{n}{d} \rfloor\lfloor \frac{m}{d} \rfloor
\end{aligned}
于是可以得到:
\begin{aligned}
ans=&\sum_{i=1}^nf(i)\sum_{i|d}\mu(\frac{d}{i})\lfloor \frac{n}{d} \rfloor\lfloor \frac{m}{d} \rfloor \\
=&\sum_{d=1}^n\lfloor \frac{n}{d} \rfloor\lfloor \frac{m}{d} \rfloor\sum_{i|d}f(i)\mu(\frac{d}{i})
\end{aligned}
我们令\(h(x)=\sum_{i|x}f(i)\mu(\frac{x}{i})\),那么我们只要有\(h(x)\)的前缀和,就可以在\(O(\sqrt{n})\)的时间内求出\(ans\)了。
然后我们来考虑\(a\)的限制。显然,只有\(f(x)\le a\)的\(f(x)\)才会对答案有影响。所以,我们把所有的询问按\(a\)从小到大排好序,并且把所有的\(f(x)\)从小到大依次加进来,使用树状数组维护前缀和,依次处理每个询问,就可以在\(O(Q\sqrt{n}\log{n}+n\log^2 n)\)的时间内解决所有询问了。
下面贴代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
#define maxn 100010 using namespace std;
typedef long long llg; struct data{
int n,m,a,b;
bool operator < (const data &h)const{return a<h.a;}
}s[maxn];
int n,m,nm,a[maxn],T,c[maxn],ans[maxn];
int mu[maxn],f[maxn],pr[maxn],lp;
const int mod=2147483647;
bool vis[maxn]; int getint(){
int w=0;bool q=0;
char c=getchar();
while((c>'9'||c<'0')&&c!='-') c=getchar();
if(c=='-') c=getchar(),q=1;
while(c>='0'&&c<='9') w=w*10+c-'0',c=getchar();
return q?-w:w;
} bool cmp(int x,int y){return f[x]<f[y];}
void add(int x,int y){if(y) while(x<=nm) c[x]+=y,x+=x&(-x);}
int sum(int x){
int t=0;
while(x) t+=c[x],x-=x&(-x);
return t;
} int main(){
File("a");
mu[1]=f[1]=a[1]=1; T=getint();
for(int i=1;s[i].b=i,i<=T;i++){
s[i].n=getint(),s[i].m=getint();
if(s[i].n>s[i].m) swap(s[i].n,s[i].m);
s[i].a=getint(); nm=max(nm,s[i].n);
}
for(int i=2;a[i]=i,i<=nm;i++){
if(!vis[i]) pr[++lp]=i,mu[i]=-1,f[i]=i+1;
for(int j=1;pr[j]*i<=nm;j++){
vis[pr[j]*i]=1;
if(i%pr[j]) mu[pr[j]*i]=-mu[i],f[pr[j]*i]=f[i]*f[pr[j]];
else{ f[pr[j]*i]=f[i]+(f[i]-f[i/pr[j]])*pr[j]; break;}
}
}
sort(s+1,s+T+1); sort(a+1,a+nm+1,cmp);
for(int i=1,now=1,b,la=0;i<=T;i++){
n=s[i].n; m=s[i].m; b=s[i].b;
while(now<=nm && f[a[now]]<=s[i].a){
for(int j=a[now];j<=nm;j+=a[now])
add(j,f[a[now]]*mu[j/a[now]]);
now++;
}
for(int j=1,nt,na;j<=n;j=nt+1){
nt=min(n/(n/j),m/(m/j)); na=sum(nt);
ans[b]+=(na-la)*(n/j)*(m/j); la=na;
}
}
for(int i=1;i<=T;i++) printf("%d\n",ans[i]&mod);
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]数表
题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3529 题解: 莫比乌斯反演. 按题目的意思,令$f(i)$表示i的所有约数的和,就是要求: ...
- 【刷题】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次 ...
随机推荐
- URL上的中文编码
参考:http://www.chinaz.com/web/2013/0226/293639.shtml 一 URL上拼接中文会进行编码 在URL上拼接中文时,比如www.abc.com?season= ...
- PHP pdf转化为图片(PNG)
/** * 将pdf文件转化为多张png图片 * @param string $pdf pdf所在路径 (/www/pdf/abc.pdf pdf所在的绝对路径) * @param string $p ...
- android麦克风自录自放demo
extends:http://blog.csdn.net/trbbadboy/article/details/7865530 是一个直接播放麦克风采集到的声音线程类: class RecordThre ...
- thinkphp --- 写入日志
在开发过程中,对于一些参数,不好直接输入或者打印调试,特别是在微信开发过程中,这个时候,通过日志来查看信息就显得格外重要. 下面是在TP3.2.3框架中,写入日志的方法: public functio ...
- Python中协程Event()函数
python线程的事件用于主线程控制其他线程的执行,事件主要提供了三个方法wait.clear.set 事件处理的机制:全局定义了一个“Flag”,如果“Flag”值为 False,那么当程序执行 e ...
- 170726、常用 Git 命令清单
,下面是我整理的常用 Git 命令清单.几个专用名词的译名如下: Workspace:工作区 Index / Stage:暂存区 Repository:仓库区(或本地仓库) Remote:远程仓库 一 ...
- Docker容器之Nginx
一,pull一个Nginx镜像 docker pull nginx 二,Nginx镜像文件说明 配置文件 /etc/nginx/nginx.conf 网站根目录 /usr/share/nginx/ht ...
- Log4j最简入门及实例
Log4j真的很简单,简单到令人发指的地步.不是要记录日志吗?那就给你一个Log,然后你用Log来写东西就行了,先来一个完整类示例: package test; import org.apache.c ...
- 括号匹配问题(区间dp)
简单的检查括号是否配对正确使用的是栈模拟,这个不必再说,现在将这个问题改变一下:如果给出一个括号序列,问需要把他补全成合法最少需要多少步? 这是一个区间dp问题,我们可以利用区间dp来解决,直接看代码 ...
- Information:java: Errors occurred while compiling module 'spring'
IntellJ Idea遇到Errors occurred while compiling module的解决方法 - sully2008的专栏 - CSDN博客 https://blog.csdn. ...