POJ 2396 有源汇的上下界可行流（好题）

（P.S.   POJ终于300题了，弱菜啊！！！）

题目链接： http://poj.org/problem?id=2396

题目大意：

　　给定一个矩阵每行，每列的和，和各个元素的限制条件(>, =, <)，求出一个满足这各种限制的矩阵。 

分析;

　　我做了前面两道sgu的入门题，发现看这题的解题报告就很好懂了，其实原理方面我还是没有lin神那么理解得透彻。

　　http://hi.baidu.com/zfy0701/blog/item/6449d82a64e15e3e5343c1ba.html

　　http://blog.csdn.net/water_glass/article/details/6823741

注意：

　　1、源S，汇T， 附加源ST（X）， 附加汇ED（Y），将有源汇的转换为无源汇的应该是在添加了附加源汇之后再连一条T->S的容量为无穷大的边。

　　2、这题我错了一天啊，究其原因非常神奇把input数据的处理随便改改就神奇地过了，也就是当输入x  y  =  t 时候，判断是否有错误的代码应该是：

　　　　if( x>c[u][v] || x<l[u][v] ) Flag= 0;就足够了，而我自己多加了一句话就wa到底了（具体见下面的代码吧），这不是说明即使一个题你会算法，可以说什么都会照样有你在某在状态下ac不了的题么？

代码：

/*2396    Accepted    924K    141MS    C++    3650B    2012-06-18 21:53:56*/
#include <cstdio>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <vector>
using namespace std;

#define mpair make_pair
#define pii pair<int,int>
#define MM(a,b) memset(a,b,sizeof(a));
typedef long long lld;
typedef unsigned long long u64;
template<class T> bool up_max(T& a,const T& b){return b>a? a=b,1 : 0;}
template<class T> bool up_min(T& a,const T& b){return b<a? a=b,1 : 0;}
#define maxn 240
const int inf= 2100000000;

int n,m,N;
int ST, ED;
int g[maxn][maxn];
int c[maxn][maxn], l[maxn][maxn];

bool Flag;

void change(int u,int v,char ch,int x){
    if( '=' == ch ){
        /// if( c[u][v]==l[u][v] ) Flag= 0;  // if have this code, it gets Wrong Answer!!!
        if( x>c[u][v] || x<l[u][v] ) Flag= 0;
        c[u][v]= l[u][v]= x;
    }
    else if( '<' == ch ){
        up_min( c[u][v], x-1 );
        if( c[u][v] < l[u][v] ) Flag= 0;
    }
    else if( '>' == ch ){
        up_max( l[u][v], x+1 );
        if( c[u][v] < l[u][v] ) Flag= 0;
    }
}

void Read(){
    char ch;
    int i,j,u,v,x,t;
    cin>>n>>m;
    N= n+m+1;
    MM( c, 0 );
    MM( l, 0 );
    for(i=1;i<=n;++i){
        scanf("%d", &t);
        c[0][i]= l[0][i]= t;
    }
    for(i=1;i<=m;++i){
        scanf("%d", &t);
        c[i+n][N]= l[i+n][N]= t;
    }
    for(i=1;i<=n;++i)for(j=1;j<=m;++j){ c[i][j+n]= inf; l[i][j+n]= 0; }
    cin>>t;
    while( t-- ){
        scanf("%d %d %c %d", &u, &v, &ch, &x );
        if( 0==u && 0==v ){
            for(i=1;i<=n;++i) for(j=1;j<=m;++j) change( i, j+n, ch, x );
        }
        else if( 0==u ){
            for(i=1;i<=n;++i) change( i, v+n, ch, x );
        }
        else if( 0==v ){
            for(i=1;i<=m;++i) change( u, i+n, ch, x );
        }
        else change( u, v+n, ch, x );
    }
}

int build_graph(){
    int sum= 0;
    ST= N+1, ED= N+2;
    MM( g, 0 );
    for(int i=0;i<=N;++i){
        int t= 0;
        for(int j=0;j<=N;++j){
            t+= l[j][i]-l[i][j];
            g[i][j]= c[i][j]-l[i][j];
        }
        if(t>0) sum+= g[ST][i]= t; /// t>0!!!;
        else g[i][ED]= -t;
    }
    g[N][0]= inf; ///!!!
    return sum;
}

bool vis[maxn];
int que[maxn], pre[maxn];
bool bfs(){
    MM( vis, 0 );
    int head=0, tail=0;
    que[tail++]= ST;
    vis[ST]= 1;
    while( head<tail ){
        int u= que[head++];
        for(int v=0;v<=ED;++v){
            if( g[u][v]>0 && !vis[v] ){
                pre[v]= u;
                if( v==ED ) return 1;
                que[tail++]= v;
                vis[v]= 1;
            }
        }
    }
    return 0;
}

int Edmond_karp(){
    int ret= 0;
    while( bfs() ){
        int t= inf;
        for(int i=ED;i!=ST;i=pre[i])
            up_min( t, g[pre[i]][i] );
        ret+= t;
        for(int i=ED;i!=ST;i=pre[i]){
            g[pre[i]][i]-= t;
            g[i][pre[i]]+= t;
        }
    }
    return ret;
}

int main()
{
    //freopen("poj2396.in","r",stdin);
    int Cas;
    cin>>Cas;
    while( Cas-- ){
        Flag= 1;
        Read();
        if( !Flag ){
            puts("IMPOSSIBLE");
            if( Cas>0 ) puts("");
            continue;
        }

        int sum= build_graph();
        int ret= Edmond_karp();
        if( ret<sum ) puts("IMPOSSIBLE");
        else{
            for(int i=1;i<=n;++i)
                for(int j=1;j<=m;++j)
                    //printf("%d%c", c[i][j+n]-g[i][j+n], j==m ? '\n' : ' ' );
                    printf("%d%c", g[j+n][i]+l[i][j+n], j==m ? '\n' : ' ' );
        }
        if( Cas>0 ) puts("");
    }
}