1. CF730I Olympiad in Programming and Sports

大意: $n$个人, 第$i$个人编程能力$a_i$, 运动能力$b_i$, 要选出$p$个组成编程队, $s$个组成运动队, 每个队的收益为队员能力和, 求最大收益.

费用流做法很显然, 开两个点$X,Y$表示编程和运动, 源点向每个人连边, 代价为$0$, 每个人向$X$连边, 代价为编程能力, 每个人向$Y$连边, 代价为运动能力, $X$向汇点连边容量为$p$, $Y$向汇点连边, 容量为$s$, 然后求出最大费用最大流即为答案.

考虑用堆模拟费用流. 每次增广只有四种情况.

  • 添加一个未被选择的人去编程队
  • 添加一个未被选择的人去运动队
  • 让一个人从运动队到编程队, 再选择一个未被选择的人去运动队
  • 让一个人从编程队到运动队, 再选择一个未被选择的人去编程队

前两种情况直接用堆维护最大即可. 后两种情况的话, 在每次决策时, 都把替换的收益扔进一个堆里即可.

#include <iostream>

#include <sstream>

#include <algorithm>

#include <cstdio>

#include <cmath>

#include <set>

#include <map>

#include <queue>

#include <string>

#include <cstring>

#include <bitset>

#include <functional>

#include <random>

#define REP(i,a,n) for(int i=a;i<=n;++i)

#define PER(i,a,n) for(int i=n;i>=a;--i)

#define hr putchar(10)

#define pb push_back

#define lc (o<<1)

#define rc (lc|1)

#define mid ((l+r)>>1)

#define ls lc,l,mid

#define rs rc,mid+1,r

#define x first

#define y second

#define io std::ios::sync_with_stdio(false)

#define endl '\n'

#define DB(a) ({REP(__i,1,n) cout<<a[__i]<<' ';hr;})

using namespace std;

typedef long long ll;

typedef pair<int,int> pii;

const int P = 1e9+, INF = 0x3f3f3f3f;

ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}

ll qpow(ll a,ll n) {ll r=%P;for (a%=P;n;a=a*a%P,n>>=)if(n&)r=r*a%P;return r;}

ll inv(ll x){return x<=?:inv(P%x)*(P-P/x)%P;}

inline int rd() {int x=;char p=getchar();while(p<''||p>'')p=getchar();while(p>=''&&p<='')x=x*+p-'',p=getchar();return x;}

//head

const int N = 1e6+;

int n, x, y;

int a[N], b[N], va[N], vb[N];

pii s[];

priority_queue<pii> q[];

ll ans;

pii get(priority_queue<pii> &q, int *v1, int *v2) {

    while (q.size()&&(v1[q.top().y]||v2[q.top().y])) q.pop();

    return q.empty()?pii(-INF,):q.top();

}

void add(int u, int tp) {

    if (va[u]) ans-=a[u],va[u]=,++x;

    if (vb[u]) ans-=b[u],vb[u]=,++y;

    if (tp==) ans+=a[u],va[u]=,--x,q[].push(pii(b[u]-a[u],u));

    else ans+=b[u],vb[u]=,--y,q[].push(pii(a[u]-b[u],u));

}

int main() {

    scanf("%d%d%d", &n, &x, &y);

    REP(i,,n) scanf("%d",a+i);

    REP(i,,n) scanf("%d",b+i);

    REP(i,,n) {

        q[].push(pii(a[i],i));

        q[].push(pii(b[i],i));

    }

    int tot = x+y;

    REP(i,,tot) {

        s[] = get(q[],va,vb);

        s[] = get(q[],va,vb);

        s[] = get(q[],va,va);

        s[] = get(q[],vb,vb);

        s[].x += s[].x;

        s[].x += s[].x;

        if (!x) s[].x=s[].x=-INF;

        if (!y) s[].x=s[].x=-INF;

        int p = max_element(s+,s+)-s;

        if (p==) add(s[].y,);

        else if (p==) add(s[].y,);

        else if (p==) add(s[].y,),add(s[].y,);

        else add(s[].y,),add(s[].y,);

    }

    printf("%lld\n", ans);

    REP(i,,n) if (va[i]) printf("%d ",i);hr;

    REP(i,,n) if (vb[i]) printf("%d ",i);hr;

}

2. NOI2019 序列

