记录编号 357558 评测结果 AAAAAAAAAA
题目名称 [HAOI 2011]问题B 最终得分 100
用户昵称 GravatarHzoi_Go灬Fire 是否通过 通过
代码语言 C++ 运行时间 5.090 s
提交时间 2016-12-11 15:16:56 内存使用 2.02 MiB
显示代码纯文本
  1. #include<cstdio>
  2. #include<cstring>
  3. #include<algorithm>
  4. using namespace std;
  5. const int maxn=55000;
  6. #define LL long long
  7. #define Begin freopen("b.in","r",stdin);freopen("b.out","w",stdout);
  8. #define End fclose(stdin);fclose(stdout);
  9. bool vis[maxn];
  10. int Prime[maxn],mu[maxn],tot,sum[maxn];
  11. int a[maxn],b[maxn],c[maxn],d[maxn],k[maxn],till;
  12. void Init();
  13. void Get_mu(){
  14. mu[1]=1;
  15. for(int i=2;i<=till;i++){
  16. if(!vis[i]){
  17. Prime[++tot]=i;mu[i]=-1;
  18. }
  19. for(int j=1;j<=tot;j++){
  20. if(i*Prime[j]>till)break;
  21. vis[i*Prime[j]]=true;
  22. if(i%Prime[j]==0)break;
  23. mu[i*Prime[j]]=-mu[i];
  24. }
  25. }
  26. }
  27. LL calc(int n,int m,int z){
  28. n/=z;m/=z;if(n>m)swap(n,m);
  29. LL ans=0;int last;
  30. for(int i=1;i<=n;i=last+1){
  31. last=min(n/(n/i),m/(m/i));
  32. ans+=1ll*(n/i)*(m/i)*(sum[last]-sum[i-1]);
  33. }
  34. return ans;
  35. }
  36. int main(){
  37. Begin;
  38. int T=0;scanf("%d",&T);
  39. for(int i=1;i<=T;i++){
  40. scanf("%d%d%d%d%d",&a[i],&b[i],&c[i],&d[i],&k[i]);
  41. till=max(till,max(b[i],d[i]));
  42. }
  43. Get_mu();
  44. for(int i=1;i<=till;i++)sum[i]=sum[i-1]+mu[i];
  45. 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]));
  46. getchar();getchar();
  47. return 0;
  48. }