Cod sursa(job #1829533)

Utilizator AlexNiuclaeNiculae Alexandru Vlad AlexNiuclae Data 15 decembrie 2016 09:45:40
Problema Ghiozdan Scor 100
Compilator cpp Status done
Runda Lista lui wefgef Marime 2.28 kb
#include <bits/stdc++.h>

using namespace std;

const int v_max = 2e2;
const int vmax = 2e2 + 10;
const int gmax = 75e3 + 10;

const int inf = 2e4 + 10;

int n, g, act, prv;
int nr[vmax];
int dp[2][gmax], w_left[2][gmax];

vector < int > ans;
int dq[gmax];

void input() {
    scanf("%d %d", &n, &g);
    for (int i = 1; i <= n; ++i) {
        int x; scanf("%d", &x);
        nr[x]++;
    }
}

int get_dp(int p, int i, int val) {
    int k = (i - p) / val;
    return dp[prv][i-k*val];
}

void update_ans(int i, int val, int p, int mid) {
    int k = (i - p) / val;

    dp[act][i] = inf;
    if (dp[prv][i-k*val] == inf) return;

    dp[act][i] = dp[prv][i-k*val] + k;
    w_left[act][i] = (val <= mid) ? i : w_left[prv][i-k*val];
}

void compute_dp(int left, int right, int g) {
    int mid = (left + right) >> 1;

    act = 0;

    for (int i = 0; i <= g; ++i) dp[act][i] = inf;
    dp[act][0] = 0;

    for (int val = left; val <= right; ++val) { //gap = val
        if (nr[val] == 0)
            continue;

        prv = act; act = prv ^ 1;
        for (int shift = 0; shift < val; ++shift) {
            int l = 1, r = 0;
            for (int i = shift; i <= g; i += val) {
                if (l <= r && (i - dq[l]) / val == nr[val] + 1)
                    l++;
                while (l <= r && dp[prv][i] <= get_dp(dq[r], i, val))
                    r--;
                dq[++r] = i;

                if (l <= r)
                    update_ans(i, val, dq[l], mid);
            }
        }
    }
}

void compute_sol(int left, int right, int g_act) {
    if (g_act == 0) return;
    if (left == right) {
        for (int i = 1; i <= g_act / left; ++i)
            printf("%d\n", left);
        return;
    }

    compute_dp(left, right, g_act);

    int mid = (left + right) >> 1;
    int w_act; for (w_act = g_act; dp[act][w_act] == inf && w_act; w_act--);

    if (right - left == v_max - 1)
        printf("%d %d\n", w_act, dp[act][w_act]);

    int w_leftACT = w_left[act][w_act];
    compute_sol(left, mid, w_leftACT);
    compute_sol(mid + 1, right, w_act - w_leftACT);
}

int main()
{
    freopen("ghiozdan.in","r",stdin);
    freopen("ghiozdan.out","w",stdout);

    input();
    compute_sol(1, v_max, g);

    return 0;
}