网络流二十四题 ————(十六）、P2754 [CTSC1999]家园 / 星际转移问题 残留网络建边，最大流+并查集

分层图，每层有n+2个点，表示n个太空站和月球地球。

第i曾表示时间i。

然后根据太阳飞船运动轨迹连边建图。跑最大流即可。

刚开始用并查集判断是否可行。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int M = 2500000+7;
const int N = 100000+7;
struct Dinic{
	//N个点，M条边 
	//add建双向边，然后D.gao ,最后输出maxflow  
	#define inf 0x3f3f3f3f
	ll maxflow;int s,t,n;
	int v[N],pre[N],d[N],now[N];
	ll incf[N];
	int head[N],cnt=1;
	struct EDGE{int to,nxt;ll w;}ee[M*2];
	inline void AD(int x,int y,ll w){ee[++cnt].nxt=head[x],ee[cnt].w=w,ee[cnt].to=y,head[x]=cnt;}
	inline void add(int x,int y,ll w){
	AD(x,y,w);AD(y,x,0);}
	inline bool bfs()//在残量网络上构造分层图 
	{
		for(int i=0;i<=n;i++)d[i]=0;
		queue<int>q;
		q.push(s);d[s]=1;
		while(q.size())
		{
			int x=q.front();q.pop();
			for(int i=head[x];i;i=ee[i].nxt)
			{
				int y=ee[i].to;ll w=ee[i].w;
				if(w&&!d[y])
				{
					
					q.push(y);
					d[y]=d[x]+1;
					if(y==t)return 1;
				}
			}
		}
		return false;
	}
	inline int dinic(int x,int flow)
	{
		if(x==t)return flow;
		ll rest = flow,k;
		for(int i=now[x];i&&rest;i=ee[i].nxt)
		{
			int y=ee[i].to;ll w= ee[i].w;
			now[x]=i;
			if(w&&d[y]==d[x]+1)
			{
				k=dinic(y,min(rest,w));
				if(!k)d[y]=0;//剪枝，去掉增广完毕的点
				ee[i].w-=k;
				ee[i^1].w+=k;
				rest-=k; 
			}
		}
		return flow - rest;
	}
	inline void gao()
	{
		int flow=0;
		while(bfs())
		{
			for(int i=0;i<=n;i++)now[i]=head[i];
			while(flow=dinic(s,inf))maxflow+=flow;
		}
	}
	inline void init(int nn,int S,int T)
	{
		cnt=1;maxflow=0;
		for(int i=0;i<=n;i++)head[i]=0;
		s=S,t=T,n=nn;
	}
}D;
//求解二分图匹配问题时的复杂度是：m*sqrt(n)的
int to[22][22];
int w[22];
int nm[22];//每个飞船的周期 
int fa[22];
int gt(int x){
	if(x==fa[x])return x;
	return fa[x]=gt(fa[x]);
}
int dq,yq;
int n,m,k,l;
int id(int x,int t){
	return (n+2)*t+x+2;//S,T占2个节点 
}
int main()
{
	
	cin>>n>>m>>k;
	dq=n+1,yq=n+2;
	for(int i=1;i<=n+2;i++)fa[i]=i;
	for(int i=1;i<=m;i++){
		cin>>w[i]>>nm[i];
		for(int j=0;j<nm[i];j++){
			cin>>to[i][j];
			if(to[i][j]==0)to[i][j]=dq;
			if(to[i][j]==-1)to[i][j]=yq;
			if(j==0)continue;
			int x=to[i][j-1],y=to[i][j];
			int gx=gt(to[i][j-1]),gy=gt(to[i][j]);
			fa[gx]=gy;
		}
	}
	if(gt(n+1)!=gt(n+2)){
		cout<<0<<endl;
		return 0;
	}
	int s=1,t=2;
	D.init(n+4,s,t);
	D.add(s,id(dq,0),inf);
	D.add(id(yq,0),t,inf);
	
	for(int T=1;;T++){
		for(int i=1;i<=m;i++){
			D.add(id(to[i][(T-1)%nm[i]],T-1),id(to[i][T%nm[i]],T),w[i]);
		}
		for(int i=1;i<=n+2;i++)D.add(id(i,T-1),id(i,T),inf);
		D.add(s,id(dq,T),inf);
		D.add(id(yq,T),t,inf);
		D.n+=n+2;
		D.gao();
		if(D.maxflow>=k){
			cout<<T<<endl;
			break;
		}
		
	}
	return 0;
}
/*
1 1 1
1 2 0 -1
*/