bzoj3309 DZY Loves Math 数论

       利用公式Σ(d|n)φ(d)=n，Σ(d|n)μ(d)=[n=1]，将所求的化为Σ(k)[a/k][b/k]Σ(d|k)f(d)μ(k/d)，令后面的为g(k)，则所求的转化为Σ(k)[a/k][b/k]g(k)，就可以O(N^0.5)查询了。

       关键是求g(k)，设k=a1^p1*a2^p2*...*an^pn，d=a1^q1*a2^q2*...*an^qn。由于μ函数的性质，如果存在pi-qi>=2，那么μ(k/d)=0，不予考虑。那么我们可以假定所有0<=pi-qi<=1。分下列两种情况讨论。

       1.如果存在pi≠pj，那么g(k)=0。首先为了便于叙述，如果pi=qi，我们称d取了pi。假设k的分解式中最大的pi为pmax，满足pi=pmax的共有t个。对于d分两种情况讨论。如果qi中最大的qmax=pmax，那么f(d)=pmax。换句话说，f(d)=pmax当且仅当d取了t中的某一些。如果对于那t个中的一种取法，剩下的必然是取n-t中的某一些。而取n-t中的某一些的所有取法中，取奇数个和偶数个的方案是一样的（可以用归纳法证明，注意0也是偶数），所以如果d取了pmax，对于t中的任何一种取法，相加以后的和为0。如果d没有取pmax，则t=0，也是同样的讨论过程。

       2.如果所有pi=pj，那么g(k)=（-1）^n。同上所有的f(d)都是可以消去的，但是其中一组对应的一个值f(a1^(p1-1)*a2^(p1-1)*...*an^(pn-1))比其他所有的都小一，所以它对应的μ值也相差了一，这也就是g(k)的值。

       那么g就可以用线性筛的方法处理了。只需要处理a[]表示最小质因数的次数，b[]表示最小质因数的a[]次方，就可以利用以前递推的值求g了。具体看代码：

AC代码如下：

#include<iostream>
#include<cstdio>
#define N 10000005
#define ll long long 
using namespace std;
 
int n,m,cnt,a[N],b[N],c[N],g[N];
void pfs(){
    int i,j,k; g[1]=0;
    for (i=2; i<=10000000; i++){
        if (!a[i]){
            c[++cnt]=i; a[i]=g[i]=1; b[i]=i;
        }
        for (j=1; j<=cnt; j++){
            k=i*c[j]; if (k>10000000) break;
            if (i%c[j]){
                a[k]=1; b[k]=c[j];
                g[k]=(a[i]==1)?-g[i]:0;
            } else{
                a[k]=a[i]+1; b[k]=b[i]*c[j];
                int tmp=i/b[i];
                if (tmp==1) g[k]=1; else g[k]=(a[tmp]==a[k])?-g[tmp]:0;
                break;
            }
        }
    }
    for (i=2; i<=10000000; i++) g[i]+=g[i-1];
}
int main(){
    int cas; scanf("%d",&cas); pfs();
    while (cas--){
        scanf("%d%d",&m,&n); int i,j; ll ans=0;
        for (i=1; i<=m && i<=n; i=j+1){
            j=min((m/(m/i)),n/(n/i));
            ans+=(ll)(g[j]-g[i-1])*(m/i)*(n/i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}