Pagini recente » Cod sursa (job #2160445) | Cod sursa (job #1477699) | Cod sursa (job #2110996) | Cod sursa (job #2240274) | Cod sursa (job #1829533)
#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;
}