一、题目

有一个 \(n\times n\) 的矩阵,每个元素可能是 .C/ 的其中一种,分别表示可以放置芯片、已经放置了芯片、不能放置芯片,你可以分别决定是否可以放置芯片的位置放置芯片。

最后需要满足 \(\forall i\),第 \(i\) 行的芯片个数等于第 \(i\) 列的芯片个数,每一行的芯片个数都不超过总芯片个数的 \(\frac{A}{B}\),问在此情况下能额为放置的芯片个数最大值,如果怎么样都不合法输出 impossible

\(n\leq 40\)

二、解法

如果直接做的话并不好入手,我们考虑调整法,也就是让所有 . 的位置都放芯片再调整。

可以把芯片看成流量,我们可以套路地建出一个二分图。想象第 \(i\) 行的点上具有 \(a_i\) 点流量,第 \(i\) 列的点上具有 \(b_i\) 点流量(分别表示初始状态下它们的芯片个数),那么个数相等的条件可以转化成第 \(i\) 行和第 \(i\) 列同时减少 \(1\) 的流量。

但是最后点上就不能有残余流量,对于 . 的位置我们可以选择不放,设它的位置是 \((i,j)\) 那么它可以让 \(i\) 行 \(j\) 列同时减少 \(1\) 的流量。在此基础上我们还要最大化总芯片数,所以可以考虑增加费用这个意义,我们把同行同列的边费用设置为 \(0\),减少芯片的边费用设置为 \(1\),跑最小费用最大流即可。

还剩下最后一个限制:每一行的芯片个数都不超过总芯片个数A/B,我们可以枚举每一行的芯片个数 \(k\),把同行同列的流量设置为 \(k\),最后可以得到总芯片数 \(sum\),所以我们只需要判断下面两点:满流;\(sum\cdot A\geq k\cdot B\)

三、总结

调整法考虑寻找一个特殊的初始状态,再考虑如何描述增加\(/\)减少的过程。网络流加调整法的应用是常见的,类似的模型还有最小链覆盖,它就是把初始设置有 \(n\) 个链来调整。

网络流是属于图论的,所以它的思维方法中也涉及到思考原问题中的元素怎么对应到图上的每个量上去,比如本题因为要最大化总芯片数我们才考虑增加费用这一维。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
using namespace std;
const int M = 105;
const int inf = 0x3f3f3f3f;
int read()
{
int x=0,f=1;char c;
while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
return x*f;
}
int n,A,B,S,T,ans,tot,f[M],a[M],b[M];
int m,dis[M],pre[M],flow[M],lst[M];
struct edge{int v,f,c,next;}e[M*M];char s[M][M];
void add(int u,int v,int F,int c)
{
e[++tot]=edge{v,F,c,f[u]},f[u]=tot;
e[++tot]=edge{u,0,-c,f[v]},f[v]=tot;
}
int bfs()
{
queue<int> q;
for(int i=0;i<=T;i++)
dis[i]=inf,flow[i]=pre[i]=lst[i]=0;
dis[S]=0;flow[S]=inf;q.push(S);
while(!q.empty())
{
int u=q.front();q.pop();
for(int i=f[u];i;i=e[i].next)
{
int v=e[i].v,c=e[i].c;
if(dis[v]>dis[u]+c && e[i].f>0)
{
dis[v]=dis[u]+c;
flow[v]=min(flow[u],e[i].f);
pre[v]=u;lst[v]=i;
q.push(v);
}
}
}
return flow[T]>0;
}
void zxy(int k)
{
int res=0,sum=0,all=0;
S=0;T=2*n+1;tot=1;
for(int i=0;i<=T;i++) f[i]=0;
for(int i=1;i<=n;i++)
{
all+=a[i];
add(S,i,a[i],0);
add(i+n,T,b[i],0);
add(i,i+n,k,0);
}
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++) if(s[i][j]=='.')
add(i,j+n,1,1);
while(bfs())
{
int t=T;
sum+=flow[t];
res+=flow[t]*dis[t];
while(t)
{
e[lst[t]].f-=flow[T];
e[lst[t]^1].f+=flow[T];
t=pre[t];
}
}
if(sum==all && k*B<=(sum-res)*A)
ans=max(ans,sum-res);
}
void work()
{
m=0;ans=-1;
for(int i=1;i<=n;i++) a[i]=b[i]=0;
for(int i=1;i<=n;i++)
{
scanf("%s",s[i]+1);
for(int j=1;j<=n;j++)
{
m+=(s[i][j]=='C');
a[i]+=(s[i][j]=='C' || s[i][j]=='.');
b[j]+=(s[i][j]=='C' || s[i][j]=='.');
}
}
for(int i=0;i<=n;i++) zxy(i);
if(ans==-1) puts("impossible");
else printf("%d\n",ans-m);
}
signed main()
{
int Case=0;
while(~scanf("%d %d %d",&n,&A,&B) && n+A+B)
{
printf("Case %d: ",++Case);
work();
}
}

