BZOJ [清华集训2013] 模积和 数论 整除分块

求${\sum }_{i=1}^n{\sum }_{j=1}^m[i̸=j](n\%i)(m\%i)$，其中$n,m\le 1e9$
直接利用上一篇的数论分块的知识。对原来式子化简实际上只要求
${\sum }_{i=1}^n(n\%i){\sum }_{i=1}^m(m\%i)-{\sum }_{i=1}^{min(n,m)}(n\%i)(m\%i)$
全部可以在$O(\sqrt{n})$的时间复杂度里面解决。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const ll INF=LONG_LONG_MAX;
const int mod=19940417;
const int inv6=3323403;
ll cal(ll n) {
	n%=mod;
	return n*(n+1)%mod*(2*n+1)%mod*inv6%mod;
}
ll solve1(ll n) {
	ll ans=0;
	for(ll r,l=1;l<=n;l=r+1) {
		r=n/(n/l);
		ans=(ans+((l+r)*(r-l+1)/2)%mod*(n/l))%mod;
	}
	return (n*n%mod-ans+mod)%mod;
}
ll solve2(ll n,ll m) {
	if(n>m) swap(n,m);
	ll ans1=0,ans2=0,ans3=0;
	for(ll r,l=1;l<=n;l=r+1) {
		r=min(n/(n/l),m/(m/l));
		ans1=(ans1+((cal(r)-cal(l-1)+mod)%mod)*(n/l)%mod*(m/l))%mod;
		ans2=(ans2+((r+l)*(r-l+1)/2%mod)*(m/l))%mod;
		ans3=(ans3+((r+l)*(r-l+1)/2%mod)*(n/l))%mod;
	}
	ans1=(ans1-ans2*n%mod+mod)%mod;
	ans1=(ans1-ans3*m%mod+mod)%mod;
	return ans1;
}
int main() {
	ll n,m;
	scanf("%lld%lld",&n,&m);
	ll ans=solve1(n);
	ans=(ans*solve1(m))%mod; 
	ans=(ans-(n*m)%mod*min(n,m)%mod+mod)%mod;
	ans=(ans-solve2(n,m)%mod+mod)%mod;
	printf("%lld\n",ans%mod);
	return 0;
}