QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#292730 | #7417. Honorable Mention | DaiRuiChen007 | WA | 1101ms | 33728kb | C++17 | 2.2kb | 2023-12-28 12:07:07 | 2023-12-28 12:07:08 |
Judging History
answer
#include<bits/stdc++.h>
#define ll long long
#define sz(v) int(v.size())
using namespace std;
const int MAXN=35005;
const ll inf=2e9;
typedef vector<ll> poly;
poly operator *(const poly &f,const poly &g) {
if(f.empty()||g.empty()) return {};
poly h(sz(f)+sz(g)-1);
int i=1,j=1,k=1; h[0]=f[0]+g[0];
while(i<sz(f)&&j<sz(g)) {
if(f[i]-f[i-1]>g[j]-g[j-1]) h[k]=h[k-1]+f[i]-f[i-1],++i,++k;
else h[k]=h[k-1]+g[j]-g[j-1],++j,++k;
}
while(i<sz(f)) h[k]=h[k-1]+f[i]-f[i-1],++i,++k;
while(j<sz(g)) h[k]=h[k-1]+g[j]-g[j-1],++j,++k;
return h;
}
int n,q,a[MAXN];
poly tr[MAXN<<2][2][2],f0,f1;
void build(int l=1,int r=n,int p=1) {
if(l==r) {
tr[p][0][0]={0,a[l]},tr[p][0][1]=tr[p][1][0]=tr[p][1][1]={-inf,a[l]};
return ;
}
int mid=(l+r)>>1;
build(l,mid,p<<1),build(mid+1,r,p<<1|1);
for(int i:{0,1}) for(int j:{0,1}) {
f0=tr[p<<1][i][0]*tr[p<<1|1][0][j];
f1=tr[p<<1][i][1]*tr[p<<1|1][1][j];
tr[p][i][j].resize(r-l+2);
for(int k=0;k<=r-l;++k) tr[p][i][j][k]=max(f0[k],f1[k+1]);
tr[p][i][j][r-l+1]=f0[r-l+1];
}
}
struct info {
ll v; int c;
friend info operator +(info x,info y) { return {x.v+y.v,x.c+y.c}; }
friend bool operator <(info x,info y) { return x.v==y.v?x.c<y.c:x.v<y.v; }
};
typedef array<array<info,2>,2> mat;
info eval(poly &f,ll d) {
int l=1,r=sz(f)-1,p=0;
while(l<=r) {
int mid=(l+r)>>1;
if(f[mid]-f[mid-1]>=d) l=mid+1,p=mid;
else r=mid-1;
}
return {f[p]-d*p,p};
}
mat qry(int ul,int ur,ll d,int l=1,int r=n,int p=1) {
if(ul<=l&&r<=ur) {
mat f;
for(int i:{0,1}) for(int j:{0,1}) f[i][j]=eval(tr[p][i][j],d);
return f;
}
int mid=(l+r)>>1;
if(ur<=mid) return qry(ul,ur,d,l,mid,p<<1);
if(mid<ul) return qry(ul,ur,d,mid+1,r,p<<1|1);
mat f=qry(ul,ur,d,l,mid,p<<1),g=qry(ul,ur,d,mid+1,r,p<<1|1),h;
for(int i:{0,1}) for(int j:{0,1}) {
h[i][j]=max(f[i][0]+g[0][j],f[i][1]+g[1][j]+info{d,-1});
}
return h;
}
signed main() {
scanf("%d%d",&n,&q);
for(int i=1;i<=n;++i) scanf("%d",&a[i]);
build();
for(int l,r,k;q--;) {
scanf("%d%d%d",&l,&r,&k);
ll L=-inf,R=inf,d=R;
while(L<=R) {
ll mid=(L+R)>>1;
if(qry(l,r,mid)[0][0].c<=k) R=mid-1,d=mid;
else L=mid+1;
}
printf("%lld\n",qry(l,r,d)[0][0].v+k*d);
}
return 0;
}
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 17012kb
input:
5 5 -1 2 -3 4 -5 1 5 1 1 5 2 1 5 3 1 5 4 1 5 5
output:
4 6 5 2 -3
result:
ok 5 number(s): "4 6 5 2 -3"
Test #2:
score: 0
Accepted
time: 0ms
memory: 16768kb
input:
5 1 7 7 7 7 7 1 5 1
output:
35
result:
ok 1 number(s): "35"
Test #3:
score: -100
Wrong Answer
time: 1101ms
memory: 33728kb
input:
25000 25000 -11889 -16468 18402 -13238 14125 12120 -3295 -12106 20152 7478 8575 9968 -13946 -13408 23715 -14670 -21837 12589 7029 20609 2868 12728 -6942 982 5884 -7561 -2058 -7892 7806 16344 14645 13904 7355 9860 15446 -539 24072 21969 19028 -24939 -21911 -11836 21045 -4727 -3882 2138 18513 -4796 -1...
output:
10342053 44514177 6687268 77971944 99605102 107726056 15473041 17383093 62018237 10703102 41214044 22949449 9407209 9022260 39351814 72311293 5178065 42027491 52701707 38135774 37045964 4761513 16326256 16812497 107987457 28306484 46710957 864270 102812861 27960495 50366764 16379719 2791377 21113338...
result:
wrong answer 1st numbers differ - expected: '10341941', found: '10342053'