记录编号 |
285942 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[HAOI 2011]问题B |
最终得分 |
100 |
用户昵称 |
TenderRun |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
9.236 s |
提交时间 |
2016-07-25 18:00:41 |
内存使用 |
0.84 MiB |
显示代码纯文本
#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxn=50010;
int prime[maxn],cnt;
int mu[maxn],sum[maxn];
bool check[maxn];
void Prepare(){
mu[1]=1;
for(int i=2;i<=50000;i++){
if(!check[i]){
prime[++cnt]=i;
mu[i]=-1;
}
for(int j=1;j<=cnt;j++){
if(prime[j]*i>50000)break;
check[prime[j]*i]=true;
if(i%prime[j]==0){
mu[prime[j]*i]=0;
break;
}
mu[prime[j]*i]=mu[i]*-1;
}
}
for(int i=1;i<=50000;i++)
sum[i]=sum[i-1]+mu[i];
}
int T,k;
int a,b,c,d;
int C(int n,int m){
int ret=0,p;
if(n>m)swap(n,m);
for(int i=1;i*k<=n;i=p+1){
p=min(n/(n/i),m/(m/i));
ret+=(sum[p]-sum[i-1])*(n/(i*k))*(m/(i*k));
}
return ret;
}
int main(){
freopen("b.in","r",stdin);
freopen("b.out","w",stdout);
Prepare();
scanf("%d",&T);
while(T--){
scanf("%d%d%d%d%d",&a,&b,&c,&d,&k);
printf("%d\n",C(b,d)-C(b,c-1)-C(a-1,d)+C(a-1,c-1));
}
return 0;
}