BZOJ2154: Crash的数字表格

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

题意&&题解：http://www.cnblogs.com/jianglangcaijin/archive/2013/11/27/3446169.html

我只是来发泄的！！！怎么每次打数论题都会被爆int，爆long long这种问题纠缠好久！！！！！！！！！！

我还以为这题是双倍经验呢

代码：

#include<cstdio>

#include<cstdlib>

#include<cmath>

#include<cstring>

#include<algorithm>

#include<iostream>

#include<vector>

#include<map>

#include<set>

#include<queue>

#include<string>

#define inf 1000000000

#define maxn 10000000+5

#define eps 1e-10

#define ll long long

#define pa pair<int,int>

#define for0(i,n) for(int i=0;i<=(n);i++)

#define for1(i,n) for(int i=1;i<=(n);i++)

#define for2(i,x,y) for(int i=(x);i<=(y);i++)

#define for3(i,x,y) for(int i=(x);i>=(y);i--)

#define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)

#define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)
#define mod 20101009

using namespace std;

inline int read()

{

    int x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}

    return x*f;

}
int p[maxn],tot,mu[maxn];
bool v[maxn];
void get(int m)
{
    mu[1]=1;
    for2(i,2,m)
    {
        if(!v[i])p[++tot]=i,mu[i]=-1;
        for1(j,tot)
        {
            int k=i*p[j];
            if(k>m)break;
            v[k]=1;
            if(i%p[j])mu[k]=-mu[i];
            else {mu[k]=0;break;}
        }
    }
    for1(i,m)mu[i]=((ll)mu[i]*i%mod*i%mod+mu[i-1])%mod;
}
inline int sum(int n,int m)
{
    return ((ll)n*(n+1)/2%mod)*((ll)m*(m+1)/2%mod)%mod;
}
inline int f(int n,int m)
{
    int ret=0;
    if(n>m)swap(n,m);
    for(int i=1,j;i<=n;i=j+1)
    {
        j=min(n/(n/i),m/(m/i));
        ret=(ret+(ll)(mu[j]-mu[i-1])*sum(n/i,m/i)%mod)%mod;
    }
    return ret;
}

int main()

{

    freopen("input.txt","r",stdin);

    freopen("output.txt","w",stdout);
        int n=read(),m=read(),ans=0;get(m);
        if(n>m)swap(n,m);
        for(int i=1,j;i<=n;i=j+1)
        {
            j=min(n/(n/i),m/(m/i));
            ans=(ans+(ll)(i+j)*(j-i+1)/2%mod*f(n/i,m/i)%mod)%mod;
        }
        printf("%d
",(ans+mod)%mod);

    return 0;

}