题目：求∑i=1n∑j=1mlcm(i,j)\sum\limits_{i=1}^n\sum\limits_{j=1}^m lcm(i,j)
#include
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#def...
分类：
其他好文 时间：
2016-05-07 06:56:07
阅读次数：
329
分析：对于给出的n个询问，每次求有多少个数对(x,y)，满足a≤x≤b，c≤y≤d，且gcd(x,y) = k，gcd(x,y)函数为x和y的最大公约数。 然后对于求这样单个的gcd(x,y)=k的，我们通常采用莫比乌斯反演 但是，时间复杂度是O(n*(n/k))的，当复杂度很坏的时候，当k=1时， ...
分类：
其他好文 时间：
2016-04-27 22:45:52
阅读次数：
334
分析：http://www.cnblogs.com/huhuuu/archive/2011/11/25/2263803.html 注：从这个题收获了两点 1,第一象限（x,y)到（0，0）的线段上整点的个数是gcd(x,y) 2,新学了一发求gcd(x,y)=k有多少对的姿势，已知0<x<=n,0< ...
分类：
其他好文 时间：
2016-04-27 22:31:19
阅读次数：
253
分析：筛素数，然后枚举，莫比乌斯反演，然后关键就是分块加速（分块加速在上一篇文章） #include<cstdio> #include<cstring> #include<queue> #include<cstdlib> #include<algorithm> #include<vector> #i ...
分类：
其他好文 时间：
2016-04-27 20:37:22
阅读次数：
300
分析： http://blog.csdn.net/acdreamers/article/details/12871643 分析参见这一篇 http://wenku.baidu.com/view/fbe263d384254b35eefd34eb.html 分块看这一篇 #include<cstdio> ...
分类：
其他好文 时间：
2016-04-27 20:28:26
阅读次数：
112
今天才知道莫比乌斯反演还可以这样：$$F(n)=\sum_{n|d}f(d) \Rightarrow f(n)=\sum_{n|d}\mu(\frac{d}{n})F(d)$$我好弱，，，对于$$F(i)=\left \lfloor \frac{n}{i} \right \rfloor\left \ ...
分类：
其他好文 时间：
2016-04-23 16:48:06
阅读次数：
265
分析：简单的莫比乌斯反演 f[i]为k=i时的答案数 然后就很简单了 #include<iostream> #include<algorithm> #include<set> #include<vector> #include<queue> #include<cstdlib> #include<cs ...
分类：
其他好文 时间：
2016-04-22 18:09:31
阅读次数：
124
分析：令f(x)为1到n的gcd(i,j)==x的个数 F(x)为1到n的x|gcd(i,j)的对数 显然F(n)=∑n|df(d) 然后由莫比乌斯反演可得f(n)=∑n|d μ(d/n)*F(d) 由题目显然可得，令cnt=n/x,当cnt<3时,F(x)为0,cnt>=3,F(x)=cnt*(c ...
分类：
其他好文 时间：
2016-04-19 19:05:53
阅读次数：
262
//容斥原理，c[i]表示i当前要算的次数，复杂度和第二层循环相关 O(nlogn~n^2) LL in_exclusion(int n,int *c) { for(int i=0;i<=n;i++) c[i]=1; //不一定是这样初始化，要算到的才初始化为1 LL ans=0; for(int ...
分类：
其他好文 时间：
2016-04-19 00:02:48
阅读次数：
407
