记录编号	367349	评测结果	AAAAAAAAAA
题目名称	[BZOJ 2693] jzptab	最终得分	100
用户昵称	AntiLeaf	是否通过	通过
代码语言	C++	运行时间	18.761 s
提交时间	2017-01-30 11:10:16	内存使用	77.51 MiB

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=10000010,p=100000009,inv_2=50000005;
void get_table(int);
bool notp[maxn]={false};
int prime[maxn]={0},f[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*((f[last]-f[i-1])%p)%p)%p;
		}
		if(ans<0)ans+=p;
		ans=(long long)ans*inv_2%p*inv_2%p;
		printf("%d\n",ans);
	}
	return 0;
}
void get_table(int n){
	f[1]=1;
	for(int i=2;i<=n;i++){
		if(!notp[i]){
			prime[++prime[0]]=i;
			f[i]=(1-i)%p;
		}
		for(int j=1;j<=prime[0]&&i*prime[j]<=n;j++){
			notp[i*prime[j]]=true;
			if(i%prime[j])f[i*prime[j]]=(long long)f[i]*f[prime[j]]%p;
			else{
				f[i*prime[j]]=f[i];
				break;
			}
		}
	}
	for(int i=2;i<=n;i++)f[i]=((long long)f[i]*i%p+f[i-1])%p;
}