工作室课题第二次

找素数

Description

给定区间[L, R] ， 请计算区间中素数的个数。

Input

两个数L和R。

Output

一行，区间中素数的个数。

Sample Input 1

2 11

Sample Output 1

5

Hint

HINT:时间限制：1.0s 内存限制：256.0MB

2 <= L <= R <= 2147483647 R-L <= 1000000

对于此题的理解：

   我首先想到的是写一个带有一个参数函数来判断这个参数是否是素数，然后在主函数里在for循环里 调用这个函数，如果函数返回值为真则sum++，最后输出这个sum ,就是我们想要找的答案。下面是我第一次的代码

#include<stdio.h>
#include<math.h> 
int fun(double n);
int main( void){
	double L,R;
	double c,a;
	int sum=0;
	scanf("%lf" "%lf",&L,&R);
	a=R-L;
	if(L>=2&&R>=L&&L<=2147483674&&a<=10000){
		for(c=L;c<=R;c++) 
	{
		fun(c);
		if (fun(c)==1)
		sum++;
	}
	printf("%d",sum);
	return 0;
	}
	else
	return 0;
}
int fun(double n) 
{
	double i;
	for(i=2;i<n;i++){
		if(fmod(n,i)==0)
		break;
	}
	if(i>=n&&i!=1)
	return 1;
	else 
	return 0;
}

在我运行的时候是正确的，调试的数据也是对的，但是我在平台提交的时候我只对了一半。反思后，我觉得这是我把问题想简单了，忽略了算法复杂度的问题于是，我学习了筛选法求素数，把从1开始的、某一范围内的正整数从小到大顺序排列， 1不是素数，首先把它筛掉。剩下的数中选择最小的 数是素数，然后去掉它的倍数。依次类推，直到“筛子”为空时结束。

学习的代码（参考的代码）

#include<bits/stdc++.h>
using namespace std;
const int N = 1e6+10;
typedef long long ll;

int l, r, ans, cnt;
int vis[N] = {1, 1}, prime[N], st[N];

int main() {
 cin>>l>>r;
 
 int n = sqrt(r);
 
 // 欧拉筛 
 for(int i = 2;i <= n;i ++) {
  if(!vis[i]) prime[++cnt] = i;
  for(int j = 1;j <= cnt;j ++) {
   if(i*prime[j] > n) break;
   vis[i*prime[j]] = 1;
   if(i%prime[j] == 0) break;
  }
 }
 
 // 埃式筛 
/* 
 for(int i = 1;i <= n;i ++) {
  if(!vis[i]) {
   prime[++cnt] = i;
   for(int j = i+i;j <= n;j += i) vis[j] = 1;
  }
 }
*/
 
 for(int i = 1;i <= cnt;i ++) {
  for(ll j = (prime[i] + l - 1)/prime[i]*prime[i];j <= r;j += prime[i]) {
   if(j == prime[i]) continue;
   else st[j-l] = 1;
  }
 }
 for(int i = 0;i <= r-l;i ++) ans += !st[i];

 cout << ans << endl; 
}