poj 2175

题意

n栋大楼，每个大楼里ni个人，m个避难所，每个避难所能容纳mi个人。每个建筑物都有坐标，每个人到达避难所的费用是距离+1。给了你一种安排，问这种安排是不是最优的。

题解

根据消圈定理
如果一个费用流的残留网络，他如果存在有一个负环
那么这个肯定不是最小费用流
因为你可以用这个负环去更新答案
反过来也是一样，如果没有负环，那么就是最小费用流
然后这题是一个模板题
我们考虑根据他所给的条件，把残留网络建出来
然后检查一下就可以了

CODE:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
const int N=105*2;
int n,m;
struct qq
{
    int x,y,z;
}a[N],b[N];
int o[N];
int s[N][N];
int dis (int x,int y)   {return abs(a[x].x-b[y].x)+abs(a[x].y-b[y].y);}
struct qt
{
    int x,y,z,last;
}e[N*N];int num,last[N];
void init (int x,int y,int z)
{
//  printf("YES:%d %d %d\n",x,y,z);
    num++;
    e[num].x=x;e[num].y=y;e[num].z=z;
    e[num].last=last[x];
    last[x]=num;
}
int cnt[N];//每一个点被经过了多少次 
int h[N];
int from[N];
bool in[N];
int SPFA ()
{
    memset(in,false,sizeof(in));
    memset(from,-1,sizeof(from));
    memset(h,63,sizeof(h));
    memset(cnt,0,sizeof(cnt));
    int lim=n+m;
    queue<int> q;q.push(0);h[0]=0;
    while (!q.empty())
    {
        int x=q.front();q.pop();
        for (int u=last[x];u!=-1;u=e[u].last)
        {
            int y=e[u].y;
            if (h[y]>h[x]+e[u].z)
            {
                h[y]=h[x]+e[u].z;
                from[y]=x;
                cnt[y]++;
                if (in[y]==false)
                {
                    in[y]=true;
                    q.push(y);
                }
                if (cnt[y]>lim) return y;
            }
        }
        in[x]=false;
    }
    return -1;
}
int main()
{
    num=0;memset(last,-1,sizeof(last));
    scanf("%d%d",&n,&m);
    for (int u=1;u<=n;u++)  scanf("%d%d%d",&a[u].x,&a[u].y,&a[u].z);
    for (int u=1;u<=m;u++)  scanf("%d%d%d",&b[u].x,&b[u].y,&b[u].z);
    for (int u=1;u<=n;u++)
        for (int i=1;i<=m;i++)
        {
            scanf("%d",&s[u][i]);
            init(u,i+n,dis(u,i));
            o[i]=o[i]+s[u][i];
            if (s[u][i]!=0)//那么肯定存在反向边
                init(i+n,u,-dis(u,i));
        }
    for (int u=1;u<=m;u++)
        if (o[u]!=0)//那么肯定存在反向边 
        {
            init(0,u+n,0);
            if (o[u]<b[u].z)//还有流量 
                init(u+n,0,0);
        }
    int x=SPFA();
    if (x==-1) printf("OPTIMAL\n");
    else
    {
        printf("SUBOPTIMAL\n");
        memset(in,false,sizeof(in));
        while (in[x]==false)    {/*printf("%d\n",x);*/in[x]=true;x=from[x];}
        int last=x;
        do
        {
            int xx=from[x];
            //printf("%d %d\n",x,xx);
            //那么xx到x这是一条边
            if (xx>n&&x<=n) s[x][xx-n]--;
            if (xx<=n&&x>n) s[xx][x-n]++;
            x=xx;
        }while (x!=last);
        for (int u=1;u<=n;u++)
        {
            for (int i=1;i<=m;i++)
                printf("%d ",s[u][i]);
            printf("\n");
        }
    }
    return 0;
}