标签：prime cnt const int ++ 3006 POJ primeQ

Code

#include <cstdio>
#include <algorithm>
using namespace std;

const int N = 1000001;
bool primeQ[N];
int prime[N / 10];
int cnt;
int a, d, n, cur;

inline int pos(int x){
    // 查询x在素数表中的位置，不存在时返回-1
    int t = lower_bound(prime + 1, prime + 1 + cnt, x) - prime;
    return prime[t] == x?t:-1;
}

int main(int argc, char const *argv[]){
    for(int i = 0; i < N; ++i)primeQ[i] = true;
    for(int i = 2; i < N; ++i){
        if(primeQ[i])prime[++cnt] = i;
        for(int j = 1; j <= cnt && prime[j] * i <= N; ++j){
            primeQ[i * prime[j]] = false;
            if(i % prime[j] == 0)break;
        }
    }
    scanf("%d%d%d", &a, &d, &n);
    while(!(a == 0 && d == 0 && n == 0)){
        cur = 0;
        for(int i = a; cur <= n; i += d){
            if(pos(i) != -1)++cur;
            if(cur == n){
                printf("%d\n", i);
                break;
            }
        }
        scanf("%d%d%d", &a, &d, &n);
    }
    return 0;
}

Review

• 简单考察一下线性筛
• 不用手写二分太好了
• 水题

标签：prime,cnt,const,int,++,3006,POJ,primeQ
来源： https://www.cnblogs.com/mojibake/p/15310129.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。