HYSBZ - 2005 莫比乌斯反演

链接

对于gcd（i，j）的位置来说，对答案的贡献是2*（gcd（i，j）-1）+1，所以答案ans

ans=Σ(1<=i<=n)(1<=j<=m)2*(gcd(i,j)-1)+1

ans=2*Σ(1<=i<=n)(1<=j<=m)gcd(i,j)-n*m

前者可以通过莫比乌斯反演来计算，便很容易得出答案

//#pragma comment(linker, "/stack:200000000")
//#pragma GCC optimize("Ofast,no-stack-protector")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
//#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define pi acos(-1.0)
#define ll long long
#define mod 1000000007
#define C 0.5772156649
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
#define pil pair<int,ll>
#define pii pair<int,int>
#define ull unsigned long long
#define base 1000000000000000000
#define fio ios::sync_with_stdio(false);cin.tie(0)

using namespace std;

const double g=10.0,eps=1e-;
const int N=+,maxn=+,inf=0x3f3f3f3f;

int mu[N],prime[N],sum[N];
bool mark[N];
int num[N];
void init()
{
    mu[]=;
    int cnt=;
    for(int i=;i<N;i++)
    {
        if(!mark[i])prime[++cnt]=i,mu[i]=-,num[i]=;
        for(int j=;j<=cnt;j++)
        {
            int t=i*prime[j];
            if(t>N)break;
            mark[t]=;
            num[t]=num[i]+;
            if(i%prime[j]==){mu[t]=;break;}
            else mu[t]=-mu[i];
        }
    }
    for(int i=;i<N;i++)sum[i]=sum[i-]+mu[i];
}
int main()
{
    init();
    int n,m;
    scanf("%d%d",&n,&m);
    ll ans=;
    for(int j=;j<=max(n,m);j++)
    {
        ll te=;
        int ten=n/j,tem=m/j;
        for(int i=,last;i<=min(ten,tem);i=last+)
        {
            last=min(ten/(ten/i),tem/(tem/i));
            te+=(ll)(sum[last]-sum[i-])*(ten/i)*(tem/i);
        }
//        printf("%lld\n",te);
        ans+=te*j;
    }
    printf("%lld\n",*ans-(ll)n*m);
    return ;
}
/********************

********************/