有一个结论就是每次操作直接取一个存在目前最左边的最小值区间即可。
但是我不会证啊......
大家感性理解。
点击查看代码
/*
Tips:
你数组开小了吗？
你MLE了吗？
你觉得是贪心，是不是该想想dp？
一个小时没调出来，是不是该考虑换题？
打 cf 不要用 umap！！！
记住，rating 是身外之物。
该冲正解时冲正解！
Problem:
算法：
思路：
*/
#include<bits/stdc++.h>
using namespace std;
//#define map unordered_map
#define forl(i,a,b) for(register long long i=a;i<=b;i++)
#define forr(i,a,b) for(register long long i=a;i>=b;i--)
#define forll(i,a,b,c) for(register long long i=a;i<=b;i+=c)
#define forrr(i,a,b,c) for(register long long i=a;i>=b;i-=c)
#define lc(x) x<<1
#define rc(x) x<<1|1
#define mid ((l+r)>>1)
#define cin(x) scanf("%lld",&x)
#define cout(x) printf("%lld",x)
#define lowbit(x) (x&-x)
#define pb push_back
#define pf push_front
#define IOS ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
#define endl '\n'
#define QwQ return 0;
#define ll long long
#define ull unsigned long long
#define lcm(x,y) x/__gcd(x,y)*y
#define Sum(x,y) 1ll*(x+y)*(y-x+1)/2
#define aty cout<<"Yes\n";
#define atn cout<<"No\n";
#define cfy cout<<"YES\n";
#define cfn cout<<"NO\n";
#define xxy cout<<"yes\n";
#define xxn cout<<"no\n";
#define printcf(x) x?cout<<"YES\n":cout<<"NO\n";
#define printat(x) x?cout<<"Yes\n":cout<<"No\n";
#define printxx(x) x?cout<<"yes\n":cout<<"no\n";
ll t;
ll n,m,k;
ll a[100010],b[100010],sq;
struct node{
ll L,R,minn,tag;
}bl[110];
void upd(ll id)
{
forl(i,bl[id].L,bl[id].R)
a[i]+=bl[id].tag;
bl[id].tag=0;
// bl[id].mp.clear();
bl[id].minn=1e18;
forl(i,bl[id].L,bl[id].R)
/*bl[id].mp[a[i]]++,*/bl[id].minn=min(bl[id].minn,a[i]);
}
void add(ll l,ll r,ll x)
{
if(b[l]==b[r])
{
forl(i,l,r)
a[i]+=x;
upd(b[l]);
return ;
}
forl(i,l,bl[b[l]].R)
a[i]+=x;
upd(b[l]);
forl(i,bl[b[r]].L,r)
a[i]+=x;
upd(b[r]);
forl(i,b[l]+1,b[r]-1)
bl[i].tag+=x;
}
ll queryminid(ll l,ll r)
{
if(b[l]==b[r])
{
upd(b[l]);
ll mi=1e18,id=0;
forl(i,l,r)
if(a[i]+bl[b[l]].tag<mi)
mi=a[i]+bl[b[l]].tag,id=i;
return id;
}
ll mi=1e18,id=0;
mi=min(mi,bl[b[l]].minn+bl[b[l]].tag);
mi=min(mi,bl[b[r]].minn+bl[b[r]].tag);
forl(i,b[l]+1,b[r]-1)
mi=min(mi,bl[i].minn+bl[i].tag);
if(bl[b[l]].minn+bl[b[l]].tag==mi)
forl(j,bl[b[l]].L,bl[b[l]].R)
if(a[j]+bl[b[l]].tag==mi)
return j;
forl(i,b[l]+1,b[r]-1)
if(bl[i].minn+bl[i].tag==mi)
forl(j,bl[i].L,bl[i].R)
if(a[j]+bl[i].tag==mi)
return j;
if(bl[b[r]].minn+bl[b[r]].tag==mi)
forl(j,bl[b[r]].L,bl[b[r]].R)
if(a[j]+bl[b[r]].tag==mi)
return j;
}
ll querymin(ll l,ll r)
{
if(b[l]==b[r])
{
upd(b[l]);
ll mi=1e18,id=0;
forl(i,l,r)
if(a[i]+bl[b[l]].tag<mi)
mi=a[i]+bl[b[l]].tag,id=i;
return mi;
}
ll mi=1e18,id=0;
mi=min(mi,bl[b[l]].minn+bl[b[l]].tag);
mi=min(mi,bl[b[r]].minn+bl[b[r]].tag);
forl(i,b[l]+1,b[r]-1)
mi=min(mi,bl[i].minn+bl[i].tag);
return mi;
}
void solve()
{
cin>>n>>m>>k;
forl(i,1,n)
cin>>a[i];
if(n*m<=1e8 && n!=1e5)
{
forl(i,1,m)
{
ll minn=1e18,id=0;
forl(j,1,n)
minn=min(minn,a[j]);
forl(j,1,n)
{
if(a[j]==minn)
{
id=j;
break;
}
}
if(id+k-1<=n)
forl(j,id,id+k-1)
a[j]++;
else
forl(j,n-k+1,n)
a[j]++;
}
ll minn=1e18;
forl(i,1,n)
minn=min(minn,a[i]);
cout<<minn<<endl;
}
else
{
sq=1000;
forl(i,1,n)
{
bl[i].L=bl[i-1].R+1;
bl[i].R=min(n,sq*i);
bl[i].minn=1e18;
forl(j,bl[i].L,bl[i].R)
bl[i].minn=min(bl[i].minn,a[j]),b[j]=i,bl[i].mp[a[j]]++;
if(bl[i].R==n)
break;
}
// cout<<querymin(1,n)<<endl;
forl(i,1,m)
{
ll id=queryminid(1,n);
if(id+k-1<=n)
add(id,id+k-1,1);
else
add(n-k+1,n,1);
}
cout<<querymin(1,n)<<endl;
// forl(i,23223,23227)
// cout<<a[i]<<' ';
// cout<<queryminid(1,n)<<endl;
// cout<<a[queryminid(1,n)]<<endl;
}
}
int main()
{
IOS;
t=1;
// cin>>t;
while(t--)
solve();
/******************/
/*while(L<q[i].l) */
/* del(a[L++]);*/
/*while(L>q[i].l) */
/* add(a[--L]);*/
/*while(R<q[i].r) */
/* add(a[++R]);*/
/*while(R>q[i].r) */
/* del(a[R--]);*/
/******************/
QwQ;
}