[数论]jb数2

题目大意：

定义jb数为：一个数分解质因数后，所有质因子的和。

例如jb(4)=4,jb(5)=5，jb(6)=5等。

现给出一个范围，求出这个范围内jb数的和。

由于这个数可能过大，只需要求出这个数模1919的值就可以了。

solution：

已知一个数a和质数b,如果a % b=0，那么Jb[a]=b+Jb[a / b]

可以类似筛法的求出每个数是由哪个质数的到的，再次循环即可求出Jb[x]

维护前缀和，输出答案。

输入输出都好大，输入输出挂都加上了  = =

#include<iostream>#include<cstdio>#include<cctype>#include<cstring>using namespace std;const int M=1000010;long long s[2*M],g[2*M],d[2*M];inline int Init(){    char c=getchar();    while(!isdigit(c))c=getchar();    int x=0;    while(isdigit(c))    {        x=x*10+c-'0';        c=getchar();    }    return x;}int buf[10];inline void Out(int i){    int p=0;    if(i==0)buf[p++]=0;    else    while(i)    {        buf[p++]=i%10;        i/=10;    }    for(int j=p-1;j>=0;--j)putchar('0'+buf[j]);}int main(){    int T=Init();    for(int i=2;i<=M;++i)g[i]=i;    for(int i=2;i<=1001;++i)    if(g[i]==i)for(int j=2;j<=M/i+1;++j)g[i*j]=i;    for(int i=2;i<=M;++i)    {        if(g[i]==i)d[i]=i;        else d[i]=d[i/g[i]]+g[i];        s[i]=s[i-1]+d[i];    }    ios::sync_with_stdio(false);    while(T--)    {        int l=Init(),r=Init();        Out((s[r]-s[l-1])%1919);        putchar('\n');    }    return 0;}