Submission #1266241

#TimeUsernameProblemLanguageResultExecution timeMemory
1266241wjqFun Tour (APIO20_fun)C++20
100 / 100
181 ms23144 KiB
#include"fun.h"
#include<bits/stdc++.h>
int n;
namespace Task1{
  const int N=25;
  int dis[N][N],dp[1<<20][20],pre[1<<20][20];
  std::vector<int> solve(int n){
    std::vector<int> ans;
    for(int i=0;i<n;i++)
      for(int j=i+1;j<n;j++){
        dis[i][j]=dis[j][i]=hoursRequired(i,j);
      }
    for(int i=0;i<n;++i)dp[1<<i][i]=1e9;
    for(int i=1;i<n;++i)
      for(int S=1;S<1<<n;++S)if(__builtin_popcount(S)==i)
        for(int u=0;u<n;++u)if((S>>u)&1)
          for(int v=0;v<n;++v)if(!((S>>v)&1))
            if(dp[S][u]>=dis[u][v]&&dp[S|(1<<v)][v]<dis[u][v]){
              dp[S|(1<<v)][v]=dis[u][v];
              pre[S|(1<<v)][v]=u;
            }
    for(int i=0;i<n;++i)if(dp[(1<<n)-1][i]<1e9){
      int j,s=(1<<n)-1;
      while(s){
        ans.push_back(i);
        j=pre[s][i],s^=1<<i,i=j;
      }
      std::reverse(ans.begin(),ans.end());
      return ans;
    }
    assert(0);
  }
}
namespace Task0{
  inline bool cmp(const std::pair<int,int> &a,const std::pair<int,int> &b){
    return a.second<b.second;
  }
  int find(int rt){
    std::vector<int> s(n);
    for(int i=0;i<n;i++)s[i]=attractionsBehind(rt,i);
    bool fl=true;
    for(int i=0;i<n;i++)fl&=s[i]<=n/2;
    if(fl)return rt;
    int min=INT_MAX,cho=-1;
    for(int i=0;i<n;i++)if(s[i]>n/2&&s[i]<min){
      min=s[i],cho=i;
    }
    return cho;
  }
  inline void sort(std::vector<std::pair<int,int>> &a,std::vector<std::pair<int,int>> &b,std::vector<std::pair<int,int>> &c){
    if(a.size()<c.size())std::swap(a,c);
    if(a.size()<b.size())std::swap(a,b);
    if(b.size()<c.size())std::swap(b,c);
  }
  std::vector<std::pair<int,int>> build(std::vector<std::pair<int,int>> a,std::vector<std::pair<int,int>> b,std::vector<std::pair<int,int>> c,int u){
    std::vector<std::pair<int,int>> ans;
    sort(a,b,c);
    if(a.size()&&b.size()){
      if(cmp(a.back(),b.back())&&b.size()+c.size()>=a.size())std::swap(a,b);
    }
    while(a.size()||b.size()||c.size()){ //a褰撳墠,b涓婁竴娆?c鍓╀笅鐨?
#ifdef local
      printf("(%d %d) (%d %d) (%d %d)\n",
        a.size()?a.back().first:-1,a.size()?a.back().second:-1,
        b.size()?b.back().first:-1,b.size()?b.back().second:-1,
        c.size()?c.back().first:-1,c.size()?c.back().second:-1);
#endif
      if(a.size()>b.size()+c.size()&&c.size()){
#ifdef local
        printf("merge !!\n");
#endif
        b.insert(b.end(),c.begin(),c.end()),c.clear();
        std::sort(b.begin(),b.end(),cmp);
        // for(const auto &x:b)printf("> %d %d\n",x.first,x.second);
        if(b.size()&&ans.size()&&cmp(ans.back(),b.back())){
          ans.push_back(b.back()),b.pop_back();
        }
      }
      if(!c.size()){ //cross
        if(!a.size())a.push_back({u,0}),u=-1;
        ans.push_back(a.back()),a.pop_back();
        std::swap(a,b);
        continue;
      }
      if(cmp(a.back(),c.back()))std::swap(a,c);
      ans.push_back(a.back()),a.pop_back();
      std::swap(a,b);
      if(c.size()&&a.size()<c.size())std::swap(a,c);
    }
    if(~u)ans.push_back({u,0});
#ifdef local
    for(const auto &x:ans)printf("> %d %d\n",x.first,x.second);
#endif
    return ans;
  }
  std::vector<int> solve(int u){
    // printf("solve %d\n",u);
    std::vector<int> ret,du(n),dv(n),G;
    for(int i=0;i<n;i++)if(i!=u&&hoursRequired(u,i)==1)G.push_back(i);
    std::vector<std::pair<int,int>> ans,S[3];
    for(int i=0;i<n;i++)if(i!=u)du[i]=hoursRequired(G[0],i);
    for(int i=0;i<n;i++)if(i!=u)dv[i]=hoursRequired(G[1],i);
    for(int i=0;i<n;i++)if(i!=u){
      if(du[i]<dv[i])S[0].push_back({i,du[i]+1});
      if(du[i]>dv[i])S[1].push_back({i,dv[i]+1});
      if(du[i]==dv[i])S[2].push_back({i,du[i]-1});
    }
    for(int fl=0;fl<G.size();fl++){
      std::sort(S[fl].begin(),S[fl].end(),cmp);
    }
#ifdef local
    for(int fl=0;fl<G.size();fl++){
      printf("[%d] => ",fl);
      for(auto x:S[fl])printf("(%d,%d) ",x.first,x.second);
      puts("");
    }
#endif
    ans=build(S[0],S[1],S[2],u);
    for(const auto &it:ans)ret.push_back(it.first);
    return ret;
  }
  std::vector<int> solve(){
    return solve(find(rand()%n));
  }
}
std::vector<int> createFunTour(int N,int Q){
  n=N;
  if(n<=17)return Task1::solve(n);
  return Task0::solve();
}






























