#pragma GCC optimize ("Ofast") #include using namespace std; void *wmem; char memarr[96000000]; template inline void walloc1d(T **arr, int x, void **mem = &wmem){ static int skip[16] = {0, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1}; (*mem) = (void*)( ((char*)(*mem)) + skip[((unsigned long long)(*mem)) & 15] ); (*arr)=(T*)(*mem); (*mem)=((*arr)+x); } void sortF_L(int N, unsigned A[], void *mem = wmem){ int i; int m; int bt; unsigned *arr; unsigned c; int *sz; if(N < 256){ sort(A, A+N); return; } bt = sizeof(unsigned) * 8; walloc1d(&arr, N, &mem); walloc1d(&sz, N, &mem); for(m=0;m>m)&255u) ]++; } for(i=(1);i<(257);i++){ sz[i] += sz[i-1]; } for(i=(0);i<(N);i++){ c = ((A[i]>>m)&255u); arr[sz[c]++] = A[i]; } swap(A, arr); } } void sortF_L(int N, int A[], void *mem = wmem){ int i; int x; int y; int z; int *arr; unsigned *send; if(N < 256){ sort(A, A+N); return; } send = (unsigned*)A; sortF_L(N, send, mem); if(A[0] < 0 || A[N-1] >= 0){ return; } x = 0; y = N; while(x < y){ z = (x+y) / 2; if(A[z] < 0){ y = z; } else{ x = z+1; } } walloc1d(&arr, N, &mem); z = 0; for(i=(x);i<(N);i++){ arr[z++] = A[i]; } for(i=(0);i<(x);i++){ arr[z++] = A[i]; } for(i=(0);i<(N);i++){ A[i] = arr[i]; } } void sortF_L(int N, unsigned long long A[], void *mem = wmem){ int i; int m; int bt; unsigned long long *arr; unsigned c; int *sz; if(N < 512){ sort(A, A+N); return; } bt = sizeof(unsigned long long) * 8; walloc1d(&arr, N, &mem); walloc1d(&sz, N, &mem); for(m=0;m>m)&255u) ]++; } for(i=(1);i<(257);i++){ sz[i] += sz[i-1]; } for(i=(0);i<(N);i++){ c = ((A[i]>>m)&255u); arr[sz[c]++] = A[i]; } swap(A, arr); } } void sortF_L(int N, long long A[], void *mem = wmem){ int i; int x; int y; int z; long long *arr; unsigned long long *send; if(N < 512){ sort(A, A+N); return; } send = (unsigned long long*)A; sortF_L(N, send, mem); if(A[0] < 0 || A[N-1] >= 0){ return; } x = 0; y = N; while(x < y){ z = (x+y) / 2; if(A[z] < 0){ y = z; } else{ x = z+1; } } walloc1d(&arr, N, &mem); z = 0; for(i=(x);i<(N);i++){ arr[z++] = A[i]; } for(i=(0);i<(x);i++){ arr[z++] = A[i]; } for(i=(0);i<(N);i++){ A[i] = arr[i]; } } inline void rd(int &x){ int k; int m=0; x=0; for(;;){ k = getchar_unlocked(); if(k=='-'){ m=1; break; } if('0'<=k&&k<='9'){ x=k-'0'; break; } } for(;;){ k = getchar_unlocked(); if(k<'0'||k>'9'){ break; } x=x*10+k-'0'; } if(m){ x=-x; } } inline void wt_L(char a){ putchar_unlocked(a); } inline void wt_L(long long x){ int s=0; int m=0; char f[20]; if(x<0){ m=1; x=-x; } while(x){ f[s++]=x%10; x/=10; } if(!s){ f[s++]=0; } if(m){ putchar_unlocked('-'); } while(s--){ putchar_unlocked(f[s]+'0'); } } template inline S chmax(S &a, T b){ if(a T qselect(int n, int k, T a[]){ int i; int j; T v; T c[64]; while(n > 256){ for(i=(0);i<(64);i++){ c[i] = a[rand()%n]; } sort(c, c+64); v = c[((k*7LL*64LL/n) + 32)/8]; i = -1; j = n; for(;;){ while(a[++i]v); if(i>=j){ break; } swap(a[i],a[j]); } if(k < i){ n = i; } else{ n -= i; k -= i; a += i; } } sort(a, a+n); return a[k]; } int N; int A[10000]; int arr[10000]; long long lf[10001]; long long rg[10001]; int main(){ int K; wmem = memarr; int mx; int s; long long res = 0; rd(N); { int Lj4PdHRW; for(Lj4PdHRW=(0);Lj4PdHRW<(N);Lj4PdHRW++){ rd(A[Lj4PdHRW]); } } int KL2GvlyY = N+1; for(K=(1);K<(KL2GvlyY);K++){ int i; mx = N / K; for(i=(0);i<(mx);i++){ int j; for(j=(0);j<(K);j++){ arr[j] = A[i*K+j]; } lf[i+1] = lf[i] + (long long) K * qselect(K, (K-1)/2, arr); } for(i=(0);i<(mx);i++){ int j; for(j=(0);j<(K);j++){ arr[j] = A[N-1-(i*K+j)]; } sortF_L(K, arr); rg[i+1] = rg[i] + (long long) K * qselect(K, (K-1)/2, arr); } for(i=(0);i<(mx);i++){ chmax(lf[i+1], lf[i]); } for(i=(0);i<(mx);i++){ chmax(rg[i+1], rg[i]); } for(i=(0);i<(mx+1);i++){ chmax(res, lf[i] + rg[mx-i]); } } wt_L(res); wt_L('\n'); return 0; } // cLay varsion 20191012-1 // --- original code --- // template // T qselect(int n, int k, T a[]){ // int i, j; // T v, c[64]; // while(n > 256){ // rep(i,64) c[i] = a[rand()%n]; // sort(c, c+64); // v = c[((k*7LL*64LL/n) + 32)/8]; // // i = -1; j = n; // for(;;){ // while(a[++i]v); // if(i>=j)break; // swap(a[i],a[j]); // } // if(k < i) n = i; // else n -= i, k -= i, a += i; // } // sort(a, a+n); // return a[k]; // } // // int N, A[1d4]; // int arr[1d4]; // ll lf[10001], rg[10001]; // { // int mx, s; // ll res = 0; // rd(N,A(N)); // // REP(K,1,N+1){ // mx = N / K; // rep(i,mx){ // rep(j,K) arr[j] = A[i*K+j]; // lf[i+1] = lf[i] + (ll) K * qselect(K, (K-1)/2, arr); // } // rep(i,mx){ // rep(j,K) arr[j] = A[N-1-(i*K+j)]; // sortF(K, arr); // rg[i+1] = rg[i] + (ll) K * qselect(K, (K-1)/2, arr); // } // rep(i,mx) lf[i+1] >?= lf[i]; // rep(i,mx) rg[i+1] >?= rg[i]; // rep(i,mx+1) res >?= lf[i] + rg[mx-i]; // } // wt(res); // }