2693: jzptab

https://www.cnblogs.com/ljh2000-jump/p/6358095.html
这个链接推得和我一样，线性筛F函数
最后算的时候不能T 和F(T)都求前缀和
我们可以把F[T]=F[T]*T;
这样就可以求了；

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const int N=10000000+1000,MOD=100000009;
int mu[N],F[N],s[N],prime[N],tot;
int n,m; bool flag[N];
LL ans;

void prepare()
{
    mu[1]=F[1]=1;
    for (int i=2;i<N;i++)
    {
        if (!flag[i]) 
        {
            prime[++tot]=i;
            mu[i]=-1;
            F[i]=((1-i)%MOD+MOD)%MOD;
        }
        for (int j=1;j<=tot;j++)
        {
            int p=prime[j]*i;
            if (p>=N) break;
            flag[p]=true;
            if (i%prime[j]==0)
            {
                F[p]=F[i];
                mu[i]=0;
                break;
            }
            F[p]=(LL)F[i]*F[prime[j]]%MOD;
            mu[p]=-mu[i];
        }
    }
    for (int i=2;i<N;i++) F[i]=(F[i-1]+(LL)F[i]*i)%MOD;
    for (int i=1;i<N;i++) s[i]=(s[i-1]+i)%MOD;
}
void solve()
{
    ans=0;
    scanf("%d%d",&n,&m);
    if (n>m) swap(n,m);
    for (int i=1,last=0;i<=n;i=last+1)
    {
        last=min(n/(n/i),m/(m/i));
        ans+=(LL)s[n/i]*s[m/i]%MOD*(F[last]-F[i-1]+MOD)%MOD;
        ans%=MOD;
    }
    printf("%lld\n",ans);
}
int main()
{
    prepare();
    int T;
    scanf("%d",&T);
    while (T--) solve();
    return 0;
}