#ifdef memset0
static void wrongAnswer(std::string message) {
  printf("WA: %s\n", message.c_str());
  exit(0);
}

namespace tree_helper {

static int N;
static int logN;
static std::vector<std::vector<int>> parent;
static std::vector<int> depth;
static std::vector<int> subtree_size;

static void dfs(
    const std::vector<std::vector<int>>& adj_list,
    int current_node, int parent_node) {
  parent[0][current_node] = parent_node;
  subtree_size[current_node] = 1;
  for (int i = 0; i < static_cast<int>(adj_list[current_node].size()); ++i) {
    const int next_node = adj_list[current_node][i];
    if (next_node != parent_node) {
      depth[next_node] = depth[current_node] + 1;
      dfs(adj_list, next_node, current_node);
      subtree_size[current_node] += subtree_size[next_node];
    }
  }
}

static void initializeTree(const std::vector<std::vector<int>>& adj_list) {
  N = static_cast<int>(adj_list.size());

  depth = std::vector<int>(N, 0);
  subtree_size = std::vector<int>(N, 0);
  for (logN = 0; (1 << logN) < N; ++logN) {}
  parent = std::vector<std::vector<int>>(logN, std::vector<int>(N, 0));

  dfs(adj_list, 0, 0);
  for (int i = 1; i < logN; ++i) {
    for (int j = 0; j < N; ++j) {
      parent[i][j] = parent[i - 1][parent[i - 1][j]];
    }
  }
}

static int getLowestCommonAncestor(int X, int Y) {
  if (depth[X] < depth[Y]) {
    std::swap(X, Y);
  }
  for (int i = logN - 1; i >= 0; --i) {
    if (depth[parent[i][X]] >= depth[Y]) {
      X = parent[i][X];
    }
  }
  if (X == Y) {
    return X;
  }
  for (int i = logN - 1; i >= 0; --i) {
    if (parent[i][X] != parent[i][Y]) {
      X = parent[i][X];
      Y = parent[i][Y];
    }
  }
  return parent[0][X];
}

static int getDistance(int X, int Y) {
  return depth[X] + depth[Y] - 2 * depth[getLowestCommonAncestor(X, Y)];
}

static int attractionsBehind(int X, int Y) {
  if (X == Y) {
    return N;
  }
  for (int i = logN - 1; i >= 0; --i) {
    if (depth[parent[i][X]] > depth[Y]) {
      X = parent[i][X];
    }
  }
  if (Y == parent[0][X]) {
    return N - subtree_size[X];
  }
  return subtree_size[Y];
}

static void checkFunTour(const std::vector<int>& fun_tour) {
  if (static_cast<int>(fun_tour.size()) != N) {
    wrongAnswer("Invalid size");
  }

  std::vector<bool> visited_attractions(N, false);
  for (int i = 0; i < N; ++i) {
    if (fun_tour[i] < 0 || fun_tour[i] >= N) {
      wrongAnswer("Invalid index");
    }
    if (visited_attractions[fun_tour[i]]) {
      wrongAnswer("Repeated index");
    }
    visited_attractions[fun_tour[i]] = true;
  }

  int last_travel_time = getDistance(fun_tour[0], fun_tour[1]);
  for (int i = 2; i < N; ++i) {
    int travel_time = getDistance(fun_tour[i - 1], fun_tour[i]);
    if (travel_time > last_travel_time) {
      wrongAnswer("Tour is not fun");
    }
    last_travel_time = travel_time;
  }
}

}  // namespace tree_helper

static int N, Q;

int hoursRequired(int X, int Y) {
  if (--Q < 0) {
    wrongAnswer("Too many queries");
  }
  if (X < 0 || X >= N || Y < 0 || Y >= N) {
    wrongAnswer("Invalid index");
  }
  return tree_helper::getDistance(X, Y);
}

int attractionsBehind(int X, int Y) {
  if (--Q < 0) {
    wrongAnswer("Too many queries");
  }
  if (X < 0 || X >= N || Y < 0 || Y >= N) {
    wrongAnswer("Invalid index");
  }
  return tree_helper::attractionsBehind(X, Y);
}

int main() {
#ifdef local
  freopen("../examples/03.in", "r", stdin);
#endif
  assert(2 == scanf("%d %d", &N, &Q));

  std::vector<std::vector<int>> adj_list(N);
  for (int i = 0; i < N - 1; ++i) {
    int A, B;
    assert(2 == scanf("%d %d", &A, &B));
    adj_list[A].push_back(B);
    adj_list[B].push_back(A);
  }
  tree_helper::initializeTree(adj_list);

  std::vector<int> fun_tour = createFunTour(N, Q);
  tree_helper::checkFunTour(fun_tour);

  for (int i = 0; i < N; ++i) {
    printf("%d%c", fun_tour[i], " \n"[i == N - 1]);
  }
  return 0;
}
#endif
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...