Crash的数字表格

题目描述:

题目分析：

即算出除以d后互质的对数，两个数都乘d的乘积就得到了两个数的乘积，在除以d就是他们的最小公倍数。
如何求f？

总复杂度 O（sqrt(n)*sqrt(n)）=O(n)

题目链接:

Luogu 1829
BZOJ 2154
COGS 1886

Ac 代码:

#include <cstdio>
#include <iostream>
#define ll long long
#define il inline
#define Sum(x,y) (1ll*x*(x+1)/2%mod*(1ll*y*(y+1)/2%mod)%mod)
const int maxm=1e7+100;
const int mod=20101009;
int sum[maxm],mu[maxm],prime[maxm],cnt,n,m;
bool vis[maxm];
il void getmu()
{
    mu[1]=1,vis[1]=1,sum[1]=1;
    for(int i=2;i<=std::min(n,m);i++)
    {
        if(!vis[i])
        {
            prime[++cnt]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=cnt&&prime[j]*i<=std::min(n,m);j++)
        {
            vis[prime[j]*i]=1;
            if(i%prime[j]==0)
            {
                mu[prime[j]*i]=0;
                break;
            }
            mu[prime[j]*i]=-mu[i];
        }
        sum[i]=(sum[i-1]+1ll*i*i*mu[i])%mod;
    }
}
il int getsum(int x,int y)
{
    int ret=0,last=0;
    for(int i=1;i<=std::min(x,y);i=last+1)
    {
        last=std::min(x/(x/i),y/(y/i));
        ret=(ret+((sum[last]-sum[i-1])%mod*Sum(x/i,y/i))%mod+mod)%mod;
    }
    return ret%mod;
} 
int main()
{
    //freopen("nt2011_table.in","r",stdin);
    //freopen("nt2011_table.out","w",stdout);
    scanf("%d%d",&n,&m);
    getmu();
    int ans=0,last=0;
    for(int i=1;i<=std::min(n,m);i=last+1)
    {
        last=std::min(n/(n/i),m/(m/i));
        ans=(ans+1ll*(i+last)*(last-i+1)/2%mod*getsum(n/i,m/i))%mod;
    }
    printf("%d\n",ans);
    return 0;
}