记录编号	600797	评测结果	AAAAAAAAAAAAAAAAAAAA
题目名称	[国家集训队 2011] Crash的数字表格	最终得分	100
用户昵称	李奇文	是否通过	通过
代码语言	C++	运行时间	1.203 s
提交时间	2025-05-16 22:34:07	内存使用	24.29 MiB

#include<bits/stdc++.h>
using namespace std;
const int N=1e7+5;
const int mod=20101009;
int n,m,mu[N],p[N],sum[N];
bool flg[N];
int tot;
void mobius(){
	mu[1]=1;
	int x=min(n,m);
	for(int i=2;i<=x;++i){
		if(!flg[i]) p[++tot]=i,mu[i]=-1;
		for(int j=1;j<=tot&&i*p[j]<=x;++j){
			flg[i*p[j]]=1;
			if(i%p[j]==0){
				mu[i*p[j]]=0;
				break;
			} 
			mu[i*p[j]]=-mu[i];
		}
	}
	for(int i=1;i<=x;++i) sum[i]=(sum[i-1]+1ll*i*i%mod*(mu[i]+mod))%mod; 
}
int g(int x,int y){
	return (1ll*x*(x+1)/2%mod)*(1ll*y*(y+1)/2%mod)%mod;
}
int Sum(int x,int y){
	int res=0;
	for(int i=1,j;i<=min(x,y);i=j+1){
		j=min(x/(x/i),y/(y/i));
		res=(res+1ll*(sum[j]-sum[i-1]+mod)*g(x/i,y/i)%mod)%mod;
	}
	return res;
}
int main(){
	freopen("nt2011_table.in","r",stdin);
	freopen("nt2011_table.out","w",stdout);
	cin>>n>>m;
	mobius();
	int ans=0;
	for(int i=1,j;i<=min(n,m);i=j+1){
		j=min(n/(n/i),m/(m/i));
		ans=(ans+1ll*(j-i+1)*(i+j)/2%mod*Sum(n/i,m/i)%mod)%mod;
	}
	cout<<ans<<endl;
	return 0;
}