大意: 给定两个序列$a,b$, 要求每个序列选出$K$个元素, 且位置相同的元素要至少选$L$个, 求选出元素的最大和.

首先可以得到一个显然的费用流做法.

源点向$a$中每个元素连边, 代价为$a_i$, $b$中每个元素向汇点连边, 代价为$b_i$.

$a$中每个点向对应位置的$b$中的点连边.

再开两个点$X,Y$, $X$向$Y$连容量为$K-L$的边, 表示可以选出$K-L$个不一样的位置.

$a$中每个点再连向$X$, $Y$连向$b$中每个点.

最后求出源点到汇点容量为$K$的最大费用流即为答案.

考虑用堆模拟费用流, 分四种情况.

  • 从$a$中选最大,$b$中也选最大 (要求X到Y仍有剩余流量)
  • $a$,$b$选出位置相同
  • $a$中选最大,$b$中选出之前选过的$a$对应位置
  • $b$中选最大,$a$中选出之前选过的$b$对应位置
#include <iostream>

#include <sstream>

#include <algorithm>

#include <cstdio>

#include <cmath>

#include <set>

#include <map>

#include <queue>

#include <string>

#include <cstring>

#include <bitset>

#include <functional>

#include <random>

#define REP(i,a,n) for(int i=a;i<=n;++i)

#define PER(i,a,n) for(int i=n;i>=a;--i)

#define hr putchar(10)

#define pb push_back

#define lc (o<<1)

#define rc (lc|1)

#define mid ((l+r)>>1)

#define ls lc,l,mid

#define rs rc,mid+1,r

#define x first

#define y second

#define io std::ios::sync_with_stdio(false)

#define endl '\n'

#define DB(a) ({REP(__i,1,n) cout<<a[__i]<<' ';hr;})

using namespace std;

typedef long long ll;

typedef pair<int,int> pii;

const int P = 1e9+, INF = 0x3f3f3f3f;

ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}

ll qpow(ll a,ll n) {ll r=%P;for (a%=P;n;a=a*a%P,n>>=)if(n&)r=r*a%P;return r;}

ll inv(ll x){return x<=?:inv(P%x)*(P-P/x)%P;}

inline int rd() {int x=;char p=getchar();while(p<''||p>'')p=getchar();while(p>=''&&p<='')x=x*+p-'',p=getchar();return x;}

//head

#ifdef ONLINE_JUDGE

const int N = 1e6+;

#else

const int N = ;

#endif

int n, k, l, a[N], b[N];

int va[N], vb[N];

priority_queue<pii> q[];

struct t3 {int x,a,b;} f[];

ll ans;

pii s[];

int cost(int x, int y) {

    int ret = ;

    if (x!=y) {

        if (vb[x]) ++ret; else --ret;

        if (va[y]) ++ret; else --ret;

    }

    return ret<?-:ret>?:;

}

void add(int x, int y) {

    l += cost(x,y);

    ans += a[x]+b[y];

    va[x] = vb[y] = ;

    if (!vb[x]) q[].push(pii(b[x],x));

    if (!va[y]) q[].push(pii(a[y],y));

}

pii get(priority_queue<pii> &q, int *v1, int *v2) {

    while (q.size()&&(v1[q.top().y]||v2[q.top().y])) q.pop();

    return q.empty()?pii(-INF,):q.top();

}

int main() {

    int t;

    scanf("%d", &t);

    while (t--) {

        REP(i,,) while (q[i].size()) q[i].pop();

        scanf("%d%d%d", &n, &k, &l);

        REP(i,,n) scanf("%d",a+i);

        REP(i,,n) scanf("%d",b+i);

        REP(i,,n) {

            va[i] = vb[i] = ;

            q[].push(pii(a[i]+b[i],i));

            q[].push(pii(a[i],i));

            q[].push(pii(b[i],i));

        }

        l = k-l, ans = ;

        REP(i,,k) {

            s[] = get(q[],va,vb);

            s[] = get(q[],va,va);

            s[] = get(q[],vb,vb);

            s[] = get(q[],va,va);

            s[] = get(q[],vb,vb);

            f[] = {s[].x+s[].x,s[].y,s[].y};

            f[] = {s[].x,s[].y,s[].y};

            f[] = {s[].x+s[].x,s[].y,s[].y};

            f[] = {s[].x+s[].x,s[].y,s[].y};

            if (l+cost(f[].a,f[].b)<) f[].x=-INF;

            int p = max_element(f,f+,[](t3 a, t3 b){return a.x<b.x;})-f;

            add(f[p].a,f[p].b);

        }

        printf("%lld\n", ans);

    }

}

