Description
幸运数是波兰数学家乌拉姆命名的。它采用与生成素数类似的“筛法”生成。
首先从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, ...
Input
输入两个正整数m n, 用空格分开 (m < n < 1000*1000)
Output
程序输出 位于m和n之间的幸运数的个数(不包含m和n)。
Sample Input 1
1 20
Sample Output 1
5
Hint
HINT:时间限制:1.0s 内存限制:256.0MB
思路:
我们只需要找n到m范围之内的幸运数就行了,也就是意味着大于m的我们可以不用管,我们把当前的被除数当作op,如果op的值比m还要大的话,我们就没有必要往下找了,因为这样他删除的都是大于m的数,跟我们所需要求的并没有关系,就可以不再往下循环。
代码:
#include <stdio.h>
int a[1000001];
int main(){
int m,n,mm=500000;//mm控制数组范围
int i,count=0;
scanf("%d%d",&m,&n);
for(int i=0;i<=500000;i++)//首先给整个数组赋值,我们知道第一次除以2也就是把偶数除去;
a[i]=i*2-1;//只给奇数赋值
int op=3;
int k;//下次循环的数组大小
int flag=2;//表示的是当前的下标
while(op<=n){//op为当前的被除数,op到n就可以截止了
flag++;
k=2;
for(int j=2;j<=mm;j++){//循环整个数组
if(j%op!=0){
a[k++]=a[j];//刷新数组
}
}
mm=k;
op=a[flag];
}
for(int i=2;i<=k;i++){
if(a[i]>m&&a[i]<n){
count++;
}
if(a[i]>=n)
break;
}
printf("%d",count);
return 0;
}