2017 Multi-University Training Contest - Team 4：1003. Counting Divisors（积性函数）

最新推荐文章于 2017-08-23 23:32:00 发布

原创最新推荐文章于 2017-08-23 23:32:00 发布 · 402 阅读

0 ·

CC 4.0 BY-SA版权

Extra 专栏收录该内容

148 篇文章

订阅专栏

本文介绍了一种通过质因数分解来求解特定数学问题的算法。该算法主要针对给定区间内的整数进行处理，通过分解小于sqrt(r)的质数，并考虑大于sqrt(r)的质数的影响，来计算特定形式的数的个数。

公式：

上面的a|b表示a是b的约数，ai是每个质因数p分解后的项数

枚举每个x(l<=x<=r)

对于每个数暴力分解小于sqrt(r)的质数就好了，因为大于sqrt(r)的质数最多只有一个，且对应的ai也只能为1

也就是说如果分解完那个数不是1的话就将d(x^k)乘上(k+1)

最后全加在一起

#include<stdio.h>
#include<math.h>
#define mod 998244353
#define LL long long
LL cnt, flag[1000005] = {1,1}, pri[1000005], val[1000005], d[1000005];
int main(void)
{
	LL T, l, r, k, i, j, ans, sum;
	for(i=2;i<=1000001;i++)
	{
		if(flag[i])
			continue;
		for(j=i*i;j<=1000001;j+=i)
			flag[j] = 1;
		pri[++cnt] = i;
	}
	scanf("%lld", &T);
	while(T--)
	{
		ans = 0;
		scanf("%lld%lld%lld", &l, &r, &k);
		for(i=l;i<=r;i++)
			val[i-l] = i, d[i-l] = 1;
		for(i=1;pri[i]*pri[i]<=r;i++)
		{
			sum = l%pri[i];
			j = (pri[i]-sum)%pri[i];
			for(;j<=r-l;j+=pri[i])
			{
				sum = 0;
				while(val[j]%pri[i]==0)
					sum++, val[j] /= pri[i];
				d[j] = d[j]*(sum*k+1)%mod;
			}
		}
		for(i=0;i<=r-l;i++)
		{
			if(val[i]>1)
				d[i] = (d[i]*(k+1))%mod;
			ans = (ans+d[i])%mod;
		}
		printf("%lld\n", ans);
	}
	return 0;
}