連結(luogu)
連結(codeforces)
有一個結論就是每次操作直接取一個存在目前最左邊的最小值區間即可。
但是我不會證啊......
大家感性理解。
點選檢視程式碼
/*
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;
}