3. CF802O April Fools' Problem (hard)

大意: $n$道题, 第$i$天可以花费$a_i$准备一道题, 花费$b_i$打印一道题, 每天最多准备一道, 最多打印一道, 准备的题可以留到以后打印, 求最少花费使得准备并打印$k$道题.

这个题好有意思, 直接用堆很难模拟, 似乎可以用线段树来模拟. 一个非常好写的做法是带权二分,带权二分具体原理其实一直都没搞懂。。。

#include <iostream>

#include <sstream>

#include <algorithm>

#include <cstdio>

#include <cmath>

#include <set>

#include <map>

#include <queue>

#include <string>

#include <cstring>

#include <bitset>

#include <functional>

#include <random>

#define REP(i,a,n) for(int i=a;i<=n;++i)

#define PER(i,a,n) for(int i=n;i>=a;--i)

#define hr putchar(10)

#define pb push_back

#define lc (o<<1)

#define rc (lc|1)

#define ls lc,l,mid

#define rs rc,mid+1,r

#define x first

#define y second

#define io std::ios::sync_with_stdio(false)

#define endl '\n'

#define DB(a) ({REP(__i,1,n) cout<<a[__i]<<' ';hr;})

using namespace std;

typedef long long ll;

typedef pair<int,int> pii;

const int P = 1e9+, INF = 0x3f3f3f3f;

ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}

ll qpow(ll a,ll n) {ll r=%P;for (a%=P;n;a=a*a%P,n>>=)if(n&)r=r*a%P;return r;}

ll inv(ll x){return x<=?:inv(P%x)*(P-P/x)%P;}

inline int rd() {int x=;char p=getchar();while(p<''||p>'')p=getchar();while(p>=''&&p<='')x=x*+p-'',p=getchar();return x;}

//head

const int N = 1e6+;

int n, k, a[N], b[N];

priority_queue<int,vector<int>,greater<int> > q[];

int main() {

    scanf("%d%d", &n, &k);

    REP(i,,n) scanf("%d", a+i);

    REP(i,,n) scanf("%d", b+i);

    int l = , r = 2e9;

    ll ans = ;

    while (l<=r) {

        REP(i,,) while (q[i].size()) q[i].pop();

        ll ret = ;

        int cnt = , mid = (ll)l+r>>;

        REP(i,,n) {

            q[].push(a[i]);

            ll t1 = q[].top()+(ll)b[i]-mid;

            ll t2 = q[].empty()?INF:b[i]+q[].top();

            if (t1<=&&t1<=t2) {

                ret += t1, ++cnt;

                q[].push(-b[i]);

                q[].pop();

            }

            else if (t2<) {

                ret += t2;

                q[].pop();

                q[].push(-b[i]);

            }

        }

        if (cnt>=k) ans=ret+(ll)k*mid,r=mid-;

        else l=mid+;

    }

    printf("%lld\n", ans);

}

4. hdu 6698 Coins

大意: $n$组硬币, 第$i$组有一个$a_i$元的和$b_i$元的, 每组要么不选, 要么选$a_i$, 要么全选, 对于$k\le 2n$, 求出拿$k$个硬币的最大钱数.

堆模拟这四种方案即可:

  • 取一个最大的$a$
  • 补上一个最大的$b$
  • 删掉一个$a$, 取一个最大的$a+b$
  • 删掉一个$b$, 取一个最大的$a+b$
#include <iostream>

#include <sstream>

#include <algorithm>

#include <cstdio>

#include <cmath>

#include <set>

#include <map>

#include <queue>

#include <string>

#include <cstring>

#include <bitset>

#include <functional>

#include <random>

#define REP(i,a,n) for(int i=a;i<=n;++i)

#define PER(i,a,n) for(int i=n;i>=a;--i)

#define hr putchar(10)

#define pb push_back

#define lc (o<<1)

#define rc (lc|1)

#define mid ((l+r)>>1)

#define ls lc,l,mid

#define rs rc,mid+1,r

#define x first

#define y second

#define io std::ios::sync_with_stdio(false)

#define endl '\n'

#define DB(a) ({REP(__i,1,n) cout<<a[__i]<<' ';hr;})

