# | Time | Username | Problem | Language | Result | Execution time | Memory |
---|---|---|---|---|---|---|---|
796449 | AbdullahMohammedAhmad | Minerals (JOI19_minerals) | C++14 | 1 ms | 464 KiB |
This submission is migrated from previous version of oj.uz, which used different machine for grading. This submission may have different result if resubmitted.
#include <bits/stdc++.h>
#include "minerals.h"
using namespace std;
int common[18][18];
int cost[18][18];
int have[18];
int nodes;
int dp[18][1<<19];
int par[18][1<<19];
int tsp(int i, int mask)
{
if(mask == 1<<i)
{
return 0;
}
if(dp[i][mask] != 0)
{
return dp[i][mask];
}
int res = INT_MAX;
for(int j = 0; j < nodes; j++)
{
if((mask&(1<<j)) && j != i)
{
int now = tsp(j, mask^(1<<i))+cost[j][i];
if(now < res)
{
res = now;
par[i][mask] = j;
}
}
}
return dp[i][mask] = res;
}
vector<int> get_best_order(int n)
{
for(int i = 0; (1<<i) < n; i++) {
for(int k = 0; k < n; k++)
{
if(k&(1<<i))
{
have[i]++;
}
}
for(int j = i+1; (1<<j) < n; j++)
{
for(int k = 0; k < n; k++)
{
if(k&(1<<i) && k&(1<<j))
{
common[i][j]++;
common[j][i]++;
}
}
cost[i][j] = cost[j][i] = have[i]+have[j]-2*common[i][j];
}
}
nodes = 31-__builtin_clz(n);
int ans = INT_MAX;
int idx = -1;
for(int i = 0; i < nodes; i++)
{
tsp(i, (1<<nodes)-1);
if(dp[i][(1<<nodes)-1]+have[i] < ans)
{
ans = dp[i][(1<<nodes)-1]+have[i];
idx = i;
}
}
assert(idx != -1);
vector<int> path;
path.push_back(idx);
int msk = (1<<nodes)-1;
while(true)
{
if(msk == 1<<idx)
{
break;
}
int nxt = par[idx][msk];
msk^=(1<<idx);
idx = nxt;
path.push_back(idx);
}
assert(path.size() == nodes);
return path;
}
mt19937 rng(66);
void Solve(int n) {
vector<int> p1, p2;
int ans = 0, prevans = -1;
vector<int> bulanea(2 * n);
iota(bulanea.begin(), bulanea.end(), 1);
shuffle(bulanea.begin(), bulanea.end(), rng);
for(int i: bulanea) {
ans = Query(i);
if(ans == prevans) {
p2.push_back(i);
} else {
p1.push_back(i);
}
prevans = ans;
}
assert(p1.size() == n && p2.size() == n);
vector<int> pr(n, 0);
vector<bool> bagat(n, 1);
vector<int> best_order = get_best_order(n);
for(int bit: best_order) {
int cntt = 0;
for(int i = 0; i < n; ++i) {
if(i >> bit & 1) {
++cntt;
if(!bagat[i]) prevans = Query(p1[i]);
bagat[i] = 1;
} else {
if(bagat[i]) prevans = Query(p1[i]);
bagat[i] = 0;
}
}
for(int i = 0; i < n && cntt; ++i) {
if(pr[i] + (1 << bit) >= n) continue;
int ans = Query(p2[i]);
if(ans == prevans) {
pr[i] += (1 << bit);
--cntt;
}
prevans = ans;
}
}
for(int i = 0; i < n; ++i) {
Answer(p2[i], p1[pr[i]]);
}
}
Compilation message (stderr)
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
Fetching results... |