习题10-4 素数间隔 UVa1644

XDU_Skyline

于 2015-02-26 23:25:15 发布

阅读量1k

点赞数

CC 4.0 BY-SA版权

分类专栏：算法竞赛入门经典（第二版）数学——数论搜索——二分查找文章标签： uva 二分查找

本文链接：https://blog.youkuaiyun.com/u014800748/article/details/43958095

算法竞赛入门经典（第二版）同时被 3 个专栏收录

136 篇文章

订阅专栏

数学——数论

33 篇文章

订阅专栏

搜索——二分查找

10 篇文章

订阅专栏

该博客主要介绍了如何解决UVa在线判题系统中的10-4题——素数间隔问题。博主通过计算1299709以内所有素数，并采用二分查找算法来寻找特定数值n附近的素数位置，以此来求解题目。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

1.题目描述：点击打开链接

2.解题思路：根据题意可知最大的素数在int范围内，可以先算出1299709以内的所有素数，随后二分查找n附近的素数的位置即可。

3.代码：

#define _CRT_SECURE_NO_WARNINGS 
#include<iostream>
#include<algorithm>
#include<string>
#include<sstream>
#include<set>
#include<vector>
#include<stack>
#include<map>
#include<queue>
#include<deque>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<functional>
using namespace std;

#define N 1300000
int vis[N];
vector<int>primes;
void init()
{
	int m = sqrt(N + 0.5);
	for (int i = 2; i <= m;i++)
	if (!vis[i])
	for (int j = i*i; j < N; j += i)
		vis[j] = 1;
	for (int i = 2; i <= N;i++)
	if (!vis[i])
		primes.push_back(i);
}
int main()
{
	//freopen("test.txt", "r", stdin);
	int n;
	init();
	while (scanf("%d", &n) != EOF&&n)
	{
		if (!vis[n])cout << 0 << endl;
		else
		{
			int L = 0, R = 100000;
			while (L < R)
			{
				int m = L + (R - L) / 2;
				if (primes[m] > n)R = m;
				else L = m + 1;
			}
			cout << primes[L] - primes[L - 1] << endl;
		}
	}
	return 0;
}