【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 100n100

显然是个最小费用可行流。

行列分别建点,对每个 iii 建一个虚点。行的点出度代表这一行的 111 的个数,列的点入度代表这一列的 111 的个数,通过虚点来达成限制。

具体而言,把这三个点连成一个环流,与虚点连接的边限对应的上下界。对于差的绝对值,相当于该虚点允许 kik_iki 范围内的流量不平衡,新建超级源汇 S,TS,TS,T,SSS 往虚点,虚点往 TTT 连流量为 kik_iki 的边。

对于矩阵中的点,如果原来是 000 就正常连边,如果是 111 就强制流 111 建反向边。

因为两条 kik_iki 的边可以不满流,所以不能直接跑。从 TTTSSS 连无穷的流量,然后新建超超级源汇 SS,STSS,STSS,ST 来达到可行流条件。

然后跑最小费用可行流即可,即 SSSSSSSTSTST 的最小费用最大流。

注意原来的 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;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314808.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

动手造轮子:实现一个简单的 EventBus

动手造轮子&#xff1a;实现一个简单的 EventBusIntroEventBus 是一种事件发布订阅模式&#xff0c;通过 EventBus 我们可以很方便的实现解耦&#xff0c;将事件的发起和事件的处理的很好的分隔开来&#xff0c;很好的实现解耦。微软官方的示例项目 EShopOnContainers 也有在使…

Codeforces Round #733 (Div. 1 + Div. 2) E. Minimax 分情况讨论 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个串&#xff0c;你可以随意安排这个串&#xff0c;使得这个串的每个前缀的kmpkmpkmp数组最大值最小&#xff0c;定义为f(a)f(a)f(a)&#xff0c;并且字典序最小&#xff0c;输出安排之后的串。 n≤1e…

【Ynoi2011】成都七中【树论】【点分树】【离线】【树状数组】

题意&#xff1a;给一棵树&#xff0c;点有颜色&#xff0c;qqq 次询问&#xff0c;每次给定 l,r,xl,r,xl,r,x &#xff0c;求只保留编号在 [l,r][l,r][l,r] 中的点时点 xxx 所在连通块的颜色数。 所有数 ≤105\leq 10^5≤105 题目背景好评 首先所有颜色不同的话就是数连通块…

通过 nginx-proxy 实现自动反向代理和 HTTPS

本章节代码已经上传至 https://github.com/siegrainwong/.NET-Core-with-Docker/tree/master/Part3系列大纲这次我们讲第三篇&#xff1a;用 docker-compose 启动 WebApi 和 SQL Server在容器中集成 Skywalking APM通过 nginx-proxy 对 Portainer、Skywalking、WebApi 实现自动…

P4781 【模板】拉格朗日插值

传送门 把公式实现一下即可&#xff1a; 当xxx连续的时候可以优化为O(N)O(N)O(N)。 // Problem: P4781 【模板】拉格朗日插值 // Contest: Luogu // URL: https://www.luogu.com.cn/problem/P4781 // Memory Limit: 125 MB // Time Limit: 1000 ms // // Powered by CP Edi…

【HNOI/AHOI2018】毒瘤【容斥】【虚树/动态dp】

题意&#xff1a;nnn 个点 mmm 条边的连通无向图的独立集个数模 998244353998244353998244353。 n≤105,m≤n10n\leq 10^5,m\leq n10n≤105,m≤n10 为什么标题要把两个算法写一起&#xff1f;因为这两个东西在这类问题上是本质相同的&#xff0c;这也是写这篇博客的原因。 显…

MediatR-进程内的消息通信框架

MediatR是一款进程内的消息订阅、发布框架&#xff0c;提供了Send方法用于发布到单个处理程序、Publish方法发布到多个处理程序&#xff0c;使用起来非常方便。目前支持 .NET Framework4.5、.NET Stardand1.3、.NET Stardand2.0等版本&#xff0c;可跨平台使用。要在项目中使用…

Codeforces Round #586 (Div. 1 + Div. 2) D. Alex and Julian 数学 + 思维

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个无限个点的坐标轴&#xff0c;一个集合BBB&#xff0c;如果存在∣i−j∣bk|i-j|b_k∣i−j∣bk​的话&#xff0c;那么i,ji,ji,j之间就连边。现在问你至少要从集合BBB中去掉多少个数才能使得连完边之…

【十二省联考2019】字符串问题【后缀自动机】【拓扑排序】

题意&#xff1a;给一个字符串 SSS&#xff0c;以子串的形式给出一些 A 类串和 B 类串以及 mmm 对 A 类串支配 B 类串的关系。求一个总长度最长的 A 类串序列&#xff0c;使得每个串都存在一个 B 类串前缀被后一个串支配。无穷输出 −1-1−1。 ∣S∣,m≤2105|S|,m\leq 2\times …

