BZOJ5251: [2018多省省队联测]劈配-网络流

最新推荐文章于 2019-11-06 13:32:17 发布
最新推荐文章于 2019-11-06 13:32:17 发布
阅读量290 收藏
点赞数
分类专栏： 网络流
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Charlie_jilei/article/details/79866149
版权

传送门

题意：

有n位选手和m位导师，每位导师战队有容量限制 bi b i 。

每位选手有一张志愿表。对于每位选手，导师之间允许并列（最多允许C位导师并列）。所有选手有一个排名，选手之间不允许并列。

导师录取的规则可以简单概括为：对每位选手，在优先满足更高位选手的前提下，尽可能满足该选手最高位志愿。

每位选手有一个理想的志愿s_i（梦想），表示他希望被s_i或更高位的志愿录取。

求：

1.每位选手的最终录取情况。

2.每位选手至少要上升几名，才能被理想的志愿录取（实现梦想）。

n,m≤200 C≤10 n , m ≤ 200   C ≤ 10

Solution：

考场上只想到了费用流，但是T掉了，自测是跑了3s，在此分享一下费用流做法：

因为我们每次需要优先满足排名在前的选手，志愿在前的导师，所以说我们S向每个选手连费用为1000*排名，流量为1的边，每个导师向T连费用为0，流量为人数上限的边，每个选手向导师连费用为志愿优先级，流量为1的边，这样对于每次跑出的费用V, V1000 V 1000 是这个选手的排名， V%1000 V % 1000 是这个选手能满足的最高志愿。

但是这样直接跑是错的（过不了大样例），我们需要每次重新建图，对于每个已知结果的选手，我们连边的时候只能连和结果志愿一样的边，理由大概就是后面的流在跑的时候可能会通过前面建的回流来影响前面的结果

这样我们需要枚举每个选手，然后每次重新建图+费用流

至于第二问只需要二分一下就好了

复杂度大概是 O(n∗n∗n∗c) O ( n ∗ n ∗ n ∗ c ) （费用流复杂度大概是 n∗n∗c n ∗ n ∗ c ）

