BZOJ.5251.[八省联考2018]劈配mentor(最大流)

于 2019-02-11 20:14:00 发布
阅读量68 收藏
点赞数
原文链接：http://www.cnblogs.com/SovietPower/p/10363057.html
版权

BZOJ
洛谷

对于每个人，每次枚举一个志愿看是否能增广即可。
对于第二问，可以保留第一问中\(n\)次增广前后的\(n\)张图，二分，在对应图上看是否能增广即可。

貌似匈牙利的某种写法比网络流优多了...懒得看...
写的常数莫名大...

增广失败的话原图流量不会变啊，不用特意删掉新加的边，把表头head改了即可，而且判断是否能增广只需要判断BFS()是否返回1。

Dinic：

//195228kb  6676ms
#include <cstdio>
#include <cctype>
#include <vector>
#include <cstring>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=408,M=(N+40005)*2,INF=0x3f3f3f3f;

int T;
struct Graph
{
    int Enum,H[N],nxt[M],to[M],cap[M],cur[N],lev[N];
    inline void Init(int n)
    {
        Enum=1, memset(H,0,n+1<<2);
    }
    inline void AE(int u,int v,int w)
    {
        to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w;
        to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0;
    }
    bool BFS()
    {
        static int q[N];
        memset(lev,0,T+1<<2);
        int h=0,t=1; q[0]=0, lev[0]=1;
        while(h<t)
        {
            int x=q[h++];
            for(int i=H[x]; i; i=nxt[i])
                if(!lev[to[i]] && cap[i])
                {
                    q[t++]=to[i], lev[to[i]]=lev[x]+1;
                    if(to[i]==T) return 1;
                }
        }
        return 0;
    }
    int DFS(int x,int flow)
    {
        if(x==T) return flow;
        int used=0;
        for(int &i=cur[x]; i; i=nxt[i])
            if(lev[to[i]]==lev[x]+1 && cap[i])
            {
                int delta=DFS(to[i],std::min(flow,cap[i]));
                if(delta)
                {
                    cap[i]-=delta, cap[i^1]+=delta, used+=delta;
                    if(used==flow) return flow;
                }
            }
        return lev[x]=0,used;
    }
    void Dinic()
    {
        while(BFS()) memcpy(cur,H,T+1<<2), DFS(0,INF);
    }
}G[203];

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-48,c=gc());
    return now;
}

int main()
{
    static int B[203],s[203],Ans1[203],tmpH[N];
    static std::vector<int> v[203][203];

    freopen("mentor.in","r",stdin);
    freopen("mentor.out","w",stdout);

    for(int Tests=read(),C=read(); Tests--; )
    {
        int n=read(),m=read(); T=n+m+1;
        G[0].Init(T);
        for(int i=1; i<=n; ++i) G[0].AE(0,i,1);
        for(int i=1; i<=m; ++i) G[0].AE(i+n,T,B[i]=read());
        for(int i=1; i<=n; ++i)
            for(int j=1; j<=m; ++j) v[i][j].clear();
        for(int i=1; i<=n; ++i)
            for(int j=1; j<=m; ++j) v[i][read()].push_back(j+n);
        for(int i=1; i<=n; ++i) s[i]=read();

        for(int i=1; i<=n; ++i)
        {
            G[i]=G[i-1];
            int j;
            for(j=1; j<=m; ++j)
            {
                std::vector<int> &vec=v[i][j];
                for(int k=0,l=vec.size(); k<l; ++k)  G[i].AE(i,vec[k],1);
                if(G[i].BFS()) {G[i].Dinic(); break;}
                if(vec.size()) G[i].Enum=G[i-1].Enum, memcpy(G[i].H,G[i-1].H,T+1<<2);
            }
            printf("%d ",Ans1[i]=j);
        }
        putchar('\n');
    
        for(int i=1; i<=n; ++i)
            if(Ans1[i]<=s[i]) putchar('0'), putchar(' ');
            else
            {
                int l=1,r=i-1,mid,ans=0,dream=s[i];
                while(l<=r)
                {
                    mid=l+r>>1;
                    int tmpEnum=G[mid-1].Enum;
                    memcpy(tmpH,G[mid-1].H,T+1<<2);
                    bool ok=0;
                    for(int j=1; j<=dream; ++j)
                    {
                        std::vector<int> &vec=v[i][j];
                        for(int k=0,l=vec.size(); k<l; ++k)  G[mid-1].AE(i,vec[k],1);
                        if(G[mid-1].BFS()) {ok=1; break;}
                    }
                    if(ok) ans=mid, l=mid+1;
                    else r=mid-1;
                    G[mid-1].Enum=tmpEnum, memcpy(G[mid-1].H,tmpH,T+1<<2);
                }
                printf("%d ",i-ans);
            }
        putchar('\n');
    }
    return 0;
}

