【NOI2008】志愿者招募(费用流)

传送门


题解:

很裸的费用流。

建图也很显然。

从源点向第一天连容量为 I N F INF INF的边,从第 n + 1 n+1 n+1天向汇点连 I N F INF INF的边。

然后从第 i i i天向第 i + 1 i+1 i+1天连容量为 I N F − A [ i ] INF-A[i] INFA[i]的边。

现在的问题就是,我们需要花一定的代价对流量进行补偿,对于一类志愿者 ( s , t , c ) (s,t,c) (s,t,c),我们从第 s s s天向第 t + 1 t+1 t+1天连容量为 I N F INF INF,费用为 c c c的边。

然后跑最小费用最大流。


代码:

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

cs int N=1e3+7;
cs int INF=1e9;

struct edge{int to,rev,cap,w;};
typedef std::vector<edge>::iterator iter;
std::vector<edge> G[N];iter cur[N];

inline void adde(int u,int v,int val,int cost){
	G[u].push_back((edge){v,G[v].size(),val,cost});
	G[v].push_back((edge){u,G[u].size()-1,0,-cost});
}

int n,m,S,T,tot;

int dis[N];bool vis[N];
inline bool SPFA(){
	memset(dis,0x3f,sizeof(int)*(tot+1));
	std::queue<int> q;q.push(S);dis[S]=0;
	while(!q.empty()){
		int u=q.front();q.pop();
		vis[u]=false;cur[u]=G[u].begin();
		for(iter e=G[u].begin();e!=G[u].end();++e)
		if(e->cap&&dis[u]+e->w<dis[e->to]){
			dis[e->to]=dis[u]+e->w;
			if(!vis[e->to])vis[e->to]=true,q.push(e->to);
		}
	}
	return dis[T]<INF;
}

int tot_cost,tot_flow;
int dfs(int u,int flow){
	if(u==T){
		tot_cost+=flow*dis[T];
		tot_flow+=flow;
		return flow;
	}
	int ans=0;vis[u]=true;
	for(iter &e=cur[u];e!=G[u].end();++e)
	if(e->cap&&dis[e->to]==dis[u]+e->w&&!vis[e->to]){
		int delta=dfs(e->to,std::min(flow-ans,e->cap));
		if(delta){
			e->cap-=delta;
			G[e->to][e->rev].cap+=delta;
			if((ans+=delta)==flow)break;
		}
	}
	vis[u]=false;return ans;
}

inline void Flow(){
	tot_cost=tot_flow=0;
	while(SPFA())dfs(S,INF);
}

signed main(){
#ifdef zxyoi
	freopen("volunteer.in","r",stdin);
#endif
	scanf("%d%d",&n,&m);S=n+2,T=tot=S+1;
	adde(S,1,INF,0);adde(n+1,T,INF,0);
	for(int re i=1;i<=n;++i){
		int w;scanf("%d",&w);
		adde(i,i+1,INF-w,0);
	}
	while(m--){
		int s,t,c;scanf("%d%d%d",&s,&t,&c);
		adde(s,t+1,INF,c);
	}
	Flow();std::cout<<tot_cost<<"\n";
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值