HDU6069 Counting Divisors 区间素数筛法

题目链接：HDU6069

题目大意：给你l,r,k，给你一个公式，如果最后结果很大，取模在输出。

代码参考：大佬的代码

思路：这题当时没想出来，结束之后对着大佬的代码看了蛮久终于有点明白了，原博主讲的蛮详细，我就写一下注释，方便大家理解。

1、首先需要知道一个结论，比如 100有多少个因子？最朴素的方法从1遍历到100 一直除对吧，现在我们稍微转换一下  100=10^2=(2*5)^2=2^2*5^2  这个时候100的因子就可以用 （2+1）*（2+1）=9 为什么能这样算呢？ 2的2次方，对5的2次方分别做乘法也就是  2^0 2^1 2^2 分别对 5^0 5^1 5^2 乘 ，就是3*3=9。这个应该非常好懂。所以我们求因子的个数就是把一个数分解成质因子的乘积的形式。那么我们怎么去求一个数的质因子呢？

2、一个定理：任何一个合数,只需要测试它的平方根以下的所有质数即可将这个数分解质因数. 

证明不会……从网上找了一段解释，方便理解。

任何一个数n,如果知道这个数有k个质因子,只需要测试它的k方根以下的所有质数即可将这个数分解质因数. 
现在将这个设为2个质因子... 
a*b = n 
若b>=a,则a<=sqrt(n). 为何? 简单不等式.  
a小于sqrt(n),测试小于sqrt(n)的所有整数也就可以测试n是否是质数. 
同理,可以扩展到k个质因子...应该不是很难想对不对? 
因为我们不知道这个数有多少个质因子...所以达到了sqrt(n)都没有发现任何质因子,说明n是质数... 
同样的方法也可以知道一个数字里面最多有多少个质因子. 
比如到了n^(1/k)都没有发现任何质因子,说明这个数最多也只能有k-1个质因子

AC代码：

/*
2017年8月4日15:11:54
HDU6069
AC 
*/
#include<stdio.h>
#include<string.h>
#include<math.h>

typedef long long ll;
const ll mod=998244353;
const ll maxn=1e6+10;
int flag[maxn],cnt;
ll prime[maxn],f[maxn],num[maxn];
ll l,r,k;
//生成1e6以内所有素数 
void isprime(){
	memset(flag,0,sizeof(flag));
	cnt=0;
	flag[1]=1;
	for(ll i=2;i<maxn;i++){
		if(!flag[i]){
			prime[cnt++]=i;
			for(ll j=i*i;j<maxn;j+=i) flag[j]=1;
		} 
	}
}

int main(){
	isprime();
	 int t;
	 scanf("%d",&t);
	 while(t--){
	 	scanf("%I64d%I64d%I64d",&l,&r,&k);
	 	ll ans=0;
	 	/*如果左端点是1，1的因子就是本身，所以ans=1，然后左端点右移一位l=2*/
	 	if(l==1) ans=1,l=2;
	 	
		 ll t=r-l;
		 /*f[]数组代表l-r中的每一个数字，初始值设为其本身 
		   num[i]数组代表l-r中第i个数字的所有因子，初始值设为1，*/ 
		 for(ll i=0;i<=t;i++)  f[i]=l+i,num[i]=1;
		 /*遍历素数表[0-cnt)，并且要求prime[i]*prime[i]<=r*/
		 for(ll i=0;i<cnt&&prime[i]*prime[i]<=r;i++){
		 	ll temp=l;
		 	/*如果temp小于第i个质数，那么就把temp赋值为该质数*/
		 	if(temp%prime[i]) temp=(temp/prime[i]+1)*prime[i];
		 	/*从该质数开始，把该区间所有这个质数的倍数分解成质因子的乘积形式*/
		 	for(ll j=temp;j<=r;j+=prime[i]){
		 		ll time=0;
		 		/*如果能整除，证明该数是质因子的倍数，一直除到不能整除*/
		 		while(f[j-l]%prime[i]==0){
		 			time++;//记录该质因子的次数 
		 			f[j-l]/=prime[i];
				 }
				 //存储这个因子对num[]的贡献 
				 num[j-l]=(num[j-l]*(time*k+1))%mod;
			 }
		 } 
		 for(ll i=0;i<=t;i++){
		 	//如果该数已经分解完 直接用num[i] 
		 	if(f[i]==1) ans=(ans+num[i])%mod;
		 	//如果该数没有被分解，证明该数是大于1e6的在l-r区间的质数
			 //质因子就是他自己 
		 	else ans=(ans+(k+1)*num[i])%mod;
		 }
		 printf("%I64d\n",ans);
	 }
	return 0;
}