题解【bzoj3529 [SDOI2014]数表】
Description
\(T\) 组询问,定义 \(F(n)=\sum\limits_{d|n}d\)。每次给出 \(n,m,a\) 求
\]
\(T \leq 20000;n,m,a\leq 10^5\)
Solution
首先 \(F\) 可以直接暴力地 \(O(n \log n)\) 筛出来。
考虑 \(a\) 的限制不是很好处理,假设没有这个 \(a\) 的限制,则所求为
\]
令 \(G(i)=\sum\limits_{x=1}^{n}\sum\limits_{y=1}^{m}[\gcd(x,y)=i]\)。这个东西是什么呢?在 这里 有它的推导过程。根据里面的过程,可以得到 \(G(i) = \sum\limits_{d=1}^{\lfloor\frac{n}{i}\rfloor}\mu(d)\lfloor\frac{n}{id}\rfloor\lfloor\frac{m}{id}\rfloor\)(默认 \(n \leq m\))
令下面的过程中\(t = id\),则所求的是
\]
观察后面的式子,正好是一个狄利克雷卷积的形式。这种样子的都可以类似于那种 \(O(n \log n)\) 地质数筛法在调和级数内求出来,再结合分块就可以做完这个没有 \(a\) 的题。
现在有了 \(a\) 的限制之后,离线。把询问按照 \(a\) 从小到大排序,然后按照 \(F(i)\) 从小到大加入。每当有一个新的 \(a\) ,就可以移动指针,将一些 \(F\) 用处理 \(\sum\limits_{i | t}F(i)\mu(\frac{t}{i})\) 的方式加入到这个里面。然后加入完之后用分块计算就行。
现在需要维护单点操作,查询前缀和,树状数组是不错的选择。
由于取模是 \(2^{32} - 1\) ,可以直接 int 自然溢出最后和 \(2147483647\) 取一个 & 就行了。
时间复杂度:\(O(n + n \log n + n \log ^ 2 (n) + T \log (n)\sqrt n)\)
Code
#include <bits/stdc++.h>
using namespace std;
const int N = 100000;
int T, cnt, flag[N + 5], p[N + 5], F[N + 5], mu[N + 5], ans[N + 5];
inline void prework() {
flag[1] = mu[1] = 1;
for(int i = 2; i <= N; i++) {
if(!flag[i]) { p[++cnt] = i, mu[i] = -1; }
for(int j = 1; j <= cnt && i * p[j] <= N; j++) {
flag[i * p[j]] = 1; if(i % p[j] == 0) {
mu[i * p[j]] = 0; break ;
} mu[i * p[j]] = mu[i] * -1;
}
}
for(int i = 1; i <= N; i++)
for(int j = i; j <= N; j += i) F[j] += i;
}
int c[N + 5];
inline int lb(int x) { return x & (-x); }
inline void add(int x, int d) {
for(int i = x; i <= N; i += lb(i)) c[i] += d;
}
inline int sum(int x) {
int ret = 0;
for(int i = x; i; i -= lb(i))
ret += c[i]; return ret;
}
inline int calc(int n, int m) {
int ret = 0;
for(int l = 1, r; l <= min(n, m); l = r + 1) {
r = min(n / (n / l), m / (m / l));
ret += (n / l) * (m / l) * (sum(r) - sum(l - 1));
} return ret;
}
struct Query {
int n, m, a, id;
inline bool operator < (const Query &x) const {
return a < x.a;
}
}Q[N + 5];
struct node {
int id, d;
inline bool operator < (const node &x) const {
return d < x.d;
}
}A[N + 5];
int main() {
prework(); scanf("%d", &T);
for(int i = 1; i <= T; i++) scanf("%d %d %d", &Q[i].n, &Q[i].m, &Q[i].a), Q[i].id = i;
for(int i = 1; i <= N; i++) A[i].d = F[i], A[i].id = i;
sort(Q + 1, Q + T + 1); sort(A + 1, A + N + 1); int pos = 0;
for(int i = 1; i <= T; i++) {
while(pos < N && A[pos + 1].d <= Q[i].a) {
++pos;
for(int j = 1; A[pos].id * j <= N; j++)
add(j * A[pos].id, A[pos].d * mu[j]);
} ans[Q[i].id] = calc(Q[i].n, Q[i].m);
}
for(int i = 1; i <= T; i++) printf("%d\n", ans[i] & 2147483647);
return 0;
}
题解【bzoj3529 [SDOI2014]数表】的更多相关文章
- BZOJ3529 [Sdoi2014]数表 【莫比乌斯反演】
3529: [Sdoi2014]数表 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 2151 Solved: 1080 [Submit][Status ...
- [bzoj3529][Sdoi2014]数表_树状数组_莫比乌斯反演
数表 bzoj-3529 Sdoi-2014 题目大意:n*m的数表,第i行第j列的数是同时整除i和j的所有自然数之和.给定a,求数表中所有不超过a的和. 注释:$1\le n,m \le 10^5$ ...
- BZOJ3529 [Sdoi2014]数表【莫比乌斯反演】
Description 有一张 n×m 的数表,其第 i 行第 j 列(1 <= i <= n, 1 <= j <= m)的数值为 能同时整除 i 和 j 的所有自然数之和.给 ...
- bzoj千题计划205:bzoj3529: [Sdoi2014]数表
http://www.lydsy.com/JudgeOnline/problem.php?id=3529 有一张n*m的数表,其第i行第j列(1 < =i < =n,1 < =j & ...
- 莫比乌斯反演套路二--(n/d)(m/d)给提出来--BZOJ3529: [Sdoi2014]数表
一个数表上第i行第j列表示能同时整除i和j的自然数,Q<=2e4个询问,每次问表上1<=x<=n,1<=y<=m区域内所有<=a的数之和.n,m<=1e5,a ...
- BZOJ3529: [Sdoi2014]数表(莫比乌斯反演,离线)
Description 有一张 n×m 的数表,其第 i 行第 j 列(1 <= i <= n, 1 <= j <= m)的数值为 能同时整除 i 和 j 的所有自然数之和.给 ...
- BZOJ3529: [Sdoi2014]数表
题目:http://www.lydsy.com/JudgeOnline/problem.php?id=3529 挺恶心的数论TAT... 设f[i]是i的约数和,这个可以nln(n)扫出来. ans= ...
- bzoj3529: [Sdoi2014]数表 莫比乌斯反演
题意:求\(\sum_{i=1}^n\sum_{j=1}^nf(gcd(i,j))(gcd(i,j)<=a),f(x)是x的因子和函数\) 先考虑没有限制的情况,考虑枚举gcd为x,那么有\(\ ...
- BZOJ3529: [Sdoi2014]数表(莫比乌斯反演 树状数组)
题意 题目链接 Sol 首先不考虑\(a\)的限制 我们要求的是 \[\sum_{i = 1}^n \sum_{j = 1}^m \sigma(gcd(i, j))\] 用常规的套路可以化到这个形式 ...
随机推荐
- cmake-index-3.11.4机翻
index next | CMake » git-stage git-master latest release 3.13 3.12 3.11.4 3.10 3.9 3.8 3.7 3.6 3.5 3 ...
- Windows和Linux系统如何退出python命令行
python命令行是新手学习python过程中必须要学的一个工具,下面我们来看一下怎么退出python命令行. 第一种方式: 使用python提供的exit()函数,linux平台和windows平台 ...
- Twitter推广消息可使品牌线下销售额增长三成
新浪科技讯 北京时间8月9日上午消息,Twitter周四宣布,该公司的推广消息(Promoted Tweet)可以让品牌的线下销售增长29%. 此外,Twitter当天还推出了一个新项目,让品牌可以追 ...
- CSS Grid布局指南
简介 CSS Grid布局 (又名"网格"),是一个基于二维网格布局的系统,主要目的是改变我们基于网格设计的用户接口方式.如我们所知,CSS 总是用于网页的样式设置,但它并没有起到 ...
- jquery选择器 直观实验
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/ ...
- 查看dll依赖项
win7 系统: 开始-->所有程序->vs2012文件夹->vs tools->对应的命令提示符 输入命令: dumpbin /dependents 你的文件(可以是exe, ...
- OpenNF tutorial复现
这篇博客记录了自己实现OpenNF官网上tutorial的过程和遇见的问题,如果有不对的地方还请批评指正! tutorial链接 实验内容 这个实验展示了如何迅速且安全地把一个TCP流从一个NF实例迁 ...
- 软工 · BETA 版冲刺前准备(团队)
软工 · BETA 版冲刺前准备(团队) 过去存在的问题 组员之间缺乏沟通,前后端缺乏沟通协作 组员积极性不高 基础知识不够扎实 手动整合代码效率过低 我们已经做了哪些调整/改进 通过会议加强组员之间 ...
- springboot+vue+element:echarts开发遇见问题---后端sql(三)
<select id="getSumRequestRankingCount" parameterType="java.lang.String" resul ...
- DNS缓存服务器的配置步骤
yum安装bind 编辑主配置文件/etc/named.conf 修改全局配置文件段 listen-on port 53 {172.16.19.45;}; //allow-query ...