51Nod1192 Gcd表中的质数

题目看这里
又到了推式子的时候了，莫比乌斯反演入门题

∑i=1n∑j=1m[prime(i,j)] ∑ i = 1 n ∑ j = 1 m [ p r i m e ( i , j ) ]

$$\sum_{i=1}^n\sum_{j=1}^m[prime(i,j)]$$
首先搞一个$f(d)$表示有多少对(i,j)的gcd就是d
那么原式$=\sum_{d=1}^n[prime(d)]*f(d)$
设$F(x)=\sum_{x|d}f(d)=[\frac{n}{x}][\frac{m}{x}]$
就得到原式$=\sum_{d=1}^n[prime(d)]*\sum_{i}F(i*d)*\mu(i)$
$=\sum_{i*d<=n}F(i*d)*\mu(i)*[prime(d)]$
$=\sum_{T<=n}F(T)*\sum_{prime(d)=1,d|T}\mu(\frac{T}{d})$
预处理后面那个$g(T)=\sum_{prime(d)=1,d|T}\mu(\frac{T}{d})$就可以分块求答案了

#pragma GCC optimize("O3")
#pragma G++ optimize("O3")
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 5000010
#define LL long long
using namespace std;
bool vis[N];
int mu[N],s[N],w[N>>2],t,n,m,T;
inline LL F(int d){
    return (LL)(n/d)*(m/d);
}
inline void cal(){
    scanf("%d%d",&n,&m);
    if(n>m) swap(n,m); LL S=0;
    for(int i=1,j;i<=n;i=j+1){
        j=min(n/(n/i),m/(m/i));
        S+=F(i)*(s[j]-s[i-1]);
    }
    printf("%lld\n",S);
}
int main(){
    mu[1]=1;
    for(int i=2;i<=5000000;++i){
        if(!vis[i]) mu[w[++t]=i]=-1;
        for(int j=1,k;j<=t && (k=i*w[j])<=5000000;++j){
            vis[k]=1;
            if(i%w[j]==0) { mu[k]=0; break; } else mu[k]=-mu[i];
        }
    }
    for(int i=1;i<=t;++i)
        for(int j=w[i],k=1;j<=5000000;++k,j+=w[i]) s[j]+=mu[k];
    for(int i=1;i<=5000000;++i) s[i]+=s[i-1];
    for(scanf("%d",&T);T--;cal());
}