幸运数是波兰数学家乌拉姆命名的。它采用与生成素数类似的“筛法”生成。
首先从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, ...
1 #include <iostream> 2 #include <cstdio> 3 #include <cstdlib> 4 #include <cstring> 5 #include <string> 6 #include <cmath> 7 #include <algorithm> 8 #define INF 0x3f3f3f3f 9 #define zero 1e-7 10 11 using namespace std; 12 typedef long long ll; 13 const ll mod=1e9+7; 14 const ll max_n=1e6+7; 15 16 int a[max_n]={0}; 17 18 int main() { 19 int m, n, ans=0; 20 cin>>m>>n; 21 int sum=n/2+1;//排除2的倍数之后还剩少于（n/2+1）个数 22 for(int i=1; i<sum; i++) a[i]=i*2-1; 23 for(int i=2; i<sum; i++) { 24 int cnt=i; 25 for(int j=i; j<sum; j++) { 26 if(j%a[i]!=0) a[cnt++]=a[j]; 27 } 28 sum=cnt;//更新sum的值为 排除幸运数a[i]的倍数下标对应的数后剩余的个数 29 } 30 for(int i=1; i<sum; i++) { 31 if(a[i]>=n) break; 32 if(a[i]>m && a[i]<n) ans++; 33 } 34 cout<<ans<<endl; 35 return 0; 36 }
原文：https://www.cnblogs.com/wwqzbl/p/13571792.html