結果
問題 | No.919 You Are A Project Manager |
ユーザー |
|
提出日時 | 2019-11-06 20:00:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2,303 ms / 3,000 ms |
コード長 | 2,511 bytes |
コンパイル時間 | 1,089 ms |
コンパイル使用メモリ | 121,888 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-06-27 19:42:00 |
合計ジャッジ時間 | 34,153 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 55 |
ソースコード
#define _USE_MATH_DEFINES#include <cstdio>#include <iostream>#include <sstream>#include <fstream>#include <iomanip>#include <algorithm>#include <cmath>#include <complex>#include <string>#include <vector>#include <array>#include <list>#include <queue>#include <stack>#include <set>#include <map>#include <bitset>#include <numeric>#include <limits>#include <climits>#include <cfloat>#include <functional>#include <iterator>#include <memory>#include <regex>using namespace std;class Xorshift{private:uint32_t x, y, z, w;public:Xorshift(uint32_t seed=88675123, int32_t loop=50){x = 123456789;y = 362436069;z = 521288629;w = seed;while(--loop >= 0){(*this)();}}uint32_t operator()(){uint32_t t=(x^(x<<11));x=y; y=z; z=w;return w=(w^(w>>19))^(t^(t>>8));}uint32_t operator()(uint32_t size){return (*this)() % size;}};Xorshift xorshift;int quickSelect(vector<int>::iterator left, vector<int>::iterator right, int index){int n = right - left;if(n == 1)return *left;swap(*left, *(left + xorshift(n)));auto it1 = left + 1;auto it2 = right - 1;int m1 = 0;int m2 = 0;while(it1 <= it2){if(*it1 < *left || (*it1 == *left && m1 < m2)){++ it1;++ m1;}else{swap(*it1, *it2);-- it2;++ m2;}}swap(*left, *it2);if(index < m1)return quickSelect(left, it2, index);else if(index == m1)return *it2;elsereturn quickSelect(it1, right, index - m1 - 1);}int main(){int n;cin >> n;vector<int> a(n);for(int i=0; i<n; ++i)cin >> a[i];long long ans = 0;for(int len=1; len<=n; ++len){int m = n / len;vector<int> a1 = a;vector<int> a2 = a;vector<long long> v1(m+1), v2(m+1);for(int i=0; i<m; ++i){v1[i+1] = v1[i] + quickSelect(a1.begin() + len * i, a1.begin() + len * (i + 1), (len - 1) / 2);v2[i+1] = v2[i] + quickSelect(a2.end() - len * (i + 1), a2.end() - len * i, (len - 1) / 2);}for(int i=0; i<m; ++i){v1[i+1] = max(v1[i+1], v1[i]);v2[i+1] = max(v2[i+1], v2[i]);}for(int i=0; i<=m; ++i)ans = max(ans, (v1[i] + v2[m-i]) * len);}cout << ans << endl;return 0;}