BZOJ1101 Mobius函数+分块

http://www.lydsy.com/JudgeOnline/problem.php?id=1101

题意很简单

意思就是给你下a,b,d;

问x<=a&&y<=b&&gcd(x,y)==k

典型的mobius函数的模板题

但是朴素的求法会超时

所以可以利用一种除数分块的思想

在某种距离之内是不会改变的

AC代码:

#include
using namespace std;
typedef long long ll;
const int N=50000;
int miu[N+10],v[N+10],P[N+10];
inline void read(int &x){char ch;bool ok;
for(ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if(ch=='-') ok=1;
for(x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());if(ok) x=-x;}
int main(){
    for(int i=1;i<=N;i++) miu[i]=1,v[i]=0;
    for(int i=2;i<=N;i++){
        if(v[i]) continue;
        miu[i] = -1;
        for(int j=2*i;j<=N;j+=i){
            v[j]=1;
            if((j/i)%i==0) miu[j]=0;
            else miu[j]*=-1;
        }
    }
    for(int i=1;i<=N;++i)
        P[i]=P[i-1]+miu[i];
    int T;
    read(T);
    while(T--){
        int a,b,k;
        read(a);
        read(b);
        read(k);
        a=a/k;
        b=b/k;
        int ans=0;
        int c;
        for(int i=1;i<=min(a,b);i=c+1){
            c=min(a/(a/i),b/(b/i));
            ans+=(P[c]-P[i-1])*(a/i)*(b/i);
        }
        printf("%d\n",ans);
    }
}

你可能感兴趣的:(数论)