Submission #239954

#	Time	Username	Problem	Language	Result	Execution time	Memory
239954		Vimmer	Priglavci (COCI18_priglavci)	C++14	96 / 160	46 ms	896 KiB

priglvaci

#include <bits/stdc++.h>
//#include <ext/pb_ds/assoc_container.hpp>
//#include <ext/pb_ds/tree_policy.hpp>

//#pragma GCC optimize("unroll-loops")
//#pragma GCC optimize("-O3")
//#pragma GCC optimize("Ofast")
//#pragma GCC optimize("fast-math")
//#pragma GCC optimize("no-stack-protector")

#define F first
#define S second
#define sz(x) int(x.size())
#define pb push_back
#define N 101001
#define M ll(1e9 + 7)
#define inf 1e9 + 1e9


using namespace std;
//using namespace __gnu_pbds;

typedef long double ld;
typedef long long ll;
typedef short int si;
typedef array <int, 3> a3;

//typedef tree <int, null_type, less_equal<int>, rb_tree_tag, tree_order_statistics_node_update> ordered_set;



int id[105], a[105];

int mn[105][105], number[105][105], n, m, c, k;

set <int> se[105];

ll limit;

int rst(int x, int y, int cx, int cy) {return (cx - x) * (cx - x) + (cy - y) * (cy - y);}

set <pair <int, int> > g[105];

bool mk[105];

bool kuna(int v)
{
    if (mk[v]) return 0;

    mk[v] = 1;

    for (auto it : g[v])
        if (it.F <= limit && (sz(se[it.S]) == 0 || kuna(*se[it.S].begin())))
            {
                if (sz(se[it.S]) != 0) se[it.S].erase(se[it.S].begin());

                se[it.S].insert(v);

                return 1;
            }

    return 0;
}

bool gd()
{
    memset(a, 0, sizeof(a));

    for (int i = 1; i <= k; i++)  se[i].clear();

    for (int i = 0; i < n; i++)
    {
        memset(mk, 0, sizeof(mk));

        kuna(i);
    }

    for (int i = 1; i <= k; i++)
        for (auto it : se[i]) a[it] = i;

    for (int i = 0; i < n; i++) if (a[i] == 0) return 0;

    return 1;
}
int main()
{
    //freopen("input.txt", "r", stdin); //freopen("output4.txt", "w", stdout);

    ios_base::sync_with_stdio(0); istream::sync_with_stdio(0); cin.tie(0); cout.tie(0);

    cin >> n >> m >> c >> k;

    if (c * k < n) {cout << -1 << endl; exit(0);}

    int cx[n], cy[n];

    for (int i = 0; i < n; i++) cin >> cx[i] >> cy[i];

    for (int i = 0; i < 105; i++)
        for (int j = 0; j < 105; j++) mn[i][j] = 1e9;

    int sx[m], sy[m];

    for (int i = 0; i < m; i++) cin >> sx[i] >> sy[i];

    for (int i = 1; i <= k; i++)
    {
        int kol;

        cin >> kol;

        for (; kol > 0; kol--) {int x; cin >> x; x--; id[x] = i;}
    }

    vector <int> values; values.clear();

    for (int i = 0; i < n; i++)
        for (int j = 0; j < m; j++)
          {ll rast = rst(cx[i], cy[i], sx[j], sy[j]); if (mn[i][id[j]] > rast) {mn[i][id[j]] = rast; number[i][id[j]] = j + 1;} values.pb(rast);}

    for (int i = 0; i < n; i++)
        for (int j = 1; j <= k; j++)
          g[i].insert({mn[i][j], j});

    sort(values.begin(), values.end());

    int l = 0, r = sz(values) - 1;

    while (l + 1 < r)
    {
        int md = (l + r) >> 1;

        limit = values[md];

        if (gd()) r = md; else l = md;
    }

    limit = values[l];

    while (gd() && l >= 0) {l--; if (l != -1) limit = values[l];};

    l++;

    limit = values[l];

    gd();

    cout << limit << endl;

    for (int i = 0; i < n; i++) cout << number[i][a[i]] << endl;
}

• Subtask 1: 96 / 160