转载: http://blog.csdn.net/axuan_k/article/details/47297395

题目描述:

现在要针对多赛区竞赛制定一个预算,该预算是一个行代表不同种类支出、列代表不同赛区支出的矩阵。组委会曾经开会讨论过各类支出的总和,以及各赛区所需支出的总和。

另外,组委会还讨论了一些特殊的约束条件:例如,有人提出计算机中心至少需要1000K 里亚尔(伊朗货币),用于购买食物;也有人提出Sharif 赛区用于购买体恤衫的费用不能超过30000K 里亚尔。

组委会的任务是制定一个满足所有约束条件且行列和满足要求的预算。

解题报告:

该题解题关键在于建图, 将行置于左边,列置于右边,然后连接所有行和列,每条边表示第x行第y列的物品的价格预算。

根据约束条件来限定这些边的上界与下界,从而将题目转化为 有上下界的可行流求解:

已有边的权为c[i][j]-b[i][j]

另外构造源点和汇点    分别连接 所有行 和 所有列

边权为 sum[x]-sum(b[x][i])        sum[y]-sum(b[i][y])

如果源流和汇流相等 则说明有可行解

相当于无源,我们只需汇点向源点连一条[0,inf]的边就行了

 #include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <queue>
const int MAXN =;
const int MAXM=;
const int INF=0x3f3f3f;
const int maxx=0x3f3f3f3f;
using namespace std;
struct Edge {
int to,cap,flow,next;
} edge[MAXM];
int head[MAXN],tot,gap[MAXN],d[MAXN],cur[MAXN],que[MAXN],p[MAXN]; void init()
{
tot=;
memset(head,-,sizeof(head));
} int isap(int source,int sink,int N)
{
memset(gap,,sizeof(gap));
memset(d,,sizeof(d));
memcpy(cur,head,sizeof(head));
int top = ,x = source,flow = ;
while(d[source] < N) {
if(x == sink) {
int Min = maxx,inser=;
for(int i = ; i < top; ++i) {
if(Min > edge[p[i]].cap - edge[p[i]].flow) {
Min = edge[p[i]].cap - edge[p[i]].flow;
inser = i;
}
}
for(int i = ; i < top; ++i) {
edge[p[i]].flow += Min;
edge[p[i]^].flow -= Min;
}
if(Min!=INF) flow += Min;
top = inser;
x = edge[p[top]^].to;
continue;
}
int ok = ;
for(int i = cur[x]; i != -; i = edge[i].next) {
int v = edge[i].to;
if(edge[i].cap > edge[i].flow && d[v]+ == d[x]) {
ok = ;
cur[x] = i;
p[top++] = i;
x = edge[i].to;
break;
}
}
if(!ok) {
int Min = N;
for(int i = head[x]; i != -; i = edge[i].next) {
if(edge[i].cap > edge[i].flow && d[edge[i].to] < Min) {
Min = d[edge[i].to];
cur[x] = i;
}
}
if(--gap[d[x]] == ) break;
gap[d[x] = Min+]++;
if(x != source) x = edge[p[--top]^].to;
}
}
return flow;
} int b[][],h[][];
int sumr[],sumc[];
int id[][]; //用来保存[i][j]在邻接表的边的编号
int s,t;
int n,m;
int sum;
int flag; void addedge(int u,int v,int c,int x,int y)
{
edge[tot]=(Edge){v,c,,head[u]};
id[x][y]= head[u] = tot++;
edge[tot]=(Edge){u,c,c,head[v]};
head[v] = tot++;
} void set_(int x,int y,char ch,int v)
{
if(ch=='=')
{
if(b[x][y]>v)
flag=;
if(h[x][y]<v)
flag=;
b[x][y]=h[x][y]=v;
}
else if(ch=='>')
{
if(h[x][y]<=v)
flag=;
b[x][y]=max(b[x][y],v+);
}
else
{
if(b[x][y]>=v)
flag=;
h[x][y]=min(h[x][y],v-);
}
} void build()
{
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
addedge(i,n+j,h[i][j]-b[i][j],i,j);
for(int i=;i<=n+m;i++)
{
int f=i<=n?sumr[i]:-sumc[i-n];
if(i<=n)
for(int j=;j<=m;j++)
f-=b[i][j];
else
for(int j=;j<=n;j++)
f+=b[j][i-n];
if(f>){
addedge(s,i,f,,);
sum+=f;
}
else
addedge(i,t,-f,,);
}
} int main()
{
// freopen("in.txt","r",stdin);
int T;
int nn,a,bb,c;
char ch;
scanf("%d",&T);
for(int Case=;Case<=T;Case++)
{
if(Case!=)
printf("\n");
scanf("%d%d",&n,&m);
init();
s=,t=n+m+;
sum=;
flag=; for(int i=;i<=n;i++)
for(int j=;j<=m;j++) //初始上下界
{
b[i][j]=;
h[i][j]=INF;
}
for(int i=;i<=n;i++)
scanf("%d",&sumr[i]);
for(int j=;j<=m;j++)
scanf("%d",&sumc[j]);
scanf("%d",&nn);
while(nn--) //限定上下界
{
cin>>a>>bb>>ch>>c;
if(!a&&!bb)
{
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
set_(i,j,ch,c);
}
else if(!a&&bb)
{
for(int i=;i<=n;i++)
set_(i,bb,ch,c);
}
else if(a&&!bb)
{
for(int i=;i<=m;i++)
set_(a,i,ch,c);
}
else
set_(a,bb,ch,c);
}
if(flag==){ //约束条件矛盾
printf("IMPOSSIBLE\n");
continue;
}
build();
int ans=isap(s,t,t+); if(sum!=ans) //源流==汇流
printf("IMPOSSIBLE\n");
else {
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
if(j==m)
printf("%d\n",b[i][j]+edge[id[i][j]].flow);
else
printf("%d ",b[i][j]+edge[id[i][j]].flow);
}
}
return ;
}

