結果
問題 |
No.919 You Are A Project Manager
|
ユーザー |
|
提出日時 | 2019-10-26 00:13:47 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
RE
|
実行時間 | - |
コード長 | 2,027 bytes |
コンパイル時間 | 1,369 ms |
コンパイル使用メモリ | 101,464 KB |
実行使用メモリ | 6,948 KB |
最終ジャッジ日時 | 2024-09-13 11:10:23 |
合計ジャッジ時間 | 10,230 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 18 WA * 2 RE * 35 |
コンパイルメッセージ
main.cpp:39:1: warning: ISO C++ forbids declaration of 'main' with no type [-Wreturn-type] 39 | main() | ^~~~
ソースコード
#include<iostream> #include<algorithm> #include<vector> #include<map> using namespace std; //1-indexed #include<vector> template<typename T> struct BIT{ int n; vector<T>bit; BIT(int n_=0):n(n_),bit(n_+1){} T sum(int i) { T ans=0; for(;i>0;i-=i&-i)ans+=bit[i]; return ans; } void add(int i,T a) { if(i==0)return; for(;i<=n;i+=i&-i)bit[i]+=a; } int lower_bound(T k)//k<=sum(ret) { if(k<=0)return 0; int ret=0,i=1; while((i<<1)<=n)i<<=1; for(;i;i>>=1) if(ret+i<=n&&bit[ret+i]<k)k-=bit[ret+=i]; return ret+1; } }; int N,A[10000]; long L[20000],R[20000]; vector<pair<int,int> >Q; vector<pair<int,int> >Ai; map<pair<int,int>,int>mid; main() { cin>>N; for(int i=0;i<N;i++)cin>>A[i]; for(int i=0;i<N;i++)Ai.push_back(make_pair(A[i],i)); for(int K=1;K<=N;K++) { for(int i=0;i+K<=N;i+=K) { Q.push_back(make_pair(i,i+K)); Q.push_back(make_pair(N-i-K,N-i)); } } for(int i=0;i<Q.size();i++)L[i]=-1e9-1,R[i]=1e9; sort(Ai.begin(),Ai.end()); for(int c=0;c<100;c++) { vector<pair<long,int> >M; for(int i=0;i<Q.size();i++) { M.push_back(make_pair((L[i]+R[i])/2,i)); } sort(M.begin(),M.end()); int id=0; BIT<int>P(N); for(int i=0;i<=Ai.size();i++) { long T=i<Ai.size()?Ai[i].first:1e18; while(id<M.size()&&M[id].first<T) { int I=M[id].second; int l=Q[I].first,r=Q[I].second; bool flag=P.sum(r)-P.sum(l)>=(r-l+1)/2; if(flag)R[I]=M[id].first; else L[I]=M[id].first; id++; } if(i<Ai.size())P.add(Ai[i].second+1,1); } } for(int i=0;i<Q.size();i++)mid[Q[i]]=R[i]; long ans=0; for(int K=1;K<=N;K++) { int t=(K-1)/2; vector<long>L(N/K),R(N/K); for(int i=0;i+K<=N;i+=K) { L[i/K]=mid[make_pair(i,i+K)]; R[i/K]=mid[make_pair(N-i-K,N-i)]; } for(int i=1;i<N/K;i++) { L[i]+=L[i-1]; R[i]+=R[i-1]; } for(int i=1;i<N/K;i++) { L[i]=max(L[i],L[i-1]); R[i]=max(R[i],R[i-1]); } for(int i=-1;i<N/K;i++) { int j=N/K-i-2; ans=max(ans,((i<0?0L:L[i])+(j<0?0L:R[j]))*K); } } cout<<ans<<endl; }