解题报告：HDU_5663 Hillan and the girl 莫比乌斯反演

题目链接

题意:

a为小于等于n的正整数，b为小于等于m的正整数，求gcd(a,b)为平方数的计数

思路：

令

设 f(x) 表示gcd(a,b)为x的组数，它的答案为

我们要求的是

可以发现若k为平方数，那么g(x*k)前的系数等同于g(x)前的系数

那么我们在O(n)的时间内预处理出对应的系数，对于每个询问在O( sqrt(n)）的时间内得到答案

总复杂度O(T*sqrt(n))，属于O(能过)

代码：

#include<bits/stdc++.h>

const int N = 1e7+10;
using namespace std;

vector<int>pr;
int mu[N];
bool Np[N];
int sum[N];



void init(){
   sum[1] = mu[1] = 1;
   for(int i=2;i<N;i++){
      if(!Np[i]){
         pr.emplace_back(i);
         mu[i] = -1;
      }for(int j=0;j<pr.size();j++){
         int k = i * pr[j];
         if(k>=N)break;
         Np[k] = true;
         if(i%pr[j]==0){
            mu[k] = mu[i/pr[j]];
            break;
         }mu[k] = -mu[i];
      }
      sum[i]=sum[i-1]+mu[i];
   }
}

int main()
{
   init();
   int T;
   scanf("%d",&T);
   while(T--){
      int n,m;
      scanf("%d%d",&n,&m);
      if(n>m)swap(n,m);
      long long ans = 1LL*n*m ;
      for(int i=1,last;i<=n;i=last+1){
         int ni = n/i , mi = m/i;
         last = min(n/ni,m/mi);
         ans += 1LL * ( sum[i-1] - sum[last] ) * ni * mi;
      }printf("%I64d\n",ans);
   }return 0;
}