社区讨论
本地过了,洛谷全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 条回复,欢迎继续交流。
正在加载回复...