test/yukicoder/1867.test.cpp

• View this file on GitHub
• Last update: 2024-09-01 14:21:07+09:00
• Problem: https://yukicoder.me/problems/no/1867

Depends on

• 関数型や関数オブジェクトに関するテンプレート (src/Internal/function_traits.hpp)
• Li-Chao-Tree (src/Optimization/LiChaoTree.hpp)
• 最大最小を指定するための列挙型 (src/Optimization/MinMaxEnum.hpp)

Code

// competitive-verifier: PROBLEM https://yukicoder.me/problems/no/1867
// competitive-verifier: TLE 2
// competitive-verifier: MLE 256
// anti monge
#include <iostream>
#include <algorithm>
#include <vector>
#include "src/Optimization/LiChaoTree.hpp"
using namespace std;
signed main() {
 cin.tie(0);
 ios::sync_with_stdio(false);
 int N;
 cin >> N;
 int P[N];
 for (int i= 0; i < N; ++i) cin >> P[i];
 static constexpr int INF= 1e8;
 vector D(N, vector(N + 1, INF));
 for (int i= 0; i < N - 1; ++i) {
  int R[N + 1];
  fill_n(R, N + 1, 0);
  for (int j= 0; j <= i; ++j) ++R[P[j]];
  for (int j= 0; j < N; ++j) R[j + 1]+= R[j];
  int sum= 0;
  for (int j= i + 1; j < N; ++j) D[i + 1][j + 1]= sum+= R[N] - R[P[j]];
 }
 int dp[N + 1];
 fill_n(dp, N + 1, 0);
 cout << 0 << '\n';
 for (int k= 1; k < N; ++k) {
  int ndp[N + 1];
  LiChaoTree lct([&](int i, int j) { return dp[j] + D[j][i]; }, k, N + 1);
  auto tree= lct.make_tree<MINIMIZE>();
  for (int i= k; i <= N; ++i) {
   ndp[i]= tree.query(i).first;
   if (i < N) tree.insert(i, i + 1);
  }
  copy(ndp + k, ndp + N + 1, dp + k);
  cout << dp[N] << '\n';
 }
 return 0;
}

