链接:http://acm.hdu.edu.cn/showproblem.php?pid=5663
题意:给定n,m,求数对(i,j)的个数{gcd(i,j)!=x,1<=i<=n,1<=j<=m,x为完全平方数}。
分析:很明显的莫比乌斯反演。首先我们将问题变为求有多少数对是等于完全平方数的,令g[x]={(i,j),gcd(i,j)=x},f[x]=∑g[d]{x|d}。显然f[x]=(n/x)*(m/x)所以g[x]=∑mu[d/x]*(n/d)*(m/d){x|d}。显然我们如果枚举每一个x(x为完全平方数)然后如果对每个x求一次g[x]的话时间复杂度会是√n*√n即O(n)但是这样还不够快因为有一个T。我们继续分析如果我们将其整体考虑即ans=∑ (k)∑(k|d)mu[d/k]*(n/d)*m/d。我们交换两个∑的位置其实就是交换k和d的枚举先后得:ans=∑(d)(n/d)*(m/d)*∑(k|d)mu[d/k]。这样的话我们只要能够预处理是∑(k|d)mu[d/k]的前缀和就能和前面那一坨一起用分块加速了。显然这样都可以在预处理莫比乌斯函数的时候一起处理好。详见代码。O(nloglogn+t*sqrt(n))。
代码:
#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=10000010;
const int MAX=151;
const int mod=100000000;
const int MOD1=100000007;
const int MOD2=100000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000009;
const ll INF=10000000010;
typedef double db;
typedef unsigned long long ull;
int a[N],q[N],g[N],mu[N];
void deal(int n) {
int i,j,k=0;
q[1]=mu[1]=1;
for (i=2;i<=n;i++) {
if (!q[i]) { a[++k]=i;mu[i]=-1; }
for (j=1;j<=k;j++) {
if (i*a[j]>n) break ;
q[i*a[j]]=1;mu[i*a[j]]=-mu[i];
if (i%a[j]==0) {
mu[i*a[j]]=0;break ;
}
}
}
for (i=1;i*i<=n;i++)
for (j=i*i;j<=n;j+=i*i) g[j]+=mu[j/i/i];
for (i=2;i<=n;i++) g[i]+=g[i-1];
}
int main()
{
int i,n,m,t,last;
ll ans;
deal(10000000);
scanf("%d", &t);
while (t--) {
scanf("%d%d", &n, &m);
if (n>m) { n^=m;m^=n;n^=m; }
ans=last=0;
for (i=1;i<=n;i=last+1) {
last=min(n/(n/i),m/(m/i));
ans+=(ll)(n/i)*(m/i)*(g[last]-g[i-1]);
}
printf("%I64d\n", (ll)n*m-ans);
}
return 0;
}