新生赛第二场 C题

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");
        }
    }
}