BZOJ2693: jzptab

题目：http://www.lydsy.com/JudgeOnline/problem.php?id=2693

题意：同2154 多组数据

题解：按2154再往后转化一下就可以把n，m放到一边儿，然后线性筛右边。

硬要看推导的话可以戳：http://blog.csdn.net/popoqqq/article/details/42078725

不会latex感觉跟不上时代了T_T

还有此题模数真是坑爹

代码：

 #include<cstdio>

 #include<cstdlib>

 #include<cmath>

 #include<cstring>

 #include<algorithm>

 #include<iostream>

 #include<vector>

 #include<map>

 #include<set>

 #include<queue>

 #include<string>

 #define inf 1000000000

 #define maxn 10000000+5

 #define maxm 10000000

 #define eps 1e-10

 #define ll long long

 #define pa pair<int,int>

 #define for0(i,n) for(int i=0;i<=(n);i++)

 #define for1(i,n) for(int i=1;i<=(n);i++)

 #define for2(i,x,y) for(int i=(x);i<=(y);i++)

 #define for3(i,x,y) for(int i=(x);i>=(y);i--)

 #define for4(i,x) for(int i=head[x],y=e[i].go;i;i=e[i].next,y=e[i].go)

 #define for5(n,m) for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)

 #define mod 100000009

 using namespace std;

 inline int read()

 {

     int x=,f=;char ch=getchar();

     while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}

     while(ch>=''&&ch<=''){x=*x+ch-'';ch=getchar();}

     return x*f;

 }
 int tot,p[maxn],mu[maxn];
 bool v[maxn];
 void get()
 {
     mu[]=;
     for2(i,,maxm)
     {
         if(!v[i])p[++tot]=i,mu[i]=(i-(ll)i*i)%mod;
         for1(j,tot)
         {
             int k=i*p[j];
             if(k>maxm)break;
             v[k]=;
             if(i%p[j])mu[k]=(ll)mu[p[j]]*mu[i]%mod;
             else {mu[k]=(ll)p[j]*mu[i]%mod;break;}
         }
     }
     for1(i,maxm)mu[i]=(mu[i]+mu[i-])%mod;
 }
 inline int sum(int n,int m)
 {
     return ((ll)n*(n+)/%mod)*((ll)m*(m+)/%mod)%mod;
 }

 int main()

 {

     freopen("input.txt","r",stdin);

     freopen("output.txt","w",stdout);
     get();

     int T=read();
     while(T--)
     {
         int n=read(),m=read(),ans=;if(n>m)swap(n,m);
         for(int i=,j;i<=n;i=j+)
         {
             j=min(n/(n/i),m/(m/i));
             ans=(ans+(ll)(mu[j]-mu[i-])*sum(n/i,m/i)%mod)%mod;
         }
         printf("%d\n",(ans+mod)%mod);
     }

     return ;

 }