p205.ans1
======================
5
1 1 3

=================
p205.cpp
======================
#include <bits/stdc++.h>

using namespace std;

template<typename T1, typename T2>
ostream& operator<<(ostream& out, const pair<T1, T2>& x) {
    return out << x.first << ' ' << x.second;
}

template<typename T1, typename T2>
istream& operator>>(istream& in, pair<T1, T2>& x) {
    return in >> x.first >> x.second;
}

template<typename T>
istream& operator>>(istream& in, vector<T>& a) {
    for(auto& x: a) {
        in >> x;
    }
    return in;
};

template<typename T>
ostream& operator<<(ostream& out, const vector<T>& a) {
    for(auto x: a) {
        out << x << ' ';
    }
    return out;
};

const int64_t inf = (int64_t)1e18 + 42;

int n, m, s;
vector<int> x;
vector<vector<int>> L;

void read() {
    cin >> n;
    x.resize(n);
    cin >> x;
    cin >> m >> s;
    L.resize(m, vector<int>(s));
    for(auto& row: L) {
        cin >> row;
    }
}

void solve() {
    // Given n values x[0..n-1], m=2^p level sets each of size s=2^q, where
    // L[i][k] is the k-th level in set i. We must replace each x[j] with some
    // level L[set][k] from the current set, minimizing sum |x[j] - L[set][k]|.
    // The first value uses set 0. The chosen index k's p LSBs
    // determine the next set: next_set = k & (m-1). So picking a suboptimal
    // level now can lead to a better set (and lower total error) later.
    //
    // We solve this with DP backwards. dp[j][set] is the min total deviation
    // for x[j..n-1] when position j uses level set `set`. For each (j, set), we
    // try all s level indices k, pay |x[j] - L[set][k]| and transition to
    // dp[j+1][k & (m-1)]. We track choices for reconstruction. The complexity
    // is O(n * m * s), with n=1000, m=128, s=128 that's ~16M ops, fast enough
    // for the 0.25s time limit.

    vector<vector<int64_t>> dp(n + 1, vector<int64_t>(m, 0));
    vector<vector<int>> choice(n, vector<int>(m, -1));

    for(int j = n - 1; j >= 0; j--) {
        for(int set = 0; set < m; set++) {
            dp[j][set] = inf;
            for(int k = 0; k < s; k++) {
                int next_set = k & (m - 1);
                int64_t cost = abs(x[j] - L[set][k]) + dp[j + 1][next_set];
                if(cost < dp[j][set]) {
                    dp[j][set] = cost;
                    choice[j][set] = k;
                }
            }
        }
    }

    cout << dp[0][0] << '\n';
    int cur_set = 0;
    for(int j = 0; j < n; j++) {
        int k = choice[j][cur_set];
        cout << k << " \n"[j == n - 1];
        cur_set = k & (m - 1);
    }
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    int T = 1;
    // cin >> T;
    for(int test = 1; test <= T; test++) {
        read();
        solve();
    }

    return 0;
}

=================
statement.txt
======================
205. Quantization Problem
time limit per test: 0.25 sec.
memory limit per test: 65536 KB
input: standard
output: standard



When entering some analog data into a computer, this information must be quantized. Quantization transforms each measured value x to some other value l(x) selected from the predefined set L of levels. Sometimes to reduce the influence of the levels set to the information, the group of levels sets Li is used. The number of levels sets is usually chosen to be the power of 2.

When using the number of levels sets, some additional information should be used to specify which set was used for each quantization. However, providing this information may be too expensive — the better solution would be to choose more levels and use one set. To avoid the specification of the quantization set, the following technique is used. Suppose that n values x1, x2, ..., xn are to be quantized and the group of m=2p levels sets Li, i=0, ..., m-1; each of size s=2q is used to quantize it. After quantization xj is replaced with some number lj = in Lf(j). Instead of sending lj, its ordinal number in Lf(j) is usually sent, let kj be the ordinal number of lj in Lf(j) (levels are numbered starting with 0). Take p least significant bits of kj and say that the number kj & (2p-1) is the number of the levels set that will be used for next quantization, that is f(j+1) = kj & (2p-1).

Since the least significant bits of kj are usually distributed quite randomly, the sets used for quantization change often and weakly depend on values of quantized data, thus this technique provides the good way to perform the quantization.

Usually to perform the quantization the closest to the value level of the levels set is chosen. However, using the technique described, sometimes it pays off to choose not the optimal level, but some other one, the ordinal number of which has other least significant bits, thus choosing another levels set for next measure and providing better approximation of quantized values in the future. Let us call the deviation of quantization the value of sum(j=1.. n, |xj - lj|). Your task is given measures and levels sets to choose quantized value for each measure in such a way, that the deviation of quantization is minimal possible.

The first value is always quantized using set L0.

Input
The first line of the input file contains n (1 ≤ n ≤ 1000). The second line contains n integer numbers xi ranging from 1 to 106. The next line contains m and s (1 ≤ m ≤ 128, m ≤ s ≤ 128). Next m lines contain s integer numbers each — levels of the quantization sets given in increasing order for each set, all levels satisfy 1 ≤ level ≤ 106.

Output
First output the minimal possible deviation of the quantization. Then output n integer numbers in range from 0 to s - 1. For each input value output the number of the level in the corresponding levels set (kj) used for this number to achieve the quantization required.

Sample test(s)

Input
3
8 8 19
2 4
5 10 15 20
3 7 13 17

Output
5
1 1 3
Author:	Andrew Stankevich
Resource:	Petrozavodsk Summer Trainings 2003
Date:	2003-08-23

=================
p205.in1
======================
3
8 8 19
2 4
5 10 15 20
3 7 13 17

=================