#line 1 "test/yukicoder/1867.test.cpp"
// competitive-verifier: PROBLEM https://yukicoder.me/problems/no/1867
// competitive-verifier: TLE 2
// competitive-verifier: MLE 256
// anti monge
#include <iostream>
#include <algorithm>
#include <vector>
#line 2 "src/Optimization/LiChaoTree.hpp"
#include <limits>
#line 5 "src/Optimization/LiChaoTree.hpp"
#include <tuple>
#include <cassert>
#line 2 "src/Internal/function_traits.hpp"
#include <type_traits>
// clang-format off
namespace function_template_internal{
template<class C>struct is_function_object{
 template<class U,int dummy=(&U::operator(),0)> static std::true_type check(U *);
 static std::false_type check(...);
 static C *m;
 static constexpr bool value= decltype(check(m))::value;
};
template<class F,bool,bool>struct function_type_impl{using type= void;};
template<class F>struct function_type_impl<F,true,false>{using type= F *;};
template<class F>struct function_type_impl<F,false,true>{using type= decltype(&F::operator());};
template<class F> using function_type_t= typename function_type_impl<F,std::is_function_v<F>,is_function_object<F>::value>::type;
template<class... Args>struct result_type_impl{using type= void;};
template<class R,class... Args>struct result_type_impl<R(*)(Args...)>{using type= R;};
template<class C,class R,class... Args>struct result_type_impl<R(C::*)(Args...)>{using type= R;};
template<class C,class R,class... Args>struct result_type_impl<R(C::*)(Args...)const>{using type= R;};
template<class F> using result_type_t= typename result_type_impl<function_type_t<F>>::type;
template<class... Args>struct argument_type_impl{using type= void;};
template<class R,class... Args>struct argument_type_impl<R(*)(Args...)>{using type= std::tuple<Args...>;};
template<class C,class R,class... Args>struct argument_type_impl<R(C::*)(Args...)>{using type= std::tuple<Args...>;};
template<class C,class R,class... Args>struct argument_type_impl<R(C::*)(Args...)const>{using type= std::tuple<Args...>;};
template<class F> using argument_type_t= typename argument_type_impl<function_type_t<F>>::type;
}
using function_template_internal::result_type_t,function_template_internal::argument_type_t;
// clang-format on
#line 2 "src/Optimization/MinMaxEnum.hpp"
enum MinMaxEnum { MAXIMIZE= -1, MINIMIZE= 1 };
#line 9 "src/Optimization/LiChaoTree.hpp"
template <class F, class T> class LiChaoTree {};
template <class F, class T, class... Prms> class LiChaoTree<F, std::tuple<T, Prms...>> {
 using R= result_type_t<F>;
 F f;
 const T LB, UB;
 std::vector<std::tuple<Prms...>> ps;
 template <MinMaxEnum sgn, bool persistent> class LiChaoTreeInterface {
  LiChaoTree *ins;
  struct Node {
   int id= -1;
   Node *ch[2]= {nullptr, nullptr};
  } *root;
  inline R eval(int id, T x) const { return std::apply(ins->f, std::tuple_cat(std::make_tuple(x), ins->ps[id])); }
  static inline bool cmp(const R &p, const R &n, int pi, int ni) {
   if constexpr (sgn == MINIMIZE) return p > n || (p == n && pi > ni);
   else return p < n || (p == n && pi > ni);
  }
  inline bool cmp(T x, int pi, int ni) const {
   if (ni == -1) return false;
   if (pi == -1) return true;
   return cmp(eval(pi, x), eval(ni, x), pi, ni);
  }
  static inline T ub(T r) {
   if constexpr (std::is_floating_point_v<T>) return r;
   else return r - 1;
  }
  inline void addl(Node *&t, int id, T xl, T xr) {
   if (!t) return t= new Node{id}, void();
   bool bl= cmp(xl, t->id, id), br= cmp(ub(xr), t->id, id);
   if (!bl && !br) return;
   if constexpr (persistent) t= new Node(*t);
   if (bl && br) return t->id= id, void();
   T xm= (xl + xr) / 2;
   if (cmp(xm, t->id, id)) std::swap(t->id, id), bl= !bl;
   if constexpr (std::is_floating_point_v<T>)
    if (xr - xl < 1e-9) return;
   bl ? addl(t->ch[0], id, xl, xm) : addl(t->ch[1], id, xm, xr);
  }
  inline void adds(Node *&t, int id, T l, T r, T xl, T xr) {
   if (r <= xl || xr <= l) return;
   if (l <= xl && xr <= r) return addl(t, id, xl, xr);
   if (!t) t= new Node;
   else if constexpr (persistent) t= new Node(*t);
   T xm= (xl + xr) / 2;
   adds(t->ch[0], id, l, r, xl, xm), adds(t->ch[1], id, l, r, xm, xr);
  }
  inline std::pair<R, int> query(const Node *t, T x, T xl, T xr) const {
   if (!t) return {R(), -1};
   T xm= (xl + xr) / 2;
   auto b= x < xm ? query(t->ch[0], x, xl, xm) : query(t->ch[1], x, xm, xr);
   if (t->id == -1) return b;
   R a= eval(t->id, x);
   return b.second != -1 && cmp(a, b.first, t->id, b.second) ? b : std::make_pair(a, t->id);
  }
 public:
  LiChaoTreeInterface()= default;
  LiChaoTreeInterface(LiChaoTree *ins): ins(ins), root(nullptr) {}
  void insert(const Prms &...args) { ins->ps.emplace_back(args...), addl(root, ins->ps.size() - 1, ins->LB, ins->UB); }
  // [l, r)
  void insert(const Prms &...args, T l, T r) {
   l= std::max(l, ins->LB), r= std::min(r, ins->UB);
   if (l < r) ins->ps.emplace_back(args...), adds(root, ins->ps.size() - 1, l, r, ins->LB, ins->UB);
  }
  // [l, UB)
  void insert(const Prms &...args, T l) { insert(args..., l, ins->UB); }
  std::pair<R, int> query(T x) const { return assert(ins->LB <= x && x < ins->UB), query(root, x, ins->LB, ins->UB); }
  const std::tuple<Prms...> &params(int id) const { return ins->ps[id]; }
 };
public:
 LiChaoTree(const F &f, T LB= -2e9, T UB= 2e9): f(f), LB(LB), UB(UB) {}
 template <MinMaxEnum sgn= MINIMIZE, bool persistent= false> LiChaoTreeInterface<sgn, persistent> make_tree() { return this; }
};
template <class F, class T, class U> LiChaoTree(F, T, U) -> LiChaoTree<F, argument_type_t<F>>;
template <class F, class T> LiChaoTree(F, T) -> LiChaoTree<F, argument_type_t<F>>;
template <class F> LiChaoTree(F) -> LiChaoTree<F, argument_type_t<F>>;
#line 9 "test/yukicoder/1867.test.cpp"
using namespace std;
signed main() {
 cin.tie(0);
 ios::sync_with_stdio(false);
 int N;
 cin >> N;
 int P[N];
 for (int i= 0; i < N; ++i) cin >> P[i];
 static constexpr int INF= 1e8;
 vector D(N, vector(N + 1, INF));
 for (int i= 0; i < N - 1; ++i) {
  int R[N + 1];
  fill_n(R, N + 1, 0);
  for (int j= 0; j <= i; ++j) ++R[P[j]];
  for (int j= 0; j < N; ++j) R[j + 1]+= R[j];
  int sum= 0;
  for (int j= i + 1; j < N; ++j) D[i + 1][j + 1]= sum+= R[N] - R[P[j]];
 }
 int dp[N + 1];
 fill_n(dp, N + 1, 0);
 cout << 0 << '\n';
 for (int k= 1; k < N; ++k) {
  int ndp[N + 1];
  LiChaoTree lct([&](int i, int j) { return dp[j] + D[j][i]; }, k, N + 1);
  auto tree= lct.make_tree<MINIMIZE>();
  for (int i= k; i <= N; ++i) {
   ndp[i]= tree.query(i).first;
   if (i < N) tree.insert(i, i + 1);
  }
  copy(ndp + k, ndp + N + 1, dp + k);
  cout << dp[N] << '\n';
 }
 return 0;
}