ISAP：

//259320kb  7584ms
#include <cstdio>
#include <cctype>
#include <vector>
#include <cstring>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=408,M=(N+40005)*2,INF=0x3f3f3f3f;

int T;
struct Graph
{
    int Enum,H[N],nxt[M],fr[M],to[M],cap[M],pre[N],lev[N];
    inline void Init(int n)
    {
        Enum=1, memset(H,0,n+1<<2);
    }
    inline void AE(int u,int v,int w)
    {
        to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum, fr[Enum]=u, cap[Enum]=w;
        to[++Enum]=u, nxt[Enum]=H[v], H[v]=Enum, fr[Enum]=v, cap[Enum]=0;
    }
    bool BFS()
    {
        static int q[N];
        const int lim=T+1;
        for(int i=0; i<T; ++i) lev[i]=lim;
        int h=0,t=1; q[0]=T, lev[T]=0;
        while(h<t)
        {
            int x=q[h++];
            for(int i=H[x]; i; i=nxt[i])
                if(cap[i^1] && lev[to[i]]==lim) q[t++]=to[i], lev[to[i]]=lev[x]+1;
        }
        return lev[0]<=T;
    }
    inline void Augment()
    {
//      int mn=INF;
//      for(int i=T; i; i=fr[pre[i]])
//          mn=std::min(mn,cap[pre[i]]);
        for(int i=T; i; i=fr[pre[i]])
            --cap[pre[i]], ++cap[pre[i]^1];
    }
    void ISAP()
    {
        static int cur[N],num[N];
        memset(num,0,T+1<<2);
        for(int i=0; i<=T; ++i) cur[i]=H[i], ++num[lev[i]];
        int x=0;
        while(lev[0]<=T)
        {
            if(x==T) x=0, Augment();
            bool can=0;
            for(int i=cur[x]; i; i=nxt[i])
                if(lev[to[i]]==lev[x]-1 && cap[i])
                {
                    can=1, cur[x]=i, pre[x=to[i]]=i;
                    break;
                }
            if(!can)
            {
                int mn=T;
                for(int i=H[x]; i; i=nxt[i])
                    if(cap[i]) mn=std::min(mn,lev[to[i]]);
                if(!--num[lev[x]]) break;
                ++num[lev[x]=mn+1], cur[x]=H[x];
                if(x) x=fr[pre[x]];
            }
        }
    }
}G[203];

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-48,c=gc());
    return now;
}

