hdu 6069 Counting Divisors

/*    1~1e6的  素数个数 约 等于 1e6/ lg1e6 =78498 ,直接对区间内每个数作因数 分析会超时， 那可以枚举 质因数，再类似于 筛法，累计每个数 当前的因数和 */ #include<iostream>#include<algorithm>#include<cstring>#include<cstdio>#define rep(i,a,b) for(int i=a;i<=b;i++)using namespace std;typedef long long ll;const int mod =998244353;const int Size = 1e6+105;int vis[Size],e[Size],cnt=0;int ans[Size];ll as[Size];void prime(){int q=1e6+3; for(int i=2;i<q;i++) if(!vis[i]) {   e[++cnt]=i;   for(int j=i+i;j<q; j += i) vis[j]=1;     }}int main(){memset(vis,0,sizeof(vis));prime();    int k,t;    scanf("%d",&t);    ll l,r;    while(t--)     {      scanf("%lld%lld%d",&l,&r,&k);      for(int i=1;i<=r-l+1;i++)     {    ans[i]=1;    as[i]=l-1+i;        }         for(int i=1;i<=cnt&&e[i]<=r;i++)          {           ll x = l/e[i];          ll p = (ll)x*e[i];          if(p<l)    {     p +=e[i];   }          for(;p<=r;p +=e[i])           {            ll nb=p;            int xx=0;            int pos=p-l+1;            while(nb%e[i]==0)             {             xx++;             nb /=e[i];             as[pos] /=e[i];  }  ll nn=ans[pos]; nn *=(xx*k+1)%mod; nn %=mod;  ans[pos] =nn; }         } ll ansans=0;for(int i=1;i<=r-l+1;i++) {  if(as[i]!=1)   {  ll nn=(k+1); nn *=ans[i]; nn %=mod;  ans[i] =nn;  }  ansans +=ans[i];  ansans %=mod;    }    printf("%lld\n",ansans);}return 0;}