首页 > 其他 > 详细

[luoguP1835] 素数密度_NOI导刊2011提高(04)(素数筛)

时间:2017-08-14 19:29:51      阅读:139      评论:0      收藏:0      [点我收藏+]

传送门

 

数据辣么大,怎么搞?(L≤R≤2147483647)

注意到R-L≤1000000

所以可以直接筛R-L区间内的数,

但是需要用已知的小的素数筛,

R-L区间内的大部分数肯定能用较小的素数筛去,但是还有一些较大的数,可能等于两个大质数的乘积,没法被筛去。

但是又注意到,数据最大才10位,也就是说我们只需要用位数<=5的素数筛就可以了,所以先预处理出来,直接筛就ok了。

 

#include <cstdio>
#define N 1000001
#define max(x, y) ((x) > (y) ? (x) : (y))

int L, R, cnt, ans;
int prime[N];
bool notprime[N], ansprime[N];

int main()
{
	int i, j;
	scanf("%d %d", &L, &R);
	notprime[0] = notprime[1] = 1;
	for(i = 2; i <= 100000; i++)
	{
		if(!notprime[i]) prime[++cnt] = i;
		for(j = i + i; j <= 100000; j += i) notprime[j] = 1;
	}
	for(i = 1; i <= cnt; i++)
		for(j = max(2, L / prime[i]); j <= R / prime[i]; j++)
			ansprime[j * prime[i] - L] = 1;
	for(i = 0; i <= R - L; i++)
		if(!ansprime[i])
			ans++;
	printf("%d\n", ans);
	return 0;
} 

  

[luoguP1835] 素数密度_NOI导刊2011提高(04)(素数筛)

原文:http://www.cnblogs.com/zhenghaotian/p/7359637.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!