# | Submission time | Handle | Problem | Language | Result | Execution time | Memory |
---|---|---|---|---|---|---|---|
796281 | 2023-07-28T08:49:32 Z | AbdullahMohammedAhmad | Minerals (JOI19_minerals) | C++14 | 1 ms | 592 KB |
#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; } } 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
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 1 ms | 592 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |
# | Verdict | Execution time | Memory | Grader output |
---|---|---|---|---|
1 | Incorrect | 0 ms | 336 KB | Wrong Answer [4] |
2 | Halted | 0 ms | 0 KB | - |