[BJWC2018]Kakuro

最新推荐文章于 2021-11-17 22:27:34 发布

StaroForgin

最新推荐文章于 2021-11-17 22:27:34 发布

阅读量284

点赞数 1

分类专栏： ------图论------ # 网络流

本文链接：https://blog.csdn.net/Tan_tan_tann/article/details/113726471

版权

------图论------ 同时被 2 个专栏收录

101 篇文章 0 订阅

订阅专栏

网络流

20 篇文章 0 订阅

订阅专栏

Kakuro

题解

~~多水的一道题呀~~

其实看到题目应该是很容易想到网络流的。
由于它要求的是每列或每行的总和一定，我们可以先找出一组可行解，再来再这组解上修改。
观察到对于一个数，它无论增加还是减少它的单位花费都是一定的，所以我们可以先找出一组最小的可行解，这样之后就处理增加多少。
而最小的可行解明显是所有空格都取1（题目要求的是正整数），所有线索都取它右侧或下侧的方格数。我们可以先预处理出这样的解它的答案是多少。
对于不能变动的点我们就把它的变动的花费设为inf，只需要最后特判一下，不会影响到最后的答案。

我们考虑再这组解上做改动来找出最优解。
由于这组解已经是最小的了，我们只考虑每个空格或线索增加的情况。
很容易发现，当一个空格增加时它会影响到它所对应的上方与左方线索。由于它的增加量与这些线索的增加量是一样的，我们可以将它想象成网络流的一条边，连接上方线索与左方线索两个点。
而对于它所对应上方线索与左方线索，他们的性质一定是不一样的，我们可以将其分列两边，分别与源点与汇点相连。
对于任意一条经过这个空格的路径，它一定经过了不超过一个的上方线索和左方线索，它的流量表示它们共同增加了多少(与源点相连的边流量表示上方线索的增加量，两个线索的连边的流量代表它们共同对应的空格的增加量，与汇点相连的边的流量代表左方线索的增加量)，故它的费用也该对应增加这么多流量的费用。

我们注意到对于这些最小解中的线索与空格，它的增加一定是先填补之前减少的，再开始从初始状态增加（也由可能没有第一个阶段）。
对于第一个阶段，它的代价为 $cost_{i,j}$ ，会增加 $A_{i,j}-B_{i,j}$ 的流量， $A_{i,j}$ 代表它的初始值， $B_{i,j}$ 代表它在我们求出的可行解的值。对于这个阶段，我们再两者直接连一条流量为 $A_{i,j}-B_{i,j}$ ，代价为 $cost_{i,j}$ 的边。
对于第二个阶段，它的代价就是 $cost_{i,j}$ ，会增加的流量是inf，因为它可以增加到无限大。我们也只需要像第一个阶段一样，在它所需要连同的两个点间连接即可。
对于那些没有对应两个的空格，我们可以直接在源点或汇点上去连边。

只要我们像上面说的一样去建图，再去跑最小费用可行流即可。
由于本题的所有边的下限都是 $0$ ，我们只需要将所有长度小于0的增广路径跑了即可，用 $E K$ 即可解决。
对于合法的图，即无论如何也找不出一组解的图，我们只需要找找有没有权值为inf的边有流量经过，或为**-inf的边无流量经过即可。这里的inf**指我们最开始将不能改的边设的权值。

时间复杂度 $\left(n^3m^3\right)$ ？反正网络流的时间复杂度一点都不准就对了。~~应该也没人故意卡你网络流的时间复杂度。~~

