[BZOJ 2693] jzptab

2693: jzptab

Time Limit: 10 SecMemory Limit: 512 MB

Description

同BZOJ 2154

Input

一个正整数T表示数据组数

接下来T行 每行两个正整数 表示N、M

Output

T行 每行一个整数 表示第i组数据的结果

 

Sample Input

1

4 5

Sample Output

122

HINT
T <= 10000
N, M<=10000000

【题解】

同BZOJ 2154，只是有多组数据，需要用到前缀和优化

注意能预处理尽量预处理，还有就是ans要清零！！

#include<bits/stdc++.h>
using namespace std;
int prime[1000010];
int tot,mz[10000010];
int f[10000010],n,m;
long long ans;
const int mod=100000009;
long long sum(int ki,int kj) {
    return ((((long long)ki*(ki+1)/2)%mod)*(((long long)kj*(kj+1)/2)%mod))%mod;
}
int main() {
    //freopen("a.out","w",stdout);
    int T;int nx=10000000;
    scanf("%d",&T);
    for (int i=2;i<=nx;++i) {
        if (!mz[i]) {
            prime[++tot]=i;
            mz[i]=i;
        }
        for (int j=1;j<=tot&&i*prime[j]<=nx;++j) {
            mz[i*prime[j]]=prime[j];
            if(i%prime[j]==0) break;
        }
    }
    //for (int i=1;i<=tot;++i) cout<<prime[i]<<endl;
    //for (int i=1;i<=n;++i) cout<<mz[i]<<endl;
    f[1]=1;
    for (int i=2;i<=nx;++i)
        if(mz[i/mz[i]]==mz[i]) f[i]=f[i/mz[i]];
        else f[i]=((long long)f[i/mz[i]]*(1-mz[i]))%mod;
    for (int i=2;i<=nx;++i)
        f[i]=((long long)f[i]*i+f[i-1])%mod;
    while(T--) {ans=0;
        scanf("%d%d",&n,&m);
        if (n>m) {
            int t=n;
            n=m;
            m=t;
        }int j=0;
        for (int i=1;i<=n;i=j+1) {
            j=min(n/(n/i),m/(m/i));
            ans=(ans+(sum(n/i,m/i)*(long long)(f[j]-f[i-1])%mod))%mod;
        }
        printf("%lld
",(ans+mod)%mod);
    }
    
    return 0;
}