HDU 1695 GCD

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1695

题目大意：给出区间（1,b）,(1,d)，正整数k,在(1,b)中取出一个数x，在(1,d)中取出一个数y,求满足gcd(x,y)==k的所有（x,y）的对数，（x,y）和(y,x)被当作相同的一对。

由欧几里得定理可以知道 gcd(x,y)==k可以写成gcd(x/k,y/k)==1，那么问题转化为在(1,b/k)和（1,d/k）中分别找出一个数x,y使得它们互质，

假设b<d，当x<y时就可以保证对数是不会重复的，在区间（1，b/k）中的对数可以通过求每个数的欧拉函数值加和得到，因此可以先预处理出1e5之内的欧拉函数值

而区间(b/k+1,d/k)中的数n与(1,b/k)中互质的数的个数求解不能用欧拉函数，那么可以用容斥定理求得，这里给出一个公式：

所有不与n互质的数的个数=1个因子倍数的个数-2个因子乘积的倍数的个数+3个因子乘积的倍数的个数...

最后所有与n互质的个数只需要用b/k减去上面式子的和。

而如果w是n的素因子，那么在区间(1,b/k)中w的倍数共有b/k/w个，对(b/k+1,d/k)中的每一个数枚举计算相加即可得到结果

 

代码：

#include<bits/stdc++.h>
using namespace s