POJ 2396 Budget (上下界网络流有源可行流)的更多相关文章

  1. hdu 4940 Destroy Transportation system( 无源汇上下界网络流的可行流推断 )

    题意:有n个点和m条有向边构成的网络.每条边有两个花费: d:毁坏这条边的花费 b:重建一条双向边的花费 寻找这样两个点集,使得点集s到点集t满足 毁坏全部S到T的路径的费用和 > 毁坏全部T到 ...

  2. sgu 176 上下界网络流最小可行流带输出方案

    算法步骤: 1. 先将原图像最大可行流那样变换,唯一不同的是不加dst->src那条边来将它变成无源无汇的网络流图.直接跑一边超级源到超级汇的最大流. 2. 加上刚才没有加上的那条边p 3. 再 ...

  3. BZOJ 3876 支线剧情 有源汇有上下界最小费用可行流

    题意: 给定一张拓扑图,每条边有边权,每次只能从第一个点出发沿着拓扑图走一条路径,求遍历所有边所需要的最小边权和 分析: 这道题乍一看,可能会想到什么最小链覆盖之类的,但是仔细一想,会发现不行,一是因 ...

  4. ZOJ2314 Reactor Cooling(无源汇流量有上下界网络的可行流)

    题目大概说一个核反应堆的冷却系统有n个结点,有m条单向的管子连接它们,管子内流量有上下界的要求,问能否使液体在整个系统中循环流动. 本质上就是求一个无源汇流量有上下界的容量网络的可行流,因为无源汇的容 ...

  5. zoj 3229 上下界网络最大可行流带输出方案

    收获: 1. 上下界网络流求最大流步骤: 1) 建出无环无汇的网络,并看是否存在可行流 2) 如果存在,那么以原来的源汇跑一次最大流 3) 流量下界加上当前网络每条边的流量就是最大可行流了. 2. 输 ...

  6. ZOJ 1314 Reactor Cooling | 上下界无源汇可行流

    ZOJ 1314 Reactor Cooling | 上下界无源汇可行流 题意 有一个网络,每条边有流量的上界和下界,求一种方案,让里面的流可以循环往复地流动起来. 题解 上下界无源汇可行流的模型: ...

  7. BZOJ 2055 80人环游世界 有上下界最小费用可行流

    题意: 现在有这么一个m人的团伙,也想来一次环游世界. 他们打算兵分多路,游遍每一个国家.    因为他们主要分布在东方,所以他们只朝西方进军.设从东方到西方的每一个国家的编号依次为1...N.假若第 ...

  8. POJ2396 Budget(有源汇流量有上下界网络的可行流)

    题目大概给一个有n×m个单元的矩阵,各单元是一个非负整数,已知其每行每列所有单元的和,还有几个约束条件描述一些单元是大于小于还是等于某个数,问矩阵可以是怎样的. 经典的流量有上下界网络流问题. 把行. ...

  9. [BZOJ2502]清理雪道 有上下界网络流(最小流)

    2502: 清理雪道 Time Limit: 10 Sec  Memory Limit: 128 MB Description        滑雪场坐落在FJ省西北部的若干座山上. 从空中鸟瞰,滑雪场 ...

随机推荐

  1. Oracle创建dblink进行跨数据库操作

    创建dblink:create public database link link_to_test connect to 用户名 IDENTIFIED by "密码大写" usin ...

  2. 【ABAP系列】SAP DOI技术中I_OI_SPREADSHEET接口的使用

    公众号:SAP Technical 本文作者:matinal 原文出处:http://www.cnblogs.com/SAPmatinal/ 原文链接:[ABAP系列]SAP DOI技术中I_OI_S ...

  3. Windows下通过GitHub+Hexo搭建个人博客的步骤

    Windows下通过GitHub+Hexo搭建个人博客的步骤  https://blog.csdn.net/namechenfl/article/details/90442312 https://bl ...

  4. ‘’‘安装PyMouse,个人日志'''

    管理员启动CMD 1.直接pip install pymouse,成功安装, 2.安装PyHook: https://www.lfd.uci.edu/~gohlke/pythonlibs/ (向下滑动 ...

  5. python学习第二十八天函数局部变量的用法

    函数局部变量是在函数里面的变量,只能在函数内部使用,如果函数没有找对应变量,函数将去函数外部找对应变量,局部变量优先级大于外部变量,详细说明一下 1,局部变量已经定义值 name='zhan san' ...

  6. CSS-06 外部JS,CSS文件的寻址问题

    如果js.css外部文件有使用到相对路径,其相对路径的基准是不一样的 当一个index.html中引入外部的JS和CSS文件时: 在index.css文件中,相对路径的写法是以css文件相对于img图 ...

  7. 2019-1-24-WPF-文字描边

    title author date CreateTime categories WPF 文字描边 lindexi 2019-01-24 19:47:18 +0800 2019-1-24 19:40:7 ...

  8. 几种激活Profile的方式

    方法一: 选择spring.profiles.active spring.profiles.active=prodspring.profiles.active=dev 方法二: 选择spring.pr ...

  9. python断言语句assert

    断言语句的格式 assert test, [msg] test是一个表达式,表达式求值为Fals时引发AssertionError异常,msg是可选的异常消息. def test_assert(a): ...

  10. MySQL--15 MHA简介

    目录 一.MHA简介 二.工作流程 三.MHA架构图 四.MHA工具介绍 五.基于GTID的主从复制 六.部署MHA 一.MHA简介 松信嘉範: MySQL/Linux专家 2001年索尼公司入职 ...