BZOJ 2820: YY的GCD | 数论

题目:

　　

---

题解:

http://hzwer.com/6142.html

#include<cstdio>
#include<algorithm>
#define N 10000005
typedef long long ll;
using namespace std;
int T,n,m,cnt;
bool mark[N];
int pri[N],mu[N];
ll f[N];
void getphi()
{
    mu[]=;
    for (int i=;i<N;i++)
    {
    if (!mark[i]) pri[++cnt]=i,mu[i]=-;
    for (int j=;j<=cnt && pri[j]*i<N;j++)
    {
        mark[i*pri[j]]=;
        if (i%pri[j]==)
        {mu[i*pri[j]]=;break;}
        else mu[i*pri[j]]=-mu[i];
    }
    }
    for (int i=;i<=cnt;i++)
    {
    int p=pri[i];
    for (int j=;j*p<N;j++)
        f[j*p]+=mu[j];
    }
    for (int i=;i<N;i++) f[i]+=f[i-];
}
int main()
{
    getphi();
    scanf("%d",&T);
    while (T--)
    {
    ll ans=;
    scanf("%d%d",&n,&m);
    if (n>m) swap(n,m);
    for (int i=,j;i<=n;i=j+)
    {
        j=min(n/(n/i),m/(m/i));
        ans+=(f[j]-f[i-])*(n/i)*(m/i);
    }
    printf("%lld\n",ans);
    }
    return ;
}