这份代码为什么RE了？

#include<bits/stdc++.h>
#define int long long
const int maxn=1e7;
const int mod=20101009;
int mu[maxn+10],h[maxn+10],pri[maxn/100],tot;
int n,m;
inline void init(){
	mu[1]=1;
	for(int i=2;i<=n;i++){
		if(!h[i])mu[i]=-1,pri[++tot]=i;
		for(int j=1;j<=tot&&pri[j]*i<=n;j++){
			h[pri[j]*i]=1;
			if(i%pri[j]==0)break;
			mu[pri[j]*i]=-mu[i];
		}
	}
	memset(h,0,sizeof(h));
	for(int i=1;i<=n;i++)
		if(mu[i]!=0)
			for(int j=1;j<=n/i;j++)
				h[i*j]=(h[i*j]+mu[i]*i)%mod;
	for(int i=1;i<=n;i++)
		h[i]=(i*h[i]+h[i-1])%mod;
}
signed main(){
	scanf("%lld%lld",&n,&m);
	if(n>m)std::swap(n,m);
	init();
	int ans=0;
	for(int l=1,r;l<=n;l=r+1){
		r=std::min(n/(n/l),m/(m/l));
		ans=(ans+((((n/l)*(n/l+1)/2)%mod)*(((m/l)*(m/l+1)/2)%mod)%mod*(h[r]-h[l-1])%mod))%mod;
	}
	printf("%lld\n",(ans%mod+mod)%mod);
	return 0;
}