问题描述
幸运数是波兰数学家乌拉姆命名的。它采用与生成素数类似的“筛法”生成
。首先从1开始写出自然数1,2,3,4,5,6,….
1 就是第一个幸运数。
我们从2这个数开始。把所有序号能被2整除的项删除,变为:
1 _ 3 _ 5 _ 7 _ 9 ….
把它们缩紧,重新记序,为:
1 3 5 7 9 …. 。这时,3为第2个幸运数,然后把所有能被3整除的序号位置的数删去。注意,是序号位置,不是那个数本身能否被3整除!! 删除的应该是5,11, 17, …
此时7为第3个幸运数,然后再删去序号位置能被7整除的(19,39,…)
最后剩下的序列类似:
1, 3, 7, 9, 13, 15, 21, 25, 31, 33, 37, 43, 49, 51, 63, 67, 69, 73, 75, 79, …
输入格式
输入两个正整数m n, 用空格分开 (m < n < 1000*1000)
输出格式
程序输出 位于m和n之间的幸运数的个数(不包含m和n)。
样例输入1
1 20
样例输出1
5
样例输入2
30 69
样例输出2
8
只在范围内查找,每次变化查询范围:
代码如下:
- #include <iostream>
- #include <cstring>
- using namespace std;
- #define MAX 1000010
- int a[MAX];
- int filter(int s,int e)
- {
- int k=s;
- int t=a[s];
- for (int i=s;i<e;i++)
- {
- if (i%t)
- {
- a[k]=a[i];
- k++;
- }
- }
- if (t<e)
- return filter(s+1,k);
- else return k;
- }
- int main()
- {
- int s,e,sum;
- while(cin>>s>>e){
- sum=0;
- memset(a,0,sizeof(a));
- int temp=e/2+1;
- for (int i=1;i<temp;i++){
- a[i]=2*i-1;
- }
- int k=filter(2,temp);
- for (int i=1;i<k;i++)
- {
- if (a[i]>s) sum++;
- else if (a[i]>=e)break;
- }
- cout<<sum<<endl;
- }
- return 0;
- }
第二种实现:
- #include<iostream>
- #include<memory.h>
- using namespace std;
- int a[500005];
- int main()
- {
- int i,k=1,temp;
- int n,m,sum=0;
- cin>>n>>m;
- for(i=1;i<=1000000;i+=2)//除去2之后
- a[k++]=i;
- temp=2;
- while(temp<=1000)
- {
- int b[500005],t=1;
- for(i=1;i<k;i++)
- {
- if(i%a[temp]!=0)
- b[t++]=a[i];
- }
- for(i=1;i<t;i++)
- a[i]=b[i];
- temp++;
- k=t;
- }
- for(i=1;i<=m;i++)
- {
- if(a[i]>n&&a[i]<m)
- sum++;
- }
- cout<<sum<<endl;
- }