#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 ;

 }