Pagini recente » Cod sursa (job #1967967) | Cod sursa (job #1306891) | Cod sursa (job #1471018) | Cod sursa (job #22576) | Cod sursa (job #19226)
Cod sursa(job #19226)
#include <cstdio>
#include <string.h>
#include <algorithm>
using namespace std;
#define inf 1000000000
#define rad 14
#define Nmax 20001
#define Gmax 75001
#define Vmax 201
int n, g, a, smax;
int ct[Vmax];
int sol[21][Gmax];
int sol2[11][Gmax];
int nr[21][Gmax];
int Q[Gmax];
inline void citire()
{
int i;
scanf("%d %d\n", &n, &g);
for (i = 1; i <= n; ++i)
{
scanf("%d\n", &a);
++ct[a];
}
}
void solve()
{
int i, j, k, l, st, dr, pas, gr = 1, num;
memset(sol, 0x3f, sizeof(sol));
memset(sol2, 0x3f, sizeof(sol2));
sol[0][0] = 0;
sol2[0][0] = 0;
for (pas = 0; pas <= 199; pas += 20, ++gr)
{
for (i = 1; i <= 20; ++i)
{
num = pas + i;
if (ct[num])
{
for (j = 0; j < num; ++j)
{
st = 1;
dr = 0;
if (j == 0)
Q[++dr] = 0;
for (k = j; k <= g; k += num)
{
if (st <= dr)
while (k - Q[st] > ct[num] * num)
++st;
sol[i][k] = sol[i-1][k];
if (st <= dr)
if (sol[i-1][Q[st]] + ((k - Q[st]) / num) < sol[i][k])
sol[i][k] = sol[i-1][Q[st]] + ((k - Q[st]) / num);
while (st <= dr && sol[i - 1][Q[dr]] > sol[i - 1][k])
--dr;
Q[++dr] = k;
}
}
}
for (j = 0; j <= g; ++j)
sol[i][j] = min(sol[i][j], sol[i-1][j]);
}
for (j = 0; j <= g; ++j)
sol2[gr][j] = sol[20][j];
memset(sol, 0x3f, sizeof(sol));
for (j = 0; j <= g; ++j)
sol[0][j] = sol2[gr][j];
}
--gr;
for (smax = g; smax >= 0; --smax)
if (sol2[gr][smax] < inf)
break;
printf("%d %d\n", smax, sol2[gr][smax]);
for (pas = 180; pas >= 0; pas -= 20, --gr)
{
//memset(sol, 0x3f, sizeof(sol));
//memset(nr, 0, sizeof(nr));
for (j = 0; j <= g; ++j)
sol[0][j] = sol2[gr - 1][j];
for (i = 1; i <= 20; ++i)
{
num = pas + i;
if (ct[num])
{
for (j = 0; j < num; ++j)
{
st = 1;
dr = 0;
if (j == 0)
Q[++dr] = 0;
for (k = j; k <= g; k += num)
{
if (st <= dr)
while (k - Q[st] > ct[num] * num)
++st;
sol[i][k] = sol[i-1][k];
nr[i][k] = 0;
if (st <= dr)
if (sol[i-1][Q[st]] + ((k - Q[st]) / num) < sol[i][k])
{
sol[i][k] = sol[i-1][Q[st]] + ((k - Q[st]) / num);
nr[i][k] = (k - Q[st]) / num;
}
while (st <= dr && sol[i - 1][Q[dr]] > sol[i - 1][k])
--dr;
Q[++dr] = k;
}
}
}
for (j = 0; j <= g; ++j)
sol[i][j] = min(sol[i][j], sol[i-1][j]);
}
for (i = 20; i > 0; --i)
if (nr[i][smax])
{
for (j = 1; j <= nr[i][smax]; ++j)
printf("%d\n", i + pas);
smax -= nr[i][smax] * (pas + i);
}
}
}
int main()
{
freopen("ghiozdan.in", "r", stdin);
freopen("ghiozdan.out", "w", stdout);
citire();
solve();
return 0;
}