首页 > 代码库 > hdu 6069
hdu 6069
思路:n=p1^x1*p2^x2....pm^xm,则p的约数个数为(x1+1)*(x2+1)....(xm+1),那么n^k=p1^(x1+k)....pm^(xm+k),约数个数为(x1*k+1)*....*(xm*k+1)。
先求出1-1e6内的质数,再对l--r之间的数求xi
1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long ll; 4 const int N=1e6+10; 5 const ll mod=998244353; 6 7 ll prime[N],vis[N]; 8 9 int init(){ 10 int x=0; 11 for(ll i=2;i<=N;i++){ 12 if(!vis[i]) { 13 // cout<<i<<endl; 14 prime[++x]=i; 15 } 16 else continue; 17 for(ll j=i*i;j<=N;j+=i){ 18 vis[j]=1; 19 } 20 } 21 return x; 22 } 23 ll a[N],b[N]; 24 int main(){ 25 int t; 26 ll l,r,k; 27 int len=init(); 28 scanf("%d",&t); 29 while(t--){ 30 scanf("%lld%lld%llld",&l,&r,&k); 31 for(int i=0;i<=N-1;i++){ 32 a[i]=b[i]=1; 33 } 34 for(int i=1;i<=len;i++){ 35 ll p=prime[i];ll L; 36 if(l%p!=0) L=(l/p+1)*p; else L=l; 37 for(ll j=L;j<=r;j+=p){ 38 ll x=0; 39 ll base=1, y=j; 40 while(y%p==0){ 41 x++; 42 y/=p; 43 base*=p; 44 } 45 a[j-l+1]=a[j-l+1]*(x*k+1)%mod; 46 b[j-l+1]*=base; 47 48 } 49 } 50 ll sum=0; 51 for(ll i=l;i<=r;i++){ 52 if(b[i-l+1]!=i) 53 a[i-l+1]=a[i-l+1]*(k+1)%mod; 54 sum=(sum+a[i-l+1])%mod; 55 } 56 printf("%lld\n",sum); 57 } 58 }
hdu 6069
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。