代码（TLE）：

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int N=410;
const int inf=1e9;
int n,m,t,C,S,T,sz,head[N],ans[N],a[300][300],b[300],dis[N],q[N*N],s[N];
bool vis[N];
int pre[N],prv[N];
struct edg{
    int to,next,v,f;
}e[200*N];
void add(int x,int y,int f,int v)
{
    sz++;e[sz].to=y;e[sz].f=f;e[sz].v=v;e[sz].next=head[x];head[x]=sz;
    sz++;e[sz].to=x;e[sz].f=0;e[sz].v=-v;e[sz].next=head[y];head[y]=sz;
}
bool SPFA()
{
    int tt=0,hh=0;
    for (int i=S;i<=T;++i) vis[i]=0,dis[i]=inf;
    vis[S]=1;dis[S]=0;q[++tt]=S;
    while (hh<tt)
    {
        int x=q[++hh];vis[x]=0;
        for (int i=head[x];i;i=e[i].next)
        {
            int y=e[i].to;
            if (e[i].f&&dis[y]>dis[x]+e[i].v)
            {
                dis[y]=dis[x]+e[i].v;
                pre[y]=i;prv[y]=x;
                if (!vis[y]) q[++tt]=y,vis[y]=1;
            }
        }
    }
    return (dis[T]!=inf); 
}
int main()
{
    freopen("mentor.in","r",stdin);
    freopen("mentor.out","w",stdout);
    scanf("%d%d",&t,&C);
    while (t--)
    {
        sz=1;memset(head,0,sizeof(head));
        scanf("%d%d",&n,&m);S=0;T=n+m+1;
        for (int i=1;i<=n;++i) ans[i]=m+1;
        for (int i=1;i<=m;++i) scanf("%d",&b[i]);
        for (int i=1;i<=n;++i)
        {
            sz=1;memset(head,0,sizeof(head));
            for (int j=1;j<=n;++j) add(S,j,1,j*1000);
            for (int j=1;j<=m;++j) add(n+j,T,b[j],0);
            for (int j=1;j<i;++j)
                for (int k=1;k<=m;++k)
                    if (a[j][k]==ans[j]) add(j,n+k,1,a[j][k]);
            for (int j=1;j<=m;++j)
            {
                scanf("%d",&a[i][j]);
                if (a[i][j]) add(i,n+j,1,a[i][j]);
            }
            while (SPFA())
            {
                ans[dis[T]/1000]=dis[T]%1000;
                for (int i=T;i!=S;i=prv[i])
                    e[pre[i]].f--,e[pre[i]^1].f++;
            }
        }
        for (int x,i=1;i<=n;++i)
        {
            scanf("%d",&x);s[i]=0;
            if (x>=ans[i]) continue;
            int l=1,r=i-1;s[i]=i;
            while (l<=r)
            {
                sz=1;memset(head,0,sizeof(head));
                for (int j=1;j<=n;++j) add(S,j,1,j*1000);
                for (int j=1;j<=m;++j) add(n+j,T,b[j],0);
                int mid=l+r>>1;
                for (int j=1;j<mid;++j)
                    for (int k=1;k<=m;++k)
                        if (a[j][k]==ans[j]) add(j,n+k,1,a[j][k]);
                for (int j=1;j<=m;++j)
                    if (a[i][j]) add(i,n+j,1,a[i][j]);
                bool p=0;
                while (SPFA())
                {
                    if (dis[T]/1000==i)
                    {
                        p=1;
                        if (dis[T]%1000<=x) {s[i]=i-mid;l=mid+1;break;}
                        else {r=mid-1;break;} 
                    }
                    //cout<<dis[T]<<endl; 
                    for (int i=T;i!=S;i=prv[i])
                        e[pre[i]].f--,e[pre[i]^1].f++;
                }
                if (!p) r=mid-1;
            }
        } 
        for (int i=1;i<=n;++i)
            printf("%d ",ans[i]);
        printf("\n");
        for (int i=1;i<=n;++i)
            printf("%d ",s[i]);
        printf("\n");
    }
}

普通网络流思路和费用流差不多，只需要对于每个选手枚举一个志愿即可了，对于每个选手只需要重新建边一次，枚举志愿时可以直接加边

第二问也是二分

复杂度 O(n∗m∗nc−−√) O ( n ∗ m ∗ n c ) ( nc−−√ n c 是二分图dinic复杂度)

