【CF1344D】Résumé Review【数学】【二分】
題意:給定長度為nnn的自然數序列aaa和整數kkk,要求構造長度為nnn的自然數序列bbb,使得0≤bi≤ai,∑bi=k0\leq b_i\leq a_i,\sum b_i=k0≤bi?≤ai?,∑bi?=k,且∑bi(ai?bi2)\sum b_i(a_i-b_i^2)∑bi?(ai??bi2?)最大化。
n≤105n\leq 10^5n≤105
要求最大化的是b(a?b2)b(a-b^2)b(a?b2),記為f(x)=x(a?x2)f(x)=x(a-x^2)f(x)=x(a?x2)
注意到xxx每增加111,f(x)f(x)f(x)會增加Δi(x)=ai?3x2?3x?1\Delta_i(x)= a_i-3x^2-3x-1Δi?(x)=ai??3x2?3x?1
注意到對于每個iii,當x≥0x\geq 0x≥0時,Δi(x)\Delta_i(x)Δi?(x)單調遞減,題目所求就是所有數中的前kkk大之和,而每個Δi\Delta_iΔi?被取的一定是前面若干個,所以可以套路性地維護當前的xix_ixi?和一個堆,每次取堆頂讓對應的xi+1x_i+1xi?+1
這樣復雜度是O(klog?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?被選了多少個,判斷∑xi<k\sum x_i<k∑xi?<k
最后∑xi\sum x_i∑xi?可能小于kkk,但每個位置最多只能+1+1+1,直接在當前基礎上暴力選k?∑xik-\sum x_ik?∑xi?個
復雜度O(nlog?2V)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; }總結
以上是生活随笔為你收集整理的【CF1344D】Résumé Review【数学】【二分】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 【CF1349D】Slime and B
- 下一篇: 八珍益母丸的功效与作用