POJ-2689 Prime Distance(区间素数筛--经典题)

23 篇文章 0 订阅
19 篇文章 2 订阅

题意:

给出一个区间[L, R],求该区间内的所有素数(1 <= L,R <= 2.1E9, R-L <= 1e6)

思路:

区间素数筛,枚举不超过sqrt(R)的所有质数p,然后再去枚举[L, R]区间内p的倍数,将之划去,最后区间内剩下的就都是素数了。注意一点,1需要特判不是素数。


代码:

#include <algorithm>
#include <string.h>
#include <cstdio>
#define LL long long
using namespace std;
const int N = 1e5;
const int maxn = 1e6+5;
int isprime[maxn], prime[maxn], len;
LL L, R;
int ans[maxn];
void prime_init()
{
	len = 0;
	memset(isprime, 0, sizeof isprime);
	for(int i = 2; i <= N; ++i)
	if(!isprime[i])
	{
		prime[++len] = i;
		if(1ll*i*i > N) continue;
		for(int j = i*i; j <= N; j+=i)
		isprime[j] = 1;
	}
}
void work()
{
	memset(ans, 0, sizeof ans);
	for(int i = 1; i <= len; ++i)
	{
		if(1ll*prime[i]*prime[i] > R) break;
		LL star = ((L-1)/prime[i]+1)*prime[i];
		for(LL j = star; j <= R; j+=prime[i])
		if(j != prime[i]) ans[j-L+1] = 1;
	}
	if(L == 1) ans[1] = 1;
}
int main()
{
	prime_init();
	while(~scanf("%lld %lld", &L, &R))
	{
		work();
		LL i, mins, maxx, ans1, ans2, ans3, ans4;
		ans1 = ans2 = -1, i = 0;
		while(i < R-L+1 && (ans1 == -1 || ans2 == -1))
		{
			++i;
			if(ans[i]) continue;
			if(ans1 == -1) ans1 = i+L-1;
			else if(ans2 == -1) ans2 = i+L-1;
		}
		if(ans1 == -1 || ans2 == -1) puts("There are no adjacent primes.");
		else
		{
			mins = maxx = ans2-ans1;
			ans3 = ans1, ans4 = ans2;
			int pre = ans2;
			for(++i; i <= R-L+1; ++i)
			{
				if(ans[i]) continue;
				if(i+L-1-pre < mins)
				{
					mins = i+L-1-pre;
					ans1 = pre, ans2 = i+L-1;
				}
				if(i+L-1-pre > maxx)
				{
					maxx = i+L-1-pre;
					ans3 = pre, ans4 = i+L-1;
				}
				pre = i+L-1;
			}
			printf("%lld,%lld are closest, %lld,%lld are most distant.\n", ans1, ans2, ans3, ans4);
		}
	}
	return 0;
}


继续加油~

POJ - 3616是一个目,目描述如下: 给定一组区间,每个区间有一个权重,要求选择一些区间,使得这些区间的右端点都小于等于k,并且权重之和最大。请问最大的权重和是多少? 解决这个问的思路是使用动态规划。首先,将区间按照左端点从小到大进行排序。然后,定义一个dp数组,dp[i]表示右端点小于等于i的所有区间所能得到的最大权重。 接下来,遍历每一个区间,对于每个区间i,将dp[i]初始化为区间i的权重。然后,再遍历i之前的每个区间j,如果区间j的右端点小于等于k,并且区间j的权重加上区间i的权重大于dp[i],则更新dp[i]为dp[j]加上区间i的权重。 最后,遍历整个dp数组,找到最大的权重和,即为所求的答案。 下面是具体的代码实现: ```cpp #include <cstdio> #include <cstring> #include <algorithm> using namespace std; struct interval{ int start, end, weight; }; interval intervals[10005]; int dp[10005]; int n, m, k; bool compare(interval a, interval b) { if (a.start == b.start) { return a.end < b.end; } else { return a.start < b.start; } } int main() { while(~scanf("%d %d %d", &n, &m, &k)) { memset(dp, 0, sizeof dp); for (int i = 0; i < m; i++) { scanf("%d %d %d", &intervals[i].start, &intervals[i].end, &intervals[i].weight); } sort(intervals, intervals + m, compare); for (int i = 0; i < m; i++) { dp[i] = intervals[i].weight; for (int j = 0; j < i; j++) { if (intervals[j].end <= k && dp[j] + intervals[i].weight > dp[i]) { dp[i] = dp[j] + intervals[i].weight; } } } int maxWeight = 0; for (int i = 0; i < m; i++) { maxWeight = max(maxWeight, dp[i]); } printf("%d\n", maxWeight); } } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值