#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
const int maxn=205;
int n,m,g[maxn][maxn],flow[maxn],pre[maxn];
bool bfs(int s,int t,int v){//搜索增广路,源点,汇点,顶点数
memset(pre,-1,sizeof(pre));
memset(flow,inf,sizeof(flow));
queue<int> q;
q.push(s);
pre[s]=0;//pre[s]==0表示s已经走过,避免重复算
while(!q.empty()){
int now=q.front();
if(now==t) break;
for(int i=1;i<=v;i++){//v标记图的顶点数
if(pre[i]==-1&&g[now][i]>0){//pre[]==-1表示未走过,g[][]>0表示正向有路
pre[i]=now;
flow[i]=min(flow[now],g[now][i]);
q.push(i);
}
}
q.pop();
}
return pre[t]!=-1;//pre[t]==-1表示t没有前驱节点即找不到增广路
}
void update(int s,int t){//更新残留网络
for(int fa,now=t;now!=s;now=fa){
fa=pre[now];
g[fa][now]-=flow[t];
g[now][fa]+=flow[t];
}//f[t]标记这条路上最小的流量,所以整条路都用t
}
int main()
{
int t,x,sum,maxf;
char str[maxn];
cin>>t;
while(t--){
memset(g,0,sizeof(g));
sum=0;
maxf=0;
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%d",&x);
sum+=x;
g[n+m+1][i]=x;
}
for(int i=1;i<=m;i++){
scanf("%d",&x);
g[n+i][n+m+2]=x;
}
for(int i=1;i<=n;i++){
scanf("%s",str+1);
for(int j=1;j<=m;j++) if(str[j]=='1') g[i][n+j]=inf;
}
while(bfs(n+m+1,n+m+2,n+m+2)){
maxf+=flow[n+m+2];
update(n+m+1,n+m+2);
}
printf("%d\n",sum-maxf);
}
return 0;
}
偷钻石(网络流)
最新推荐文章于 2024-10-30 15:52:39 发布