using namespace std;

typedef long long ll;

typedef pair<int,int> pii;

const int P = 1e9+, INF = 0x3f3f3f3f;

ll gcd(ll a,ll b) {return b?gcd(b,a%b):a;}

ll qpow(ll a,ll n) {ll r=%P;for (a%=P;n;a=a*a%P,n>>=)if(n&)r=r*a%P;return r;}

ll inv(ll x){return x<=?:inv(P%x)*(P-P/x)%P;}

inline int rd() {int x=;char p=getchar();while(p<''||p>'')p=getchar();while(p>=''&&p<='')x=x*+p-'',p=getchar();return x;}

//head

const int N = 1e6+;

int n, a[N], b[N];

int va[N], vb[N];

priority_queue<pii> q[];

priority_queue<pii,vector<pii>,greater<pii> > qq[];

pii s[];

void add(int u, int tp) {

    if (tp==) va[u]=,q[].push(pii(b[u],u)),qq[].push(pii(a[u],u));

    else vb[u]=,qq[].push(pii(b[u],u));

}

void work() {

    scanf("%d", &n);

    REP(i,,) while (q[i].size()) q[i].pop();

    REP(i,,) while (qq[i].size()) qq[i].pop();

    REP(i,,n) {

        scanf("%d%d",a+i,b+i);

        q[].push(pii(a[i],i));

        q[].push(pii(a[i]+b[i],i));

        va[i] = vb[i] = ;

    }

    ll ans = ;

    REP(i,,*n) {

        while (q[].size()&&va[q[].top().y]) q[].pop();

        s[] = q[].empty()?pii(-INF,):q[].top();

        while (q[].size()&&(!va[q[].top().y]||vb[q[].top().y])) q[].pop();

        s[] = q[].empty()?pii(-INF,):q[].top();

        while (q[].size()&&va[q[].top().y]) q[].pop();

        s[] = s[] = q[].empty()?pii(-INF,):q[].top();

        if (i==) s[].x = s[].x = -INF;

        else {

            while (qq[].size()&&(!va[qq[].top().y]||vb[qq[].top().y])) qq[].pop();

            if (qq[].size()) s[].x -= qq[].top().x;

            else s[].x = -INF;

            while (qq[].size()&&!vb[qq[].top().y]) qq[].pop();

            if (qq[].size()) s[].x -= qq[].top().x;

            else s[].x = -INF;

        }

        int p = max_element(s,s+)-s;

        ans += s[p].x;

        if (p<=) add(s[p].y,p);

        else {

            //要注意先删除在添加, 不然qq会变

            if (p==) va[qq[].top().y]=,q[].push(qq[].top());

            else vb[qq[].top().y]=,q[].push(qq[].top());

            add(s[p].y,),add(s[p].y,);

        }

        printf("%lld%c",ans," \n"[i==*n]);

    }

}

int main() {

    int t;

    scanf("%d", &t);

    while (t--) work();

}

