Project-Euler/Problem37.cpp at master · krishnaanaril/Project-Euler · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
/*
    Krishna Mohan
    Project Euler - Problem 37
    Date: 3/22/2015
*/
#include <bits/stdc++.h>

using namespace std;

#define fo(i, n) for(int (i)=0;(i)<(n); (i)++)
#define MAX 1000001

bitset<MAX> isPrime;

void init()
{
            isPrime.set();
            isPrime[0]=isPrime[1]=0;
            for(int i=4;i<MAX;i+=2)
            {
                        isPrime[i]=0;
            }
            for(int i=3;i<MAX;i++)
            {
                        if(isPrime[i])
                        {
                                    for(int j=2*i;j<MAX;j+=i)
                                    {
                                                isPrime[j]=0;
                                    }
                        }
            }
}

int main() {
            init();
            int num[25], len, n;
            long long sum=0;
            for(int i=11;i<MAX;i++)
            {
                        len=0;
                        n=i;
                        while(n)
                        {
                                    num[len++]=n%10;
                                    n/=10;
                        }
                        /*if(i==51)
                        {
                                    for(int j=0;j<len;j++)
                                                cout<<num[j];
                                    cout<<endl;
                        }*/
                        int f=0;
                        for(int j=len-1;j>=0;j--)
                        {
                                    n=0;
                                    for(int k=j;k>=0;k--)
                                    {
                                                n=n*10+num[k];
                                    }
                                    /*if(i==51)
                                                cout<<"Digs: "<<n<<endl;*/
                                    if(!isPrime[n])
                                    {
                                                f=1;
                                                break;
                                    }
                        }
                        int l=0;
                        for(int j=0;j<len;j++)
                        {
                                    n=0;
                                    for(int k=len-1;k>=j;k--)
                                    {
                                                n=n*10+num[k];
                                    }
                                    if(!isPrime[n])
                                    {
                                                l=1;
                                                break;
                                    }
                        }
                        if(!f && !l)
                        {
                                    sum+=i;
                                    //cout<<i<<endl;
                        }
            }
            cout<<"Ans: "<<sum<<endl;
            return 0;
}