【BZOJ】【1101】【POI2007】Zap

莫比乌斯反演

　　PoPoQQQ的讲义例一的一半……好吧这题是那题的基础部分= =很水……

WA了一次：因为没强制类型转换LL

 /**************************************************************
     Problem: 1101
     User: Tunix
     Language: C++
     Result: Accepted
     Time:6896 ms
     Memory:1908 kb
 ****************************************************************/

 //BZOJ 1101
 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #define rep(i,n) for(int i=0;i<n;++i)
 #define F(i,j,n) for(int i=j;i<=n;++i)
 #define D(i,j,n) for(int i=j;i>=n;--i)
 using namespace std;

 int getint(){
     int v=,sign=; char ch=getchar();
     while(ch<''||ch>'') {if (ch=='-') sign=-; ch=getchar();}
     while(ch>=''&&ch<='') {v=v*+ch-''; ch=getchar();}
     return v*=sign;
 }
 /*******************tamplate********************/
 const int N=;
 typedef long long LL;
 int prime[N],mu[N],sum[N];
 bool check[N];
 void getmu(){
     mu[]=;
     int tot=;
     for(int i=;i<N;++i){
         if (!check[i]){
             prime[tot++]=i;
             mu[i]=-;
         }
         rep(j,tot){
             if (i*prime[j]>N) break;
             check[i*prime[j]]=;
             if (i%prime[j])
                 mu[i*prime[j]]=-mu[i];
             else{
                 mu[i*prime[j]]=;
                 break;
             }
         }
     }
     F(i,,N-) sum[i]=sum[i-]+mu[i];
 }
 LL calc(int m,int n,int k){
     int i,last;
     LL re=;
     n/=k; m/=k;
     for(i=;i<=m && i<=n;i=last+){
         last=min(n/(n/i),m/(m/i));
         re+=(LL)(sum[last]-sum[i-])*(m/i)*(n/i);
     }
     return re;
 }
 int main(){
     int n=getint(),a,b,d;
     getmu();
     F(i,,n){
         a=getint(); b=getint(); d=getint();
         printf("%lld\n",calc(a,b,d));
     }
     return ;
 }

1101: [POI2007]Zap

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1559  Solved: 534
[Submit][Status][Discuss]

Description

FGD正在破解一段密码，他需要回答很多类似的问题：对于给定的整数a,b和d，有多少正整数对x,y，满足x<=a，y<=b，并且gcd(x,y)=d。作为FGD的同学，FGD希望得到你的帮助。

Input

第一行包含一个正整数n，表示一共有n组询问。（1<=n<= 50000）接下来n行，每行表示一个询问，每行三个正整数，分别为a,b,d。（1<=d<=a,b<=50000）

Output

对于每组询问，输出到输出文件zap.out一个正整数，表示满足条件的整数对数。

Sample Input

2
4 5 2
6 4 3

Sample Output

3
2

HINT

对于第一组询问，满足条件的整数对有(2,2)，（2,4），（4，2）。对于第二组询问，满足条件的整数对有（6,3），（3,3）。

Source

[Submit][Status][Discuss]