贪心构造专题

CF1530D Secret Santa

给你一个序列 \(a\),满足 \(a_i\neq i\),每次操作可以将任意位置上的数改为 \([1,n]\) 上的数,求将此序列变为一个 \(1\sim n\) 的排列且 \(a_i\neq i\) 的最小操作次数,输出一种方案。

\(n\le 2\times 10^5\)

贪心构造,首先次数很简单,就把所有不合法的位置搞掉即可。对于不合法的位置,将可以填上去的数集 \(S\) 倒序放在每个位置上,然后依次检查每个位置,若出现 \(a_i=i\) 的情况,就和同类型的,第一个未修改的位置交换,可以保证合法。时间复杂度 \(O(n\log n)\)

code
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+3;
int T,n,a[maxn],siz[maxn],vis[maxn];
vector<int>v[maxn];
vector<int>t1,s;
int r1[maxn];
void solve(){
    cin>>n;
    for(int i=1;i<=n;i++){
        vis[i]=0;
        cin>>a[i];
        siz[a[i]]++;
        v[a[i]].emplace_back(i);
    }
    for(int i=1;i<=n;i++){
        if(siz[i]){
            vis[i]=1;
            if(siz[i]>1){
                for(int j=1;j<siz[i];j++){
                    t1.emplace_back(v[i][j]);
                }
                r1[v[i][0]]=i;
            }else{
                if(v[i][0]==i){
                    t1.emplace_back(v[i][0]);
                }else{
                    r1[v[i][0]]=i;
                }
            }
        }
    }
    sort(t1.begin(),t1.end());
    for(int i=1;i<=n;i++){
        if(!vis[i]){
            s.emplace_back(i);
        }
    }
    int len=s.size();
    for(int i=0;i<len;i++){
        r1[t1[i]]=s[len-i-1];
    }
    int f1=0;
    for(int i=1;i<=n;i++){
        if(r1[i]==i){
            swap(r1[i],r1[v[a[i]][0]]);
        }
    }
    if(!f1){
        cout<<n-len<<'\n';
        for(int i=1;i<=n;i++){
            cout<<r1[i]<<" \n"[i==n];
        }
    }
    for(int i=1;i<=n;i++){
        v[i].clear();
        siz[i]=0;
        r1[i]=0;
    }
    s.clear();
    t1.clear();
}
signed main(){
    cin>>T;
    while(T--){
        solve();
    }
    return 0;
}

ABC137D Summer Vacation

给你 \(n\) 个工作,\((a_i,b_i)\) 表示工作 \(i\) 需要花费 \(a_i\) 天,收益为 \(b_i\),每天只能做一个,求 \(m\) 天内得到的最大收益。

\(n,m\le 10^5\)

考虑带悔贪心,先按 \(a_i\) 为第一关键字,\(b_i\) 为第二关键字排序,若 \(m-cnt\ge a_i\) 则把 \(i\) 扔进堆里,否则就考虑把堆里最小的元素扔掉直到满足条件,再把 \(i\) 放进堆里,每步求最优答案,时间复杂度 \(O(n\log n)\)

code
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+3;
int n,m;
struct node{
    int a,b;
    bool operator<(const node o)const{
        if(a!=o.a) return a>o.a;
        return b<o.b;
    }
}a[maxn];
struct nodee{
    int a,b;
    bool operator<(const nodee o)const{
        if(b!=o.b) return b>o.b;
        return a>o.a;
    }
};
priority_queue<nodee>q;
int ans,sum,cnt;
signed main(){
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        cin>>a[i].a>>a[i].b;
    }
    sort(a+1,a+n+1);
    for(int i=1;i<=n;i++){
        int tot=0;
        if(m-cnt>=a[i].a){
            q.push({a[i].a,a[i].b});
            cnt++;
            sum+=a[i].b;
            ans=max(ans,sum);
        }else{
            while(!q.empty()&&m-cnt<a[i].a&&tot+q.top().b<=a[i].b){
                tot+=q.top().b;
                sum-=q.top().b;
                cnt--;
                q.pop();
            }
            if(m-cnt>=a[i].a){
                q.push({a[i].a,a[i].b});
                cnt++;
                sum+=a[i].b;
            }
            ans=max(ans,sum);
        }
    }

    cout<<ans<<'\n';
    return 0;
}

P4377 [USACO18OPEN] Talent Show G

给你一序列 \((w_i,t_i)\),求一种 \(c_i\in\{0,1\}\) 序列,满足 \(R=\frac{\sum c_it_i}{\sum c_iw_i}\) 最大且分母不小于 \(W\),输出 \(1000\times R_{\max}\)

\(n\le 250,W\le 1000\)

分数规划板子题。

分数规划具有单调性,于是我们可以二分一个 \(x\),判断 \(R\ge x\) 即合法。拆式子

\[\begin{aligned}\sum c_it_i & \ge \sum c_ixw_i\\\sum c_i(t_i-xw_i) & \ge 0\end{aligned} \]

满足上式即可,记 \(d_i=t_i-xw_i\)。对于没有 \(W\) 的限制,可以直接贪心选取 \(d_i>0\) 的位置;而有了 \(W\) 的限制,就使用背包,\(w_i\) 为重量,\(d_i\) 为价值,\(f_W>0\) 即合法,时间复杂度 \(O(nW\log nW)\)

code
#include<bits/stdc++.h>
using namespace std;
const int maxn=253;
const double eps=1e-6;
int n,W,a[maxn],b[maxn];
double f[maxn<<3];
bool check(double x){
    for(int i=1;i<=W;i++) f[i]=-1e9;
    for(int i=1;i<=n;i++){
        for(int j=W;j>=0;j--){
            int k=min(W,j+a[i]);
            f[k]=max(f[k],f[j]+b[i]-x*a[i]);
        }
    }
    return f[W]>0;
}
signed main(){
    cin>>n>>W;
    for(int i=1;i<=n;i++){
        cin>>a[i]>>b[i];
    }
    double l=0,r=250000,ans;
    while(l+eps<=r){
        double mid=(l+r)/2;
        if(check(mid)){
            l=mid;
        }else{
            r=mid;
        }
    }
    cout<<int(l*1000);
    return 0;
}
posted @ 2024-09-26 17:35  view3937  阅读(6)  评论(0编辑  收藏  举报
Title