题目：求∑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
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5584 给一个坐标(ex, ey)，问是由哪几个点走过来的。走的规则是x或者y加上他们的最小公倍数lcm(x, y)。 考虑(ex, ey)是由其他点走过来的，不妨设当走到(x,y)时候，gcd(x, y ...
分类：
其他好文 时间：
2016-04-24 15:32:31
阅读次数：
173
这个算法称为欧几里得算法。不会溢出，因为<!--EndF-->gcd函数的递归层数不超过4.785lgN + 1.6723，其中N=max{a,b}。 让gcd递归层数最多的是gcd(Fn,Fn-1)。利用gcd还可以求出两个整数a和b的最小公倍数lcm(a,b)。 这个结论很容易由唯一分解定理得到 ...
分类：
编程语言 时间：
2016-04-21 20:09:34
阅读次数：
222
设$f(x)=\sum_{x|d}p(d)$。 则$ans=\sum_{i=1}^n\sum_{j=1}^n\sum_{k=1}^n\mu(i)\mu(j)\mu(k)f(lcm(i,j))f(lcm(i,k))f(lcm(j,k))$。 转化成图论模型，$i$到$j$有边的条件是$\mu(i)\n ...
分类：
其他好文 时间：
2016-04-21 13:47:02
阅读次数：
137
每次合并两个同余模方程，然后用exgcd解即可。 ll LCM(ll a,ll b) { return a/__gcd(a,b)*b; } void exgcd(ll a,ll b,ll &d,ll &x,ll &y) { if(b==0){ x=1;y=0;d=a; return; } exgcd ...
分类：
其他好文 时间：
2016-04-18 13:24:51
阅读次数：
154
1、在调试“E:\htwang\smart2200v201\ARM嵌入式系统实验教程(二)\开发板出厂编程程序\液晶显示程序\LCM_Disp”的程序时，想使用外部RAM进行仿真调试，在将ADS1.2中的“DebugInExram Settings…->Arm Linker -> Output -> ...
分类：
其他好文 时间：
2016-04-12 14:18:20
阅读次数：
161
题目链接：http://gdutcode.sinaapp.com/problem.php?cid=1031&pid=5 这个题目一看就是一道数论题，应该考虑使用容斥原理，这里对lcm进行容斥。 不过直接上去是T，考虑到序列中同时存在i和ki的话，其实只需要考虑i，所以先对序列中为倍数的对进行处理。 ...
分类：
其他好文 时间：
2016-04-10 21:30:34
阅读次数：
405
题意： 给n,m,和m个数（k1~km）。求1~n中有多少个数不是(k1~km)中任意一数的倍数。 题解： 容斥模板题。a的倍数有n/a个；既是a,也是b的倍数，即为lcm(a,b)的倍数。是a,b,c的倍数，即为lcm(a,b,c)的倍数。 #include #include #include u... ...
分类：
其他好文 时间：
2016-04-09 23:13:12
阅读次数：
180
设答案为ans1,ans2 ans1=a1*gcd，ans2=a2*gcd，a1,a2互质 gcd*a1*b1=lcm，gcd*a2*b2=lcm a1*b1=lcm=(ans1*ans2)/gcd=a1*a2 综上所诉，a1=b2,a2=b1。 也就是说，ans1=gcd*k1，ans2=gcd* ...
分类：
其他好文 时间：
2016-04-08 21:34:59
阅读次数：
242
挺简单的题目，但是有一堆恶心的边界 在刨去恶心的边界之后： 假定我们知道两边的循环节为b1，b2 其中h第一次到达目标的时间为a1，a2 又知道对于答案t t=a1+b1*t1=a2+b2*t2 不妨枚举t1，判断是否存在可行解即可 又因为LCM(b1,b2)就开始循环了 且b1*b2<=b1*mo ...
分类：
其他好文 时间：
2016-04-04 22:37:04
阅读次数：
160