模拟费用流 & 可撤销贪心的更多相关文章

  1. [NOI2019]序列（模拟费用流）

    题意: 有两个长度为n的序列,要求从每个序列中选k个,并且满足至少有l个位置都被选,问总和最大是多少. \(1\leq l\leq k\leq n\leq 2*10^5\). 首先,记录当前考虑到的位 ...

  2. UOJ #455 [UER #8]雪灾与外卖 (贪心、模拟费用流)

    题目链接 http://uoj.ac/contest/47/problem/455 题解 模拟费用流,一个非常神奇的东西. 本题即为WC2019 laofu的讲课中的Problem 8,经典的老鼠进洞 ...

  3. luogu P5470 [NOI2019]序列 dp 贪心 费用流 模拟费用流

    LINK:序列 考虑前20分 容易想到爆搜. 考虑dp 容易设\(f_{i,j,k,l}\)表示前i个位置 选了j对 且此时A选择了k个 B选择了l个的最大值.期望得分28. code //#incl ...

  4. 贪心（模拟费用流）：NOIP2011 观光公交

    [问题描述] 风景迷人的小城Y 市,拥有n 个美丽的景点.由于慕名而来的游客越来越多,Y 市特意安排了一辆观光公交车,为游客提供更便捷的交通服务.观光公交车在第0 分钟出现在1号景点,随后依次前往2. ...

  5. BZOJ4977[Lydsy1708月赛]跳伞求生——贪心+堆+模拟费用流

    题目链接: 跳伞求生 可以将题目转化成数轴上有$n$个人和$m$个房子,坐标分别为$a_{i}$和$b_{i}$,每个人可以进一个他左边的房子,每个房子只能进一个人.每个房子有一个收益$c_{i}$, ...

  6. 【BZOJ3502/2288】PA2012 Tanie linie/【POJ Challenge】生日礼物 堆+链表（模拟费用流）

    [BZOJ3502]PA2012 Tanie linie Description n个数字,求不相交的总和最大的最多k个连续子序列. 1<= k<= N<= 1000000. Sam ...

  7. 【bzoj1150】[CTSC2007]数据备份Backup 模拟费用流+链表+堆

    题目描述 你在一家 IT 公司为大型写字楼或办公楼(offices)的计算机数据做备份.然而数据备份的工作是枯燥乏味的,因此你想设计一个系统让不同的办公楼彼此之间互相备份,而你则坐在家中尽享计算机游戏 ...

  8. [UOJ455][UER #8]雪灾与外卖——堆+模拟费用流

    题目链接: [UOJ455]雪灾与外卖 题目描述:有$n$个送餐员(坐标为$x_{i}$)及$m$个餐厅(坐标为$y_{i}$,权值为$w_{i}$),每个送餐员需要前往一个餐厅,每个餐厅只能容纳$c ...

  9. Codeforces 280D k-Maximum Subsequence Sum [模拟费用流，线段树]

    洛谷 Codeforces bzoj1,bzoj2 这可真是一道n倍经验题呢-- 思路 我首先想到了DP,然后矩阵,然后线段树,然后T飞-- 搜了题解之后发现是模拟费用流. 直接维护选k个子段时的最优 ...

随机推荐

  1. Spark（四十九）：Spark On YARN启动流程源码分析（一）

    引导: 该篇章主要讲解执行spark-submit.sh提交到将任务提交给Yarn阶段代码分析. spark-submit的入口函数 一般提交一个spark作业的方式采用spark-submit来提交 ...

  2. R获取指定GO term和KEGG pathway的gene list基因集

    clusterProfiler没有显性的接口,但是可以直接扣取clusterProfiler里的函数. 核心函数就是get_GO_data GO_DATA <- get_GO_data(&quo ...

  3. Tosca case status PLANNED,IN-WORK,COMPLETED 对应的图标

    #PLANNED #IN-WORK #COMPLETED

  4. C#多线程和异步——Task和async/await详解

    阅读目录 一.什么是异步 二.Task介绍 1 Task创建和运行 2 Task的阻塞方法(Wait/WaitAll/WaitAny) 3 Task的延续操作(WhenAny/WhenAll/Cont ...

  5. ubuntu系统TCP连接参数优化-TIME_WAIT过多解决办法

    状态:描述CLOSED:无连接是活动的或正在进行LISTEN:服务器在等待进入呼叫SYN_RECV:一个连接请求已经到达,等待确认SYN_SENT:应用已经开始,打开一个连接ESTABLISHED:正 ...

  6. Spring Cloud微服务安全实战- 2-1 环境安装

    下面这些.后续随着讲课逐步再去安装. 2019年1月开始 jdk是收费的 找jdk最后一个免费版本 8u192这是jdk1.8最后的一个免费版本 STS spring提供的ide可以方便的开发spri ...

  7. 【429】关于ADT的访问权限

    在看老师代码的时候,发现ADT中的 struct 有时候写到了 adt.c 里面,有时候写到了 adt.h 里面,其实有些困惑,经过仔细研究,发现写在 adt.h 中的 struct 可以在 test ...

  8. java 时间的原生操作和工具类操作

    package com.xc.test.dateoperation; import org.apache.commons.lang3.time.DateFormatUtils; import org. ...

  9. 【源码解析】Flink 滑动窗口数据分配到多个窗口

    之前一直用翻滚窗口,每条数据都只属于一个窗口,所有不需要考虑数据需要在多个窗口存的事情. 刚好有个需求,要用到滑动窗口,来翻翻 flink 在滑动窗口中,数据是怎么分配到多个窗口的 一段简单的测试代码 ...

  10. js文档系统-jsdoc-docdash

    一.参考文档 模版:https://github.com/clenemt/docdash 例子:http://clenemt.github.io/docdash/index.html jsdoc:ht ...