Codeforces 237C Primes on Interval(素数统计)

C. Primes on Interval
time limit per test
1 second
memory limit per test
256 megabytes
input
standard input
output
standard output

You've decided to carry out a survey in the theory of prime numbers. Let us remind you that a prime number is a positive integer that has exactly two distinct positive integer divisors.

Consider positive integers aa + 1...b (a ≤ b). You want to find the minimum integer l (1 ≤ l ≤ b - a + 1) such that for any integer x(a ≤ x ≤ b - l + 1) among l integers xx + 1...x + l - 1 there are at least k prime numbers.

Find and print the required minimum l. If no value l meets the described limitations, print -1.

Input

A single line contains three space-separated integers a, b, k (1 ≤ a, b, k ≤ 106a ≤ b).

Output

In a single line print a single integer — the required minimum l. If there's no solution, print -1.

Examples
input
2 4 2
output
3
input
6 13 1
output
4
input
1 4 3
output
-1

题意:

在给定的区间  [a, b] 之内,问是否存在一个最小的整数 l ,使得任意一个在区间内(a ≤ x ≤ b - l + 1)x ,至少在 x, x + 1, ..., x + l - 1 k 个素数。

思路:

统计每个数字之前的素数个数,二分即可。

就这题,四级怎么办?!哭

AC CODE:

#include<stdio.h>
#include<cstring>
#include<algorithm>
#define HardBoy main()
#define ForMyLove return 0;
using namespace std;
const int MYDD = 1103+1e6;

bool isprime[MYDD];
int cntisp[MYDD];
void init() {
	memset(isprime, 1, sizeof(isprime));
	isprime[0] = 0, isprime[1] = 0;
	for(int j = 2; j < MYDD; j++) {
		if(isprime[j]) {
			for(int k = 2*j; k < MYDD; k = k+j)
				isprime[k] = 0;
		}
	}
	
	cntisp[0] = 0, cntisp[1] = 0;//当前数字之前素数个数统计
	for(int j = 2; j < MYDD; j++) {
		cntisp[j] = cntisp[j-1];
		if(isprime[j])  cntisp[j]++;
	}
}

int HardBoy {
	init();
	int a, b, k;
	scanf("%d %d %d", &a, &b, &k);

	int l = k, r = b-a+1;
	while(l <= r) {
		int j, mid = (l+r)/2;
		for(j = a; j+mid-1 <= b; j++) {
			if(cntisp[j+mid-1] - cntisp[j-1] < k) break;
		}
		if(j+mid-1 <= b) l = mid+1;
		else r = mid-1;
	}

	if(l > b-a+1) puts("-1");
	else
		printf("%d\n", l);
		
	ForMyLove
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值