Test cases

Env	Name	Status	Elapsed	Memory
g++-13	N15fil1.txt	AC	6 ms	4 MB
g++-13	N15fil2.txt	AC	5 ms	4 MB
g++-13	N15fil3.txt	AC	5 ms	4 MB
g++-13	N15fil4.txt	AC	5 ms	3 MB
g++-13	N15fil5.txt	AC	5 ms	4 MB
g++-13	N2000fil1.txt	AC	399 ms	108 MB
g++-13	N2000fil10.txt	AC	402 ms	108 MB
g++-13	N2000fil11.txt	AC	406 ms	108 MB
g++-13	N2000fil12.txt	AC	405 ms	108 MB
g++-13	N2000fil13.txt	AC	397 ms	108 MB
g++-13	N2000fil14.txt	AC	404 ms	108 MB
g++-13	N2000fil15.txt	AC	403 ms	108 MB
g++-13	N2000fil16.txt	AC	402 ms	108 MB
g++-13	N2000fil17.txt	AC	399 ms	108 MB
g++-13	N2000fil18.txt	AC	401 ms	108 MB
g++-13	N2000fil19.txt	AC	401 ms	108 MB
g++-13	N2000fil2.txt	AC	404 ms	108 MB
g++-13	N2000fil20.txt	AC	405 ms	108 MB
g++-13	N2000fil3.txt	AC	405 ms	108 MB
g++-13	N2000fil4.txt	AC	399 ms	108 MB
g++-13	N2000fil5.txt	AC	402 ms	108 MB
g++-13	N2000fil6.txt	AC	403 ms	108 MB
g++-13	N2000fil7.txt	AC	400 ms	108 MB
g++-13	N2000fil8.txt	AC	401 ms	108 MB
g++-13	N2000fil9.txt	AC	405 ms	108 MB
g++-13	N500fil1.txt	AC	24 ms	10 MB
g++-13	N500fil2.txt	AC	24 ms	10 MB
g++-13	N500fil3.txt	AC	23 ms	10 MB
g++-13	N500fil4.txt	AC	23 ms	10 MB
g++-13	N500fil5.txt	AC	24 ms	10 MB
g++-13	corner1.txt	AC	6 ms	3 MB
g++-13	corner2.txt	AC	6 ms	4 MB
g++-13	maximum1.txt	AC	968 ms	238 MB
g++-13	maximum10.txt	AC	959 ms	237 MB
g++-13	maximum11.txt	AC	966 ms	238 MB
g++-13	maximum12.txt	AC	954 ms	238 MB
g++-13	maximum13.txt	AC	954 ms	238 MB
g++-13	maximum14.txt	AC	959 ms	238 MB
g++-13	maximum15.txt	AC	956 ms	238 MB
g++-13	maximum16.txt	AC	960 ms	237 MB
g++-13	maximum17.txt	AC	955 ms	237 MB
g++-13	maximum18.txt	AC	963 ms	238 MB
g++-13	maximum19.txt	AC	962 ms	238 MB
g++-13	maximum2.txt	AC	959 ms	238 MB
g++-13	maximum20.txt	AC	955 ms	238 MB
g++-13	maximum21.txt	AC	950 ms	238 MB
g++-13	maximum22.txt	AC	951 ms	238 MB
g++-13	maximum23.txt	AC	958 ms	237 MB
g++-13	maximum24.txt	AC	960 ms	237 MB
g++-13	maximum25.txt	AC	962 ms	237 MB
g++-13	maximum26.txt	AC	956 ms	237 MB
g++-13	maximum27.txt	AC	956 ms	238 MB
g++-13	maximum28.txt	AC	961 ms	237 MB
g++-13	maximum29.txt	AC	961 ms	238 MB
g++-13	maximum3.txt	AC	957 ms	238 MB
g++-13	maximum30.txt	AC	956 ms	237 MB
g++-13	maximum4.txt	AC	962 ms	237 MB
g++-13	maximum5.txt	AC	951 ms	238 MB
g++-13	maximum6.txt	AC	959 ms	237 MB
g++-13	maximum7.txt	AC	959 ms	238 MB
g++-13	maximum8.txt	AC	958 ms	237 MB
g++-13	maximum9.txt	AC	950 ms	238 MB
g++-13	minimum.txt	AC	6 ms	4 MB
g++-13	revsortedperm.txt	AC	892 ms	238 MB
g++-13	sample1.txt	AC	6 ms	4 MB
g++-13	sample2.txt	AC	5 ms	4 MB
g++-13	sample3.txt	AC	5 ms	4 MB
g++-13	sortedperm.txt	AC	867 ms	238 MB
clang++-18	N15fil1.txt	AC	9 ms	4 MB
clang++-18	N15fil2.txt	AC	6 ms	4 MB
clang++-18	N15fil3.txt	AC	6 ms	4 MB
clang++-18	N15fil4.txt	AC	5 ms	4 MB
clang++-18	N15fil5.txt	AC	5 ms	4 MB
clang++-18	N2000fil1.txt	AC	429 ms	108 MB
clang++-18	N2000fil10.txt	AC	430 ms	108 MB
clang++-18	N2000fil11.txt	AC	431 ms	108 MB
clang++-18	N2000fil12.txt	AC	431 ms	108 MB
clang++-18	N2000fil13.txt	AC	430 ms	108 MB
clang++-18	N2000fil14.txt	AC	427 ms	108 MB
clang++-18	N2000fil15.txt	AC	427 ms	108 MB
clang++-18	N2000fil16.txt	AC	425 ms	108 MB
clang++-18	N2000fil17.txt	AC	427 ms	108 MB
clang++-18	N2000fil18.txt	AC	434 ms	108 MB
clang++-18	N2000fil19.txt	AC	429 ms	108 MB
clang++-18	N2000fil2.txt	AC	426 ms	108 MB
clang++-18	N2000fil20.txt	AC	436 ms	108 MB
clang++-18	N2000fil3.txt	AC	429 ms	108 MB
clang++-18	N2000fil4.txt	AC	428 ms	108 MB
clang++-18	N2000fil5.txt	AC	423 ms	108 MB
clang++-18	N2000fil6.txt	AC	430 ms	108 MB
clang++-18	N2000fil7.txt	AC	424 ms	108 MB
clang++-18	N2000fil8.txt	AC	429 ms	108 MB
clang++-18	N2000fil9.txt	AC	434 ms	108 MB
clang++-18	N500fil1.txt	AC	24 ms	10 MB
clang++-18	N500fil2.txt	AC	24 ms	10 MB
clang++-18	N500fil3.txt	AC	23 ms	10 MB
clang++-18	N500fil4.txt	AC	23 ms	10 MB
clang++-18	N500fil5.txt	AC	24 ms	10 MB
clang++-18	corner1.txt	AC	5 ms	4 MB
clang++-18	corner2.txt	AC	5 ms	4 MB
clang++-18	maximum1.txt	AC	1037 ms	238 MB
clang++-18	maximum10.txt	AC	1036 ms	238 MB
clang++-18	maximum11.txt	AC	1043 ms	237 MB
clang++-18	maximum12.txt	AC	1040 ms	238 MB
clang++-18	maximum13.txt	AC	1033 ms	238 MB
clang++-18	maximum14.txt	AC	1033 ms	237 MB
clang++-18	maximum15.txt	AC	1036 ms	238 MB
clang++-18	maximum16.txt	AC	1031 ms	237 MB
clang++-18	maximum17.txt	AC	1027 ms	238 MB
clang++-18	maximum18.txt	AC	1043 ms	237 MB
clang++-18	maximum19.txt	AC	1037 ms	238 MB
clang++-18	maximum2.txt	AC	1040 ms	238 MB
clang++-18	maximum20.txt	AC	1049 ms	238 MB
clang++-18	maximum21.txt	AC	1031 ms	238 MB
clang++-18	maximum22.txt	AC	1025 ms	237 MB
clang++-18	maximum23.txt	AC	1063 ms	238 MB
clang++-18	maximum24.txt	AC	1029 ms	237 MB
clang++-18	maximum25.txt	AC	1028 ms	237 MB
clang++-18	maximum26.txt	AC	1036 ms	237 MB
clang++-18	maximum27.txt	AC	1035 ms	237 MB
clang++-18	maximum28.txt	AC	1034 ms	237 MB
clang++-18	maximum29.txt	AC	1042 ms	237 MB
clang++-18	maximum3.txt	AC	1042 ms	238 MB
clang++-18	maximum30.txt	AC	1043 ms	238 MB
clang++-18	maximum4.txt	AC	1034 ms	237 MB
clang++-18	maximum5.txt	AC	1048 ms	237 MB
clang++-18	maximum6.txt	AC	1032 ms	238 MB
clang++-18	maximum7.txt	AC	1030 ms	237 MB
clang++-18	maximum8.txt	AC	1043 ms	237 MB
clang++-18	maximum9.txt	AC	1044 ms	238 MB
clang++-18	minimum.txt	AC	6 ms	4 MB
clang++-18	revsortedperm.txt	AC	1000 ms	237 MB
clang++-18	sample1.txt	AC	6 ms	4 MB
clang++-18	sample2.txt	AC	5 ms	4 MB
clang++-18	sample3.txt	AC	5 ms	4 MB
clang++-18	sortedperm.txt	AC	975 ms	238 MB

Back to top page