一.题目 #124. 除数函数求和 二.分析 比较好的一题,首先我们要对题目和样例进行分析,明白题目的意思. 由于对于每一个$d$,它所能整除的数其实都是定的,且数量是$ \lfloor \frac{n}{d} \rfloor $ 最终推导出这个公式 $$  ans =   \sum_{d=1}^{n} \lfloor \frac{n}{d} \rfloor d^{k}$$ 对于$n <= 10^{7}$其实复杂度是可以接受的.但是对于求$d^{k}$这个复杂度如果直接用快速幂预处理肯定会T.…

题目大意 给出正整数 n 和 k 计算 \(G(n, k)=k\ \bmod\ 1 + k\ \bmod\ 2 + k\ \bmod\ 3 + \cdots + k\ \bmod\ n\) 的值 其中 \(k\ \bmod\ i\) 表示 k 除以 i 的余数. 解析 整除分块的一个典型例子. 整除分块解决的是形如 \[ \sum^n_{i=1} ~ \lfloor\frac{n}{i}\rfloor \] 的问题,其复杂度为\(O(\sqrt{n})\). 实际上是规律性的一类问题,打表可以发…

\(\sum_{i=1}^n\;k\;mod\;i\) Solution \(\sum_{i=1}^n\;k\;mod\;i\\=\sum_{i=1}^n(k-i\lfloor{\frac{k}{i}}\rfloor)\\=k\times n-\sum_{i=1}^ni\lfloor{\frac{k}{i}}\rfloor\) 至于后面那项,整除分块即可 #include <bits/stdc++.h> using namespace std; #define int long long si…

题目描述 $\sigma_k(n) = \sum_{d | n} d ^ k$​ 求 $\sum_{i=1}^n\sigma_k(i)$ 的值对 109 取模的结果. 输入格式 第一行两个正整数 n,k . 输出格式 第一行输出答案. 样例 输入样例 5 2 输出样例 63 数据范围与提示 对于 100% 的数据,1≤n,k≤107​7​​ . Solution: 本题ZYYS... 直接枚举显然不行,我们考虑改为求$n$的某一因子$d$在整个函数中的贡献是多少. 套上数论分块的思想,一个因子$…

链接:https://loj.ac/problem/124 就是筛一下积性函数. #include<bits/stdc++.h> #define ll long long #define maxn 10000000 #define ha 1000000007 using namespace std; int zs[maxn/10],t=0; int low[maxn+5]; int f[maxn+5]; int mik[maxn+5]; bool v[maxn+5]; int n,k; inl…

上模板题例题: [CQOI2007]余数求和 洛谷 BZOJ 题目大意:求 $\sum^n_{i=1}k\ mod\ i$ 的值. 等等……这题就学了三天C++的都会吧? $1\leq n,k\leq 10^9$.(一口老血喷到屏幕上) $O(n)$ 行不通了,考虑别的做法. 我们来看一下 $\lfloor\frac{x}{i}\rfloor$ 的值. $x=9$:(不包括0,只有4种取值?) i 1 2 3 4 5 6 7 8 9 10 x/i 9 4 3 2 1 1 1 1 1 0 $x=1…

题目链接 https://www.luogu.org/problemnew/show/P2261 分析 显然\(k\) \(mod\) \(i=k-\lfloor {k/i}\rfloor\) \(\times\) \(i\),于是我们只需要求\(N * k-\sum_{i=1}^N {\lfloor {k/i}\rfloor\times i}\) 这里就需要数论分块,也称作整除分块的知识 结论: \(\forall{i} \in [x,\lfloor {k/{\lfloor {k/x}\rfl…

一.题面 P2261 [CQOI2007]余数求和 二.分析 参考文章:click here 对于整除分块,最重要的是弄清楚怎样求的分得的每个块的范围. 假设$ n = 10 ,k = 5 $ $$   i : 1 \  2 \ 3 \ 4 \ 5 \ 6 \ 7 \ 8 \ 9 \ 10  \\  \lfloor \frac{k}{i} \rfloor :  5 \ 2 \ 1 \ 1 \ 1 \ 0 \ 0 \ 0 \ 0 \ 0   $$ 我们推导出假设$ L = i $,那么,对应的 $…

洛谷传送门 bzoj传送门 这道题要用到学习莫比乌斯反演时掌握的整除分块算法,也就是对于一个数n" role="presentation" style="position: relative;">nn,n" role="presentation" style="position: relative;">nn除以1" role="presentation" style…

参考:题解 令f(i)=k%i,[p]表示不大于p的最大整数f(i)=k%i=k-[k/i]*i令q=[k/i]f(i)=k-qi如果k/(i+1)=k/i=qf(i+1)=k-q(i+1)=k-qi-q=f(i)-q于是,对于区间[l,r],使其之内任意两个整数i,j,都满足k/i=k/j,则f(l)到f(r)是一个递减的等差数列,公差为[k/i].现在就是要把1到n分成这样的一些区间,设某个区间的商(公差)为p设区间内某数为x,则现在要做的是解方程[k/x]=p显然px<=k,因此x<=k…