-
Notifications
You must be signed in to change notification settings - Fork 59
/
Copy pathcopaci2.cpp
76 lines (69 loc) · 1.81 KB
/
copaci2.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
# include <fstream>
# include <algorithm>
# include <deque>
# define NR 1050
# define inf 999999999
using namespace std;
ifstream f("copaci.in");
ofstream g("copaci.out");
deque <int> d;
int i,j,n,m,ci,cs,mij,maxx,bani,COST,minn,F;
int C[NR], H[NR], a[NR], b[NR];
void init ()
{
for (int i=0; i<=maxx; ++i)
a[i]=inf;
}
int main ()
{
f>>n>>bani;
for (i=1; i<=n; ++i)
{
f>>H[i]>>C[i];
cs=max(cs, H[i]);
}
ci=0; cs=cs; maxx=cs;
while (ci<=cs)
{
mij=(ci + cs)/2;
//init
init();
for (j=0; j<=maxx; ++j)
a[j]=abs(H[1]-j)*C[1];
//facem costurile
for (i=2; i<=n; ++i) //fiecare copac si incepem de la inaltimea 0
{
//facem deque-ul
for (j=0; j<=mij; ++j)
if (a[j]!=inf)
{
while (! d.empty() && a[j]<=a[d.back()])
d.pop_back();
d.push_back(j);
}
for (j=0; j<=maxx; ++j) //fiecare inaltime
{
if (d.empty()) b[j]=inf;
else b[j]=abs(H[i]-j)*C[i]+a[d.front()];
if (! d.empty() && d.front()==j-mij) d.pop_front();
if (j+mij+1<=maxx && a[j+mij+1]!=inf)
{
while (! d.empty() && a[j+mij+1]<=a[d.back()])
d.pop_back();
d.push_back(j+mij+1);
}
}
for (j=0; j<=maxx; ++j)
a[j]=b[j], b[j]=inf;
while (! d.empty()) d.pop_back();
}
COST=inf;
for (i=0; i<=maxx; ++i)
COST=min(COST, a[i]);
if (COST<=bani) cs=mij-1, minn=mij, F=COST;
else ci=mij+1;
}
g<<minn<<"\n";
//g<<F<<"\n";
return 0;
}