bzoj千题计划253：bzoj2154: Crash的数字表格

http://www.lydsy.com/JudgeOnline/problem.php?id=2154

#include<cstdio>
#include<algorithm>

using namespace std;

#define N 10000001

const int mod=;

int cnt,p[N];
int miu[N];
bool vis[N];

int sum[N];

void pre()
{
    miu[]=;
    for(int i=;i<N;++i)
    {
        if(!vis[i])
        {
            p[++cnt]=i;
            miu[i]=-;
        }
        for(int j=;j<=cnt;++j)
        {
            if(p[j]*i>=N) break;
            vis[p[j]*i]=true;
            if(!(i%p[j])) break;
            miu[i*p[j]]=-miu[i];
        }
    }
    for(int i=;i<N;++i) sum[i]=(sum[i-]+1LL*i*i*miu[i]%mod)%mod;
}

int get_S(int n,int m)
{
    return 1LL*n*(n+)/%mod*(1LL*m*(m+)/%mod)%mod;
}

int get_F(int n,int m)
{
    if(n>m) swap(n,m);
    int ans=,j;
    for(int i=;i<=n;i=j+)
    {
        j=min(n/(n/i),m/(m/i));
        ans=(ans+1LL*(sum[j]-sum[i-]+mod)%mod*get_S(n/i,m/i)%mod)%mod;
    }
    return ans;
}

int main()
{
    pre();
    int n,m;
    scanf("%d%d",&n,&m);
    if(n>m) swap(n,m);
    int ans=,j;
    for(int i=;i<=n;i=j+)
    {
        j=min(n/(n/i),m/(m/i));
        ans=(ans+(1LL*(j+i)*(j-i+)/%mod*get_F(n/i,m/i)%mod))%mod;
    }
    printf("%d",ans);
}