省选专练「网络流 24 题」星际转移

分层图上网络流。

以时间为键值开点。

对于每一个仓库,(i,time-1)->(i,time)(INF)表示可以放永远。

对于动态的飞船(a[i][time-1],time-1)->(a[i][time],time)(w[i])表示这一天可以搬这么多

那么,枚举时间,在参量网络上继续跑最大流,直到跑满就完了。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#define C c
using namespace std;
const int N=20000;
const int INF=1e9+7;
struct Front_star{
	int u,v,w,nxt;
}e[N*4];
int cnt=1;
int first[N]={0};
void addedge(int u,int v,int w){
	cnt++;
	e[cnt].u=u;
	e[cnt].w=w;
	e[cnt].v=v;
	e[cnt].nxt=first[u];
	first[u]=cnt;
}
void add(int u,int v,int w){
	addedge(u,v,w);
	addedge(v,u,0);
}
int n,m,k;
int S=9201;
int T=1;
int w[N]={0};
int a[101][101]={0};
int fa[101]={0};
int C[101]={0};//record a new start
int getfa(int x){
	if(x==fa[x]){
		return x;
	}
	else{
		return fa[x]=getfa(fa[x]);
	}
}
void uni(int x,int y){
	int dx=getfa(x);
	int dy=getfa(y);
	fa[dx]=dy;
}
int d[50001]={0};
int bfs(){
	memset(d,-1,sizeof(d));
	queue<int> q;
	q.push(S);
	d[S]=1;
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=first[x];i;i=e[i].nxt){
			int v=e[i].v;
			if(e[i].w>0&&d[v]==-1){
				d[v]=d[x]+1;
				q.push(v);
			}
		}
	}
	return d[T]!=-1;
}
int dfs(int st,int ed,int nowdat){
	int dat=0;
	if(st==ed)
		return nowdat;
	for(int i=first[st];i;i=e[i].nxt){
		int v=e[i].v;
		if(e[i].w>0&&d[v]==d[st]+1){
			int now=min(nowdat-dat,e[i].w);
			now=dfs(v,ed,now);
			dat+=now;
			e[i].w-=now;
			e[i^1].w+=now;
		}
	}
	if(!dat){
		d[st]=-2;
	}
	return dat;
}
int Dinic(){
	int ans=0;
	while(bfs()){
//		cout<<"#$@#"<<endl;
		ans+=dfs(S,T,0x3f3f3f3f);
	}
	return ans;
}
int main(){
	scanf("%d%d%d",&n,&m,&k);n+=2;
	for(int i=1;i<=n;i++){
		fa[i]=i;
	}
	for(int i=1;i<=m;i++){
		scanf("%d",&w[i]);
		scanf("%d",&a[i][0]);
		for(int j=1;j<=a[i][0];j++){
			scanf("%d",&a[i][j]);
			a[i][j]+=2;
			if(j!=1){
				uni(a[i][j],a[i][j-1]);
			}
		}
	}
	if(getfa(1)!=getfa(2)){
		cout<<0<<endl;
		return 0;
	}
	add(S,2,k);
	for(int i=1;i<=m;i++){
		C[i]=1;
	}
	int nowdat=0;
	int stream=0;
	int ans=0;
	while(nowdat<k){
//		cout<<ans<<" "<<nowdat<<endl;
		ans++;
		for(int i=1;i<=n;i++){
			add(i+(ans-1)*n,i+ans*n,INF);
		}
		for(int i=1;i<=m;i++){
			int p=C[i]+1;
			if(p>a[i][0]){
				p=1;
			}
			add(a[i][C[i]]+n*(ans-1),a[i][p]+n*ans,w[i]);
			C[i]=p;
		}
		T=ans*n+1;
		nowdat+=Dinic();
//		cout<<ans<<" "<<nowdat<<endl;
	}
	cout<<ans<<endl;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值