代码：

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int N=410;
const int inf=1e9;
int n,m,t,C,S,T,sz,head[N],v[N],ans[N],a[300][300],b[300],dis[N],q[N],s[N];
int dep[N];
int pre[N],prv[N];
struct edg{
    int to,next,f;
}e[200*N];
void add(int x,int y,int f)
{
    sz++;e[sz].to=y;e[sz].f=f;e[sz].next=head[x];head[x]=sz;
    sz++;e[sz].to=x;e[sz].f=0;e[sz].next=head[y];head[y]=sz;
}
bool bfs()
{
    int hh=0,tt=0;
    for (int i=S;i<=T;i++) dep[i]=0;
    dep[S]=1;q[++tt]=S;
    while (hh<tt)
    {
        int x=q[++hh];
        for (int i=head[x];i;i=e[i].next)
        {
            int y=e[i].to;
            if (e[i].f&&!dep[y]) dep[y]=dep[x]+1,q[++tt]=y;
        }
    }
    return (dep[T]!=0);
}
int dfs(int x,int flow)
{
    if (x==T) return flow;
    int used=0,now=0;
    for (int i=v[x];i;i=e[i].next)
    {
        int y=e[i].to;
        if (dep[y]==dep[x]+1&&e[i].f)
        {
            now=dfs(y,min(flow-used,e[i].f));
            e[i].f-=now;
            e[i^1].f+=now;
            used+=now;
            v[x]=i;
            if (flow==used) break;
        }
    }
    if (used==0) dep[x]=-1;
    return used;
}
int main()
{
//  freopen("mentor.in","r",stdin);
//  freopen("mentor.out","w",stdout);
    scanf("%d%d",&t,&C);
    while (t--)
    {
        sz=1;memset(head,0,sizeof(head));
        scanf("%d%d",&n,&m);S=0;T=n+m+1;
        for (int i=1;i<=n;++i) ans[i]=m+1;
        for (int i=1;i<=m;++i) scanf("%d",&b[i]);
        int tag=0;
        for (int i=1;i<=n;++i)
        {
            sz=1;memset(head,0,sizeof(head));
            for (int j=1;j<=i;++j) {if (j==i) tag=sz+1;add(S,j,1);}
            for (int j=1;j<=m;++j) add(n+j,T,b[j]);
            for (int j=1;j<i;++j)
                for (int k=1;k<=m;++k)
                    if (a[j][k]==ans[j]) add(j,n+k,1);
            for (int j=1;j<=m;++j)
                scanf("%d",&a[i][j]);
            while (bfs())
            {
                for (int j=S;j<=T;++j) v[j]=head[j]; 
                dfs(S,inf);
            }
            for (int l=1;l<=m;++l)
            {
                for (int j=1;j<=m;++j)
                    if (a[i][j]==l) add(i,n+j,1);
                while (bfs())
                {
                    for (int j=S;j<=T;++j) v[j]=head[j]; 
                    dfs(S,inf);
                }
                if (e[tag].f==0) {ans[i]=l;break;}
            }
        }
        for (int x,i=1;i<=n;++i)
        {
            scanf("%d",&x);s[i]=0;
            if (x>=ans[i]) continue;
            int l=1,r=i-1;s[i]=i;
            while (l<=r)
            {
                sz=1;memset(head,0,sizeof(head));
                for (int j=1;j<=n;++j) {if (j==i) tag=sz+1;add(S,j,1);}
                for (int j=1;j<=m;++j) add(n+j,T,b[j]);
                int mid=l+r>>1;
                for (int j=1;j<mid;++j)
                    for (int k=1;k<=m;++k)
                        if (a[j][k]==ans[j]) add(j,n+k,1);
                while (bfs())
                {
                    for (int j=S;j<=T;++j) v[j]=head[j]; 
                    dfs(S,inf);
                }
                for (int j=1;j<=m;++j)
                    if (a[i][j]>0&&a[i][j]<=x) add(i,n+j,1);

                while (bfs())
                {
                    for (int j=S;j<=T;++j) v[j]=head[j]; 
                    dfs(S,inf);
                }
                if (e[tag].f==0) {l=mid+1,s[i]=i-mid;}
                else r=mid-1;
            }
        } 
        for (int i=1;i<=n;++i)
            printf("%d ",ans[i]);
        printf("\n");
        for (int i=1;i<=n;++i)
            printf("%d ",s[i]);
        printf("\n");
    }
}
Fizzmy
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
wangjunrui666#bzoj-problem##5218. [Lydsy2017省队十连测]友好城市1
07-25
对于两座不同的城市a 和b，如果a 能通过这些单向道路直接或间接到达b，且b 也能如此到达a，那么它们就会被认为是一对友好城市。eland 的交通系统十分特殊，
【BZOJ5251】【2018多省省队联测】劈配
cz_xuyixuan的博客
04-09 560
【题目链接】点击打开链接【思路要点】按照题意建图，跑网络流即可。由于\(C\)的限制，图中点数为\(O(N+M)\)级别，边数为\(O(N*C)\)级别的。并且，每一次运行Dinic算法只会增广一条路径，因此，时间复杂度为\(O(T*C*N^3)\)(\(N\)、\(M\)同阶)。【代码】#include&lt;bits/stdc++.h&gt; using namespace std; cons...
Bzoj5251: [2018多省省队联测]劈配
Cyhlnj
04-09 281
题面 传送门 Sol 网络流辣 枚举每个点的最优匹配，然后只建最优匹配的边跑网络流 下一个点的图中，保证上一个点只连了最优匹配的边 直接整个图复制过来 然后二分答案一个点向上跳多少 和上面一样建图checkcheckcheck # include &lt;bits/stdc++.h&gt; # define RG register # define IL inline #...
bzoj5251 [2018多省省队联测]劈配（网络流+贪心+二分）
Icefox的博客
04-09 219
我们就贪心地，对于每个人，按第一志愿第二志愿。。。试图匹配，直到能满流。 然后对于第二问我们二分一下。
5251: [2018多省省队联测]劈配
Tgotp的博客
04-11 219
简单观察即可发现是个网络流。 然后建边很好弄。 考虑第一个询问实际上就只用枚举一下当前放第几志愿能跑出流量。 第二个询问实际上就是重新建图跑流量。 发现直接不是很好求出答案， 那么发现可以在跑第i个人之前，处理一下哪些点还可以放人，然后更新答案即可 c++代码如下： #include&lt;bits/stdc++.h&gt; #define rep(i,x,y) for(regist...
[八省联考2018]劈配【动态加边/删边 网络流】【二分答案】
gigo_64 sink in my world
11-06 324
传送门 md神题，， 大致题意： n个人，有m个志愿档，有m个老师。每个老师有最大限额学生。每个学生的每个志愿档可以指定0~c个老师。 从第一个学生开始，选择合理的，最高志愿。 对于第i个人来说，目标仅仅是指1~i个人全都被可能的最高志愿档次录取，跟后面的人无关。 原题是这样说的：如果一种方案满足“前n名的录取结果最优”，那么我们可以简称这种方案是最优的。 两个子问题。 子问题...
wangjunrui666#bzoj-problem##1596. [Usaco2008 Jan]电话网络1
07-25
第1行: 1个整数，N第2..N行: 每行为2个用空格隔开的整数A、B，为两块相邻草地的编号第1行: 输出1个整数，即FJ最少建立无线电通讯塔的数目输入说明:F
wangjunrui666#bzoj-problem##5007. TCP协议1
07-25
1.将所有之前执行的时间戳大于t的push和pop指令全部撤销 2.执行当前这条指令 3.按时间戳顺序重新执行在第1步被撤销的指令
Node.js-ZOJCH是BZOJ题库的离线版
08-09
ZOJCH是BZOJ题库的离线版
BZOJ_problem_rank:BZOJ问题等级按交流比排序
05-25
#BZOJ Problem Rankrank.cpp 程序文件data.dat bzoj题库数据done.dat AC过的题，初始可以把所有A过的题粘进去，正常退出的话自动维护。black.dat 黑名单。选题时会跳过。错题、神题、没题面、不想做等等。//Thank ...
[BZOJ5251][九省联考2018]劈配（网络流， Dinic 动态加边）
女装的你，如此好看！
10-20 351
Address 洛谷P4382 BZOJ5251 LOJ#2477 Solution 好评！感觉做完这题之后对网络流有了更深的理解。 首先要知道， Dinic 是可以动态加边的。 知道了这一点之后，就非常好做了。 建图方案：源点和汇点 SSS 和 TTT ，点分成三层，第一层为选手，第二层为志愿，第三层为导师。 （1）如果第 iii 名选手的第 jjj 档志愿填写了导师 kkk ，则由第 iii ...
九省联考2018总结
zhouyuheng2003的博客
04-15 1619
前言 我作为一名ZJOI2018day1凉凉选手，为获得更多的省选体验，去上海参加了九省联考2018，当然，是属于外省选手参加的（所以实际并测的时候出现了一些问题也得不到解决），做完两天的题，深有感触，在此写一篇博客进行总结（姗姗来迟~） （另：为啥ZJOI不写博客呢，因为day2还没考，逃~） day0 从杭州出发，乘坐高铁来到上海，高铁的速度还挺快，实际上就花了一个小时就到了上海虹...
[八省联考2018] 劈配 mentor
weixin_30794491的博客
04-09 79
Description 一年一度的综艺节目《中国新代码》又开始了。Zayid 从小就梦想成为一名程序员，他觉得这是一个展示自己的舞台，于是他毫不犹豫地报名了。 Input 轻车熟路的Zayid 顺利地通过了海选，接下来的环节是导师盲选，这一阶段的规则是这样的： 总共n 名参赛选手（编号从1 至n）每人写出一份代码并介绍自己的梦想。接着 由所有导师对这些选手进行排名。为了避免后续的麻烦，...
bzoj5251 【九省联考 2018】劈配 【网络流】
cdsszjj的博客
05-02 285
解题思路： 一眼可看出是网络流，但如何建图才能满足优先性呢？ 注意dinic算法是可以边加边边增广的，如果知道这一点就比较好想了。 首先我们建图，将所有导师向汇点连一条流量为招收学员上限的边，我们按次序枚举每一个学员的志愿，从原点向学员连流量为1的边，按当前枚举到的志愿由学员向导师连边，流量为1，然后开始增广，若失败，则删去这一志愿连的边，继续枚举下一志愿，若成功则这就是第一问答案，这样的话...
[2018多省省队联测]劈配
weixin_30569033的博客
03-27 55
题目 第一问非常裸，就是一个动态加边的\(dinic\)，我们从\(1\)到\(n\)考虑每一个人的每一个志愿，每次把这个人这一志愿的所有边都连出来，跑最大流发现流量增加就说明这个志愿匹配上了，退出做下一个人；否则就做这个人的下一个志愿 在\(dinic\)的时候删掉没有用的边会让代码快很多 第二问看起来非常二分，但是二分需要我们每次都重构出残量网络或者是暴力存下所有残量网络，这有点过于...
[BZOJ5251][多省联测2018]劈配
abcyan1235的博客
04-09 57
bzoj luogu sol 从前往后依次加边，每次对一个人做完劈配后就把当前这个残余网络存下来。这样第二问就可以二分排在第几名然后check一下在对应排名的残余网络上还能不能再增广。 给网络流开结构体直接赋值美滋滋。 然后这样就可以根据评测机的实际情况获得\(60-80\)不等的好成绩。 一个显而易见的优化就是：对于一个人，如果他已经确认了选择某个志愿，那么不仅更靠后的志愿是没用的，更...
bzoj5251: [2018多省省队联测]劈配
weixin_30642267的博客
04-08 107
bzoj5251: [2018多省省队联测]劈配 Description 一年一度的综艺节目《中国新代码》又开始了。 Zayid从小就梦想成为一名程序员，他觉得这是一个展示自己的舞台，于是他毫不犹豫地报名了。 题目描述 轻车熟路的Zayid顺利地通过了海选，接下来的环节是导师盲选，这一阶段的规则是这样的： 总共n名参赛选手（编号从1至n）每人写出一份代码并介绍自己的梦想。接着由所有导师对...
[八省联考2018]劈配（洛谷p4382）题解（网络流）
GKxx
09-26 223
传送门 发题解祭一下我调了两个晚上的结果。 首先可以明显看出（标题也暗示了）这题是一个二分图的问题，左边选手右边导师。 但是导师可以匹配多个点，如果用匈牙利的话就比较麻烦其实是我不会匈牙利，官方题解里好像有一种做法是把每个导师都拆成b[i]b[i]b[i]个点非常暴力啊，这难道就是所谓的“劈配”？（逃 所以我们当然是选择网络流了！只要把导师到汇点的容量设为b[i]b[i]b[i]就行了 再看看第二...
luogu4382 && bzoj5251[八省联考2018]劈配
maxtir的博客
01-20 239
[八省联考2018]劈配 luogu bzoj 分析 头一回见到Dinic可以动态加边的。 残量网络真是一个神奇的东西。 如果想到了动态加边，那就是一个简单的匹配问题了。 源点连导师，学员连汇点。加边的时候按志愿加，把每个志愿加进去之后跑DinicDinicDinic的BfsBfsBfs分层部分。如果有增广路就增广一下然后breakbreakbreak，否则的话把边删掉。 这样第一问就搞定了。 第...
BZOJ2839:集合计数(容斥,组合数学)
最新发布
06-08
这个问题是一个比较经典的容斥原理的应用题目。我们可以考虑使用容斥原理来求解这个问题。具体来说，我们可以考虑求出所有集合的个数，然后减去至少包含一个特定元素的集合的个数，再加上至少包含两个特定元素的集合...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值