又没推出来…… 不过通过这道题还是学到好多东西呢，比如积性函数，线筛什么的。 $\sum \limits _{i=1}^{n} \sum \limits _{j=1}^{m} lcm(i,j)$ =$\sum \limits _{d=1}^{min(n,m)} \sum \limits _{i=1} ...
分类：
其他好文 时间：
2019-10-05 18:48:45
阅读次数：
79
最终的柿子：sigema(min)T K(T)*S(T) K(T)=sigema(n/T)i sigema(m/T)j i*j S(T)=sigema(d|T)u(d)*d^d*T 收获：写柿子的时候不能省掉范围。。。要分块找向下整除 ...
分类：
其他好文 时间：
2019-01-02 17:27:24
阅读次数：
138
https://www.luogu.org/problemnew/show/P3768 题面来自洛谷，因为没用markdown所以直接截的图。 剩余的图是我用markdown写完然后截的图。 参考洛谷第一篇题解。 这个式子直观感受就需要莫比乌斯反演，大致的过程参考：BZOJ2693：jzptab 那 ...
分类：
其他好文 时间：
2018-03-20 20:47:31
阅读次数：
193
【题目描述】 为了研究最小公倍数，他画了一张N*M的表格。每个格子里写了一个数字，其中的i行的j列写着数lcm（i，j）。crash表格里所有数的和mod100000009的值。 【输入格式】 第一行输入T，表示数据组数。接下来的T行，每行输入N和M 【输出格式】 对于每个询问，输出表格中所有数的和 ...
分类：
其他好文 时间：
2018-01-24 12:36:03
阅读次数：
125
线性筛+莫比乌斯反演 盗波图 来自candy?大神 反演很重要的一条公式就是[gcd(i,j)==1]= 线性筛怎么推呢？ 我们分4个步骤，1.先推出f[1]，2.推出f[p],p是一个质数，3.由于线性筛筛的是积性函数，那么当gcd(i,p[j])==1的时候，f[i*p[j]]=f[i]*f[p ...
分类：
其他好文 时间：
2017-09-15 13:50:51
阅读次数：
157
Time Limit: 10 Sec Memory Limit: 512 MB Description Input 一个正整数T表示数据组数 接下来T行 每行两个正整数 表示N、M 一个正整数T表示数据组数 接下来T行 每行两个正整数 表示N、M Output T行 每行一个整数 表示第i组数据的结 ...
分类：
其他好文 时间：
2017-04-15 01:30:42
阅读次数：
226
题意：给定N,M，求$\sum\limits_{i = 1}^N {\sum\limits_{j = 1}^M {lcm(i,j)} }$，多组询问 题解： 设$F(x,y) = \sum\limits_{i = 1}^x {\sum\limits_{j = 1,\gcd (i,j) = 1}^y ...
分类：
其他好文 时间：
2017-02-27 01:06:50
阅读次数：
202
题目：http://www.lydsy.com/JudgeOnline/problem.php?id=2693题意：同2154 多组数据题解：按2154再往后转化一下就可以把n，m放到一边儿，然后线性筛右边。 硬要看推导的话可以戳：http://blog.csdn.net/popoqqq/artic...
分类：
其他好文 时间：
2015-01-09 08:03:39
阅读次数：
467
题目大意：同2154 多组数据
后面那坨东西 由于积性函数的约数和仍是积性函数 因此只需要线性筛一下就行
i%prime[j]==0那部分由于多出来的因数都不是无平方因子数因此μ值都为0 增加的只有原先的D/i
#include
#include
#include
#include
#define M 10001000
#define MOD 100000009
using ...
分类：
其他好文 时间：
2014-12-22 14:34:59
阅读次数：
126
