#include<bits/stdc++.h>
#define int __int128
using namespace std;
const int mod=19940417;
int n,m,s1,s2,s3;
int read()
{
	int x=0;
	char c=getchar();
	while(c<'0'||c>'9') c=getchar();
	while(c>='0'&&c<='9')
	{
		x=(x<<3)+(x<<1)+(c&15);
		c=getchar();
	}
	return x;
}
void print(int x)
{
	if(x>9) print(x/10);
	putchar(x%10+'0');
}
int po(int x,int y)
{
	int res=1ll,r=x;
	for(int i=y;i;i>>=1)
	{
		if(i&1) res=res*r%mod;
		r=r*r%mod;
	}
	return res;
}
signed main()
{
	n=read(),m=read();
	for(int l=1,r=0;l<=n;l=r+1)
	{
		r=n/(n/l);
		s1=(s1+(r-l+1)*n%mod-(n/l)*((l+r)*(r-l+1)/2%mod)%mod+mod)%mod;
	}
	for(int l=1,r=0;l<=m;l=r+1)
	{
		r=m/(m/l);
		s2=(s2+(r-l+1)*m%mod-(m/l)*((l+r)*(r-l+1)/2%mod)%mod+mod)%mod;
	}
	for(int l=1,r=0;l<=min(n,m);l=r+1)
	{
		r=min(n/(n/l),m/(m/l));
		s3=(s3+(r-l+1)*n%mod*m%mod-n*(m/l)%mod*((l+r)*(r-l+1)/2%mod)%mod-m*(n/l)%mod*((l+r)*(r-l+1)/2%mod)%mod+(n/l)*(m/l)%mod*((r*(r+1)*(r*2%mod+1)/6%mod)%mod-(l*(l-1)*(l*2%mod-1)/6%mod)+mod)%mod+mod*3)%mod;
	}
	print((s1*s2%mod-s3+mod)%mod);
}