记录编号	359984	评测结果	AAAAAAAAAA
题目名称	2031.[BZOJ 2693] jzptab	最终得分	100
用户昵称	AntiLeaf	是否通过	通过
代码语言	C++	运行时间	16.366 s
提交时间	2016-12-25 21:22:28	内存使用	86.12 MiB

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=10000010,p=100000009,inv=75000007;
void get_table(int);
bool isp[maxn];
int prime[maxn]={0},h[maxn]={0};
int T,n,m,ans;
int main(){
	freopen("bzoj_2693.in","r",stdin);
	freopen("bzoj_2693.out","w",stdout);
	get_table(10000000);
	scanf("%d",&T);
	while(T--){
		scanf("%d%d",&n,&m);
		if(n>m)swap(n,m);
		ans=0;
		for(int i=1,last;i<=n;i=last+1){
			last=min(n/(n/i),m/(m/i));
			ans=(ans+(long long)(n/i)*(n/i+1)%p*(m/i)%p*(m/i+1)%p*((h[last]-h[i-1]+p)%p)%p)%p;
		}
		ans=(long long)ans*inv%p;
		printf("%d\n",ans);
	}
	return 0;
}
void get_table(int n){
	fill(isp+1,isp+n+1,true);
	h[1]=1;
	for(int i=2;i<=n;i++){
		if(isp[i]){
			prime[++prime[0]]=i;
			h[i]=(1-i+p)%p;
		}
		for(int j=1;j<=prime[0]&&i*prime[j]<=n;j++){
			isp[i*prime[j]]=false;
			if(i%prime[j])h[i*prime[j]]=(long long)h[i]*h[prime[j]]%p;
			else{
				h[i*prime[j]]=h[i];
				break;
			}
		}
	}
	for(int i=1;i<=n;i++)h[i]=((long long)h[i]*i%p+h[i-1]+p)%p;
}