【UOJ575】光伏元件【网络流建图】【上下界网络流】【费用流】
题意: n×nn\times nn×n 的 01 矩阵,对于 i∈[1,n]i \in [1,n]i∈[1,n] 有三个参数 li,li,kil_i,l_i,k_ili,li,ki,表示第 iii 行,第 iii 列的 111 的个数分别在 [li,ri][l_i,r_i][li,ri] 中,且差的绝对值不超过 kik_iki。每个元素可以花费 ci,jc_{i,j}ci,j 修改,也有不可修改的元素。求最少代价满足限制。保证有解。
n≤100n\leq 100n≤100
显然是个最小费用可行流。
行列分别建点,对每个 iii 建一个虚点。行的点出度代表这一行的 111 的个数,列的点入度代表这一列的 111 的个数,通过虚点来达成限制。
具体而言,把这三个点连成一个环流,与虚点连接的边限对应的上下界。对于差的绝对值,相当于该虚点允许 kik_iki 范围内的流量不平衡,新建超级源汇 S,TS,TS,T,SSS 往虚点,虚点往 TTT 连流量为 kik_iki 的边。
对于矩阵中的点,如果原来是 000 就正常连边,如果是 111 就强制流 111 建反向边。
因为两条 kik_iki 的边可以不满流,所以不能直接跑。从 TTT 往 SSS 连无穷的流量,然后新建超超级源汇 SS,STSS,STSS,ST 来达到可行流条件。
然后跑最小费用可行流即可,即 SSSSSS 到 STSTST 的最小费用最大流。
注意原来的 S,TS,TS,T 就只是来调整平衡的,不要在上面做些奇怪的事情。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cctype>
#include <queue>
#define MAXN 405
#define MAXM 2000005
using namespace std;
const int INF=0x7fffffff;
typedef long long ll;
inline int read()
{int ans=0,f=1;char c=getchar();while (!isdigit(c)) (c=='-')&&(f=-1),c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return f*ans;
}
struct edge{int u,v,c,w;}e[MAXM];
int head[MAXN],cur[MAXN],nxt[MAXM],cnt=1;
inline void insert(int u,int v,int c,int w)
{e[++cnt]=(edge){u,v,c,w};nxt[cnt]=head[u];head[u]=cnt;
}
inline void addnode(int u,int v,int c,int w){insert(u,v,c,w),insert(v,u,0,-w);}
int dis[MAXN],vis[MAXN],S,T;
bool spfa()
{memset(dis,0x3f,sizeof(dis));memset(vis,0,sizeof(vis));dis[S]=0,vis[S]=1;queue<int> q;q.push(S);while (!q.empty()){int u=q.front();q.pop();vis[u]=0;for (int i=head[u];i;i=nxt[i])if (e[i].c&&dis[u]+e[i].w<dis[e[i].v]){dis[e[i].v]=dis[u]+e[i].w;if (!vis[e[i].v]) q.push(e[i].v),vis[e[i].v]=1;}}return dis[T]<dis[0];
}
ll cost;
int dfs(int u,int f)
{if (u==T||!f) return f;int used=0;vis[u]=1;for (int& i=cur[u];i;i=nxt[i])if (e[i].c&&!vis[e[i].v]&&dis[u]+e[i].w==dis[e[i].v]){int w=dfs(e[i].v,min(f,e[i].c));if (!w) continue;e[i^1].c+=w,e[i].c-=w;used+=w,f-=w;cost+=(ll)w*e[i].w;if (!f) break;}return used;
}
inline int dinic()
{int mflow=0;while (spfa()){memcpy(cur,head,sizeof(cur));mflow+=dfs(S,INF);}return mflow;
}
int A[MAXN][MAXN],C[MAXN][MAXN],s[MAXN];
int main()
{int n=read();for (int i=1;i<=n;i++)for (int j=1;j<=n;j++)A[i][j]=read();for (int i=1;i<=n;i++)for (int j=1;j<=n;j++)C[i][j]=read();int SS=3*n+1,ST=3*n+2;for (int i=1;i<=n;i++){int l,r,k;l=read(),r=read(),k=read();s[i]+=l,s[i+2*n]-=l,addnode(i+2*n,i,r-l,0);s[i+2*n]+=l,s[i+n]-=l,addnode(i+n,i+2*n,r-l,0);addnode(SS,i+2*n,k,0),addnode(i+2*n,ST,k,0);}addnode(ST,SS,INF,0);int cur=cnt+1;for (int i=1;i<=n;i++)for (int j=1;j<=n;j++){if (A[i][j]) ++s[j+n],--s[i];if (~C[i][j]){if (A[i][j]) addnode(j+n,i,1,C[i][j]);else addnode(i,j+n,1,C[i][j]);}}S=3*n+3,T=3*n+4;int sum=0;for (int i=1;i<=3*n;i++){if (s[i]>0) sum+=s[i],addnode(S,i,s[i],0);else if (s[i]<0) addnode(i,T,-s[i],0);}cerr<<dinic()<<' '<<sum<<'\n';cout<<cost<<'\n';for (int i=1;i<=n;i++,puts(""))for (int j=1;j<=n;j++){if (~C[i][j]){if (!e[cur].c) A[i][j]^=1;cur+=2;}printf("%d ",A[i][j]);}return 0;
}
【UOJ575】光伏元件【网络流建图】【上下界网络流】【费用流】相关推荐
- 【bzoj3876】【AHOI2014】【支线剧情】【有上下界的费用流】
Description [故事背景] 宅男JYY非常喜欢玩RPG游戏,比如仙剑,轩辕剑等等.不过JYY喜欢的并不是战斗场景,而是类似电视剧一般的充满恩怨情仇的剧情.这些游戏往往 都有很多的支线剧情,现 ...
- HDU4411 Arrest(负边权或上下界网络流建图)
牵涉到必须经过每个点至少一次,所以拆点 把点iii拆分为iii和i+ni+ni+n iii到i+ni+ni+n可以走多次,但至少需要经过一次(代表经过这个点) 所以iii到i+ni+ni+n首先连一条 ...
- 有上下界网络流 ---- P4843 清理雪道(DAG图上最小路径重复边覆盖)【模板】有源汇上下界最小流
题目链接 题目大意: 解题思路: 首先我们发现对于每条边至少要覆盖一次,最多覆盖无数次 那么就有点像上下界网络流了[1,INF][1,INF][1,INF]的限制关系 跑一边最小流就可以了!! #in ...
- 有上下界的网络流1-无源汇带上下界网络流SGU194
有上下界的网络流1-无源汇带上下界网络流SGU194 今天开始啃网络流了.对于求解无源汇带上下界的网络流,我们可以这样建图: 建图模型: 以前写的最大流默认的下界为0,而这里的下界却 ...
- 上下界网络流-无源汇可行流与有源汇最大流
上下界网络流 2021.9.3 无源汇上下界可行流 之前的最大流讨论一般为有源无下届情况,那么无源汇有上下界可行流应如何求解? 首先要做的是消除下边界,应如何做?在有下届情形下,流网络中的任意一条边的 ...
- 有上下界网络流问题汇总
无源汇有上下界可行流 法一(据说适合点少边多的图): 建图方法 首先建立附加源点ss和附加汇点tt 对于原图中的边x->y,若限制为[b,c],那么连边x->y,流量为c-b 对于原图中的 ...
- [学习笔记]上下界网络流
有的时候,网络流建模要考虑某些边必须选择若干次,又不能多于若干次,而且不太容易转化成比较好的限制模型, 就简单粗暴地给每条边定一个流量的上下界,求在满足上下界的基础上的一些问题. 大概有以下几种. 基 ...
- 【学习笔记】上下界网络流
零.前言 昨天刚刚被神秘建图方式打败了.今天又发现连知识点都没学完-- 学习知料:"人菜无水平" 国集学长太阳神.看上去像法语的家伙.网络流建模基础(只是挂个链接)."简 ...
- 学习上下界网络流小记
前言 这个上下界网络流是一个以前我这个巨弱弱想都不敢想的一个东西. 然而,最近一次比赛居然考了这个东东. 于是整个机房掀起了破烂学上下界网络流的热. 那么我也来学学. 预备知识 要懂得很多很多的网络流 ...
- HDU 4940 Destroy Transportation system(无源汇上下界网络流)
Problem Description Tom is a commander, his task is destroying his enemy's transportation system. Le ...
最新文章
- C++ 如何获取数组的长度
- Linux那些事儿 之 戏说USB(28)设备的生命线(十一)
- 5、图书类别查询功能
- 怎样提高团队管理能力5
- superMap添加marker及连线
- Array.prototype.slice.call(arguments)
- oracle 加密怎么解密,oracle加密encrypt,解密decrypt,
- 在Linux中模拟击键和鼠标移动
- Succinctly 中文系列教程(二) 20220109 更新
- android 系统(99)---MTK 平台系统重启分类
- MySql基础汇总-DDL DML DQL DCL,存储过程,事务,触发器,函数等
- 前端基于Canvas生成等值面的方案
- Xcode安装及卸载
- 挑战程序设计竞赛-小笔记
- 百度云网盘一直显示“下载请求中”,一个 解决办法
- 应用程序无法正常启动(0x000007b)或者找不到dll文件(以vcruntime140d.dll为例)的原因原理分析和解决方法(亲测已解决)
- 一些dalao的博客
- Cubase Pro 12 v12.0.0 Full R2R WiN 音乐制作软件完整版
- VC6 SP6下载地址
- 谈谈未来 APP 产品设计的革命性趋势
热门文章
- mathwin.top 联系我们_设计:“有目的”的“因果联系”
- 打印文件前,千万记得把弹窗叉掉!!!
- 世界上最奇异的10种树,你都见过吗?
- 那些你从未见过的神奇物理化学实验,今天全给你整理了!
- 孩子不是笨,他和“最强大脑”差的是这个!
- 三次握手和四次挥手图解_三次握手和四次挥手简单理解
- 如何查询服务器是否安装系统时间,如何查看系统当前的NTP配置?
- java 进程同步代码_java 实现进程间的同步(源代码)
- aqs java 简书,Java AQS源码解读
- github可以传java吗_如何在github上传本地项目代码(新手使用)----亲测使用