Cod sursa(job #18876)

Utilizator victorsbVictor Rusu victorsb Data 18 februarie 2007 14:38:15
Problema Ghiozdan Scor 10
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.3 kb
#include <cstdio>
#include <deque>

using namespace std;

#define inf 1000000000
#define Nmax 20001
#define Gmax 75001
#define Vmax 40
#define x first
#define pos second
#define mp make_pair
#define sz size()

int n, g, a, smax, s, maxv, size;
int ct[Vmax];
int sol[Vmax][Gmax];
char nr[Vmax][Gmax];
int Q[Gmax];

void citire()
{
    int i;
    scanf("%d %d\n", &n, &g);
    for (i = 1; i <= n; ++i)
    {
        scanf("%d\n", &a);
        ++ct[a];
        if (a > maxv)
            maxv = a;
    }
}

void scrie(int num, int suma)
{
    if (num == 0)
        return;
    ct[num] = nr[num][suma];
    suma -= num * ct[num];
    scrie(num - 1, suma);
}

void solve()
{
    int i, j, k, l, st, dr;
    memset(sol, 0x3f, sizeof(sol));
    sol[0][0] = 0;
    for (i = 1; i <= maxv; ++i)
    {
    if (ct[i])
    {
        for (j = 0; j < i; ++j)
        {
            st = 1;
            dr = 0;
            if (j == 0)
                Q[++dr] = 0;
                
            for (k = j; k <= g; k += i)
            {
                if (st <= dr)
                while (k - Q[st] > ct[i] * i)
                    ++st;
                sol[i][k] = sol[i-1][k];
                nr[i][k] = 0;
                
                if (st <= dr)
                    if (sol[i-1][Q[st]] + ((k - Q[st]) / i) < sol[i][k])
                    {
                        sol[i][k] = sol[i-1][Q[st]] + ((k - Q[st]) / i);
                        nr[i][k] = (k - Q[st]) / i;
                    }
                
                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 (smax = g; smax >= 0; --smax)
        if (sol[maxv][smax] < inf)
            break;

    memset(ct, 0, sizeof(ct));
    scrie(maxv, smax);

    for (i = 1; i <= Vmax; ++i)
        s += ct[i];
    printf("%d %d\n", smax, s);
    for (i = 1; i <= Vmax; ++i)
        if (ct[i])
            for (j = 1; j <= ct[i]; ++j)
                printf("%d\n", i);
}


int main()
{
    freopen("ghiozdan.in", "r", stdin);
    freopen("ghiozdan.out", "w", stdout);
    citire();
    solve();
    return 0;
}