区间内素数的个数
来源:互联网 发布:大学有linux课程吗 编辑:程序博客网 时间:2024/05/17 01:19
给定整数a和b,请问区间[a,b)内有多少个素数?
限制条件:①a<b<=10^12, ②b-a<=10^6
输入:22 37
输出:3(23、23、31共3个素数)
输入:a=22801763489, b=22801787297
输出:1000
分析:
b以内的合数的最小质因数一定不超过√b。如果有√b以内的素数表的话,就可以把埃氏筛法运用在[a,b)上了。也就是说,先分别做好[2,√b)的表和[a,b)的表,然后从[2,√b)的表中筛得素数的同时,也将其倍数从[a,b)的表中划去,最后剩下的就是区间[a,b)内的素数了。
#include<stdio.h>#include<algorithm>using namespace std;typedef long long ll;bool is_prime[1000010];bool is_prime_small[1000010];void segment_sieve(ll a, ll b){ for(int i=0; (ll)i*i<b; i++) is_prime_small[i] = true;//i是素数 for(int i=0; i<b-a; i++) is_prime[i] = true; //利用0~len代表a~b的数 for(int i=2; (ll)i*i<b; i++) { if(is_prime_small[i]) { for(int j=i+i; (ll)j*j<b; j+=i) is_prime_small[j] = false; //筛[2,√b) for(ll j=max(2LL, (a+i-1)/i)*i; j<b; j+=i) is_prime[j-a] = false; //筛[a,b) //j代表素数,j-a是将a~b变为0~b-a以便数组好存储 //2LL是2的长整形形式,与其比较意思是j最少是i的两倍 //((a+i-1)/i)*i得出的是(>=a && %i==0)离a最近的数,其实 //也可以写成a%i==0 ? a : (a/i+1)*i } }}int main(){ ll a, b; scanf("%lld%lld", &a,&b); segment_sieve(a,b); int ans = 0; for(ll i=0; i<b-a; i++) if(is_prime[i]) printf("%lld ", i+a), ans++; printf("\nans=%d\n", ans); return 0;}
0 0
- 区间内素数的个数
- 求一段区间内素数的个数
- 判断区间内的素数个数
- 判断区间内素数个数
- 筛选素数,打表,区间内特定元素的个数
- 《挑战程序设计竞赛》 大区间内素数的个数
- LightOJ 1197 Help Hanzo 求区间内素数的个数
- 每日一题 No.47 区间内素数的个数
- 循环三 统计给定区间内素数的个数
- [tips]筛选特定区间内素数个数
- 求大区间内素数的个数(区间筛法)
- 求区间内的素数
- 统计给定整数M和N区间内素数的个数
- LightOJ 1197(大数范围筛区间内素数个数)
- Help Hanzo 区间的素数的个数
- 求解1-10000区间内的素数
- 求解1-10000区间内的素数
- 计算一段区间内的所有素数
- jQuery Easing 动画效果扩展
- SQL Sever资源整合
- POJ 2586 Y2K Accounting Bug
- java 获取字体宽度
- struct ifreq 结构体
- 区间内素数的个数
- PHP +fopen 显示图片
- Json系列之五 揭开JsonConfig的神秘面纱 json to java
- Maven依赖
- C# 生成语音文件
- mysql安装和使用
- JAVA知识点总结(八)calendar类及Math-random类
- Android4.4 wpa_supplicant深入分析之初识wpa_supplicant
- 多线程开发——线程安全