BZOJ2693: jzptab

沃日 昨天过了那个单次询问的。。。具体内幕我还是不说了 说多了都是泪
我居然是Rank1？

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;

#define ll long long
const
   int  maxn=10000011;
const
  int  Mod=100000009;
ll  Pre[maxn],tp;
bool check[maxn];
int  prime[maxn],tot;
int  sum[maxn];

int main()
{
    int  i,j,k;
    Pre[1]=1;
    sum[1]=1;
    int  n,m,T;
    for(i=2;i<maxn;i++)
      {
        if(!check[i])prime[++tot]=i,Pre[i]=i-i*1ll*i;
        tp=Pre[i];
        Pre[i]+=Pre[i-1];
        if(Pre[i]>=Mod||Pre[i]<=-Mod)Pre[i]%=Mod;
        for(j=1;j<=tot;j++)
         {
            k=prime[j]*i;
            if(k>maxn)break;
            else if(i%prime[j]==0){check[k]=true,Pre[k]=prime[j]*tp;if(Pre[k]>=Mod||Pre[k]<=-Mod)Pre[k]%=Mod;;break;}
            check[k]=true,Pre[k]=((-prime[j]+1ll)*prime[j]%Mod)*tp;
            if(Pre[k]>=Mod||Pre[k]<=-Mod)Pre[k]%=Mod;
         }
      }
    scanf("%d",&T);
    int t=0;
    for(i=2;i<maxn;i++)sum[i]=(sum[i-1]+i)%Mod;
    while(T--)
    {
    scanf("%d%d",&n,&m);
    if(n>m)swap(n,m);
      int  G=1,next;
      int  ans=0;
      while(G<=n)
      {
        next=min(m/(m/G)+1,n/(n/G)+1);
        ans+=(sum[n/G]*1ll*sum[m/G]%Mod)*(Pre[next-1]-Pre[G-1])%Mod;
        ans%=Mod;
        G=next;
        t++;
      }
      printf("%d\n",(ans+Mod)%Mod);
 }
}