源码

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
#define MAXM 50005
#define MAXN 35
#define reg register
typedef long long LL;
typedef pair<int,LL> pii;
const int INF=0x7f7f7f7f;
const LL inf=0x7f7f7f7f7f7f;
template<typename _T>
inline _T Fabs(_T x){return x<0?-x:x;}
template<typename _T>
inline void read(_T &x){
	_T f=1;x=0;char s=getchar();
	while('0'>s||'9'<s){if(s=='-')f=-1;s=getchar();}
	while('0'<=s&&s<='9'){x=(x<<3)+(x<<1)+(s^48);s=getchar();}
	x*=f;
}
int n,m,head[MAXM],tot,pre[MAXM],pw[MAXM],S,T,idx,cnt;LL dis[MAXM];
int l1[MAXN][MAXN],l2[MAXN][MAXN],lg[MAXN][MAXN],rg[MAXN][MAXN];
struct ming{int typ,val,vL,vR,cL,cR,cost,idl,idr,id;}s[MAXN][MAXN];
queue<int> q;bool vis[MAXM];
struct edge{int to,nxt,flow,paid,op;}e[MAXM*10];
inline void addEdge(const int u,const int v,const int w,const int f){e[++tot]=(edge){v,head[u],f,w};head[u]=tot;}
inline void addedge(const int u,const int v,const int w,const int f){addEdge(u,v,w,f);e[tot].op=tot+1;addEdge(v,u,-w,0);e[tot].op=tot-1;}
inline pii sakura(){
	int mf=0;LL res=0;
	while(1){
		while(!q.empty())q.pop();
		for(reg int i=1;i<=cnt;++i)dis[i]=inf;
		q.push(S);dis[S]=0;vis[S]=1;
		while(!q.empty()){
			int u=q.front();q.pop();vis[u]=0;
			for(int i=head[u];i;i=e[i].nxt){
				int v=e[i].to;
				if(dis[v]>dis[u]+1ll*e[i].paid&&e[i].flow>0){
					dis[v]=dis[u]+1ll*e[i].paid;pre[v]=u;pw[v]=i;
					if(!vis[v])q.push(v),vis[v]=1;
				}
			}
		}
		if(dis[T]>=0)break;int minn=INF;
		for(reg int i=T;i^S;i=pre[i])minn=min(minn,e[pw[i]].flow);mf+=minn;res+=1ll*minn*dis[T];
		for(reg int i=T;i^S;i=pre[i])e[pw[i]].flow-=minn,e[e[pw[i]].op].flow+=minn;
	}
	return make_pair(mf,res);
}
signed main(){
	read(n);read(m);LL ans=0;bool flag=0;
	for(reg int i=1;i<=n;++i)for(reg int j=1;j<=m;++j)read(s[i][j].typ);
	for(reg int i=1;i<=n;++i)
		for(reg int j=1;j<=m;++j){
			if(s[i][j].typ&1)read(s[i][j].vL),s[i][j].idl=++idx;
			if(s[i][j].typ&2)read(s[i][j].vR),s[i][j].idr=++idx;
			if(s[i][j].typ&4)read(s[i][j].val),s[i][j].id=++idx;	
		}
	S=++idx;cnt=T=++idx;
	for(reg int i=n;i>0;--i)
		for(reg int j=m;j>0;--j)
			if(s[i][j].typ&4)l1[i][j]=l1[i+1][j]+1,l2[i][j]=l2[i][j+1]+1;
	for(reg int i=1;i<=n;++i)
		for(reg int j=1;j<=m;++j){
			if(s[i][j].typ&1)read(s[i][j].cL);if(s[i][j].cL<0)s[i][j].cL=INF;
			if(s[i][j].typ&2)read(s[i][j].cR);if(s[i][j].cR<0)s[i][j].cR=INF;
			if(s[i][j].typ&4)read(s[i][j].cost);if(s[i][j].cost<0)s[i][j].cost=INF;
			if(s[i][j].typ&1)ans+=1ll*s[i][j].cL*Fabs(s[i][j].vL-l1[i+1][j]);
			if(s[i][j].typ&2)ans+=1ll*s[i][j].cR*Fabs(s[i][j].vR-l2[i][j+1]);
			if(s[i][j].typ&4)ans+=1ll*s[i][j].cost*Fabs(s[i][j].val-1);
			if(s[i][j].typ&1&&s[i][j].cL>INF-1&&s[i][j].vL<l1[i+1][j])flag=1;
			if(s[i][j].typ&2&&s[i][j].cR>INF-1&&s[i][j].vR<l2[i][j+1])flag=1;
		}
	if(flag){puts("-1");return 0;}
	for(reg int i=1;i<=n;++i)
		for(reg int j=1;j<=m;++j){
			if(s[i][j].typ&1)addedge(S,s[i][j].idl,s[i][j].cL,INF),lg[i][j]=s[i][j].idl,
				addedge(S,s[i][j].idl,-s[i][j].cL,s[i][j].vL-l1[i+1][j]);
			if(s[i][j].typ&2)addedge(s[i][j].idr,T,s[i][j].cR,INF),rg[i][j]=s[i][j].idr,
				addedge(s[i][j].idr,T,-s[i][j].cR,s[i][j].vR-l2[i][j+1]);
		}
	for(reg int i=1;i<=n;++i)
		for(reg int j=1;j<=m;++j){
			if(!(s[i][j].typ&4))continue;lg[i][j]=lg[i-1][j];rg[i][j]=rg[i][j-1];
			if(lg[i][j]&&!rg[i][j])addedge(lg[i][j],T,-s[i][j].cost,s[i][j].val-1),addedge(lg[i][j],T,s[i][j].cost,INF); 
			if(!lg[i][j]&&rg[i][j])addedge(S,rg[i][j],-s[i][j].cost,s[i][j].val-1),addedge(S,rg[i][j],s[i][j].cost,INF);
			if(lg[i][j]&&rg[i][j])addedge(lg[i][j],rg[i][j],-s[i][j].cost,s[i][j].val-1),addedge(lg[i][j],rg[i][j],s[i][j].cost,INF);
		}
	ans=ans+sakura().second;
	for(reg int i=1;i<=tot;i+=2){
		if(e[i].paid>INF-1&&e[e[i].op].flow){flag=1;break;}
		if(e[i].paid<-INF+1&&e[i].flow){flag=1;break;}
	}
	if(flag){puts("-1");return 0;}printf("%lld\n",ans);
	return 0;
}

谢谢！！！

StaroForgin

关注

1
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
[BJWC2018]Kakuro

Kakuro题解多水的一道题呀其实看到题目应该是很容易想到网络流的。由于它要求的是每列或每行的总和一定，我们可以先找出一组可行解，再来再这组解上修改。观察到对于一个数，它无论增加还是减少它的单位花费都是一定的，所以我们可以先找出一组最小的可行解，这样之后就处理增加多少。而最小的可行解明显是所有空格都取1（题目要求的是正整数），所有线索都取它右侧或下侧的方格数。我们可以先预处理出这样的解它的答案是多少。对于不能变动的点我们就把它的变动的花费设为inf，只需要最后特判一下，不会影响到最后的答案。
复制链接

扫一扫