社区讨论

本地过了,洛谷全WA

P4015运输问题参与者 1已保存回复 1

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
1 条
当前快照
1 份
快照标识符
@m0g1tos2
此快照首次捕获于
2024/08/30 09:42
2 年前
此快照最后确认于
2024/08/30 13:36
2 年前
查看原帖
RT,蒟蒻求调
CPP
#include <bits/stdc++.h>
#define min(a,b) a<b?a:b
using namespace std;
typedef long long ll;
typedef pair<ll,ll> pll;
const int MAXN=3e5+10;
const int inf=2e9;
namespace flow {
	int S,T,n,m,cnt,head[MAXN],nxt[MAXN],to[MAXN],use[MAXN],pre[MAXN],vis[MAXN];
	ll ansmx,ansmn,dis[MAXN],flow[MAXN],cap[MAXN],val[MAXN];
	void add(int x,int y,int w,int c) {
		to[++cnt]=y,cap[cnt]=w,val[cnt]=c;
		nxt[cnt]=head[x];
		head[x]=cnt;
	}
	void add_edge(int x,int y,int w,int c) {
		add(x,y,w,c);
		add(y,x,0,-c);
	}
	void init() {
		for(int i=0; i<=MAXN; i++) head[i]=use[i]=flow[i]=0;
		cnt=1,S=0,T=m+n+1;
		ansmx=ansmn=0;
	}
	bool bfs() {
		queue<ll> q;
		for(int i=S; i<=MAXN; i++) pre[i]=-1,vis[i]=0,flow[i]=1e15,dis[i]=1e15;
		q.push(S),vis[S]=1,dis[S]=0;
		while(!q.empty()) {
			ll x=q.front();
			q.pop();
			vis[x]=0;
			for(int i=head[x]; i; i=nxt[i]) {
				ll y=to[i];
				if(dis[y]<=dis[x]+val[i]||cap[i]==0) continue;
				dis[y]=dis[x]+val[i];
				flow[y]=min(cap[i],flow[x]);
				pre[y]=x,use[y]=i;
				if(!vis[y]) vis[y]=1,q.push(y);
			}
		}
		return pre[T]!=-1;
	}
	pll dinic() {
		while (bfs()) {
			ll x=T;
			ansmx+=flow[T];
			ansmn+=flow[T]*dis[T];
			while(x) {
				cap[use[x]]-=flow[T];
				cap[use[x]^1]+=flow[T];
				x=pre[x];
			}
		}
		return make_pair(ansmx,ansmn);
	}
}
int a[MAXN],b[MAXN],c[2000][2000];
int main() {
	ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
	cin>>flow::m>>flow::n;
	flow::cnt=1,flow::S=0,flow::T=flow::m+flow::n+1;
	for(int i=1; i<=flow::m; i++) {
		cin>>a[i];
		flow::add_edge(flow::S,i,a[i],0);
	}
	for(int i=1; i<=flow::n; i++) {
		cin>>b[i];
		flow::add_edge(i+flow::m,flow::T,b[i],0);
	}
	for(int i=1; i<=flow::m; i++) {
		for(int j=1; j<=flow::n; j++) {
			cin>>c[i][j];
			flow::add_edge(i,j+flow::m,inf,c[i][j]);
		}
	}
	cout<<flow::dinic().second<<'\n';
	flow::init();
	for(int i=1;i<=flow::m;i++) flow::add_edge(flow::S,i,a[i],0);
	for(int i=1;i<=flow::n;i++) flow::add_edge(i+flow::m,flow::T,b[i],0);
	for(int i=1;i<=flow::m;i++) 
	   for(int j=1;j<=flow::n;j++) 
	       flow::add_edge(i,j+flow::m,inf,-c[i][j]);
	cout<<-flow::dinic().second;
	return 0;
}

回复

1 条回复,欢迎继续交流。

正在加载回复...