题意:给定长度为 n n n的自然数序列 a a a和整数 k k k,要求构造长度为 n n n的自然数序列 b b b,使得 0 ≤ b i ≤ a i , ∑ b i = k 0\leq b_i\leq a_i,\sum b_i=k 0≤bi≤ai,∑bi=k,且 ∑ b i ( a i − b i 2 ) \sum b_i(a_i-b_i^2) ∑bi(ai−bi2)最大化。
n ≤ 1 0 5 n\leq 10^5 n≤105
要求最大化的是 b ( a − b 2 ) b(a-b^2) b(a−b2),记为 f ( x ) = x ( a − x 2 ) f(x)=x(a-x^2) f(x)=x(a−x2)
注意到 x x x每增加 1 1 1, f ( x ) f(x) f(x)会增加 Δ i ( x ) = a i − 3 x 2 − 3 x − 1 \Delta_i(x)= a_i-3x^2-3x-1 Δi(x)=ai−3x2−3x−1
注意到对于每个 i i i,当 x ≥ 0 x\geq 0 x≥0时, Δ i ( x ) \Delta_i(x) Δi(x)单调递减,题目所求就是所有数中的前 k k k大之和,而每个 Δ i \Delta_i Δi被取的一定是前面若干个,所以可以套路性地维护当前的 x i x_i xi和一个堆,每次取堆顶让对应的 x i + 1 x_i+1 xi+1
这样复杂度是 O ( k log n ) O(k\log n) O(klogn)无法通过
但是每次加上的 Δ i ( x ) \Delta_i(x) Δi(x)是单调递减的(废话),所以可以二分最小的 d = Δ i ( x ) d=\Delta_i(x) d=Δi(x),再二分出每个 Δ i \Delta_i Δi被选了多少个,判断 ∑ x i < k \sum x_i<k ∑xi<k
最后 ∑ x i \sum x_i ∑xi可能小于 k k k,但每个位置最多只能 + 1 +1 +1,直接在当前基础上暴力选 k − ∑ x i k-\sum x_i k−∑xi个
复杂度 O ( n log 2 V ) O(n\log^2V) O(nlog2V)
#include <iostream>
#include <cctype>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define MAXN 100005
using namespace std;
typedef long long ll;
const ll INF=4e18;
int a[MAXN],b[MAXN],x[MAXN],p[MAXN],n;
ll k,sum;
inline ll delta(const int& a,const int& x){return x==a? -INF:a-3ll*x*x-3ll*x-1;}
inline bool cmp(int A,int B){return delta(a[A],x[A])>delta(a[B],x[B]);}
inline bool check(ll v)
{
for (int i=1;i<=n;i++)
{
if (delta(a[i],0)<v) x[i]=0;
else
{
int l=0,r=a[i]-1,mid;
while (l<r)
{
mid=(l+r+1)>>1;
if (delta(a[i],mid)<v) r=mid-1;
else l=mid;
}
x[i]=l+1;
}
}
sum=0;
for (int i=1;i<=n;i++) sum+=x[i];
if (sum>k) return false;
for (int i=1;i<=n;i++) b[i]=x[i];
return true;
}
int main()
{
cin>>n>>k;
for (int i=1;i<=n;i++) scanf("%d",&a[i]);
ll l=-INF,r=INF,mid,ans;
while (l<r)
{
mid=(l+r)/2;
if (check(mid)) ans=mid,r=mid-1;
else l=mid+1;
}
check(ans);
for (int i=1;i<=n;i++) p[i]=i;
sort(p+1,p+n+1,cmp);
for (int i=1;i<=k-sum;i++) ++b[p[i]];
for (int i=1;i<=n;i++) printf("%d ",b[i]);
return 0;
}