记录编号 |
483531 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[HAOI 2011]问题B |
最终得分 |
100 |
用户昵称 |
Shirry |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
10.292 s |
提交时间 |
2018-01-17 17:13:20 |
内存使用 |
0.91 MiB |
显示代码纯文本
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=50010;
int n,a,b,c,d,k,cnt,pri[maxn],sum[maxn],mobius[maxn];
bool vis[maxn];
void GetMobius(){
mobius[1]=1;
for(int i=2;i<=maxn;i++){
if(!vis[i])mobius[i]=-1,pri[++cnt]=i;
for(int j=1;j<=cnt&&i*pri[j]<=maxn;j++){
vis[i*pri[j]]=1;
if(i%pri[j])mobius[i*pri[j]]=-mobius[i];
else{mobius[i*pri[j]]=0;break;}
}
}
for(int i=1;i<=maxn;i++)sum[i]=sum[i-1]+mobius[i];
}
long long work(long long x,long long y){
x/=k,y/=k;
if(x>y)swap(x,y);
if(x==0)return 0;
long long tot=0,tmp;
for(long long i=1;i<=x;i=tmp){
tmp=min(x/(x/i),y/(y/i));
tot+=(x/i)*(y/i)*(sum[tmp]-sum[i-1]);
tmp++;
}
return tot;
}
int main(){
freopen("b.in","r",stdin);
freopen("b.out","w",stdout);
GetMobius();
scanf("%d",&n);
while(n--){
scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
long long ans=work(b,d)+work(a-1,c-1)-work(a-1,d)-work(b,c-1);
printf("%lld\n",ans);
}
return 0;
}