bzoj 3518 Dirichlet卷积

详情见代码，回头再填坑。。。

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #define int long long
 #define p 1000000007
 using namespace std;
 int n,m;
 int phi[],su[],pr[],cnt;
 void shai()
 {
     phi[]=;
     for(int j=;j<=;j++)
     {
         if(!phi[j])phi[j]=j-,su[++cnt]=j,pr[j]=j;
         for(int i=;su[i]<=pr[j]&&i<=cnt&&su[i]*j<=;i++)
         {
             pr[su[i]*j]=su[i];
             if(!phi[su[i]*j])phi[su[i]*j]=su[i]*j;
             if(su[i]==pr[j])phi[su[i]*j]=phi[j]*su[i];
             else phi[su[i]*j]=phi[j]*(su[i]-);
         }
     }
 }
 signed main()
 {
     shai();
     scanf("%lld%lld",&n,&m);
     int ans=;
     ans+=n*(n-)*(n-)*m/;ans%=p;
     ans+=m*(m-)*(m-)*n/;ans%=p;
     int tmp=;
     for(int i=;i<=min(n-,m-);i++)tmp=(tmp+m*n%p*phi[i]*((n-)/i)%p*((m-)/i))%p;
     for(int i=;i<=min(n-,m-);i++)tmp=(tmp+phi[i]*i%p*i%p*(((n-)/i)*(+(n-)/i)/)%p*(((m-)/i)*(+(m-)/i)/)%p)%p;
     for(int i=;i<=min(n-,m-);i++)tmp=(tmp-n*phi[i]%p*i%p*((n-)/i)%p*((((m-)/i)*(+(m-)/i)/)%p)%p+p)%p;
     for(int i=;i<=min(n-,m-);i++)tmp=(tmp-m*phi[i]%p*i%p*((m-)/i)%p*((((n-)/i)*(+(n-)/i)/)%p)%p+p)%p;
     tmp-=((+n-)*(n-)/)%p*((+m-)*(m-)/)%p;tmp=(tmp+p)%p;
     ans=(ans+tmp*)%p;
     printf("%lld\n",ans);
     return ;
 }