记录编号 |
357558 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[HAOI 2011]问题B |
最终得分 |
100 |
用户昵称 |
Hzoi_Go灬Fire |
是否通过 |
通过 |
代码语言 |
C++ |
运行时间 |
5.090 s |
提交时间 |
2016-12-11 15:16:56 |
内存使用 |
2.02 MiB |
显示代码纯文本
- #include<cstdio>
- #include<cstring>
- #include<algorithm>
- using namespace std;
- const int maxn=55000;
- #define LL long long
- #define Begin freopen("b.in","r",stdin);freopen("b.out","w",stdout);
- #define End fclose(stdin);fclose(stdout);
- bool vis[maxn];
- int Prime[maxn],mu[maxn],tot,sum[maxn];
- int a[maxn],b[maxn],c[maxn],d[maxn],k[maxn],till;
- void Init();
- void Get_mu(){
- mu[1]=1;
- for(int i=2;i<=till;i++){
- if(!vis[i]){
- Prime[++tot]=i;mu[i]=-1;
- }
- for(int j=1;j<=tot;j++){
- if(i*Prime[j]>till)break;
- vis[i*Prime[j]]=true;
- if(i%Prime[j]==0)break;
- mu[i*Prime[j]]=-mu[i];
- }
- }
- }
- LL calc(int n,int m,int z){
- n/=z;m/=z;if(n>m)swap(n,m);
- LL ans=0;int last;
- for(int i=1;i<=n;i=last+1){
- last=min(n/(n/i),m/(m/i));
- ans+=1ll*(n/i)*(m/i)*(sum[last]-sum[i-1]);
- }
- return ans;
- }
- int main(){
- Begin;
- int T=0;scanf("%d",&T);
- for(int i=1;i<=T;i++){
- scanf("%d%d%d%d%d",&a[i],&b[i],&c[i],&d[i],&k[i]);
- till=max(till,max(b[i],d[i]));
- }
- Get_mu();
- for(int i=1;i<=till;i++)sum[i]=sum[i-1]+mu[i];
- for(int i=1;i<=T;i++)printf("%lld\n",calc(b[i],d[i],k[i])-calc(a[i]-1,d[i],k[i])-calc(b[i],c[i]-1,k[i])+calc(a[i]-1,c[i]-1,k[i]));
- getchar();getchar();
- return 0;
- }