Prime k-tuple UVA - 1404

就是大区间求素数  参考

LightOJ - 1197

https://www.cnblogs.com/WTSRUVF/p/9190660.html

直接套那个代码就好了

#include <iostream>
#include <cstdio>
#include <sstream>
#include <cstring>
#include <map>
#include <set>
#include <vector>
#include <stack>
#include <queue>
#include <algorithm>
#include <cmath>
#define MOD 2018
#define LL long long
#define ULL unsigned long long
#define Pair pair<int, int>
#define mem(a, b) memset(a, b, sizeof(a))
#define _  ios_base::sync_with_stdio(0),cin.tie(0)
using namespace std;
const int maxn =  * 1e9;
int primes[];
bool vis[maxn+];
int ans = ;
vector<int> G;

void init()
{
    mem(vis,);
    vis[] = ;
    for(int i=; i<=; i++)
        if(!vis[i])
        {
            primes[ans++] = i;
            for(LL j=(LL)i*i; j<=; j+=i)
                vis[j] = ;
        }
}

int main()
{
    init();
    int T;
    int kase = ;
    LL a, b, u, s;
    cin>> T;
    while(T--)
    {
        G.clear();
        int res = ;
        mem(vis,);
        cin>> a >> b >> u >>s;
       // if(a <= 2) a = 2;
        int len = b - a;
        for(int i=; i<ans && primes[i] * primes[i] <= b; i++)
        {
            int j = a/primes[i];
            if(j*primes[i] < a) j++;
            if(j <= ) j++;
            while(j*primes[i] <= b)
            {
                vis[j*primes[i] - a] = ;
                j++;
            }

        }
        if(a == ) vis[] = ;
        for(int k=; k<=len; k++)
            if(!vis[k] && a + k > )
                G.push_back(a+k);
        for(int i=;  i + u - <G.size(); i++)
            if(G[i+u-] - G[i] == s)
                res++;
        printf("%d\n",res);

    }

    return ;
}

我觉着人家写的代码比较好  参考一下

#include <cstdio>
#include <cstring>
#include <cmath>
#include <vector>
#include <algorithm>

using namespace std;
typedef long long ll;
const int sqrt_inf = ;
const int maxn =  * 1e9;

int np, pri[sqrt_inf];
bool vis[maxn+];
vector<int> vec;

void prime_table (int n) {
    np = ;
    memset(vis, , sizeof(vis));

    for (int i = ; i <= n; i++) {
        if (vis[i])
            continue;

        pri[np++] = i;
        for (int j = i * i; j <= n; j += i)
            vis[j] = ;
    }
}

int solve () {
    int ret = ;
    int a, b, s, k;
    vec.clear();
    memset(vis, , sizeof(vis));

    scanf("%d%d%d%d", &a, &b, &k, &s);

    for (int i = ; i < np && pri[i] * pri[i] <= b; i++) {
        int u = pri[i], d = (u - a % u) % u;

        if (u == a + d)
            d += u;

        while (d <= b - a) {
            vis[d] = ;
            d += u;
        }
    }

    for (int i = ; i <= b-a; i++) {
        if (vis[i] ==  && a + i > )
            vec.push_back(a+i);
    }

    for (int i = ; i + k -  < vec.size(); i++) {
        if (vec[i+k-] - vec[i] == s)
            ret++;
    }

    return ret;
}

int main () {
    prime_table(sqrt_inf);

    int cas;
    scanf("%d", &cas);

    while (cas--) {
        printf("%d\n", solve());
    }
    return ;
}