POJ2175

n 个建筑物 m 个避难所，每个建筑物有一定数量的人，每个避难所有一定的容量限制，每个人要跑到避难所里，代价为曼哈顿距离+1。给出一种疏散的方案，问是不是最优的。如果不是最优的，给出一个更优的方案。
$n,m\le 100$

按照给定的方案构建残余网络，如果网络上有负环，就说明该方案不是最优的。沿着负环增广就可以获得一个更优的方案。

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#define ll long long
#define pb push_back
#define fir first
#define sec second
#define ld long double
const int N=210,inf=1e9;
struct edge {
	int to,next,flow,cost;
}ed[N<<10];
int head[N],sz=1,x[N],y[N],c[N],vis[N],q[N<<13],deep[N],dis[N],pre[N],n,m,S,T,e[N][N],sum[N],vv[N];
void add_edge(int from,int to,int flow,int cost)
{
	ed[++sz]=(edge){to,head[from],flow,cost};
	head[from]=sz;
}
int gg(int a) {return a>0?a:-a;}
int read()
{
	int x=0;char c=getchar(),flag='+';
	while(!isdigit(c)) flag=c,c=getchar();
	while(isdigit(c)) x=x*10+c-'0',c=getchar();
	return flag=='-'?-x:x;
}
void solve(int x)
{
	memset(vis,0,sizeof(vis));
	while(!vis[x]) 
	{
		vis[x]=1;
		x=ed[pre[x]^1].to;
	}
	memset(vis,0,sizeof(vis));
	while(!vis[x])
	{
		vis[x]=1;
		int j=pre[x];
		ed[j].flow--,ed[j^1].flow++;
		x=ed[j^1].to;
	}
}
bool spfa(int s)
{
	memset(dis,0x3f,sizeof(dis));
	memset(vis,0,sizeof(vis));
	memset(deep,0,sizeof(deep));
	int h=1,t=0;
	q[++t]=s;
	dis[s]=0;
	vis[s]=1;
	while(h<=t)
	{
		int u=q[h++];
		vis[u]=0;
		vv[u]=1;
		for(int i=head[u];i;i=ed[i].next)
		{
			int v=ed[i].to;
			if(ed[i].flow&&dis[u]+ed[i].cost<dis[v])
			{
				dis[v]=dis[u]+ed[i].cost;
				pre[v]=i;
				deep[v]=deep[u]+1;
				if(deep[v]>2*(n+m+2))
				{
					solve(v);
					return true;
				}
				if(!vis[v]) vis[v]=1,q[++t]=v;
			}
		}
	}
	return false;
}
int main()
{
	n=read(),m=read(),S=0,T=n+m+1;
	for(int i=1;i<=n;i++) x[i]=read(),y[i]=read(),c[i]=read();
	for(int i=1;i<=m;i++) x[i+n]=read(),y[i+n]=read(),c[i+n]=read();
	for(int i=1;i<=n;i++)
	{
		add_edge(S,i,0,0);
		add_edge(i,S,c[i],0);
		for(int j=1;j<=m;j++)
		{
			int w=read();
			add_edge(i,j+n,inf,gg(x[i]-x[j+n])+gg(y[i]-y[j+n])+1);
			add_edge(j+n,i,w,-gg(x[i]-x[j+n])-gg(y[i]-y[j+n])-1);
			sum[j]+=w;
			e[i][j]=sz;
		}
	}
	for(int i=1;i<=m;i++) add_edge(i+n,T,c[i+n]-sum[i],0),add_edge(T,i+n,sum[i],0);
	for(int i=0;i<=n+m+1;i++) if(!vv[i]&&spfa(i))
	{
		puts("SUBOPTIMAL");
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=m;j++) std::cout<<ed[e[i][j]].flow<<' ';
			std::cout<<'\n';
		}
		return 0;
	}
	puts("OPTIMAL");
	return 0;
}