欧拉计划37题

Truncatable primes

The number 3797 has an interesting property. Being prime itself, it is possible to continuously remove digits from left to right, and remain prime at each stage: 3797, 797, 97, and 7. Similarly we can work from right to left: 3797, 379, 37, and 3.

Find the sum of the only eleven primes that are both truncatable from left to right and right to left.

NOTE: 2, 3, 5, and 7 are not considered to be truncatable primes.

可截素数

3797有着奇特的性质。它本身是一个素数；如果从左往右逐一截去数字，剩下的仍然都是素数：3797、797、97和7；如果从右往左逐一截去数字，剩下的也仍然都是素数：3797、379、37和3。

如果一个素数满足，无论从左往右还是从右往左逐一截去数字，剩下的仍然都是素数，则称之为可截素数。已知总共有十一个可截素数，求这些数的和。

注意：2、3、5和7不被视为可截素数。

         这个题目主要问题就是如何去从左往右逐一截取数字；其实方法很简单，先求到这个数的位数，然后用10的位数的次饭去求得最高位然后再减去最高位就可以实现从左往右逐一截去数字；而从右往左就很简单除10就可以了；下面是代码实现：

        

#include <stdio.h>
#include <math.h>
#define MAX_N 1000000

int prime[MAX_N + 5];
int isprime[MAX_N + 5];

void init() {
    prime[1] = prime[0] = 1;
    for (int i = 2; 2 * i <= MAX_N; i++) {
        if (!prime[i]) isprime[++isprime[0]] = i;
        for (int j = 1; i * isprime[j] <= MAX_N; j++) {
            prime[i * isprime[j]] = 1;
            if (i % isprime[j] == 0) break;
        }
    }
    return ;
}

int get_len(int x) {
    return log10(x);
}

int is_Truncatable_prime(int x) {
    int k = x;
    while (k) {//从左往右逐一删
        if (prime[k]) return 0; 
        int len = get_len(k);
        k -= k / (int)pow(10, len) * (int)pow(10, len);//必须强转，不强转会有小数，造成结果不对
    }
    k = x;
    while (k) {//从右往左逐一删
        if (prime[k]) return 0;
        k /= 10;
    }
    return 1;
}

int main() {
    init();
    int cnt = 0, ans = 0;
    for (int i = 10; cnt < 11; i++) {
        if (prime[i] || !is_Truncatable_prime(i)) continue;
        ans += i;
        cnt++;
    }
    printf("%d\n", ans);
    return 0;
}

 最终答案：748317