Parco最近迷上了探險。第三站,他來到了阿帕帕烏克咵哇部落。
在這個部落,有一個很流行的遊戲叫哇咵克烏帕帕阿。
這個遊戲是使用一塊名叫蟶磣瀚槮的石塊,如下圖,
簡單來說其實就是一個填數遊戲。蟶磣瀚槮石塊分成兩部分組成,主體部分和紅色部分。主體部分是一個N×MN \times MN×M個格子的矩陣,每個格子裏面或許已有數字,也有可能爲空,空的則需要填數,每個格子能填入0∽90 \backsim 90∽9之間任意一個數字。紅色部分,頂上那行,表示每一列的和,左邊那列,表示每一行的和。這個遊戲就是給出紅色部分,以及主體部分(部分爲空或者完全爲空),然後根據紅色部分的要求,把主體部分完成。下圖就是上圖的一種解。
阿帕帕烏克咵哇部落有一個傳統,就是每位來賓都必須通過哇咵克烏帕帕阿遊戲,才能進入部落。
Parco很想進入部落,但是他不知道如何填數,於是他把手上的蟶磣瀚槮石塊拍了張照片給你,希望你能幫他解決。
輸入格式
輸入數據包含多組數據。
對於每組數據:
第一行爲兩個整數NNN、MMM,表示主體部分的大小。
接下來NNN行,每行爲一個長度爲MMM的字符串。‘0’-‘9’表示主體部分原有的數字,‘#’表示主體部分需要填數的位置。(具體可看樣例)
接下來一行,包含MMM個整數,表示紅色部分頂上那行。
接下來一行,包含NNN個整數,表示紅色部分左邊那列。
(1≤N,M≤10001 \leq N,M \leq 10001≤N,M≤1000)
輸出格式
對於每組數據:
輸出NNN行,每行爲長度爲MMM的字符串,表示填好的主體部分。(具體可看樣例)。
數據保證有解,若有多種解,則輸出任意一種即可。
樣例數據
輸入
4 5
3#786
#68#0
52##8
#049#
18 13 24 24 21
29 25 23 23
輸出
35786
76840
52538
30497
也是isap,只不過是/*ISAP鄰接表形式
#include<bits/stdc++.h>
using namespace std;
struct node
{
int to,cap,rev;
};
const int MAXN=0x3f3f3f3f;
vector<node>edge[1100000];
queue<int>que;
int level[1100000],iter[1100000];
int tot;
void add_edge(int u,int v,int f)
{
edge[u].push_back((node){v,f,(int)edge[v].size()});
edge[v].push_back((node){u,0,(int)edge[u].size()-1});
}
int dfs(int s,int t,int f)
{
if (s==t) return f;
for (int &i=iter[s];i<edge[s].size();i++)
{
node tmp=edge[s][i];
if (tmp.cap>0&&level[tmp.to]>level[s])
{
int d=dfs(tmp.to,t,min(f,tmp.cap));
if (d>0)
{
edge[s][i].cap-=d;
edge[tmp.to][tmp.rev].cap+=d;
return d;
}
}
}
return 0;
}
void bfs(int s)
{
level[s]=0;
que.push(s);
while (!que.empty())
{
s=que.front();
que.pop();
for (int i=0;i<edge[s].size();i++)
{
node tmp=edge[s][i];
if (level[tmp.to]==-1&&tmp.cap>0)
{
level[tmp.to]=level[s]+1;
que.push(tmp.to);
}
}
}
}
int max_flow(int s,int t)
{
int flow=0;
while (true)
{
memset(level,-1,sizeof(level));
memset(iter,0,sizeof(iter));
bfs(s);
if (level[t]==-1) break;
int f;
f=dfs(s,t,MAXN);
while (f>0)
{
flow+=f;
f=dfs(s,t,MAXN);
}
}
return flow;
}
char s[1010][1010];
int id[1010][1010];
int sum1[1010],sum2[1010];
map<pair<int,int>,int>mp;
int main()
{
int n,m;
while (~scanf("%d%d",&n,&m))
{
for (int i=0;i<1100000;i++) edge[i].clear();
for (int i=1;i<=n;i++) scanf("%s",s[i]+1);
for (int i=1;i<=m;i++) scanf("%d",&sum1[i]);
for (int i=1;i<=n;i++) scanf("%d",&sum2[i]);
int tot=0;
for (int i=1;i<=n;i++)
{
for (int j=1;j<=m;j++)
{
if (s[i][j]=='#')
{
add_edge(i,j+n,9);
}
else
{
sum1[j]-=s[i][j]-'0';
sum2[i]-=s[i][j]-'0';
}
}
}
int S=0,T=n+m+2;
for (int i=1;i<=n;i++) add_edge(S,i,sum2[i]);
for (int i=1;i<=m;i++) add_edge(i+n,T,sum1[i]);
int flow=max_flow(S,T);
for (int i=1;i<=n;i++)
{
for(int k=0;k<edge[i].size();++k){
int v=edge[i][k].to;
int cost=edge[i][k].cap;
//printf("%d %d\n",v,cost);
if(v>n)s[i][v-n]='0'+9-cost;
}
}
for (int i=1;i<=n;i++)
{
for (int j=1;j<=m;j++) printf("%c",s[i][j]);
printf("\n");
}
}
}