试题 历届试题 幸运数(二分)

试题 历届试题 幸运数

资源限制
时间限制：1.0s 内存限制：256.0MB

---

题描述
幸运数是波兰数学家乌拉姆命名的。它采用与生成素数类似的“筛法”生成

。
首先从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>
const int N = 1e6+5;
using namespace std;
int a[N];
bool vis[N];
//预先筛选
int solve()
{
    int k = 0;
    for(int i = 1; i < 1e6; i+=2)
    {
        a[++k] = i;
    }
    int id = 1;
    while(1)
    {
        bool flag = false;
        int siz = k;
        for(int i = a[++id]; i <= siz; i+=a[id])
        {
            if(i % a[id] == 0)
            {
                flag = true;
                vis[a[i]] = true;
            }
        }
        if(!flag) return k;
        k = 0;
        for(int i = 1; i <= siz; ++i)
        {
            if(!vis[a[i]]) a[++k] = a[i];
        }
    }
}
//第一个大于v
int getPos1(int k,int v)
{
    int l = 1,r = k;
    int res = 0;
    while(l <= r)
    {
        int mid = (l + r) >> 1;
        if(a[mid] > v)
        {
            res = mid;
            r = mid - 1;
        }
        else l = mid + 1;
    }
    return res;
}

//最后一个小于v
int getPos2(int k,int v)
{
    int l = 1,r = k;
    int res = 0;
    while(l <= r)
    {
        int mid = (l + r) >> 1;
        if(a[mid] < v)
        {
            res = mid;
            l = mid + 1;
        }
        else r = mid - 1;
    }
    return res;
}

int main()
{
    int n,m;
    int k = solve();
    cin>>m>>n;
    int ans = getPos2(k,n)-getPos1(k,m) + 1;
    cout<<ans<<endl;
    return 0;
}