正题

---

题目大意

求$$\sum _{i=a}^b\sum _{j=c}^d(gcd(i,j)==k)$$

---

解题思路

定义
$$f(i)=\sum _{i=1}^n\sum _{j=1}^m(gcd(i,j)==i)$$
然后计算f利用容斥计算答案

之后我们考虑如何计算
$$F(i)=\sum _d^{d|i}f(i)$$
显然可得出
$$F(i)=\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{i}\rfloor$$
然后莫比乌斯反演一下
$$f(i)=\sum _{d|i}\mu (\frac{d}{k})+F(d)$$
$$f(i)=\sum _{d|i}\mu (\frac{d}{k})+\lfloor \frac{n}{d}\rfloor \lfloor \frac{m}{d}\rfloor$$
时间复杂度降低到$O(n)$
之后我们可以发现$\lfloor \frac{n}{d}\rfloor$只有$2\sqrt{n}$种取值，那么$\lfloor \frac{n}{i}\rfloor \lfloor \frac{m}{i}\rfloor$最多就有$2(\sqrt{n}+\sqrt{m})$种取值，我们可以直接计算这个范围内莫比乌斯函数的前缀和然后直接O(min{n,m})O(\sqrt {min\{n,m\}})O(min{n,m}​)计算答案

---

code

// luogu-judger-enable-o2
#include<cstdio>
#include<algorithm>
#define N 50000
#define ll long long
using namespace std;
ll miu[N+10],c,d,a,b,k,t;
bool v[N+10];
void sumul()
{for(ll i=1;i<=N;i++) miu[i]=1,v[i]=0;for(ll i=2;i<=N;i++){if(v[i]) continue;miu[i]=-1;for(ll j=2*i;j<=N;j+=i){v[j]=1;if((j/i)%i==0) miu[j]=0;else miu[j]*=-1;}}for(ll i=1;i<=N;i++)miu[i]+=miu[i-1];
}
ll ask(ll n,ll m)
{if(n>m) swap(n,m);ll last,re=0;n/=k;m/=k;for(ll i=1;i<=n;i=last+1){last=min(n/(n/i),m/(m/i));re+=(n/i)*(m/i)*(miu[last]-miu[i-1]);}return re;
}
int main()
{sumul();scanf("%lld",&t);for(ll i=1;i<=t;i++){scanf("%lld%lld%lld%lld%lld",&a,&b,&c,&d,&k);printf("%lld\n",ask(b,d)+ask(a-1,c-1)-ask(a-1,d)-ask(b,c-1));}
}