Codeforces Round #586 (Div. 1 + Div. 2) B. Multiplication Table 思维 + 公式

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个n∗nn*nn∗n的矩阵&#xff0c;每个位置由ai∗aja_i*a_jai​∗aj​得来&#xff0c;主对角线为000&#xff0c;让你求出来aia_iai​。 n≤1e3n\le1e3n≤1e3 思路&#xff1a; 由公式ai,j∗ai,kaj,…

不好意思,这么久没有更新《从零开始掌握ASP.NET Core 》

点击上方蓝字&#xff0c;关注「我们」等了快个月了&#xff0c;终于开始更新了。因为感冒&#xff0c;弄的嗓子有点沙哑。所以停了半个月才是更新&#xff0c;目前一口气更新了12个章节&#xff0c;大家可以耐心观看内容了。《从零开始学ASP.NET Core 》-- 更新通知视频课程更…

【NOI2018】你的名字【后缀自动机】【可持久化线段树合并】【乱搞】

题意&#xff1a;给一个串 SSS&#xff0c;qqq 次询问&#xff0c;每次给定串 TTT 和 l,rl,rl,r &#xff0c;求有多少个本质不同的串是 TTT 的子串而不是 Sl…rS_{l\dots r}Sl…r​ 的子串。 ∣S∣≤5105,q≤105,∑∣T∣≤106|S|\leq 5\times 10^5,q\leq 10^5,\sum|T|\leq 10^…

2021牛客暑期多校训练营1 A.Alice and Bob 博弈 SG函数

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 有两堆石子&#xff0c;两个人每次可以进行如下操作&#xff1a;从某一堆狮子中拿出x(x>0)x(x>0)x(x>0)个&#xff0c;从另一堆石子中拿出s∗x(s>0)s*x(s>0)s∗x(s>0)个。谁不能操作谁输&…

【CTSC2010】珠宝商【后缀自动机】【点分治】【根号分治】

题意&#xff1a;给一棵 nnn 个点的树&#xff0c;每个点有个字符&#xff0c;另给一个长度为 mmm 的特征串&#xff0c;求树上 n2n^2n2 条有向路径在特征串中出现的次数之和。 n,m≤5104n,m\leq 5\times 10^4n,m≤5104 看到母串先建 SAM &#xff08;bushi 树上路径统计问题…

使用Azure云原生构建博客是怎样一种体验?(下篇)

点击上方蓝字关注“汪宇杰博客”接上篇《使用Azure云原生构建博客是怎样一种体验&#xff1f;&#xff08;上篇&#xff09;》DNSAzure DNS 是一套分布全球的域名解析服务。具有超高可用性和接近实时的记录更新及生效速度。我的博客也使用了这项服务。Azure 现在可以提供域名注…

2021牛客暑期多校训练营1 G Game of Swapping Numbers 思维 + 巧妙的转换

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你两个数组A,BA,BA,B&#xff0c;你可以选择AAA的两个位置i,j,i<ji,j,i<ji,j,i<j交换Ai,AjA_i,A_jAi​,Aj​&#xff0c;需要交换正好kkk次&#xff0c;问你最大的∑i1n∣Ai−Bi∣\sum_{i1}^n|A_…

.NET Core 3.0之深入源码理解HttpClientFactory(一)

写在前面创建HttpClient实例的时候&#xff0c;在内部会创建HttpMessageHandler链&#xff0c;我们知道HttpMessageHandler是负责建立连接的抽象处理程序&#xff0c;所以HttpClient的维护实际上就是维护HttpMessageHandler的使用&#xff0c;释放HttpClient并不会及时释放连接…

WTM 构建DotNetCore开源生态,坐而论道不如起而行之

作为一个8岁开始学习编程&#xff0c;至今40岁的老程序员&#xff0c;这辈子使用过无数种语言&#xff0c;从basic开始&#xff0c;到pascal, C, C&#xff0c;到后来的 java, c#,perl,php,再到现在流行的python。小时候的我总觉得多掌握一门语言&#xff0c;我的技术能力就又前…

【WC2014】时空穿梭【组合数】【莫比乌斯反演】【整除分块】【暴力多项式】

题意&#xff1a;TTT 组数据&#xff0c;给一个 nnn 维空间&#xff0c;第 iii 维大小为 [1,mi]∩Z[1,m_i]\cap \Z[1,mi​]∩Z&#xff0c;求大小为 ccc 的严格偏序上升的共线点集个数。答案模 100071000710007。 T≤100,n≤11,m≤105,c≤20T\leq 100,n\leq 11,m\leq 10^5,c\le…

2021牛客暑期多校训练营1 H Hash Function FFT\NTT

传送门 文章目录题意&#xff1a;思路&#xff1a;题意&#xff1a; 给你一个数组aaa&#xff0c;你需要找一个最小的模数xxx&#xff0c;使得aaa中每个数都模上xxx之后互不相同。 n≤5e5,ai≤5e5,ai!ajn\le5e5,a_i\le5e5,a_i!a_jn≤5e5,ai​≤5e5,ai​!aj​ 思路&#xff1a…