int main()
{
    static int B[203],s[203],Ans1[203],tmpH[N];
    static std::vector<int> v[203][203];

    freopen("mentor.in","r",stdin);
    freopen("mentor.out","w",stdout);

    for(int Tests=read(),C=read(); Tests--; )
    {
        int n=read(),m=read(); T=n+m+1;
        G[0].Init(T);
        for(int i=1; i<=n; ++i) G[0].AE(0,i,1);
        for(int i=1; i<=m; ++i) G[0].AE(i+n,T,B[i]=read());
        for(int i=1; i<=n; ++i)
            for(int j=1; j<=m; ++j) v[i][j].clear();
        for(int i=1; i<=n; ++i)
            for(int j=1; j<=m; ++j) v[i][read()].push_back(j+n);
        for(int i=1; i<=n; ++i) s[i]=read();

        for(int i=1; i<=n; ++i)
        {
            G[i]=G[i-1];
            int j;
            for(j=1; j<=m; ++j)
            {
                std::vector<int> &vec=v[i][j];
                for(int k=0,l=vec.size(); k<l; ++k)  G[i].AE(i,vec[k],1);
                if(G[i].BFS()) {G[i].ISAP(); break;}
                if(vec.size()) G[i].Enum=G[i-1].Enum, memcpy(G[i].H,G[i-1].H,T+1<<2);
            }
            printf("%d ",Ans1[i]=j);
        }
        putchar('\n');
    
        for(int i=1; i<=n; ++i)
            if(Ans1[i]<=s[i]) putchar('0'), putchar(' ');
            else
            {
                int l=1,r=i-1,mid,ans=0,dream=s[i];
                while(l<=r)
                {
                    mid=l+r>>1;
                    int tmpEnum=G[mid-1].Enum;
                    memcpy(tmpH,G[mid-1].H,T+1<<2);
                    bool ok=0;
                    for(int j=1; j<=dream; ++j)
                    {
                        std::vector<int> &vec=v[i][j];
                        for(int k=0,l=vec.size(); k<l; ++k)  G[mid-1].AE(i,vec[k],1);
                        if(G[mid-1].BFS()) {ok=1; break;}
                    }
                    if(ok) ans=mid, l=mid+1;
                    else r=mid-1;
                    G[mid-1].Enum=tmpEnum, memcpy(G[mid-1].H,tmpH,T+1<<2);
                }
                printf("%d ",i-ans);
            }
        putchar('\n');
    }
    return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/10363057.html

weixin_30539625
关注
Node.js-ZOJCH是BZOJ题库的离线版
08-09
【标题】"Node.js-ZOJCH是BZOJ题库的离线版"揭示了这个项目的核心功能，即提供一个基于Node.js平台的离线版本的BZOJ（BestCoder Judging System）题库。BZOJ是一个在线编程竞赛平台，而ZOJCH则允许用户在没有网络...
bzoj 3730. 震波（动态点分治 + vector树状数组）
01-06
（自己写的动态点分治巨垃圾，常数是别人的两倍） 用动态开点线段树死活过不去，学了一波大佬用 vector 开树状数组立马就卡过去了 考虑点分树的做法，在点分树上每个点以距离为下标建一棵线段树，每次询问查询子树的...
BZOJ第一部分
11-03
【BZOJ第一部分】 BZOJ，全称为“北京邮电大学在线评测系统”（Beijing Zhiyuan Online Judge），是中国最早的一批在线编程竞赛平台之一，它为编程爱好者提供了一个练习和检验编程技能的环境。在这个系统中，用户...
wangjunrui666#bzoj-problem##5335. [TJOI2018]智力竞赛1
07-25
比赛规则如下:一共有m道题目，每个入都有1次答题机会，每次答题为选择一道题目回答，在回答正确后，可以从这个题目的后续题目，直达题目答错题目或者没有后续题目。其中
CTSC 2011 无穷图的桥(BZOJ 2307) 题解.ppt
03-17
CTSC 2011 无穷图的桥(BZOJ 2307) 题解.ppt
命令手册 Linux常用命令
11-07
命令手册 Linux常用命令
【超强组合】基于VMD-雪融优化算法SAO-Transformer-GRU的光伏预测算研究Matlab实现.rar
11-07
1.版本：matlab2014/2019a/2024a 2.附赠案例数据可直接运行matlab程序。 3.代码特点：参数化编程、参数可方便更改、代码编程思路清晰、注释明细。 4.适用对象：计算机，电子信息工程、数学等专业的大学生课程设计、期末大作业和毕业设计。 替换数据可以直接使用，注释清楚，适合新手
【超强组合】基于VMD-花朵授粉优化算法FPA-Transformer-BiLSTM的光伏预测算研究Matlab实现.rar
11-07
1.版本：matlab2014/2019a/2024a 2.附赠案例数据可直接运行matlab程序。 3.代码特点：参数化编程、参数可方便更改、代码编程思路清晰、注释明细。 4.适用对象：计算机，电子信息工程、数学等专业的大学生课程设计、期末大作业和毕业设计。 替换数据可以直接使用，注释清楚，适合新手
基于SpringBoot+Shiro+mysql实现的个人博客前后台管理系统 【完整源码+数据库】
11-07
一、功能描述 文章管理 分类管理 评论管理 数据库监控 通用页面 后台首页 友链管理 在线用户 权限管理 角色管理 站点管理 标签管理 主题管理 上传管理 用户管理 踢出页面 登录页面 注册页面 主题目录 默认主题 二、技术栈 Spring Boot、Apache Shiro、MyBatis-Plus、Alibaba Druid、Redis、MySQL、Thymeleaf 三、安装 将本项目源码导入本地开发工具(如 IntelliJ IDEA )，本地开发工具需要安装 lombok 插件 安装Mysql数据库：Mysql版本最低支持5.7，新建 database CREATE DATABASE pb_cms_base; 初始化数据库：找到项目数据库文件:docs/db/pb_cms_base.sql，执行 pb_cms_base.sql 前台首页，浏览器访问http://localhost:8080 后台首页，浏览器访问http://localhost:8080/admin使用账号密码admin,123456登录系统后台。
暴风电视刷机 T55FUA 通用ECHO 屏ST5461D07-2 机编60000AM6400 AM6700 V1.0.03版本
11-07
务必确认机身编号与文件名机编一致，如不一致，请勿下载 机身编号一般在机子背面的贴纸上 适配屏参：30164505 MSD6A838平台升级步骤 强制升级（不开机强制升级）： 1、下载数据，压缩包解压，将BFTVUpdate838_xx.bin拷贝至U盘根目录下，插入电视USB端口 2、插拔下电源，按一下遥控器待机键后快速不停的点按遥控器上键触发主板识别U盘软件进行升级 3、升级成功需运行至100%，此时耐心等待电视自动操作，切勿断电或拔掉U盘 4、系统升级后第一次重启系统需要3分钟左右，属于正常现象，切勿断电 升级完成后可以在系统设置——本机信息——查询软件版本更新状态 注意： 1、U盘要求使用FAT32格式，建议4G-8G的品牌U盘，刷机成功率会高 2、升级到结束，大约需要8-30分钟，中途绝对不能断电 3、升级重启第一次进入系统，请等完全正常进入开机桌面之后，才能拨下U盘 4、如无法升级，将BFTVUpdate838_xx.bin文件重命名为BFTVUpdate838.bin，再尝试刷机
wryh+pico12num.ttf
最新发布
11-07
wryh+pico12num.ttf
性能优化与加载时间控制.docx
11-07
性能优化与加载时间控制.docx
【数据分析】基于matlab数据流分类的双模型半监督自组织模糊推理系统【含Matlab源码 9067期】.zip
11-07
CSDN海神之光上传的代码均可运行，亲测可用，直接替换数据即可，适合小白； 1、代码压缩包内容 主函数：main.m； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2019b；若运行有误，根据提示修改；若不会，私信博主； 3、运行操作步骤 步骤一：将所有文件放到Matlab的当前文件夹中； 步骤二：双击打开main.m文件； 步骤三：点击运行，等程序运行完得到结果； 4、仿真咨询 如需其他服务，可私信博主或扫描博客文章底部QQ名片； 4.1 博客或资源的完整代码提供 4.2 期刊或参考文献复现 4.3 Matlab程序定制 4.4 科研合作
【超强组合】基于VMD-黑猩猩优化算法Chimp-Transformer-LSTM的光伏预测算研究Matlab实现.rar
11-07
1.版本：matlab2014/2019a/2024a 2.附赠案例数据可直接运行matlab程序。 3.代码特点：参数化编程、参数可方便更改、代码编程思路清晰、注释明细。 4.适用对象：计算机，电子信息工程、数学等专业的大学生课程设计、期末大作业和毕业设计。 替换数据可以直接使用，注释清楚，适合新手
ubuntu系统下，很好用的多线程下载工程 类似于windows系统下的迅雷
11-07
使用时，稍研究即可。
基于java+swing+applet实现的家庭理财系统(含源码+数据库+答辩PPT)
11-07
基于java Swing+applet实现的家庭理财系统(含源码+数据库+答辩PPT)
HCIE-IPv6技术
11-07
HCIE-IPv6技术
【超强组合】基于VMD-粒子群优化算法PSO-Transformer-GRU的光伏预测算研究Matlab实现.rar
11-07
1.版本：matlab2014/2019a/2024a 2.附赠案例数据可直接运行matlab程序。 3.代码特点：参数化编程、参数可方便更改、代码编程思路清晰、注释明细。 4.适用对象：计算机，电子信息工程、数学等专业的大学生课程设计、期末大作业和毕业设计。 替换数据可以直接使用，注释清楚，适合新手
数据可视化驾驶舱，包含地图，页面可以直接运行
11-07
数据可视化驾驶舱，包含地图，页面可以直接运行
bzoj.org俄罗斯方块
05-06
bzoj.org是一个在线的编程题库，其中包含了很多经典算法问题和实际应用中的编程问题。其中包括了很多有趣的游戏算法问题，比如Puzzle，Tetris等。 而bzoj.org俄罗斯方块就是在该平台上的一个经典的俄罗斯方块游戏。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值