Description
對於給出的n個詢問,每次求有多少個數對(x,y),滿足a≤x≤b,c≤y≤d,且gcd(x,y) = k,gcd(x,y)函數爲x和y的最大公約數。
Input
第一行一個整數n,接下來n行每行五個整數,分別表示a、b、c、d、k
Output
共n行,每行一個整數表示滿足要求的數對(x,y)的個數
Sample Input
2
2 5 1 5 1
1 5 1 5 2
Sample Output
14
3
HINT
100%的數據滿足:1≤n≤50000,1≤a≤b≤50000,1≤c≤d≤50000,1≤k≤50000
Source
思維難度:省選
代碼難度:NOIP+
算法:莫比烏斯反演
// luogu-judger-enable-o2
#include<cstdio>
#include<iostream>
#define ll long long
using namespace std;
const int Maxn=10000005;
ll pri[Maxn/10],u[Maxn],sum[Maxn],k,t,a,b,c,d;
bool vis[Maxn];
inline ll mn(ll x,ll y){
return x<y?x:y;
}
inline ll read(){//快速讀入,卡常必備
ll x=0;char c=getchar();
while(c<'0'||c>'9'){c=getchar();}
while(c<='9'&&c>='0'){x=x*10;x=x+c-'0';c=getchar();}
return x;
}
void prep(){//預處理
u[1]=1;
for(int i=2;i<=Maxn;++i){//篩素數+求莫比烏斯函數
if(!vis[i]){pri[++k]=i;u[i]=-1;}
for(int j=1;j<=k&&pri[j]*i<=Maxn;++j){
vis[pri[j]*i]=1;
if(i%pri[j]==0)break;
u[pri[j]*i]=-u[i];
}
}
for(int i=1;i<=Maxn;++i){
sum[i]=sum[i-1]+u[i];//前綴和
}
}
inline ll work(ll n,ll m,ll d){
ll x,ans=0,lim;n/=d;m/=d;
lim=mn(n,m);
for(int i=1;i<=lim;i=x+1){
x=mn(n/(n/i),m/(m/i));//分塊
ans+=(n/i)*(m/i)*(sum[x]-sum[i-1]);
}
return ans;
}
int main(){
prep();
scanf("%lld",&t);
while(t--){
a=read();b=read();c=read();d=read();k=read();
printf("%lld\n",work(a-1,c-1,k)+work(b,d,k)-work(a-1,d,k)-work(b,c-1,k));
}
return 0;
}