#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], mt[105], number[105][105], n, m, c, k;
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 && (mt[it.S] == -1 || kuna(mt[it.S])))
{
mt[it.S] = v;
return 1;
}
return 0;
}
bool gd()
{
for (int i = 1; i <= k; i++) mt[i] = -1;
for (int i = 0; i < n; i++)
{
memset(mk, 0, sizeof(mk));
kuna(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];
limit = values[++l];
while (!gd() && l + 1 < sz(values)) {l++; limit = values[l];}
limit = values[l - 1];
int ans = 0;
for (int i = 1; i <= k; i++) mt[i] = -1;
for (int i = 0; i < n; i++)
{
memset(mk, 0, sizeof(mk));
kuna(i);
}
for (int i = 1; i <= k; i++) a[mt[i]] = i;
for (int i = 0; i < n; i++) ans = max(ans,rst(cx[i], cy[i], sx[number[i][a[i]] - 1], sy[number[i][a[i]] - 1]));
cout << ans << endl;
for (int i = 0; i < n; i++) cout << number[i][a[i]] << endl;
}
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Incorrect |
6 ms |
640 KB |
Output isn't correct |
2 |
Correct |
5 ms |
640 KB |
Output is correct |
3 |
Correct |
5 ms |
384 KB |
Output is correct |
4 |
Incorrect |
7 ms |
896 KB |
Output isn't correct |
5 |
Incorrect |
8 ms |
768 KB |
Output isn't correct |
6 |
Incorrect |
45 ms |
896 KB |
Integer 0 violates the range [1, 50] |
7 |
Incorrect |
8 ms |
512 KB |
Integer 0 violates the range [1, 25] |
8 |
Incorrect |
6 ms |
512 KB |
Output isn't correct |
9 |
Incorrect |
39 ms |
768 KB |
Integer 0 violates the range [1, 50] |
10 |
Incorrect |
9 ms |
512 KB |
Integer 0 violates the range [1, 100] |