题目大意:
ENLIGHTENED总部有NN个Portal,编号为11~NN,编号为ii的Portal初始能量值为A[i]A[i],在Portal之间有MM条LINK,每条LINK着连接着两个不同Portal,被连接着的两个Portal可以相互传输能量,每个Portal最多总共只能向其连接着的Portal传输A[i]A[i]点能量,现在ENLIGHTENED行动指挥想让每第ii个Portal的能量值变为B[i]B[i],但他不知道着可不可行,所以找到了你,如果可行的话,需要你求出一种可行的能量传输方案。
能量只能之间传输不能间接传输
思路:
这题很明显的网络流,考虑怎么建图,我们每个点自身建一次点,然后分别连向原点和汇点,然后如果X,Y可以传输能量就X-Y*,Y-X*连边,然后直接最大流就好了,最后如果满流就是可以,然后每条边的回流就是传输了多少流量。
程序:
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
const int N=505;
const int INF=1e9;
int x,y,cnt,sum,ans,s,t,n,m;
int a[N],b[N],last[N],dep[N],cur[N],f[N][N];
queue<int>q;
struct data{int to,w,next;}e[N*100];
void add(int x,int y,int w){
e[++cnt].to=y; e[cnt].next=last[x]; e[cnt].w=w; last[x]=cnt;
e[++cnt].to=x; e[cnt].next=last[y]; e[cnt].w=0; last[y]=cnt;
}
bool bfs(){
while (!q.empty()) q.pop();
q.push(s);
memset(dep,0,sizeof(dep));
dep[s]=1;
while (!q.empty()){
int x=q.front();
q.pop();
for (int i=last[x];i;i=e[i].next)
if (e[i].w&&!dep[e[i].to]){
dep[e[i].to]=dep[x]+1;
q.push(e[i].to);
if (e[i].to==t) return 1;
}
}
return 0;
}
int dfs(int x,int maxf){
if (x==t||!maxf) return maxf;
int ret=0;
for (int &i=cur[x];i;i=e[i].next)
if (e[i].w&&dep[e[i].to]==dep[x]+1){
int f=dfs(e[i].to,min(maxf-ret,e[i].w));
e[i].w-=f;
e[i^1].w+=f;
ret+=f;
if (ret==maxf) break;
}
return ret;
}
void dinic(){
while (bfs()){
for (int i=s;i<=t;i++) cur[i]=last[i];
ans+=dfs(s,INF);
}
}
int main(){
scanf("%d%d",&n,&m);
s=0; t=2*n+1; cnt=1;
for (int i=1;i<=n;i++){
scanf("%d",&a[i]);
add(s,i,a[i]);
}
for (int i=1;i<=n;i++){
scanf("%d",&b[i]);
add(i+n,t,b[i]);
sum=sum+b[i];
}
for (int i=1;i<=n;i++) add(i,n+i,INF);
for (int i=1;i<=m;i++){
scanf("%d%d",&x,&y);
add(x,n+y,INF);
add(y,n+x,INF);
}
dinic();
for (int i=1;i<=n;i++)
for (int j=last[i];j;j=e[j].next)
if (e[j].to>n&&e[j].to<=2*n) f[i][e[j].to-n]=e[j^1].w;
if (ans!=sum) printf("NO\n");
else {
printf("YES\n");
for (int i=1;i<=n;i++){
for (int j=1;j<=n;j++)
printf("%d ",f[i][j]);
printf("\n");
}
}
}