HDU 6069 数学题,区间素数筛
2017-08-03 20:09
316 查看
Counting Divisors
Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others)Total Submission(s): 691 Accepted Submission(s): 227
Problem Description
In mathematics, the function d(n) denotes
the number of divisors of positive integer n.
For example, d(12)=6 because 1,2,3,4,6,12 are
all 12's
divisors.
In this problem, given l,r and k,
your task is to calculate the following thing :
(∑i=lrd(ik))mod998244353
Input
The first line of the input contains an integer T(1≤T≤15),
denoting the number of test cases.
In each test case, there are 3 integers l,r,k(1≤l≤r≤1012,r−l≤106,1≤k≤107).
Output
For each test case, print a single line containing an integer, denoting the answer.
Sample Input
3
1 5 1
1 10 2
1 100 3
Sample Output
10
48
2302
题意 要你输入 一个区间 l r k 要你算上式,对于求因子数 d(n) , n = (p1^c1*p2^c2...) 其中pi 为质因子,则 d(n) = (c1 + 1) * (c2 + 1) +... + (ci + 1)
而题目要求 d(n^k) = (k * c1 + 1) + (k * c2 + 1) + ... +(k * ci + 1) 所以只要求出这个 n 的质因子就可以了
由于数开的非常大,开到了 10e12,所以一般正常地打素数表是会炸的,于是用到了 区间素数筛
区间素数筛: 对于 n 以内的素数 ,只要筛出了 根号(x)以内的素数就可以筛出x以内的素数 所以对于10e12的n只需要打1e6就可以了,时间复杂度就降下来了
然后题目说 l 和 r 之间的差值不会10e6 所以数组开10e6就可以了,然后在使用的时候用 pos - l 就可以了
然后从离l最近的素数开始去计算值就可以了
#include<stdio.h> #define ll long long #define mod 998244353 #define MAXN 1000500 bool p[MAXN]; ll prm[MAXN]; ll n; ll arr[MAXN],d[MAXN]; void init(){ //区间素数筛 for(ll i = 2;i < MAXN;i++){ if(!p[i]){ prm[n++] = i; } for(int j = 0;j < n;j++){ ll t = i * prm[j]; if(t >= MAXN){ break; } p[t] = true; if(i % prm[j] == 0) break; } } } ll work(ll l,ll r,ll k){ for(int i = 0;i < n;i++){ ll pos = (l + prm[i] - 1) / prm[i] * prm[i]; while(pos <= r){ ll cnt = 0; while(arr[pos - l] % prm[i] == 0){ cnt++; arr[pos - l] /= prm[i]; } d[pos - l] *= (k * cnt + 1); d[pos - l] %= mod; pos += prm[i]; } } ll sum = 0; for(ll pos = l;pos <= r;pos++){ if(arr[pos - l] == 1){ sum += d[pos - l]; }else{ sum += d[pos - l] * (k + 1); } sum %= mod; } return sum; } int main(){ init(); int t; ll l,r,k; scanf("%d",&t); while(t--){ scanf("%lld %lld %lld",&l,&r,&k); for(ll i = l;i <= r;i++){ d[i - l] = 1; arr[i - l] = i; } ll ans = work(l,r,k); printf("%lld\n",ans); } return 0; }
相关文章推荐
- Hdu 6069 Counting Divisors【素数区间筛+预处理素因子分解】
- [hdu 6069]素数筛+区间质因数分解
- HDU 6069 Counting Divisors(枚举区间)(素数筛模版)
- HDU 6069 数论 区间素数筛(+赛后反思
- HDU 5396 Expression(组合数学+区间DP)
- HDU 6069 Counting Divisors (素数+筛法)
- HDU 5396 Expression (数学期望+区间DP)
- hdu 6069 Counting divisors 公式+区间筛
- hdu 5481(数学期望+区间合并)
- hdu 6069 素数
- HDU 4777 - Rabbit Kingdom (树状数组 区间互素数)
- hdu 6069 素数
- Hdu 6069 - Counting Divisors(区间筛质因子)
- hdu 6069 素数
- HDU 5396 区间DP 数学 Expression
- HDU-1006 Tick and Tick(模拟 + 区间交集 + 数学)
- 数学 ( 容斥原理&&求区间内互质数的个数)——Co-prime ( HDU 4135 )
- hdu 6069 素数
- HDU 4611 - Balls Rearrangement(2013MUTC2-1001)(数学,区间压缩)
- hdu 6069 素数