2693: jzptab

吐血！！！

其实理论上我本能A掉这题，但由于中途少算了个x^2，又烧了好几个小时，最后还是在连击哥的帮助下、看了答案才查出错来。

表示已经失去自己解决问题的能力了。。。

讲道理以前考虑过报数竞，但由于怕算错，果断报了信竞，现在这些数学题让我情何以堪？

不管怎样，还是感谢Xrzl的耐心教导。

恶心的代码：

#pragma GCC optimize("O3")
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define ll long long
#define pli pair<ll,int>
#define mkp make_pair
#define X first
#define Y second
#define N 10000005
#define MOD 100000009
ll _,n,m,k,ans;
inline ll add(ll a,ll b){
	return (a+b)%MOD;
}
inline ll sub(ll a,ll b){
	return (a-b)%MOD;
}
inline ll mul(ll a,ll b){
	return a*b%MOD;
}
inline ll cal(ll a){
	return a*(a+1)/2%MOD;
}
ll pri[N],cnt,mu[N],Sf[N],f[N];bool vis[N];
void shai(){
	ll i,j,v,rg=N-5;
	mu[1]=f[1]=Sf[1]=1;
	rep(i,2,rg){
		if(!vis[i])mu[pri[++cnt]=i]=-1,f[i]=1-i;
		for(j=1;j<=cnt&&(v=pri[j]*i)<=rg;j++){
			vis[v]=1;
			if(i%pri[j]==0){
				f[v]=f[i];
				break;
			}
			mu[v]=-mu[i];
			f[v]=mul(f[i],1-pri[j]);
		}
		Sf[i]=add(f[i]=mul(f[i],i),Sf[i-1]);
	}
}
int main(){
	shai();
	scanf("%lld",&_);
	ll i,j,rg;
	while(_--){
		scanf("%lld%lld",&n,&m);
		rg=min(n,m);ans=0;
		for(i=1;i<=rg;i=j+1){
			j=min(n/(n/i),m/(m/i));
			ans=add(ans,mul(mul(cal(n/i),cal(m/i)),sub(Sf[j],Sf[i-1])));
		}
		printf("%lld\n",add(ans,MOD));
	}
	return 0;
}