HDU 6069 Counting Divisors(数论)
2017-08-07 09:43
337 查看
Description
Input
第一行一整数T表示用例组数,每组用例输入三个整数l,r,k
(1<=T<=15,1<=l<=r<=1e12,r-l<=1e6,1<=k<=1e7)
Output
对于每组用例,输出答案
Sample Input
3
1 5 1
1 10 2
1 100 3
Sample Output
10
48
2302
Solution
Code
Input
第一行一整数T表示用例组数,每组用例输入三个整数l,r,k
(1<=T<=15,1<=l<=r<=1e12,r-l<=1e6,1<=k<=1e7)
Output
对于每组用例,输出答案
Sample Input
3
1 5 1
1 10 2
1 100 3
Sample Output
10
48
2302
Solution
Code
#include<cstdio> #include<cstring> using namespace std; typedef long long ll; const int maxn=1000011,mod=998244353; int mark[maxn],prime[maxn],res=0; void get_prime(int n=1000000) { for(int i=2;i<=n;i++) { if(!mark[i])mark[i]=prime[res++]=i; for(int j=0;j<res&&prime[j]*i<=n;j++) { mark[i*prime[j]]=prime[j]; if(i%prime[j]==0) break; } } } int T,k,n,ans[maxn]; ll l,r,a[maxn]; int Solve() { for(int i=0;i<res&&prime[i]*prime[i]<=r;i++) { int p=prime[i]; ll t=((l+p-1)/p)*p; for(int j=t-l+1;j<=n;j+=p) { int num=0; while(a[j]%p==0)a[j]/=p,num++; ans[j]=(ll)ans[j]*((ll)k*num%mod+1)%mod; } } int sum=0; for(int i=1;i<=n;i++) { if(a[i]>1)ans[i]=(ll)ans[i]*(k+1)%mod; sum=(sum+ans[i])%mod; } return sum; } int main() { get_prime(); scanf("%d",&T); while(T--) { scanf("%I64d%I64d%d",&l,&r,&k); n=r-l+1; for(int i=1;i<=n;i++)a[i]=l+i-1,ans[i]=1; printf("%d\n",Solve()); } return 0; }
相关文章推荐
- 【数论】hdu 6069 Counting Divisors
- HDU_6069 Counting Divisors 【数论】
- 【hdu 6069】 Counting Divisors 【思维+数论】
- HDU 6069 数论 区间素数筛(+赛后反思
- HDU 6069 Counting Divisors 【数论】
- 【原创】【数论】HDU 2462 The Luckiest number
- hdu 1019 Least Common Multiple(数论:求最小公倍数)
- hdu 1124 Factorial(数论:推理)
- hdu 5812 - 数论
- HDU 5778 abs(数论)
- hdu 1452 Happy 2004(数论:积性函数+快速幂+同余方程+扩展欧几里得算法)
- hdu 6069 区间筛
- HDU --- 5446 Unknown Treasure 数论综合【Lucas定理 + 中国剩余定理 + 快速乘 + 思维】
- hdu 6069 素数
- 【数论】HDU 4869 Turn the pokers 费马小定理
- POJ 1423 && HDU 1018 Big Number(数论)
- hdu 4002 Find the maximum——数论
- hdu 6069 Counting Divisors(约数个数)(2017 Multi-University Training Contest - Team 4 )
- hdu 4365 瞎涂颜色 数论
- hdu 2879 HeHe 数论 积性函数 优化技巧