Soldier and Traveling

B. Soldier and Traveling

Time Limit: 1000ms

Memory Limit: 262144KB

64-bit integer IO format: %I64d      Java class name: (Any)

Submit Status

In the country there are n cities and m bidirectional roads between them. Each city has an army. Army of the i-th city consists of aisoldiers. Now soldiers roam. After roaming each soldier has to either stay in his city or to go to the one of neighboring cities by atmoving along at most one road.

Check if is it possible that after roaming there will be exactly bi soldiers in the i-th city.

 

Input

First line of input consists of two integers n and m (1 ≤ n ≤ 100, 0 ≤ m ≤ 200).

Next line contains n integers a1, a2, ..., an (0 ≤ ai ≤ 100).

Next line contains n integers b1, b2, ..., bn (0 ≤ bi ≤ 100).

Then m lines follow, each of them consists of two integers p and q (1 ≤ p, q ≤ n, p ≠ q) denoting that there is an undirected road between cities p and q.

It is guaranteed that there is at most one road between each pair of cities.

 

Output

If the conditions can not be met output single word "NO".

Otherwise output word "YES" and then n lines, each of them consisting of n integers. Number in the i-th line in the j-th column should denote how many soldiers should road from city i to city j (if i ≠ j) or how many soldiers should stay in city i (if i = j).

If there are several possible answers you may output any of them.

 

Sample Input

Input

4 4
1 2 6 3
3 5 3 1
1 2
2 3
3 4
4 2

Output

YES
1 0 0 0 
2 0 0 0 
0 5 1 0 
0 0 2 1 

Input

2 0
1 2
2 1

Output

NO
思路：最大流；
首先判断变前和变后的和是否相等，如果不等则直接输出NO，否则，转换为最大流求解，原点和原来的点连边权值为原来的人口，然后每个新的状态和汇点连边，权值为后来的人口，然后
按给的边连边，权值为原来的人口，然后跑最大流，判断最大流量是否为sum。最后的矩阵由反边得来，表示从上个点有人口转移而来，反边的值就是转移人口。

#include<stdio.h>
#include<algorithm>
#include<iostream>
#include<stdlib.h>
#include<queue>
#include<string.h>
#include<map>
#include<vector>
using namespace std;
typedef long long LL;
typedef struct pp
{
    int to;
    int cap;
    int rev;
} aa;
vector<pp>vec[205];
int level[205];
int iter[205];
int ans[105];
int bns[105];
void add(int from,int to,int cap);
void bfs(int s);
int dfs(int s,int t,int f);
int max_flow(int s,int t);
const int N = 1e9;
int ma[200][200];
int main(void)
{
    int n,m;
    scanf("%d %d",&n,&m);
    int i,j;
    int sum1 = 0;
    int sum2 = 0;
    for(i = 1; i <= n; i++)
    {
        scanf("%d",&ans[i]);
        sum1 += ans[i];
    }
    for(i = 1; i <= n; i++)
    {
        scanf("%d",&bns[i]);
        sum2 += bns[i];
    }
    if(sum1 != sum2)
        printf("NO
");
    else
    {
      for(i = 1;i <= n; i++)
      {
          add(0,i,ans[i]);
          add(i+n,2*n+1,bns[i]);
          add(i,i+n,ans[i]);
      }
      while(m--)
      {
          int x,y;
          scanf("%d %d",&x,&y);
          add(x,y+n,ans[x]);
          add(y,x+n,ans[y]);
      }
      int ask = max_flow(0,2*n+1);
      if(ask != sum1)
        printf("NO
");
      else
      {
          for(i = 1;i <= n;i++)
          {
              int x = i+n;
              for(j = 0;j < vec[x].size(); j++)
              {
                  aa no = vec[x][j];
                  ma[no.to][i] = no.cap;
              }
          }printf("YES
");
          for(i = 1;i <= n; i++)
          {
              for(j = 1;j <= n; j++)
              {
                  if(j == 1)
                    printf("%d",ma[i][j]);
                  else printf(" %d",ma[i][j]);
              }
              printf("
");
          }
      }
    }return 0;
}
void add(int from,int to,int cap)
{
    pp  nn;
    nn.to = to;
    nn.cap = cap;
    nn.rev = vec[to].size();
    vec[from].push_back(nn);
    nn.to = from;
    nn.cap=0;
    nn.rev = vec[from].size()-1;
    vec[to].push_back(nn);
}
void bfs(int s)
{
    queue<int>que;
    memset(level,-1,sizeof(level));
    level[s]=0;
    que.push(s);
    while(!que.empty())
    {
        int v=que.front();
        que.pop();
        int i;
        for(i=0; i<vec[v].size(); i++)
        {
            pp e=vec[v][i];
            if(level[e.to]==-1&&e.cap>0)
            {
                level[e.to]=level[v]+1;
                que.push(e.to);
            }
        }
    }
}
int dfs(int s,int t,int f)
{
    if(s==t)
        return f;
    for(int &i=iter[s]; i<vec[s].size(); i++)
    {
        pp &e=vec[s][i];
        if(level[e.to]>level[s]&&e.cap>0)
        {
            int r=dfs(e.to,t,min(e.cap,f));
            if(r>0)
            {
                e.cap-=r;
                vec[e.to][e.rev].cap+=r;
                return r;
            }
        }
    }
    return 0;
}
int max_flow(int s,int t)
{
    int flow=0;
    for(;;)
    {
        bfs(s);
        if(level[t]<0)return flow;
        memset(iter,0,sizeof(iter));
        int f;
        while((f=dfs(s,t,N)) >0)
        {
            flow += f;
        }
    }
}