bzoj 2154 莫比乌斯反演求lcm的和

题目大意：

表格中每一个位置(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 ;
 }