Cod sursa(job #19226)

Utilizator victorsbVictor Rusu victorsb Data 18 februarie 2007 22:33:18
Problema Ghiozdan Scor 68
Compilator cpp Status done
Runda Arhiva de probleme Marime 4.12 kb
#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;
}