打卡信奥刷题（1635）用C++信奥 P7588 双重素数（2021 CoE-II A）-CSDN博客

P7588 双重素数（2021 CoE-II A）

题目描述

素数（质数）是指在大于 $1$ 的自然数中，除了 $1$ 和它本身以外不再有其他因数的自然数。定义双重素数为这样的素数：它的各位数字之和也是一个素数。给定一个闭区间，试确定在该区间内双重素数的个数。

输入格式

输入包含多组测试数据。

输入第一行包含一个整数 $T$ ，表示测试数据的组数。接下来每行一组测试数据，每组测试数据包含以空格分隔的两个整数 $L$ 和 $R$ 。

输出格式

每组测试数据输出一行，包含一个整数，表示在闭区间 $R][L,\ R]$ 内双重素数的个数。

输入输出样例 #1

输入 #1

输出 #1

说明/提示

样例说明

从 $1$ 到 $15$ 共有 $6$ 个素数： $2$ ， $3$ ， $5$ ， $7$ ， $11$ ， $13$ 。前五个素数各自的数字之和也是素数，因此都是双重素数。素数 $13$ 的各位数字之和为 $4$ ，不是素数，故 $13$ 不是双重素数。

数据范围

Subtask $1$ ： $\le L \le R \le 10^2$ ， $10$ 分。
Subtask $2$ ： $\le L \le R \le 10^4$ ， $20$ 分。
Subtask $3$ ： $\le L \le R \le 10^6$ ， $60$ 分。
Subtask $4$ ： $\le L \le R \le 10^8$ ， $10$ 分。

对于 $100%100\%$ 的数据， $\le T \le 100$ 。

提示（数据已经加强）

最后一个子任务要求你的程序必须具有较高的空间使用效率和时间效率，否则容易超出内存限制或时间限制。

C++实现

#include <bits/stdc++.h>
using namespace std;
#define GET(x) (B[x >> 5] & (1 << (x & 0x1F)))
#define SET(x) (B[x >> 5] |= (1 << (x & 0x1F)))
const int MAXB = 100000001, MAXN = 6000000;
int B[MAXB >> 5], dp[MAXN], pc = 0, dpc = 0;
int magic[25] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97};
inline bool isDoublePrime(int n){
    int ds = 0;
    while (n) ds += (n % 10), n /= 10;
    for (int i = 0; i < 25; i++)
        if (ds == magic[i])
            return true;
    return false;
}
int main(int argc, char *argv[]){
    for (int i = 2; i < MAXB; i++)    {
        if (!GET(i)) dp[++pc] = i;
        for (int j = 1; j <= pc && i * dp[j] < MAXB; j++)        {
            SET(i * dp[j]);
            if (i % dp[j] == 0) break;
        }
    }
    for (int i = 1; i <= pc; i++) if (isDoublePrime(dp[i])) dp[++dpc] = dp[i];
    int cases;
    cin >> cases;
    for (int cs = 1; cs <= cases; cs++)    {
        int L, R;
        cin >> L >> R;
        L = upper_bound(dp + 1, dp + dpc + 1, L - 1) - dp - 1;
        R = upper_bound(dp + 1, dp + dpc + 1, R) - dp - 1;
        cout << R - L << '\n';
    }
    return 0;
}