UVA1104 Chips Challenge的更多相关文章

  1. 【题解】uva1104 chips challenge

    原题传送门 题目分析 给定一张n*n的芯片. '.'表示该格子可以放一个零件. 'C'表示该格子已经放了一个零件(不能拆下). '/'表示该格子不能放零件. 要求在芯片的现有基础上,放置尽可能多的零件 ...

  2. [2011WorldFinal]Chips Challenge[流量平衡]

    Chips Challenge Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) ...

  3. UVA1104 芯片难题 Chips Challenge

    题目链接 题意 网格上放点,有些强制放,有些不能放,有些可以放可以不放.要求: 第 \(i\) 行的点数 = 第 \(i\) 列的点数 每一行每一列的点数不超过总点数的 \(k\) 倍(\(k\) 已 ...

  4. 【UVALive - 5131】Chips Challenge(上下界循环费用流)

    Description A prominent microprocessor company has enlisted your help to lay out some interchangeabl ...

  5. 【BZOJ 2673】[Wf2011]Chips Challenge

    题目大意: 传送门 $n*n$的棋盘,有一些位置可以放棋子,有一些已经放了棋子,有一些什么都没有,也不能放,要求放置以后满足:第i行和第i列的棋子数相同,同时每行的棋子数占总数比例小于$\frac{A ...

  6. Bzoj2673 3961: [WF2011]Chips Challenge 费用流

    国际惯例题面:如果我们枚举放几个零件的话,第二个限制很容易解决,但是第一个怎么办?(好的,这么建图不可做)考虑我们枚举每行每列最多放几个零件t,然后计算零件总数sum.这样如果可行的话,则有t*B&l ...

  7. 解题:BZOJ 2673 World Final 2011 Chips Challenge

    题面 数据范围看起来很像网络流诶(滚那 因为限制多而且强,数据范围也不大,我们考虑不直接求答案,而是转化为判定问题 可以发现第二个限制相对好满足,我们直接枚举这个限制就可以.具体来说是枚举所有行中的最 ...

  8. BZOJ2673 [Wf2011]Chips Challenge 费用流 zkw费用流 网络流

    https://darkbzoj.cf/problem/2673 有一个芯片,芯片上有N*N(1≤N≤40)个插槽,可以在里面装零件. 有些插槽不能装零件,有些插槽必须装零件,剩下的插槽随意. 要求装 ...

  9. bzoj3961[WF2011]Chips Challenge

    题意 给出一个n*n的网格,有些格子必须染成黑色,有些格子必须染成白色,其他格子可以染成黑色或者白色.要求最后第i行的黑格子数目等于第i列的黑格子数目,且某一行/列的格子数目不能超过格子总数的A/B. ...

随机推荐

  1. html5 sessionStorage util

      /** * Created by 13352 on 2018/7/5. */ define(function() { var session = { _id: null, _cookieCache ...

  2. 基于AM335X,如何搭建优良的Linux开发环境(下)

    接着上一篇文章的Linux开发环境搭建,文章中详细讲解了 VMware14.1.1虚拟机安装.基于虚拟机安装Ubuntu14.04.3操作系统.安装Ubuntu14.04.3操作系统.安装虚拟机工具. ...

  3. k8s负载资源StatefulSet工作解析

    在k8s中工作负载资源StatefulSet用于管理有状态应用. 什么是无状态? 组成一个应用的pod是对等的,它们之前没有关联和依赖关系,不依赖外部存储. 即我们上篇小作文中deployment创建 ...

  4. Springboot --- Bug集

    一. 启动springboot报错:找不到或无法加载主类 解决:直接选中项目,在ecplise选中"project",点击clean 清理项目再运行 问题解决. 二.报错:Fail ...

  5. Skywalking-12:Skywalking SPI机制

    SPI机制 基本概述 SPI 全称 Service Provider Interface ,是一种服务发现机制.通过提供接口.预定义的加载器( Loader )以及约定俗称的配置(一般在 META-I ...

  6. gin 集成 consul

    "github.com/hashicorp/consul/api" package initialize import ( "fmt" "github ...

  7. 关于web项目中的资源跳转

    1.跳转包括两种方式: 转发 forward 重定向 redirect 2.两种方式的代码: AServlet类: //向request范围中存储数据 request.setAttribute(&qu ...

  8. mysql join语句的执行流程是怎么样的

    mysql join语句的执行流程是怎么样的 join语句是使用十分频繁的sql语句,同样结果的join语句,写法不同会有非常大的性能差距. select * from t1 straight_joi ...

  9. 【UE4 C++】 UnrealPak 与 Pak 的制作、挂载、加载

    简介 通过 UnrealPak,可以将资源打包成 Pak 文件 Pak文件是UE4游戏生成的数据包文件. Pak 之前一般先有 Cooked 步骤,将资源烘焙为对应平台支持的资源 一般打包后的项目使用 ...

  10. 【UE4 设计模式】简单工厂模式 Simple Factory Pattern

    概述 描述 又称为静态工厂方法 一般使用静态方法,根据参数的不同创建不同类的实例 套路 创建抽象产品类 : 创建具体产品类,继承抽象产品类: 创建工厂类,通过静态方法根据传入不同参数从而创建不同具体产 ...