forked from sitz/UVa-Online-Judge
-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy path10856.cpp
109 lines (107 loc) · 1.28 KB
/
10856.cpp
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
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
#include <bits/stdc++.h>
using namespace std;
/*
10856
*/
#define mx 3000000
int P[mx + 2], tp;
char Fg[mx + 2];
typedef long long ss;
void Table()
{
ss i, j, k;
for (i = 2; i <= mx;)
{
P[i]++;
for (j = i + i; j <= mx; j += i)
{
Fg[j] = 1;
P[j]++;
}
for (j = i * i; j <= mx; j *= i)
{
P[j]++;
for (k = j + j; k <= mx; k += j)
{
P[k]++;
}
}
for (++i; Fg[i]; i++)
;
}
for (i = 2; i <= mx; i++)
{
P[i] += P[i - 1];
}
}
int Count(int n)
{
int i, p = 0;
if (Fg[n] > 0)
{
return Fg[n];
}
for (i = 0; P[i] <= n && i < tp; i++)
{
for (int j = P[i]; j <= n; j *= P[i])
{
p += n / j;
}
}
Fg[n] = p;
return p;
}
int BinarySearch(int key)
{
int lo = 2, up = mx, mid, fac, lfac, rfac;
mid = (lo + up) / 2;
if (key == 0)
{
return 0;
}
if (key == 1)
{
return 2;
}
while (lo < up)
{
fac = P[mid];
if (fac > key)
{
up = mid - 1;
}
else if (fac < key)
{
lo = mid + 1;
}
else
{
return mid;
}
mid = (lo + up) / 2;
}
fac = P[mid];
if (fac != key)
{
return -1;
}
return mid;
}
int main()
{
int n, k = 1, d;
Table();
while (scanf("%d", &n) && n >= 0)
{
d = BinarySearch(n);
if (d < 0)
{
printf("Case %d: Not possible.\n", k++);
}
else
{
printf("Case %d: %d!\n", k++, d);
}
}
return 0;
}