http://codeforces.com/contest/569/problem/C
と定義する．比が与えられるので，を満たす最大のnを求めたい．
考察
まず，と，を求める．これをどこまで必要かを自分で判断しなければならない．とあるので，多くでもがの42倍になっているnまででよいと分かる．実際に値を試した所，n = 1500000で十分だとわかった．後は，後ろから見ていき，条件を満たす最大のnを答えた．
Code
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
93
| #include <iostream>
#include <vector>
#include <string>
#include <cstring>
#include <algorithm>
#include <sstream>
#include <map>
#include <set>
#include <cmath>
#define REP(i,k,n) for(int i=k;i<n;i++)
#define rep(i,n) for(int i=0;i<n;i++)
#define INF 1<<30
#define pb push_back
#define mp make_pair
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
bool prime[10000000];
void Eratosthenes(int n) {
rep(i,n) prime[i] = true;
prime[0] = false;
prime[1] = false;
REP(i,2,(int)sqrt(n)) {
if(prime[i]) {
for(int j=0;i*(j+2)<n;j++) {
prime[i*(j+2)] = 0;
}
}
}
}
bool check(string s) {
rep(i,s.size()/2) {
if(s[i] != s[s.size()-1-i]) return false;
}
return true;
}
int main() {
double p,q;
cin >> p >> q;
int N = 1500000;
Eratosthenes(N);
vector<int> pi(N);
int cnt = 0;
rep(i,N+5) {
if(prime[i]) {
cnt++;
}
pi[i] = cnt;
}
vector<int> rub(N);
cnt = 0;
REP(i,1,N+5) {
stringstream ss;
ss << i;
if(check(ss.str())) {
cnt++;
}
rub[i] = cnt;
}
double A = p/q;
int ans = 0;
for(int i=N;i>=1;i--) {
ll a = pi[i];
ll b = rub[i];
if(a <= A * b) {
cout << i << endl;
break;
}
else {
continue;
}
}
return 0;
}
|