U41570 War1（网络流）

题目大意：

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