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. SpringBoot访问不了JSP但却能进入后台

    博主在使用sring-boot跳转HTML页面后,由于好奇心就想跳转到JSP页面,就在网上搜相关信息,结果不是跳转500错误就是下载JSP文件.各种坑啊,在博主跳了N多坑后,终于跳转JSP页面成功.故 ...

  2. Excel 如何自动调整列宽？

      excel如何自动调整列宽 1.打开Excel表格,选中要调整的表格. 2.点击"格式",选择"自动调整列宽",右键点击"设置单元格格式" ...

  3. 单例模式：Java单例模式的几种写法及它们的优缺点

    总结下Java单例模式的几种写法: 1. 饿汉式 public class Singleton { private static Singleton instance = new Singleton( ...

  4. Oracle 查询表分区相关信息

    Oracle 查询表分区相关信息 --表分区 --1,分区表信息 -- (1)显示数据库所有分区表的信息 select * from DBA_PART_TABLES a where a.owner=u ...

  5. 小D课堂 - 零基础入门SpringBoot2.X到实战_第8节 数据库操作之整合Mybaties和事务讲解_32..SpringBoot2.x持久化数据方式介绍

    笔记 1.SpringBoot2.x持久化数据方式介绍          简介:介绍近几年常用的访问数据库的方式和优缺点 1.原始java访问数据库             开发流程麻烦        ...

  6. mysql的备份与还原，安装(window)

    之前一直使用的navicat的自动的备份功能.但是在使用navicat的还原功能时,贼慢. 今天要做的是window服务器上将mysql单个数据库备份.然后还原到我的mac笔记本上. 本来以为很快的, ...

  7. linux安装6.5.3版本elastic search

    到官网https://www.elastic.co/cn/downloads/elasticsearch下载压缩包,目前最新的版本是7.3.2,我想下6.5.3,点击下面的past release链接 ...

  8. Qt widget中使用QML自定义电池

    1.效果 2.QML 在资源里新建Mybattery.qml: import QtQuick 2.0 import QtQuick 2.12 Item { id: root property colo ...

  9. 转 linux下ClamAV使用

    linux下ClamAV使用 第一步:Clamav下载http://www.clamav.net/downloads#yuminstall wget –y第二步:创建clamav用户和组groupad ...

  10. Apache使用内置插件mod_php解析php的配置

    apache安装完毕之后,修改httpd.conf配置文件,添加代码如下: LoadModule php5_module modules/libphp5.so <FilesMatch \.php ...