#define _CRT_NONSTDC_NO_WARNINGS #define _SILENCE_CXX17_ITERATOR_BASE_CLASS_DEPRECATION_WARNING #include #include #include #include #include #ifdef _MSC_VER #include #include #include #include #include #include #include #include #include /* g++ functions */ int __builtin_clz(unsigned int n) { unsigned long index; _BitScanReverse(&index, n); return 31 - index; } int __builtin_ctz(unsigned int n) { unsigned long index; _BitScanForward(&index, n); return index; } namespace std { inline int __lg(int __n) { return sizeof(int) * 8 - 1 - __builtin_clz(__n); } } int __builtin_popcount(int bits) { bits = (bits & 0x55555555) + (bits >> 1 & 0x55555555); bits = (bits & 0x33333333) + (bits >> 2 & 0x33333333); bits = (bits & 0x0f0f0f0f) + (bits >> 4 & 0x0f0f0f0f); bits = (bits & 0x00ff00ff) + (bits >> 8 & 0x00ff00ff); return (bits & 0x0000ffff) + (bits >> 16 & 0x0000ffff); } /* enable __uint128_t in MSVC */ //#include //using __uint128_t = boost::multiprecision::uint128_t; #else #pragma GCC target("avx2") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #endif /** compro io **/ namespace aux { template struct tp { static void output(std::ostream& os, const T& v) { os << std::get(v) << ", "; tp::output(os, v); } }; template struct tp { static void output(std::ostream& os, const T& v) { os << std::get(v); } }; } template std::ostream& operator<<(std::ostream& os, const std::tuple& t) { os << '{'; aux::tp, 0, sizeof...(Ts) - 1>::output(os, t); return os << '}'; } // tuple out template std::basic_ostream& operator<<(std::basic_ostream& os, const Container& x); // container out (fwd decl) template std::ostream& operator<<(std::ostream& os, const std::pair& p) { return os << '{' << p.first << ", " << p.second << '}'; } // pair out template std::istream& operator>>(std::istream& is, std::pair& p) { return is >> p.first >> p.second; } // pair in std::ostream& operator<<(std::ostream& os, const std::vector::reference& v) { os << (v ? '1' : '0'); return os; } // bool (vector) out std::ostream& operator<<(std::ostream& os, const std::vector& v) { bool f = true; os << '{'; for (const auto& x : v) { os << (f ? "" : ", ") << x; f = false; } os << '}'; return os; } // vector out template std::basic_ostream& operator<<(std::basic_ostream& os, const Container& x) { bool f = true; os << '{'; for (auto& y : x) { os << (f ? "" : ", ") << y; f = false; } return os << '}'; } // container out template())), class = typename std::enable_if::value>::type> std::istream& operator>>(std::istream& is, T& a) { for (auto& x : a) is >> x; return is; } // container in template auto operator<<(std::ostream& out, const T& t) -> decltype(out << t.stringify()) { out << t.stringify(); return out; } // struct (has stringify() func) out /** io setup **/ struct IOSetup { IOSetup(bool f) { if (f) { std::cin.tie(nullptr); std::ios::sync_with_stdio(false); } std::cout << std::fixed << std::setprecision(15); } } iosetup(true); // set false when solving interective problems /** string formatter **/ template std::string format(const std::string& f, Ts... t) { size_t l = std::snprintf(nullptr, 0, f.c_str(), t...); std::vector b(l + 1); std::snprintf(&b[0], l + 1, f.c_str(), t...); return std::string(&b[0], &b[0] + l); } /** dump **/ #define DUMPOUT std::cerr std::ostringstream DUMPBUF; #define dump(...) do{DUMPBUF<<" ";DUMPBUF<<#__VA_ARGS__<<" :[DUMP - "<<__LINE__<<":"<<__FUNCTION__<<']'< void dump_func(Head&& head, Tail&&... tail) { DUMPBUF << head; if (sizeof...(Tail) == 0) { DUMPBUF << " "; } else { DUMPBUF << ", "; } dump_func(std::move(tail)...); } /** timer **/ class Timer { double t = 0, paused = 0, tmp; public: Timer() { reset(); } static double time() { #ifdef _MSC_VER return __rdtsc() / 2.9e9; #else unsigned long long a, d; __asm__ volatile("rdtsc" : "=a"(a), "=d"(d)); return (d << 32 | a) / 2.9e9; #endif } void reset() { t = time(); } void pause() { tmp = time(); } void restart() { paused += time() - tmp; } double elapsed_ms() const { return (time() - t - paused) * 1000.0; } }; /** rand **/ struct Xorshift { Xorshift() {} Xorshift(uint64_t seed) { reseed(seed); } inline void reseed(uint64_t seed) { x = 0x498b3bc5 ^ seed; for (int i = 0; i < 20; i++) next_u64(); } inline uint64_t next_u64() { x ^= x << 7; return x ^= x >> 9; } inline uint32_t next_u32() { return next_u64() >> 32; } inline uint32_t next_u32(uint32_t mod) { return ((uint64_t)next_u32() * mod) >> 32; } inline uint32_t next_u32(uint32_t l, uint32_t r) { return l + next_u32(r - l + 1); } inline double next_double() { return next_u32() * e; } inline double next_double(double c) { return next_double() * c; } inline double next_double(double l, double r) { return next_double(r - l) + l; } private: static constexpr uint32_t M = UINT_MAX; static constexpr double e = 1.0 / M; uint64_t x = 88172645463325252LL; }; /** shuffle **/ template void shuffle_vector(std::vector& v, Xorshift& rnd) { int n = v.size(); for (int i = n - 1; i >= 1; i--) { int r = rnd.next_u32(i); std::swap(v[i], v[r]); } } /** split **/ std::vector split(const std::string& str, const std::string& delim) { std::vector res; std::string buf; for (const auto& c : str) { if (delim.find(c) != std::string::npos) { if (!buf.empty()) res.push_back(buf); buf.clear(); } else buf += c; } if (!buf.empty()) res.push_back(buf); return res; } std::string join(const std::string& delim, const std::vector& elems) { if (elems.empty()) return ""; std::string res = elems[0]; for (int i = 1; i < (int)elems.size(); i++) { res += delim + elems[i]; } return res; } /** misc **/ template inline void Fill(A(&array)[N], const T& val) { std::fill((T*)array, (T*)(array + N), val); } // fill array template auto make_vector(T x, int arg, Args ...args) { if constexpr (sizeof...(args) == 0)return std::vector(arg, x); else return std::vector(arg, make_vector(x, args...)); } template bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } return false; } template bool chmin(T& a, const T& b) { if (a > b) { a = b; return true; } return false; } /** hash **/ namespace aux { template inline void hash(std::size_t& s, const T& v) { s ^= std::hash()(v) + 0x9e3779b9 + (s << 6) + (s >> 2); } } namespace std { template struct hash> { size_t operator()(const std::pair& s) const noexcept { size_t seed = 0; aux::hash(seed, s.first); aux::hash(seed, s.second); return seed; } }; } /* fast queue */ class FastQueue { int front = 0; int back = 0; int v[4096]; public: inline bool empty() { return front == back; } inline void push(int x) { v[front++] = x; } inline int pop() { return v[back++]; } inline void reset() { front = back = 0; } inline int size() { return front - back; } }; class RandomQueue { int sz = 0; int v[4096]; public: inline bool empty() const { return !sz; } inline int size() const { return sz; } inline void push(int x) { v[sz++] = x; } inline void reset() { sz = 0; } inline int pop(int i) { std::swap(v[i], v[sz - 1]); return v[--sz]; } inline int pop(Xorshift& rnd) { return pop(rnd.next_u32(sz)); } }; #if 1 inline double get_temp(double stemp, double etemp, double t, double T) { return etemp + (stemp - etemp) * (T - t) / T; }; #else inline double get_temp(double stemp, double etemp, double t, double T) { return stemp * pow(etemp / stemp, t / T); }; #endif struct LogTable { static constexpr int M = 65536; static constexpr int mask = M - 1; double l[M]; LogTable() : l() { unsigned long long x = 88172645463325252ULL; double log_u64max = log(2) * 64; for (int i = 0; i < M; i++) { x = x ^ (x << 7); x = x ^ (x >> 9); l[i] = log(double(x)) - log_u64max; } } inline double operator[](int i) const { return l[i & mask]; } } log_table; constexpr int T = 52; constexpr int N = 10; constexpr int INITIAL_MONEY = 2000000; struct HiddenInput { const std::array D; const std::array, T> Z; static HiddenInput create(std::istream& in) { std::array D; std::array, T> Z; in >> D[0] >> D[0]; // ignore N, T in >> D >> Z; return HiddenInput(D, Z); } private: HiddenInput( const std::array D_, const std::array, T> Z_ ) : D(D_), Z(Z_) {} }; struct Result { int money; std::array S; std::array P; std::array R; }; struct Tester; using TesterPtr = std::shared_ptr; struct Tester { virtual Result query(const std::array& L) = 0; virtual Result query(const int x) = 0; }; struct ServerTester; using ServerTesterPtr = std::shared_ptr; struct ServerTester : Tester { Result receive_result() { Result res; in >> res.money; assert(res.money != -1); in >> res.S >> res.P >> res.R; return res; } Result query(const std::array& L) override { out << 1; for (int i = 0; i < N; i++) { out << ' ' << L[i]; } out << std::endl; return receive_result(); } Result query(const int x) override { out << 2 << ' ' << x << std::endl; return receive_result(); } static ServerTesterPtr create(std::istream& in, std::ostream& out) { return std::make_shared(ServerTester(in, out)); } private: std::istream& in; std::ostream& out; ServerTester(std::istream& in_, std::ostream& out_) : in(in_), out(out_) { int buf; in >> buf >> buf >> buf; // ignore T, N, Money } }; struct FileTester; using FileTesterPtr = std::shared_ptr; struct FileTester : Tester { const HiddenInput hinput; Result query(const std::array& L) override { return {}; } Result query(const int x) override { return {}; } static FileTesterPtr create(std::istream& in, std::ostream& out) { const auto hinput = HiddenInput::create(in); return std::make_shared(FileTester(hinput, out)); } private: std::ostream& out; std::string ans; FileTester(const HiddenInput& hinput_, std::ostream& out_) : hinput(hinput_), out(out_) {} }; int main([[maybe_unused]] int argc, [[maybe_unused]] char** argv) { Timer timer; #ifdef HAVE_OPENCV_HIGHGUI cv::utils::logging::setLogLevel(cv::utils::logging::LogLevel::LOG_LEVEL_SILENT); #endif #if 0 std::ifstream ifs("../../tools/in/seed_0001.txt"); std::istream& in = ifs; std::ofstream ofs("../../tools/out/seed_0001.txt"); std::ostream& out = ofs; //auto tester = FileTester::create(in, out); #else std::istream& in = std::cin; std::ostream& out = std::cout; auto tester = ServerTester::create(in, out); #endif int money = INITIAL_MONEY; std::array L{}; std::array S{}; std::array P{}; std::array R{}; std::array Dmin; std::array Dmax; for (int i = 0; i < 10; i++) { Dmin[i] = 0.5; Dmax[i] = 1.5; } using namespace std; Xorshift rnd; // インタラクティブ開始 for (int t = 1; t <= T; t++) { for (int i = 0; i < N; i++) { L[i] = (int)floor(pow(1.05, 2 * P[i]) * Dmax[i] * Dmax[i] * 0.75 * 0.75 / 0.09) - R[i]; chmax(L[i], 0); } auto res = tester->query(L); auto nmoney = res.money; auto NS = res.S; auto NP = res.P; auto NR = res.R; for (int i = 0; i < N; i++) { if (!NR[i]) { // 全部売れた chmax(Dmin[i], sqrt(NS[i]) / (pow(1.05, P[i]) * 1.25)); } else { double s = sqrt(NR[i] + NS[i]) * pow(1.05, P[i]); chmax(Dmin[i], NS[i] / (s * 1.25)); chmin(Dmax[i], NS[i] / (s * 0.75)); } } money = nmoney; S = NS; P = NP; R = NR; } return 0; }