题目大意：
表格中每一个位置(i,j)填的值是lcm(i,j) , 求n*m的表格值有多大
论文贾志鹏线性筛中过程讲的很好
最后的逆元我利用的是欧拉定理求解的
我这个最后线性扫了一遍，勉强过了，效率不是很高。。。
/*bzoj 2154*/
#include <bits/stdc++.h> using namespace std;
#define ll long long
#define N 10000000
const int MOD = ;
int mu[N+] , prime[N+] , f[N+] , tot;
bool check[N+]; void init()
{
mu[] = , f[] = ;
for(int i= ; i<=N ; i++){
if(!check[i]){
prime[tot++] = i;
mu[i] = -;
f[i] = -i;
}
for(int j= ; j<tot ; j++){
if((ll)prime[j]*i>N) break;
check[prime[j]*i] = true;
if(i%prime[j]){
mu[prime[j]*i] = -mu[i];
f[prime[j]*i] = ((ll)f[prime[j]]*f[i])%MOD;
}else{
f[prime[j]*i] = f[i];
break;
}
}
}
} int q_pow(int b)
{
ll ans = , a=;
while(b)
{
if(b&) ans = (ans*a)%MOD;
a = (a*a)%MOD;
b>>=;
}
return (int)ans;
} int solve(int n , int m)
{
if(n>m) swap(n , m);
ll ret = ;
for(int k= ; k<=n ; k++){
ll a = n/k , b = m/k;
ret = (ret+(k*(+a)%MOD)*a%MOD*b%MOD*(+b)%MOD*f[k]%MOD)%MOD;
}
ret = ((ll)ret*q_pow(MOD-)%MOD+MOD)%MOD;
return (int)ret;
} int main()
{
// freopen("in.txt" , "r" , stdin);
init();
int n , m;
/* for(int i=2 ; i<=100000 ; i++)
if(!check[i] && MOD%i==0) cout<<i<<" yes "<<endl;*/
while(~scanf("%d%d" , &n , &m)){
printf("%d\n" , solve(n , m));
}
return ;
}